Abstract visibility
Abstract text size

Papers for Friday, Jul 15 2022

Laser heterodyne interferometry plays a key role in the proof mass's monitor and control by measuring its multiple degrees of freedom motions in the Space Gravitational Wave Detection. Laboratory development of polarization-multiplexing heterodyne interferometer (PMHI) using quadrant photodetectors (QPD) is presented in this paper, intended for measuring the translation and tilt of a proof mass. The system is of symmetric design, which can expand to five degrees of freedom measurements based on polarization-multiplexing and differential wavefront sensing (DWS). The ground-simulated experimental results demonstrate that a measurement noise of 3 pm/Hz$^{1/2}$ and 2 nrad/Hz$^{1/2}$ at 1 Hz have been achieved respectively. The tilt-to-length error is dominated by geometric misalignment for the current system, the coupling of which is at micrometer level within a tilt range of 1000 {\mu}rad.

Hengxiao Guo, Aaron J. Barth, Shu Wang

25 pages, 11 figures, 4 tables. Submitted to ApJ, comments welcome

We perform a systematic survey of active galactic nucleus (AGN) continuum lags using $\sim$3 day cadence $gri$-band light curves from the Zwicky Transient Facility. We select a sample of 94 type~1 AGN at $z<0.8$ with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the lamp-post reprocessing model, our findings are: 1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin disk model; 2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks; 3) The CE size scales with continuum luminosity as $R_{\mathrm{CE}} \propto L^{0.48\pm0.04}$ with a scatter of 0.2 dex, analogous to the well-known radius-luminosity relation of broad H$\mathrm{\beta}$. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the $R_{\mathrm{CE}}-L$ relation can be explained by a photoionization model that assumes $\sim$23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anti-correlates with the continuum luminosity, indicative of a potential non-disk BLR lag contribution evolving with luminosity. Finally, a robust positive correlation between CE size and black hole mass is detected.

Stephen R. Kane, Zhexing Li

12 pages, 5 figures, accepted for publication in the Planetary Science Journal

Planetary rings are often speculated as being a relatively common attribute of giant planets, partly based on their prevalence within the Solar System. However, their formation and sustainability remain a topic of open discussion, and the most massive planet within our planetary system harbors a very modest ring system. Here, we present the results of a N-body simulation that explores dynamical constraints on the presence of substantial ring material for Jupiter. Our simulations extend from within the rigid satellite Roche limit to 10\% of the Jupiter Hill radius, and include outcomes from $10^6$ and $10^7$ year integrations. The results show possible regions of a sustained dense ring material presence around Jupiter that may comprise the foundation for moon formation. The results largely demonstrate the truncation of stable orbits imposed by the Galilean satellites, and dynamical desiccation of dense ring material within the range $\sim$3--29 Jupiter radii. We discuss the implications of these results for exoplanets, and the complex relationship between the simultaneous presence of rings and massive moon systems.

A. Vigan, M. Lopez, M. El Morsy, E. Muslimov, A. Viret, G. Zins, G. Murray, A. Costille, G. P. P. L. Otten, U. Seemann, H. Anwand-Heerwart, K. Dohlen, P. Blanchard, J. Garcia, Y. Charles, N. Tchoubaklian, T. Ely, M. Phillips, J. Paufique, J.-L. Beuzit, M. Houllé, J. Costes, R. Pourcelot, I. Baraffe, R. Dorn, M. Jaquet, M. Kasper, A. Reiners, A. Smette, L. Blanco, L. Pallanca, A. Carlotti, É. Choquet, D. Mouillet, M. N'Diaye

Submitted to SPIE Astronomical Telescopes + Instrumentation 2022, Adaptive Optics Systems VIII, Paper 12185-27

New generation exoplanet imagers on large ground-based telescopes are highly optimised for the detection of young giant exoplanets in the near-infrared, but they are intrinsically limited for their characterisation by the low spectral resolution of their integral field spectrographs ($R<100$). High-dispersion spectroscopy at $R \gg 10^4$ would be a powerful tool for the characterisation of these planets, but there is currently no high-resolution spectrograph with extreme adaptive optics and coronagraphy that would enable such characterisation. With project HiRISE we propose to use fiber coupling to combine the capabilities of two flagship instruments at the Very Large Telescope in Chile: the exoplanet imager SPHERE and the high-resolution spectrograph CRIRES+. The coupling will be implemented at the telescope in early 2023. We provide a general overview of the implementation of HiRISE, of its assembly, integration and testing (AIT) phase in Europe, and a brief assessment of its expected performance based on the final hardware.

Marcos Pellejero-Ibanez, Raul E. Angulo, Matteo Zennaro, Jens Stuecker, Sergio Contreras, Giovanni Arico, Francisco Maion

16 pages, 9 figures

We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We build the emulator using a forward neural network trained with the simulations of the BACCO project, which covers an 8-dimensional parameter space including massive neutrinos and dynamical dark energy. We show that our emulator provides unbiased cosmological constraints from the monopole, quadrupole, and hexadecapole of a mock galaxy catalogue that mimics the BOSS-CMASS sample down to nonlinear scales ($k\sim0.6$[$h/$Mpc]$^{3}$). This work opens up the possibility of robustly extracting cosmological information from small scales using observations of the large-scale structure of the Universe.

S. H. Price, H. Übler, N. M. Förster Schreiber, P. T. de Zeeuw, A. Burkert, R. Genzel, L. J. Tacconi, R. I. Davies, C. P. Price

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

Using kinematics to decompose galaxies' mass profiles, including the dark matter contribution, often requires parameterization of the baryonic mass distribution based on ancillary information. One such model choice is a deprojected S\'ersic profile with an assumed intrinsic geometry. The case of flattened, deprojected S\'ersic models has previously been applied to flattened bulges in local star-forming galaxies (SFGs), but can also be used to describe the thick, turbulent disks in distant SFGs. Here we extend this previous work that derived density ($\rho$) and circular velocity ($v_{\rm circ}$) curves by additionally calculating the spherically-enclosed 3D mass profiles ($M_{\rm sph}$). Using these profiles, we compare the projected and 3D mass distributions, quantify the differences between the projected and 3D half-mass radii ($R_{\rm e}; r_{\rm 1/2,mass,3D}$), and present virial coefficients relating $v_{\rm circ}(R)$ and $M_{\rm sph}(<r=R)$ or $M_{\rm tot}$. We then quantify differences between mass fraction estimators for multi-component systems, particularly for dark matter fractions, and consider the compound effects of measuring dark matter fractions at the projected versus 3D half-mass radii. While the fraction estimators produce only minor differences, using different aperture radius definitions can strongly impact the inferred dark matter fraction. As pressure support is important in analysis of gas kinematics (particularly at high redshifts), we also calculate the self-consistent pressure support correction profiles, which generally predict less pressure support than for the self-gravitating disk case. These results have implications for comparisons between simulation and observational measurements, and for the interpretation of SFG kinematics at high redshifts. A set of precomputed tables and the code to calculate the profiles are made publicly available. [Abridged]

M.A. Fernandez, Ming-Feng Ho, Simeon Bird

Submitted to MNRAS. 12 pages, 10 figures

In this work we extend our recently developed multi-fidelity emulation technique to the simulated Lyman-$\alpha$ forest flux power spectrum. Multi-fidelity emulation allows interpolation of simulation outputs between cosmological parameters using many cheap low-fidelity simulations and a few expensive high-fidelity simulations. Using a test suite of small box (30 Mpc/h) simulations, we show that multi-fidelity emulation is able to reproduce the Lyman-$\alpha$ forest flux power spectrum well, achieving an average accuracy when compared to a test suite of $0.8\%$. We further show that it has a substantially increased accuracy over single-fidelity emulators, constructed using either the high or low-fidelity simulations only. In particular, it allows the extension of an existing simulation suite to smaller scales and higher redshifts.

Jaikhomba Singha, Bhal Chandra Joshi, Debades Bandyopadhyay, Himanshu Grover, Shantanu Desai, P. Arumugam, Sarmistha Banik

7 figure. Accepted for publication in the Journal of Astrophysics and Astronomy : SKA special issue

There are two types of timing irregularities seen in pulsars: glitches and timing noise. Both of these phenomena can help us probe the interior of such exotic objects. This article presents a brief overview of the observational and theoretical aspects of pulsar timing irregularities and the main results from the investigations of these phenomena in India. The relevance of such Indian programs for monitoring of young pulsars with the Square Kilometer Array (SKA) is presented, highlighting possible contributions of the Indian neutron star community to the upcoming SKA endeavour.

Toshiki Saito, Shuro Takano, Nanase Harada, Taku Nakajima, Eva Schinnerer, Daizhong Liu, Akio Taniguchi, Takuma Izumi, Yumi Watanabe, Kazuharu Bamba, Kotaro Kohno, Yuri Nishimura, Sophia Stuber, Tomoka Tosaki

15 pages, 6 figures, 2 tables, accepted for publication in ApJ

Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16 well-detected molecular line intensity maps convolved to a common 150 pc resolution. In addition, we include the [SIII]/[SII] line ratio and [CI] $^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show remarkable resemblance with that of a kpc-scale biconical outflow from the central AGN. We identify two prominent features: (1) central concentration at the circumnuclear disk (CND) and (2) two peaks across the center that coincide with the biconical outflow peaks. The concentrated molecular lines in the CND are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line emissions from molecules known to be enhanced in irradiated interstellar medium, CN, C$_2$H, and HNC, show similar concentrations and extended components along the bicone, suggesting that molecule dissociation is a dominant chemical effect of the cold molecular outflow of this galaxy. Although further investigation should be made, this scenario is consistent with the faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons and electrons.

Margaret R. Rybak, Penina Axelrad, Catie LeDesma, Dana Z. Anderson, Todd Ely

30 pages, 8 figures

High-sensitivity shaken lattice interferometry (SLI) based sensors have the potential to provide deep space missions with the ability to precisely measure non-gravitational perturbing forces. This work considers the simulation of the OSIRIS-REx mission navigation in the vicinity of Bennu with the addition of measurements from onboard SLI-based accelerometers. The simulation is performed in the Jet Propulsion Laboratory's (JPL) Mission Analysis, Operations and Navigation Toolkit (MONTE) and incorporates OSIRIS-REx reconstructed trajectory and attitude data from the Navigation and Ancillary Information Facility (NAIF) database. The use of the reconstructed data from NAIF provides realistic true dynamical errors and JPL's MONTE software allows for a high-fidelity simulation of a nominal reference for the filter. The navigation performance and reduction of tracking and complex modeling enabled by the onboard SLI-based sensor are presented for two orbital phases of the OSIRIS-REx mission. Overall, the results show that the addition of SLI-based accelerometer measurements improves navigation performance, when compared to a radiometric tracking only configuration. In addition, results demonstrate that highly-precise accelerometer measurements can effectively replace at least one day of DSN passes over a three-day period, thereby reducing tracking requirements. Furthermore, it is shown that lower-fidelity surface force modeling and parameter estimation is required when using onboard SLI-based accelerometers.

Bhal Chandra Joshi (1), Achamveedu Gopakumar (2), Arul Pandian (3), Thiagaraj Prabu (3), Lankeswar Dey (2), Manjari Bagchi (4, 5), Shantanu Desai (6), Pratik Tarafdar (4), Prerna Rana (2), Yogesh Maan (1), Neelam Dhanda Batra (7), Raghav Girgaonkar (8), Nikita Agarwal (9), Paramasivan Arumugam (10), Sarmistha Banik (21) Avishek Basu (11), Adarsh Bathula (12), Subhajit Dandapat (2), Yashwant Gupta (1), Shinnosuke Hisano (13), Ryo Kato (14, 15), Divyansh Kharbanda (6), Tomonosuke Kikunaga (13), Neel Kolhe (16), M. A. Krishnakumar (17, 18), P. K. Manoharan (19), Piyush Marmat (10), Arun Naidu (20), K. Nobleson (21), Avinash Kumar Paladi (22), Dhruv Pathak (23), Jaikhomba Singha (10), Aman Srivastava (6), Mayuresh Surnis (11), Sai Chaitanya Susarla (24), Abhimanyu Susobhanan (25, 1), Keitaro Takahashi (26, 27), ((1) National Centre for Radio Astrophysics India, (2) Tata Institute of Fundamental Research India, (3) Raman Research Institute India, (4) The Institute of Mathetical Sciences India, (5) Homi Bhabha National Institute India, (6) IIT Hyderabad India, (7) University of Delhi India, (8) Amity University India, (9) Manipal Institute of Technology India, (10) Indian Institute of Technology Roorkee India, (11) Jodrell Bank Centre for Astrophysics UK, (12) The Indian Institute of Science Education and Research India, (13) Kumamoto University Japan, (14) Faculty of Advanced Science and Technology Kumamoto University Japan, (15) Osaka City University Advanced Mathematical Institute Japan, (16) St Xavier's College India, (17) Max-Planck-Institut for Radioastronomie Germany, (18) Universitat Bielefeld Germany, (19) Arecibo Observatory University of Central Florida USA, (20) University of Oxford UK, (21) BITS Pilani Hyderabad Campus India, (22) Indian Institute of Space Science and Technology Thiruvananthapuram India, (23) Inter-University Centre for Astronomy and Astrophysics India, (24) National University of Ireland Galway Ireland, (25) National Astronomical Observatories Chinese Academy of Sciences China, (26) Faculty of Advanced Science and Technology Kumamoto University Japan, (27) International Research Organization for Advanced Science and Technology Kumamoto University Japan)

Accepted for publication in Journal of Astronomy and Astrophysics for Special Issue on Indian Participation in the SKA (Editors : Abhirup Datta, Nirupam Roy, Preeti Kharb and Tirthankar Roy Choudhury)

Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsars for this purpose. We highlight InPTA's observation strategies and analysis methods, which are relevant for a future PTA experiment with the more sensitive Square Kilometer Array (SKA) telescope. We show that the unique multi-sub-array multi-band wide-bandwidth frequency coverage of the InPTA provides Dispersion Measure estimates with unprecedented precision for PTA pulsars, e.g., ~ 2 x 10{-5} pc-cm{-3} for PSR J1909-3744. Configuring the SKA-low and SKA-mid as two and four sub-arrays respectively, it is shown that comparable precision is achievable, using observation strategies similar to those pursued by the InPTA, for a larger sample of 62 pulsars requiring about 26 and 7 hours per epoch for the SKA-mid and the SKA-low telescopes respectively. We also review the ongoing efforts to develop PTA-relevant general relativistic constructs that will be required to search for nanohertz gravitational waves from isolated super-massive black hole binary systems like blazar OJ 287. These efforts should be relevant to pursue persistent multi-messenger gravitational wave astronomy during the forthcoming era of the SKA telescope, the Thirty Meter Telescope, and the next-generation Event Horizon Telescope.

E.V. Arbuzova, A.D. Dolgov

21 pages, 4 figures, 1 table. Extended version of the invited talk presented at International Workshop on Elementary Particles and Nuclear Physics, 24 -- 30 April 2022, Institute of Nuclear Physics of the Ministry of Energy of the Republic of Kazakhstan, Almaty

Universe heating in $R^2$-modified gravity is considered. The rates of particle production by the scalaron are calculated for different decay channels. Freezing of massive stable relics with the interaction strength typical for supersymmetry is studied. It is shown that the bounds on masses of supersymmetry-kind particles allowing them to form the cosmological dark matter (DM)depend upon the dominant decay mode of the scalaron. In any case the results presented open much wider mass window for DM with the interaction strength typical for supersymmetry.

L. Bonne, N. Schneider, P. García, A. Bij, P. Broos, L. Fissel, R. Guesten, J. Jackson, R. Simon, L. Townsley, A. Zavagno, R. Aladro, C. Buchbender, C. Guevara, R. Higgins, A. M. Jacob, S. Kabanovic, R. Karim, A. Soam, J. Stutzki, M. Tiwari, F. Wyrowski, A. G. G. M. Tielens

38 pages, 27 figures, 8 tables, accepted in ApJ

We present [CII] 158 $\mu$m and [OI] 63 $\mu$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE).

Malena Rice, Songhu Wang, Xian-Yu Wang, Gudmundur Stefansson, Howard Isaacson, Andrew W. Howard, Sarah E. Logsdon, Heidi Schweiker, Fei Dai, Casey Brinkman, Steven Giacalone, Rae Holcomb

Accepted for publication in AJ. 20 pages, 7 figures, 4 tables

The distribution of spin-orbit angles for systems with wide-separation, tidally detached exoplanets offers a unique constraint on the prevalence of dynamically violent planetary evolution histories. Tidally detached planets provide a relatively unbiased view of the primordial stellar obliquity distribution, since they cannot tidally realign within the system lifetime. We present the third result from our Stellar Obliquities in Long-period Exoplanet Systems (SOLES) survey: a measurement of the Rossiter-McLaughlin effect across two transits of the tidally detached warm Jupiter TOI-1478 b with the WIYN/NEID and Keck/HIRES spectrographs, revealing a sky-projected spin-orbit angle $\lambda=6.2^{+5.9}_{-5.5}$ degrees. Combining this new measurement with the full set of archival obliquity measurements, including two previous constraints from the SOLES survey, we demonstrate that, in single-star systems, tidally detached warm Jupiters are preferentially more aligned than closer-orbiting hot Jupiters. This finding has two key implications: (1) planets in single-star systems tend to form within aligned protoplanetary disks, and (2) warm Jupiters form more quiescently than hot Jupiters, which, in single-star systems, are likely perturbed into a misaligned state through planet-planet interactions in the post-disk-dispersal phase. We also find that lower-mass Saturns span a wide range of spin-orbit angles, suggesting a prevalence of planet-planet scattering and/or secular mechanisms in these systems.

Erik Petigura, Andrew Howard, Jacob Bean, Charles Beichman, Debra Fischer, BJ Fulton, Howard Isaacson, John O'Meara, Carolyn Jordan, Daniel Huber, Paul Robertson, Arpita Roy, Johanna Teske, Josh Walawender

31 pages, 10 figures, prepared for the Keck Science Steering Committee, report delivered on Nov 7, 2021

The arrival of the Keck Planet Finder (KPF) in 2022 represents a major advance in the precision radial velocity (PRV) capabilities of the W. M. Keck Observatory. In preparation for KPF science, our committee of PRV experts and WMKO staff studied the current implementation of cadence observing at Keck and other PRV facilities. We find that many of KPF's major science cases are not feasible through Keck's standard allocations of full or half nights to individual PIs. Pooling time among several PIs as is currently done by the California Planet Search (CPS) collaboration with HIRES results in lower quality science results than is possible when KPF is available at higher observational cadence. This strategy also creates barriers to entry, particularly for researchers wishing to lead small proposals. This report makes recommendations for optimizing PRV cadence at Keck subject to the following constraints: preservation of clear boundaries between cadence observations and classically scheduled time; and ensuring fairness and scientific independence of different Keck TACs and different KPF PIs. We recommend establishing a new category of Keck time allocation, "KPF Community Cadence" (KPF-CC). In many ways, KPF-CC will formalize observing strategies provided by CPS, but with higher observational cadence appropriate for KPF science and with universal access to the program for all Keck users. We recommend that KPF-CC time be scheduled classically into blocks as small as a quarter night subject to considerations of bright/dark time, variations in proposal pressure with the seasons, and the needs of non-KPF observing programs. Within KPF-CC time, the Keck Observing Assistants would execute observations generated by a dynamic scheduler. We recommend that Keck staff and a board of PRV experts design and maintain the scheduling software.

S. B. Yun, J. M. Miller, D. Barret, D. Stern, W. N. Brandt, L. Brenneman, P. Draghis, A. C. Fabian, J. Raymond, A. Zoghbi

NGC 5033 is an intriguing Seyfert galaxy because its sub-classification may change with time, and because optical and sub-mm observations find that the massive black hole does not sit at the dynamical center of the galaxy, pointing to a past merger. We obtained a new optical spectrum of NGC 5033 using the 200'' Hale telescope at Palomar that clearly reveals a broad H$\beta$ line (FWHM$=5400\pm 300~{\rm km}~{\rm s}^{-1}$). This signals a clear view of the optical broad line region (BLR) and requires Seyfert-1.5 designation. Some spectra obtained in the past suggest a Seyfert-1.9 classification, potentially signaling a variable or "changing-look" geometry. Our analysis of a 2019 Chandra spectrum of the massive black hole reveals very little obscuration, also suggesting a clean view of the central engine. However, the narrow Fe~K$\alpha$ emission line is measured to have an equivalent with of EW$=460^{+100}_{-90}$~eV. This value is extremely high compared to typical values in unobscured AGN. Indeed, the line is persistently strong in NGC 5033: the line equivalent width in a 2002 XMM-Newton snapshot is EW$=250^{+40}_{-40}$~eV, similar to the EW$=290^{+100}_{-100}$~eV equivalent width measured using ASCA in 1999. These results can likely be explained through a combination of an unusually high covering factor for reflection, and fluxes that are seen out of phase owing to light travel times. We examine the possibility that NGC 5033 may strengthen evidence for the X-ray Baldwin effect.

Yifan Zhou, Aniket Sanghi, Brendan P. Bowler, Ya-Lin Wu, Laird M. Close, Feng Long, Kimberly Ward-Duong, Zhaohuan Zhu, Adam L. Kraus, Katherine B. Follette, Jaehan Bae

Accepted for publication in ApJL

Accreting protoplanets enable the direct characterization of planet formation. As part of a high-contrast imaging search for accreting planets with the Hubble Space Telescope (HST) Wide Field Camera 3, we present H$\alpha$ images of AB Aurigae (AB Aur), a Herbig Ae/Be star harboring a transition disk. The data were collected in two epochs of direct-imaging observations using the F656N narrow-band filter. After subtracting the point spread function of the primary star, we identify a point-like source located at a P.A. of $182.5^{\circ}\pm1.4^{\circ}$ and a separation of $600\pm22$~mas relative to the host star. The position is consistent with the recently identified protoplanet candidate AB Aur b. The source is visible in two individual epochs separated by ${\sim}50$ days and the H$\alpha$ intensities in the two epochs agree. The H$\alpha$ flux density is $F_{\nu}=1.5\pm0.4$~mJy, $3.2\pm0.9$ times of the optical continuum determined by published HST/STIS photometry. In comparison to PDS 70 b and c, the H$\alpha$ excess emission is weak. The central star is accreting and the stellar H$\alpha$ emission has a similar line-to-continuum ratio as seen in AB Aur b. We conclude that both planetary accretion and scattered stellar light are possible sources of the H$\alpha$ emission, and the H$\alpha$ detection alone does not validate AB Aur b as an accreting protoplanet. Disentangling the origin of the emission will be crucial for probing planet formation in the AB Aur disk.

Stefano Anselmi, Matthew F. Carney, John T. Giblin Jr, Saurabh Kumar, James B. Mertens, Marcio ODwyer, Glenn D. Starkman, Chi Tian

The Universe is neither homogeneous nor isotropic, but it is close enough that we can reasonably approximate it as such on suitably large scales. The inflationary-$\Lambda$-Cold Dark Matter ($\Lambda$CDM) concordance cosmology builds on these assumptions to describe the origin and evolution of fluctuations. With standard assumptions about stress-energy sources, this system is specified by just seven phenomenological parameters, whose precise relations to underlying fundamental theories are complicated and may depend on details of those fields. Nevertheless, it is common practice to set the parameter that characterizes the spatial curvature, $\Omega_K$, exactly to zero. This parameter-fixed $\Lambda$CDM is awarded distinguished status as separate model, "flat $\Lambda$CDM.'' Ipso facto this places the onus on proponents of "curved $\Lambda$CDM'' to present sufficient evidence that $\Omega_K\neq0$, and is needed as a parameter. While certain inflationary model Lagrangians, with certain values of their parameters, and certain initial conditions, will lead to a present-day universe well-described as containing zero curvature, this does not justify distinguishing that subset of Lagrangians, parameters and initial conditions into a separate model. Absent any theoretical arguments, we cannot use observations that suggest small $\Omega_K$ to enforce $\Omega_K=0$. Our track record in picking inflationary models and their parameters a priori makes such a choice dubious, and concerns about tensions in cosmological parameters and large-angle cosmic-microwave-background anomalies strengthens arguments against this choice. We argue that $\Omega_K$ must not be set to zero, and that $\Lambda$CDM remains a phenomenological model with at least 7 parameters.

Kelsi N. Singer, Oliver L. White, Bernard Schmitt, Erika L. Rader, Silvia Protopapa, William M. Grundy, Dale P. Cruikshank, Tanguy Bertrand, Paul M. Schenk, William B. McKinnon, S. Alan Stern, Rajani D. Dhingra, Kirby D. Runyon, Ross A. Beyer, Veronica J. Bray, Cristina Dalle Ore, John R. Spencer, Jeffrey M. Moore, Francis Nimmo, James T. Keane, Leslie A. Young, Catherine B. Olkin, Tod R. Lauer, Harold A. Weaver, Kimberly Ennico-Smith

22 pages, including both main paper and supplement as one pdf

The New Horizons spacecraft returned images and compositional data showing that terrains on Pluto span a variety of ages, ranging from relatively ancient, heavily cratered areas to very young surfaces with few-to-no impact craters. One of the regions with very few impact craters is dominated by enormous rises with hummocky flanks. Similar features do not exist anywhere else in the imaged solar system. Here we analyze the geomorphology and composition of the features and conclude this region was resurfaced by cryovolcanic processes, of a type and scale so far unique to Pluto. Creation of this terrain requires multiple eruption sites and a large volume of material (>104 km^3) to form what we propose are multiple, several-km-high domes, some of which merge to form more complex planforms. The existence of these massive features suggests Pluto's interior structure and evolution allows for either enhanced retention of heat or more heat overall than was anticipated before New Horizons, which permitted mobilization of water-ice-rich materials late in Pluto's history.

J. S. Halekas, P. Whittlesey, D. E. Larson, M. Maksimovic, R. Livi, M. Berthomier, J. C. Kasper, A. W. Case, M. L. Stevens, S. D. Bale, R. J. MacDowall, M. P. Pulupa

Submitted to the Astrophysical Journal

We utilize observations from the Parker Solar Probe (PSP) to study the radial evolution of the solar wind in the inner heliosphere. We analyze electron velocity distribution functions observed by the Solar Wind Electrons, Alphas, and Protons suite to estimate the coronal electron temperature and the local electric potential in the solar wind. From the latter value and the local flow speed, we compute the asymptotic solar wind speed. We group the PSP observations by asymptotic speed, and characterize the radial evolution of the wind speed, electron temperature, and electric potential within each group. In agreement with previous work, we find that the electron temperature (both local and coronal) and the electric potential are anti-correlated with wind speed. This implies that the electron thermal pressure and the associated electric field can provide more net acceleration in the slow wind than in the fast wind. We then utilize the inferred coronal temperature and the extrapolated electric + gravitational potential to show that both electric field driven exospheric models and the equivalent thermally driven hydrodynamic models can explain the entire observed speed of the slowest solar wind streams. On the other hand, neither class of model can explain the observed speed of the faster solar wind streams, which thus require additional acceleration mechanisms.

Riley M. T. Connors, Javier A. Garcia, John Tomsick, Guglielmo Mastroserio, Victoria Grinberg, James F. Steiner, Jiachen Jiang, Andrew C. Fabian, Michael L. Parker, Fiona Harrison, Jeremy Hare, Labani Mallick, Hadar Lazar

18 pages, 9 figures, accepted by ApJ

The degree to which the thin accretion disks of black hole X-ray binaries are truncated during hard spectral states remains a contentious open question in black hole astrophysics. During its singular observed outburst in $2009\mbox{--}2010$, the black hole X-ray binary XTE J1752-223 spent $\sim1$~month in a long-stable hard spectral state at a luminosity of $\sim0.02\mbox{--}0.1~L_{\rm Edd}$. It was observed with 56 RXTE pointings during this period, with simultaneous Swift-XRT daily coverage during the first 10 days of the RXTE observations. Whilst reflection modeling has been extensively explored in the analysis of these data, there is a disagreement surrounding the geometry of the accretion disk and corona implied by the reflection features. We re-examine the combined, high signal-to-noise, simultaneous Swift and RXTE observations, and perform extensive reflection modeling with the latest relxill suite of reflection models, including newer high disk density models. We show that reflection modeling requires that the disk be within $\sim5~R_{\rm ISCO}$ during the hard spectral state, whilst weaker constraints from the thermal disk emission imply higher truncation ($R_{\rm in}=6\mbox{--}80~R_{\rm ISCO}$). We also explore more complex coronal continuum models, allowing for two Comptonization components instead of one, and show that the reflection features still require only a mildly truncated disk. Finally we present a full comparison of our results to previous constraints found from analyses of the same dataset.

Terrestrial planets currently in the habitable zones around M dwarfs likely experienced a long-term runaway greenhouse condition because of a slow decline in host-stellar luminosity in its pre-main sequence phase. Accordingly, they might have lost significant portions of their atmospheres including water vapor at high concentration by hydrodynamic escape induced by the strong stellar XUV irradiation. However, the atmospheric escape rates remain highly uncertain due partly to a lack of understanding of the effect of radiative cooling in the escape outflows. Here we carry out 1-D hydrodynamic escape simulations for an H$_{2}$-H$_{2}$O atmosphere on a planet with mass of $1M_{\oplus}$ considering radiative and chemical processes to estimate the atmospheric escape rate and follow the atmospheric evolution during the early runaway greenhouse phase. We find that the atmospheric escape rate decreases with the basal H$_{2}$O/H$_{2}$ ratio due to the energy loss by the radiative cooling of H$_{2}$O and chemical products such as OH and H$_{3}^{+}$: the escape rate of H$_{2}$ becomes one order of magnitude smaller when the basal H$_{2}$O/H$_{2}=0.1$ than that of the pure hydrogen atmosphere. The timescale for H$_{2}$ escape exceeds the duration of the early runaway greenhouse phase, depending on the initial atmospheric amount and composition, indicating that H$_{2}$ and H$_{2}$O could be left behind after the end of the runaway greenhouse phase. Our results suggest that temperate and reducing environments with oceans could be formed on some terrestrial planets around M dwarfs.

The bubble nebula surrounding NGC 1313 X-2 is believed to be powered by high velocity winds from the central ultraluminous X-ray source (ULX) as a result of supercritical accretion. With the Multi-Unit Spectroscopic Explorer (MUSE) observation of the nebula, we find enhanced OIII emission at locations spatially coincident with clusters of stars and the central X-ray source, suggesting that photoionization in addition to shock-ionization plays an important role in powering the nebula. The X-ray luminosity of the ULX and the number of massive stars in the nebula region can account for the required ionizing luminosity derived with MAPPINGS V, which also confirms that pure shocks cannot explain the observed emission line ratios.

Kai Wang, Hai-Ming Zhang, Ruo-Yu Liu, Xiang-Yu Wang

11 pages, 5 figures and 1 table. Accepted for publication in ApJ

Isotopic and elemental abundances seen in Galactic cosmic rays imply that $\sim20\%$ of the cosmic-ray (CR) nuclei are probably synthesized by massive Wolf-Rayet (WR) stars. Massive star clusters hosting WR and OB-type stars have been proposed as potential Galactic cosmic-ray accelerators for decades, in particular via diffusive shock acceleration at wind termination shocks. Here we report the analysis of {\em Fermi} Large Area Telescope's data towards the direction of Masgomas-6a, a young massive star cluster candidate hosting two WR stars. We detect an extended $\gamma$-ray source with $\rm{TS}=183$ in the vicinity of Masgomas-6a, spatially coincident with two unassociated {\em Fermi} 4FGL sources. We also present the CO observational results of molecular clouds in this region, using the data from the Milky Way Imaging Scroll Painting project. The $\gamma$-ray emission intensity correlates well with the distribution of molecular gas at the distance of Masgomas-6a, indicating that these gamma rays may be produced by CRs accelerated by massive stars in Masgomas-6a. At the distance of $3.9{\rm \ kpc}$ of Masgomas-6a, the luminosity of the extended source is $(1.81\pm0.02)\times 10^{35}{\rm \ erg \ s^{-1}}$. With a kinetic luminosity of $\sim 10^{37}{\rm erg \ s^{-1}}$ in the stellar winds, the WR stars are capable of powering the $\gamma$-ray emission via neutral pion decay resulted from cosmic ray $pp$ interactions. The size of the GeV source and the energetic requirement suggests a CR diffusion coefficient smaller than that in the Galactic interstellar medium, indicating strong suppression of CR diffusion in the molecular cloud.

Lei Yang, Xinwen Shu, Fabao Zhang, Yogesh Chandola, Daizhong Liu, Yi Liu, Minfeng Gu, Margherita Giustini, Ning Jiang, Ya-Ping Li, Di Li, David Elbaz, Stephanie Juneau, Maurilio Pannella, Luming Sun, Ningyu Tang, Tinggui Wang, Hongyan Zhou

11 pages, 3 figures, 2 tables. Accepted for publication in ApJ

RX J1301.9+2747 is a unique active galaxy with supersoft X-ray spectrum that lacks significant emission at energies above 2 keV. In addition, it is one of few galaxies displaying quasi-periodic X-ray eruptions that recur on a timescale of 13-20 ks. We present multi-epoch radio observations of RX J1301.9+2747 using GMRT, VLA and VLBA. The VLBA imaging at 1.6 GHz reveals a compact radio emission unresolved at a scale of <0.7 pc, with a brightness temperature of T_b>5x10^7 K. The radio emission is variable by more than a factor of 2.5 over a few days, based on the data taken from VLA monitoring campaigns. The short-term radio variability suggests that the radio emitting region has a size as small as 8x10^{-4} pc, resulting in an even higher brightness temperature of T_b ~10^{12} K. A similar limit on the source size can be obtained if the observed flux variability is not intrinsic and caused by the interstellar scintillation effect. The overall radio spectrum is steep with a time-averaged spectral index alpha=-0.78+/-0.03 between 0.89 GHz and 14 GHz. These observational properties rule out a thermal or star-formation origin of the radio emission, and appear to be consistent with the scenario of episodic jet ejections driven by magnetohydrodynamic process. Simultaneous radio and X-ray monitoring observations down to a cadence of hours are required to test whether the compact and variable radio emission is correlated with the quasi-periodic X-ray eruptions.

Thorold Tronrud, Patricia B. Tissera, Facundo A. Gómez, Robert J. J. Grand, Ruediger Pakmor, Federico Marinacci, Christine M. Simpson

19 pages, 12 figures

We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that this approach allows the detection of accreted stars that are spatially mixed into the disc. Two performance measures are defined - recovery fraction of accreted stars, and the probability that a star with a positive (accreted) classification is a true-positive result, P(TP). As the NNM output is akin to an assigned probability, we are able to determine positivity based on flexible threshold values that can be adjusted easily to refine the selection of presumed-accreted stars. We find that \maga identifies accreted disc stars within simulated galaxies, with high recovery fraction and/or high P(TP). We also find that stars in Gaia-Enceladus-Sausage (GES) mass systems are over 50% recovered by our NNMs in the majority (18/24) of cases. Additionally, nearly every individual source of accreted stars is detected at 10% or more of its peak stellar mass in the disc. We also demonstrate that a conglomerated NNM, trained on the halo and satellite stars from all of the Auriga galaxies provides the most consistent results, and could prove to be an intriguing future approach as our observational capabilities expand.

D. Villaraos, G. Gonzalez-Juarez, R. Lizardo-Castro, A. Herrera-Aguilar, U. Nucamendi

7 pages, 1 figure, 3 tables

In this work we perform a Bayesian statistical fit to estimate the mass-to-distance ratio and the recessional redshift of 10 different black holes hosted at the centre of active galactic nuclei, namely the galaxies NGC 5765b, NGC 6323, UGC 3789, CGCG 074-064, ESO 558-G009, NGC 2960, NGC 6264, NGC 4388, J0437+2456 and NGC 2273. Our general relativistic method makes use of the positions in the sky and frequency shift observations of water megamasers circularly orbiting the central black hole on their accretion disks. This approach also allows us to quantify the gravitational redshift which is not considered in a Newtonian analysis. The gravitational redshift of the megamasers closest to the black hole is found to be within the range 1-8 km/s. The order of the fitted black hole masses corresponds to supermassive black holes and lies on the range $10^6 - 10^7$ M_{sun}

To understand the formation of quiescent solar prominences, the origin of their magnetic field structures, i.e., magnetic flux ropes (MFRs), must be revealed. We use three-dimensional magnetofriction simulations in a spherical subdomain to investigate the role of typical supergranular motions in the long-term formation of a prominence magnetic field. Time-dependant horizontal supergranular motions with and without the effect of Coriolis force are simulated on the solar surface via Voronoi tessellation. The vortical motions by the Coriolis effect at boundaries of supergranules inject magnetic helicity into the corona. The helicity is transferred and accumulated along the polarity inversion line (PIL) as strongly sheared magnetic field via helicity condensation. The diverging motions of supergranules converge opposite magnetic polarities at the PIL and drive the magnetic reconnection between footpoints of the sheared magnetic arcades to form an MFR. The magnetic network, negative-helicity MFR in the northern hemisphere, and fragmented-to-continuous formation process of magnetic dip regions are in agreement with observations. Although diverging supergranulations, differential rotation, and meridional flows are included, the simulation without the Coriolis effect can not produce an MFR or sheared arcades to host a prominence. Therefore Coriolis force is a key factor for helicity injection and the formation of magnetic structures of quiescent solar prominences.

Li Yingjie, Xu Ye, Xu Jin-Long, Liu Dejian, Li Jingjing, Lin Zehao, Jiang Peng, Bian Shuaibo, Hao Chaojie, Chen Xiuhui

20 pages, 6 figures

Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) 19-beam tracking observational mode, high sensitivity and high-velocity resolution HI spectral lines have been observed toward the high-mass star-forming region G176.51+00.20. This is a pilot study of searching for HI narrow-line self-absorption (HINSA) toward high-mass star-forming regions where bipolar molecular outflows have been detected. This work is confined to the central seven beams of FAST. Two HINSA components are detected in all seven beams, which correspond to a strong CO emission region (SCER; with a velocity of $\sim$ $-$18 km s$^{-1}$) and a weak CO emission region (WCER; with a velocity of $\sim$ $-$3 km s$^{-1}$). The SCER detected in Beam 3 is probably more suitably classified as a WCER. In the SCER, the HINSA is probably associated with the molecular material traced by the CO. The fractional abundance of HINSA ranges from $\sim 1.1 \times 10^{-3}$ to $\sim 2.6 \times 10^{-2}$. Moreover, the abundance of HINSA in Beam 1 is lower than that in the surrounding beams (i.e., Beams 2 and 4--7). This possible ring could be caused by ionization of HI or relatively rapid conversion from HI to H$_2$ in the higher-density inner region. In the WCER (including Beam 3 in the SCER), the HINSA is probably not associated with CO clouds, but with CO-dark or CO-faint gas.

Li Yingjie, Xu Ye, Xu Jin-Long, Liu Dejian, Li Jingjing, Lin Zehao, Jiang Peng, Bian Shuaibo, Hao Chaojie, Chen Xiuhui

21 pages, 6 figures, 4 tables, accepted to ApJ

We observed the high-mass star-forming region G176.51+00.20 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) with the 19-beam tracking observational mode. This is a pilot work of searching for neutral stellar winds traced by atomic hydrogen (i.e., HI winds) using the high sensitivity HI line toward high-mass star-forming regions where bipolar molecular outflows have been detected with high sensitivity by Liu et al. HI wind was detected in this work only in Beam 1. We find here that, similar to low-mass star formation, no matter how large the inclination is, the HI wind is likely sufficiently strong to drive a molecular outflow. We also find that the abundance of HI in the HI wind is consistent with that of the HI narrow-line self-absorption (HINSA) in the same beam (i.e., Beam 1). This implies that there is probably an internal relationship between HI winds and HINSA. This result also reinforces the assertion that HI winds and detected molecular outflows are associated with each other.

Simon Petrus, Gaël Chauvin, Mickaël Bonnefoy, Pascal Tremblin, Benjamin Charnay, Philippe Delorme, Gabriel-Dominique Marleau, Amelia Bayo, Elena Manjavacas, Anne-Marie Lagrange, Paul Mollière, Paulina Palma-Bifani, Beth Biller James-S. Jenkins

9 pages, 4 figures, 4 tables

We present simultaneous 0.65-2.5 microns medium resolution (3300 < R < 8100) VLT/X-Shooter spectra of the young low-mass (19+/-5MJup) L-T transition object VHS 1256-1257 b, a known spectroscopic analogue of HR8799d. The companion is a prime target for the JWST Early Release Science (ERS) and one of the highest-amplitude variable brown-dwarf known to date. We compare the spectrum to the custom grids of cloudless ATMO models exploring different atmospheric composition with the Bayesian inference tool ForMoSA. We also re-analyze low-resolution HST/WFC3 1.10-1.67 microns spectra at minimum and maximum variability to contextualize the X-Shooter data interpretation. The models reproduce the slope and most molecular absorption from 1.10 to 2.48 microns self-consistently but fail to provide a radius consistent with evolutionary model predictions. They do not reproduce consistently the optical spectrum and the depth of the K I doublets in the J-band. We derive Teff = 1380+/-54 K, log(g) = 3.97+/-0.48 dex, [M/H] = 0.21+/-0.29, and C/O > 0.63. Our inversion of the HST/WFC3 spectra suggests a relative change of 27+6-5 K of the disk-integrated Teff correlated with the near-infrared brightness. Our data anchor the characterization of that object in the near-infrared and could be used jointly to the ERS mid-infrared data to provide the most detailed characterization of an ultracool dwarf to date.

D. V. Lal, N. Lyskova, C. Zhang, T. Venturi, W. R. Forman, C. Jones, E. M. Churazov, R. J. van Weeren, A. Bonafede, N. A. Miller, I. D. Roberts, A. M. Bykov, L. Di Mascolo, M. Brüggen, G. Brunetti

35 pages (includes Appendices), 15 figures (includes 2 figures in Appendix, some Figures are of lower quality), 6 tables (includes 1 table in Appendix), accepted for publication in the Astrophysical Journal

We present high-resolution, high-sensitivity upgraded Giant Metrewave Radio Telescope observations of the Coma cluster (A1656) at 250-500 MHz and 550-850 MHz. At 250-500 MHz, 135 sources have extensions $>$ 0.45 arcmin (with peak-to-local-noise ratio $> 4$). Of these, 24 sources are associated with Coma-member galaxies. In addition, we supplement this sample of 24 galaxies with 20 ram pressure stripped galaxies from Chen et al. (2020, eight are included in the original extended radio source sample) and an additional five are detected and extended. We present radio morphologies, radio spectra, spectral index maps, and equipartition properties for these two samples. In general, we find the equipartition properties lie within a narrow range (e.g., $P_{\rm min}$ = 1-3 dynes cm$^{-2}$). Only NGC 4874, one of the two brightest central Coma cluster galaxies, has a central energy density and pressure about five times higher and a radio source age about 50 % lower than that of the other Coma galaxies. We find a diffuse tail of radio emission trailing the dominant galaxy of the merging NGC 4839 group that coincides with the "slingshot" tail, seen in X-rays. The southwestern radio relic, B1253$+$275, has a large extent $\approx$ 32$^\prime$ $\times$ 10$^\prime$ ($\simeq$ 1.08 $\times$ 0.34 Mpc$^2$). For NGC 4789, whose long radio tails merge into the relic and may be a source of its relativistic seed electrons, and we find a transverse radio spectral gradient, a steepening from southwest to northeast across the width of the radio source. Finally, radio morphologies of the extended and RPS samples suggest that these galaxies are on their first infall into Coma on (predominantly) radial orbits.

Steven R. Ehlert, Riccardo Ferrazzoli, Andrea Marinucci, Herman L. Marshall, Riccardo Middei, Luigi Pacciani, Matteo Perri, Pierre-Olivier Petrucci, Simonetta Puccetti, Thibault Barnouin, Stefano Bianchi, Ioannis Liodakis, Grzegorz Madejski, Fréderic Marin, Alan P. Marscher, Giorgio Matt, Juri Poutanen, Kinwah Wu, Iván Agudo, Lucio A. Antonelli, Matteo Bachetti, Luca Baldini, Wayne H. Baumgartner, Ronaldo Bellazzini, Stephen D. Bongiorno, Raffaella Bonino, Alessandro Brez, Niccolò Bucciantini, Fiamma Capitanio, Simone Castellano, Elisabetta Cavazzuti, Stefano Ciprini, Enrico Costa, Alessandra De Rosa, Ettore Del Monte, Laura Di Gesu, Niccolò Di Lalla, Alessandro Di Marco, Immacolata Donnarumma, Victor Doroshenko, Michal Dovčiak, Teruaki Enoto, Yuri Evangelista, Sergio Fabiani, et al. (54 additional authors not shown)

Accepted for Publication in ApJ, 15 pages

We present measurements of the polarization of X-rays in the $2-8 \thinspace \mathrm{keV}$ band from the nucleus of the radio galaxy Centaurus A (Cen A), using a 100ks observation from the Imaging X-ray Polarimetry Explorer (IXPE). Nearly simultaneous observations of Cen A were also taken with the Swift, NuSTAR, and INTEGRAL observatories. No statistically significant degree of polarization is detected with IXPE. These observations have a minimum detectable polarization at $99 \%$ percent confidence (MDP$_{99}$) of $6.5 \%$ using a weighted, spectral model-independent calculation in the $2-8 \thinspace \mathrm{keV}$ band. The polarization angle $\psi$ is consequently unconstrained. Spectral fitting across three orders of magnitude in X-ray energy ($0.3-400 \thinspace \mathrm{keV}$) demonstrates that the SED of Cen A is well described by a simple power law with moderate intrinsic absorption ($N_H \sim 10^{23} \thinspace \mathrm{cm}^{-2}$) and a Fe K$\alpha$ emission line, although a second unabsorbed power law is required to account for the observed spectrum at energies below $2 \thinspace \mathrm{keV}$. This spectrum suggests that the reprocessing material responsible for this emission line is optically thin and distant from the central black hole. Our upper limits on the X-ray polarization are consistent with the predictions of Compton scattering, although the specific seed photon population responsible for production of the X-rays cannot be identified. The low polarization degree, variability in the core emission, and the relative lack of variability in the Fe K$\alpha$ emission line support a picture where electrons are accelerated in a region of highly disordered magnetic fields surrounding the innermost jet.

Yuzo Ishikawa, Takahiro Morishita, Massimo Stiavelli, Nicha Leethochawalit, Harry Ferguson, Roberto Gilli, Charlotte Mason, Michele Trenti, Tommaso Treu, Colin Norman

21 pages. 10 figures

The distribution and properties of the first galaxies and quasars are critical pieces of the puzzle in understanding galaxy evolution and cosmic reionization. Previous studies have often excluded unresolved sources as potential low redshift interlopers. We combine broadband color and photometric redshift analysis with morphological selections to identify a robust sample of candidates consistent with unresolved point sources at redshift $z\sim8$ using deep Hubble Space Telescope images. We also examine G141 grism spectroscopic data to identify and eliminate dwarf star contaminants. From these analyses, we identify three, bright ($M_{UV}\lesssim-22$ ABmag) dropout point sources at $7.5<z<8.1$. Spectral energy distribution analyses suggest that these sources are either quasars or compact star-forming galaxies. The flux captured by the IRAC 4.5 $\mu$m channel suggests that they have moderate $H\beta$+$[OIII]$ equivalent widths. We calculate the number density of point sources at $z\sim7$-8, and find that a double powerlaw model well describes the point source distribution. We then extend our analysis to estimate the combined point source + galaxy luminosity function and find that the point sources have a non-negligible contribution to the bright-end excess. The fact that the point sources dominate only at $M_{UV}\lesssim-22$ suggests that their contribution to cosmic reionization is likely limited. While spectroscopic follow-up is needed to confirm the nature of these point sources, this work demonstrates that the inclusion of Lyman dropout point sources is necessary for a complete census of the early galaxies at the epoch of cosmic reionization.

V. V. Bobylev, A. T. Bajkova

14 pages, 6 figures, 2 tables

We study the kinematics of a sample of classical Cepheids younger than 120 Myr. For these stars, the estimates of distances taken from Skowron et al., which are based on the period-luminosity relation, and the line-of-sight velocities and the proper motions from the Gaia catalog are available. There are also distance estimates derived from the trigonometric parallaxes contained in the Gaia ERD3 catalog. A method, which relies on comparison of the first-order derivative of the Galactic rotation angular velocity, showed the need to lengthen the distance scales determined by Skowron et al. by about 10%. This conclusion was confirmed by direct comparison to the distances predicted on using the trigonometric parallaxes. With taking into account this result, we obtained new estimates of the Galactic rotation parameters and the parameters of a spiral density wave.

Sayanti Dasgupta, Ahmed Atteya, Pralay Kumar Karmakar

25 pages, 20 figures, 1 table

We analyze the propagatory nucleus-acoustic wave (NAW) modes excitable in the completely degenerate (CD) core and in its nearly degenerate (ND) ambience of the ONe and CO white dwarfs (WDs). It is based on three-component spherical hydrodynamic quantum plasma consisting of tiny non-thermal quantum electrons, classical thermal light nuclear species (LNS), and classical thermal heavy nuclear species (HNS). The inner concentric layer-wise electronic pressures are judiciously modelled. The electronic energy distribution governed by the Fermi-Dirac (FD) thermostatistical distribution law involves both the thermodynamical temperature and chemical potential. Our exploration emphasizes on the transition state between the thermodynamical temperature and the Fermi temperature for the borderline regions of intermediate degeneracy. A normal spherical mode analysis procedurally yields a sextic generalized linear dispersion relation highlighting the plasma multiparametric dependency of the NAW-features. A numerical illustrative platform is constructed to investigate the full NAW propagatory and dispersive behaviours. We demonstrate that the NAW in ONe (CO) WDs exhibits sensible growth characteristics at near the transcritical (supercritical) wave zone. The temperature-sensitivity of the NAW-growth is more (less) prominent in ONe (CO) WDs. It could be hopefully useful to see the internal structure of compact astroobjects from the asteroseismic probe-perspective of collective quantum interaction processes.

Yifei Jin, Lisa J. Kewley, Ralph S. Sutherland

Accepted for Publication in ApJL

We create a photoionization model embedded in the turbulent ISM by using the state-of-the-art Messenger Monte-Carlo MAPPINGS~V code (M$^3$) in conjunction with the CMFGEN stellar atmosphere model. We show that the turbulent ISM causes the inhomogeneity of electron temperature and density within the nebula. The fluctuation in the turbulent ISM creates complex ionization structures seen in nearby nebulae. The inhomogeneous density distribution within the nebula creates a significant scatter on the spatially-resolved standard optical diagnostic diagrams, which cannot be represented by the spherical constant density photoionization model. We analyze the dependence of different optical emission lines on the complexity of nebular geometry, finding that the emission-lines residing on the nebular boundary are highly sensitive to the complexity of nebular geometry, while the emission-lines produced throughout the nebula are sensitive to the density distribution of the ISM within the nebula. Our fractal photoionization model demonstrates that a complex nebular geometry is required for accurate modeling of HII regions and emission-line galaxies, especially for the high-redshift galaxies, where the ISM is highly turbulent based on the increasing observational evidence.

Nhung, P.T., Hoai, D.T., Darriulat, P., Diep, P.N., Ngoc, N.T.B., Thai, T.T., Tuan-Anh, P

8 pages, 8 figures

We use archival ALMA observations of the CO(2-1) and SiO(5-4) molecular line emissions of AGB star R Hya to illustrate the relative contributions of rotation, expansion and line broadening to the morpho-kinematics of the circumstellar envelope (CSE) within some ~0.5 arcsec from the centre of the star. We give evidence for rotation and important line broadening to dominate the inner region, within ~100 mas from the centre of the star. The former is about an axis that projects a few degrees west of north and has a projected rotation velocity of a few km/s. The latter occurs within some 50-100 mas from the centre of the star, the line width reaching two to three times its value outside this region. We suggest that it is caused by shocks induced by stellar pulsations and convective cell ejections. We show the importance of properly accounting for the observed line broadening when discussing rotation and evaluating the radial dependence of the rotation velocity.

O. Adriani, E. Berti, P. Betti, G. Bigongiari, L. Bonechi, M. Bongi, S. Bottai, P. Brogi, G. Castellini, C. Checchia, R. D Alessandro, S. Detti, N. Finetti, P. Maestro, P. S. Marrocchesi, N. Mori, M. Olmi, L. Pacini, P. Papini, C. Poggiali, S. Ricciarini, P. Spillantini, O. Starodubtsev, F. Stolzi, A. Tiberio, E. Vannuccini

to be published in JINST

The multi-TeV energy region of the cosmic-ray spectra has been recently explored by direct detection experiments that used calorimetric techniques to measure the energy of the cosmic particles. Interesting spectral features have been observed in both all-electron and nuclei spectra. However, the interpretation of the results is compromised by the disagreements between the data obtained from the various experiments, that are not reconcilable with the quoted experimental uncertainties. Understanding the reason for the discrepancy among the measurements is of fundamental importance in view of the forthcoming high-energy cosmic-ray experiments planned for space, as well as for the correct interpretation of the available results. The purpose of this work is to investigate the possibility that a systematic effect may derive from the non-proportionality of the light response of inorganic crystals, typically used in high-energy calorimetry due to their excellent energy-resolution performance. The main reason for the non-proportionality of the crystals is that scintillation light yield depends on ionisation density. Experimental data obtained with ion beams were used to characterize the light response of various scintillator materials. The obtained luminous efficiencies were used as input of a Monte Carlo simulation to perform a comparative study of the effect of the light-yield non-proportionality on the detection of high-energy electromagnetic and hadronic showers. The result of this study indicates that, if the calorimeter response is calibrated by using the energy deposit of minimum ionizing particles, the measured shower energy might be affected by a significant systematic shift, at the level of few percent, whose sign and magnitude depend specifically on the type of scintillator material used.

Hye-Jin Park, Se-Heon Oh, Jing Wang, Yun Zheng, Hong-Xin Zhang, W. J. G. de Blok

20 pages, 11 figures, 4 tables, accepted for publication in AJ

We present H I gas kinematics and star formation activities of NGC 6822, a dwarf galaxy located in the Local Group at a distance of ~ 490 kpc. We perform profile decomposition of line-of-sight velocity profiles of the H I data cube (42.4" x 12.0" spatial, corresponding to ~ 100 pc; 1.6 km s$^{-1}$ spectral) taken with the Australia Telescope Compact Array (ATCA). For this, we use a new tool, the so-called BAYGAUD which is based on Bayesian analysis techniques, allowing us to decompose a line-of-sight velocity profile into an optimal number of Gaussian components in a quantitative manner. We classify the decomposed H I gas components of NGC 6822 into cool-bulk, warm-bulk, cool-non-bulk and warm-non-bulk motions with respect to their centroid velocities and velocity dispersions. We correlate their gas surface densities with corresponding star formation rate densities derived using both the GALEX far-ultraviolet and WISE 22 $\mu$m data to examine the resolved Kennicutt-Schmidt (K-S) law for NGC 6822. Of the decomposed H I gas components, the cool-bulk component is likely to better follow the linear extension of the K-S law for molecular hydrogen (H$_2$) at low gas surface densities where H I is not saturated.

Anwesha Maharana, Alexey Isavnin, Camilla Scolini, Nicolas Wijsen, Luciano Rodriguez, Marilena Mierla, Jasmina Magdalenic, Stefaan Poedts

22 pages, 14 figures

The Flux Rope in 3D (FRi3D, Isavnin, 2016), a coronal mass ejection (CME) model with global three-dimensional (3D) geometry, has been implemented in the space weather forecasting tool EUHFORIA (Pomoell and Poedts, 2018). By incorporating this advanced flux rope model in EUHFORIA, we aim to improve the modelling of CME flank encounters and, most importantly, the magnetic field predictions at Earth. After using synthetic events to showcase FRi3D's capabilities of modelling CME flanks, we optimize the model to run robust simulations of real events and test its predictive capabilities. We perform observation-based modelling of the halo CME event that erupted on 12 July 2012. The geometrical input parameters are constrained using the forward modelling tool included in FRi3D with additional flux rope geometry flexibilities as compared to the pre-existing models. The magnetic field input parameters are derived using the differential evolution algorithm to fit FRi3D parameters to the in situ data at 1 AU. An observation-based approach to constrain the density of CMEs is adopted, in order to achieve a better estimation of mass corresponding to the FRi3D geometry. The CME is evolved in EUHFORIA's heliospheric domain and a comparison of FRi3D's predictive performance with the previously implemented spheromak CME in EUHFORIA is presented. For this event, FRi3D improves the modelling of the total magnetic field magnitude and Bz at Earth by ~30% and ~70%, respectively. Moreover, we compute the expected geoeffectiveness of the storm at Earth using an empirical Dst model and find that the FRi3D model improves the predictions of minimum Dst by ~20% as compared to the spheromak CME model. Finally, we discuss the limitations of the current implementation of FRi3D in EUHFORIA and propose possible improvements.

Velibor Velović, M. D. Filipović, L. Barnes, R. P. Norris, C. D. Tremblay, G. Heald, L. Rudnick, S. S. Shabala, T. G. Pannuti, H. Andernach, O. Titov, S. G. H. Waddell, B. S. Koribalski, D. Grupe, T. Jarrett, R. Z. E. Alsaberi, E. Carretti, J. D. Collier, S. Einecke, T. J. Galvin, A. Hotan, P. Manojlović, J. Marvil, K. Nandra, T. H. Reiprich, G. Rowell, M. Salvato, M. Whiting

16 pages, 14 figures, to be published in MNRAS journal

We present the discovery of highly-collimated radio jets spanning a total of 355 kpc around the nearby elliptical galaxy NGC 2663, and the possible first detection of recollimation on kiloparsec scales. The small distance to the galaxy (~28.5 Mpc) allows us to resolve portions of the jets to examine their structure. We combine multiwavelength data: radio observations by the Murchison Widefield Array (MWA), the Australian Square Kilometre Array Pathfinder (ASKAP) and the Australia Telescope Compact Array (ATCA), and X-ray data from Chandra, Swift and SRG/eROSITA. We present intensity, rotation measure, polarisation, spectral index and X-ray environment maps. Regions of the southern jet show simultaneous narrowing and brightening, which can be interpreted as a signature of the recollimation of the jet by external, environmental pressure, though it is also consistent with an intermittent Active Galactic Nuclei (AGN) or complex internal jet structure. X-ray data suggest that the environment is extremely poor; if the jet is indeed recollimating, the large recollimation scale (40 kpc) is consistent with a slow jet in a low-density environment.

Chin-Ping Hu, Lupin Chun-Che Lin, Kuo-Chuan Pan, Kwan-Lok Li, Chien-Chang Yen, Albert K. H. Kong, C. Y. Hui

24 pages, 26 figures, 1 table, accepted for publication in ApJ

We analyze the gravitational wave signals with a model-independent time-frequency analysis, which is improved from the Hilbert-Huang transform (HHT) and optimized for characterizing the frequency variability on the time-frequency map. Except for the regular HHT algorithm, i.e., obtaining intrinsic mode functions with ensemble empirical mode decomposition and yielding the instantaneous frequencies, we propose an alternative algorithm that operates the ensemble mean on the time-frequency map. We systematically analyze the known gravitational wave events of the compact binary coalescence observed in LIGO O1 and O2, and in the simulated gravitational wave signals from core-collapse supernovae (CCSNe) with our method. The time-frequency maps of the binary black hole coalescence cases show much better details compared to those wavelet spectra. Moreover, the oscillation in the instantaneous frequency caused by mode-mixing could be reduced with our algorithm. For the CCSNe data, the oscillation from the proto-neutron star and the radiation from the standing accretion shock instability can be precisely determined with the HHT in great detail. More importantly, the initial stage of different modes of oscillations can be clearly separated. These results provide new hints for further establishment of the detecting algorithm, and new probes to investigate the underlying physical mechanisms.

Pablo Rivière-Marichalar, Asunción Fuente, Gisela Esplugues, Valentine Wakelam, Romane le Gal, Clément Baruteau, Álvaro Ribas, Enrique Macías, Roberto Neri, David Navarro-Almaida

11 pages, 14 figures, 4 tables

The sulfur abundance is poorly known in most environments. Yet, deriving the sulfur abundance is key to understanding the evolution of the chemistry from molecular clouds to planetary atmospheres. We present observations of H$_2$S 110-101 at 168.763 GHz toward the Herbig Ae star AB Aur. We aim to study the abundance of sulfuretted species toward AB Aur and to constrain how different species and phases contribute to the sulfur budget. We present new NOrthern Extended Millimeter Array (NOEMA) interferometric observations of the continuum and H$_2$S 110-101 line at 168.763 GHz toward AB Aur. We derived radial and azimuthal profiles and used them to compare the geometrical distribution of different species in the disk. Assuming local thermodynamical equilibrium (LTE), we derived column density and abundance maps for H$_2$S, and we further used Nautilus to produce a more detailed model of the chemical abundances at different heights over the mid-plane at a distance of r=200 au. We have resolved H$_2$S emission in the AB Aur protoplanetary disk. The emission comes from a ring extending from 0.67 (109 au) to 1.69 (275 au). Under simple assumptions, we derived an abundance of (3.1$\pm$0.8)$\times$10$\rm ^{-10}$ with respect to H nuclei, which we compare with Nautilus models to deepen our understanding of the sulfur chemistry in protoplanetary disks. Chemical models indicate that H$-2$S is an important sulfur carrier in the solid and gas phase. We also find an important transition at a height of 12 au, where the sulfur budget moves from being dominated by ice species to being dominated by gas species. Studying sulfuretted species in detail in the different phases of the interstellar medium is key to solving the issue.

In cosmology, the fine-structure constant can affect the whole ionization history. However, the previous works confine themselves to the recombination epoch and give various strong constraints on the fine-structure constant. In this paper, we also take the reionization epoch into consideration and do a consistency test of the fine-structure constant from the whole ionization history. We find that the constraints from these two epochs are not consistent with each other by $4.64\sigma$. The straightforward explanation for this inconsistency is a varying fine-structure constant. Another explanation is that the data we used are not good enough. We prefer the latter explanation and a need for calibration of some important parameters involved in reconstructing the reionization history.

Yuki Okoda, Yoko Oya, Muneaki Imai, Nami Sakai, Yoshimasa Watanabe, Ana López-Sepulcre, Kazuya Saigo, Satoshi Yamamoto

Resolving physical and chemical structures in the vicinity of a protostar is of fundamental importance for elucidating their evolution to a planetary system. In this context, we have conducted 1.2 mm observations toward the low-mass protostellar source B335 at a resolution of 0."03 with ALMA. More than 20 molecular species including HCOOH, NH2 CHO, HNCO, CH3 OH, CH2 DOH, CHD2 OH, and CH3 OD are detected within a few 10 au around the continuum peak. We find a systematic chemical differentiation between oxygen-bearing and nitrogen-bearing organic molecules by using the principal component analysis for the image cube data. The distributions of the nitrogen-bearing molecules are more compact than those of the oxygen-bearing ones except for HCOOH. The temperature distribution of the disk/envelope system is revealed by a multi-line analysis for each of HCOOH, NH2 CHO, CH3 OH, and CH2 DOH. The rotation temperatures at the radius of 0."06 along the envelope direction of CH3OH and CH2DOH are derived to be 150-165 K. On the other hand, those of HCOOH and NH2CHO, which have a smaller distribution, are 75-112 K, and are significantly lower than those for CH3OH and CH2DOH. This means that the outer envelope traced by CH3OH and CH2DOH is heated by additional mechanisms rather than the protostellar heating. We here propose the accretion shock as the heating mechanism. The chemical differentiation and the temperature structure on a few au scale provide us with key information to further understand chemical processes in protostellar sources.

Stefano Marchesi, Xiurui Zhao, Núria Torres-Albà, Marco Ajello, Massimo Gaspari, Andrealuna Pizzetti, Johannes Buchner, Elena Bertola, Andrea Comastri, Anna Feltre, Roberto Gilli, Giorgio Lanzuisi, Gabriele Matzeu, Francesca Pozzi, Francesco Salvestrini, Dhrubojyoti Sengupta, Ross Silver, Francesco Tombesi, Alberto Traina, Cristian Vignali, Luca Zappacosta

24 pages, 13 figures: main results in Figures 4, 7, 8 and 9. Accepted for publication in the Astrophysical Journal

We present the multi-epoch monitoring with NuSTAR and XMM-Newton of NGC 1358, a nearby Seyfert 2 galaxy whose properties made it a promising candidate X-ray changing look AGN, i.e., a source whose column density could transition from its 2017 Compton-thick (CT-, having line-of-sight Hydrogen column density NH,los>10^24 cm^-2) state to a Compton-thin (NH,los<10^24 cm^-2) one. The multi-epoch X-ray monitoring confirmed the presence of significant NH,los variability over time-scales as short as weeks, and allowed us to confirm the "changing look" nature of NGC 1358, which has most recently been observed in a Compton-thin status. Multi-epoch monitoring with NuSTAR and XMM-Newton is demonstrated to be highly effective in simultaneously constraining three otherwise highly degenerate parameters: the torus average column density and covering factor, and the inclination angle between the torus axis and the observer. We find a tentative anti-correlation between column density and luminosity, which can be understood in the framework of Chaotic Cold Accretion clouds driving recursive AGN feedback. The monitoring campaign of NGC 1358 has proven the efficiency of our newly developed method to select candidate NH,los-variable, heavily obscured AGN, which we plan to soon extend to a larger sample to better characterize the properties of the obscuring material surrounding accreting supermassive black holes, as well as constrain AGN feeding models.

Shuai Fu, Zheyi Ding, Yongjie Zhang, Xiaoping Zhang, Cunhui Li, Gang Li, Shuwen Tang, Haiyan Zhang, Yi Xu, Yuming Wang, Jingnan Guo, Lingling Zhao, Yi Wang, Xiangyu Hu, Pengwei Luo, Zhiyu Sun, Yuhong Yu, Lianghai Xie

This manuscript has been accepted for publication in The Astrophysical Journal Letters

Solar energetic particles (SEPs) associated with flares and/or coronal mass ejection (CME)-driven shocks can impose acute radiation hazards to space explorations. To measure energetic particles in near-Mars space, the Mars Energetic Particle Analyzer (MEPA) instrument onboard China's Tianwen-1 (TW-1) mission was designed. Here, we report the first MEPA measurements of the widespread SEP event occurring on 29 November 2020 when TW-1 was in transit to Mars. This event occurred when TW-1 and Earth were magnetically well connected, known as the Hohmann-Parker effect, thus offering a rare opportunity to understand the underlying particle acceleration and transport process. Measurements from TW-1 and near-Earth spacecraft show similar double-power-law spectra and a radial dependence of the SEP peak intensities. Moreover, the decay phases of the time-intensity profiles at different locations clearly show the reservoir effect. We conclude that the double-power-law spectrum is likely generated at the acceleration site, and that a small but finite cross-field diffusion is crucial to understand the formation of the SEP reservoir phenomenon. These results provide insight into particle acceleration and transport associated with CME-driven shocks, which may contribute to the improvement of relevant physical models.

Zhang XueGuang (NNU), Zhao YuanBo

6 pages, 5 figures, accepted to be published in MNRAS

In the Letter, an interesting method is proposed to estimate size of narrow emission lines regions (NLRs) of a Type-2 AGN SDSS J083823.91+490241.1 (=SDSS J0838) at a redshift of 0.101, by comparing spectroscopic properties through the SDSS fiber (MJD=51873) (diameter of 3 arcseconds) and through the eBOSS fiber (MJD=55277) (diameter of 2 arcseconds). After subtractions of pPXF method determined host galaxy contributions, the narrow emission lines of SDSS J0838 in the SDSS spectrum and in the eBOSS spectrum can be well measured by Gaussian functions, leading more than 90\% of [O~{\sc iii}] emissions to be covered by the eBOSS fiber with diameter of 2 arcseconds. Meanwhile, both none broad emission components and none-variabilities of ZTF 3years-long g/r-band light curves can be applied to confirm SDSS J0838 as a Type-2 AGN, indicating few orientation effects on the projected NLRs size in SDSS J0838. Therefore, upper limit about 1arcsecond (2250pc) of the NLRs size can be reasonably accepted in SDSS J0838. Combining with the intrinsic reddening corrected [O~{\sc iii}] line luminosity, the upper limit of NLRs size in SDSS J0838 well lies within the 99.9999\% confidence bands of the R-L empirical relation for NLRs in AGN.

Tomek Kamiński, Mirek Schmidt, Marcin Hajduk, Aleksandra Kiljan, Inna Izviekova, Adam Frankowski

submitted to A&A; comments welcome

Context. Galactic red novae are thought to be produced in stellar mergers between non-compact stars, such as main-sequence stars and cool giants. They are hoped to help in explaining physical processes involved in common envelope evolution and stellar binary collisions. Aims. We investigate the presence of lithium in three best-observed Galactic red nova remnants. Explaining the origin of lithium may point to mixing mechanism present before, during, or after the merger. Methods. The lithium line at 6707.81 A was compared to a feature of [Ca I at 6572.78 A to derive relative abundances in circumstellar gas. Absolute abundances were next calculated assuming the Solar calcium to lithium abundance ratio. Results. Lithium abundances were measured in the merger remnants of V838 Mon with A(Li)=2.3, CK Vul with A(Li)=2.5, and V1309 Sco with A(Li)=1.8. Conclusions. Lithium is overabundant in red novae, suggesting that at least some merger products activate the Cameron-Fowler mechanism whereby convective mixing can reach the deep stellar interior. Whether deep convection and associated diffusion alone or some other processes (e.g. spin down) can be responsible for driving the Cameron-Fowler mechanism in the remnants requires further studies. Early observations of lithium in V838 Mon hint that these mechanisms can activate early, perhaps already in the common envelope phase. These observations should be taken into account in modelling these complex systems.

Marco Barbera, Ugo Lo Cicero, Luisa Sciortino

45 pages, 11 figures Contribution to the Handbook of X-ray and Gamma-ray Astrophysics in the chapter "Detectors for X-ray Astrophysics", edited by Jan-Willem den Herder, Norbert Meidinger, Marco Feroci

Thin filters and gas tight windows are used in Space to protect sensitive X-ray detectors from out-of-band electromagnetic radiation, low-energy particles, and molecular contamination. Though very thin and made of light materials, filters are not fully transparent to X-rays. For this reason, they ultimately define the detector quantum efficiency at low energies. In this chapter, we initially provide a brief overview of filter materials and specific designs adopted on space experiments with main focus on detectors operating at the focal plane of grazing incidence X-ray telescopes. We then provide a series of inputs driving the design and development of filters for high-energy astrophysics space missions. We begin with the identification of the main functional goals and requirements driving the preliminary design, and identify modeling tools and experimental characterization techniques needed to prove the technology and consolidate the design. Finally, we describe the calibration activities required to derive the filter response with high accuracy.We conclude with some hints on materials and technologies presently under investigation for future X-ray missions.

Laurent Legrand, Pratika Dayal, Anne Hutter, Stefan Gottlöber, Gustavo Yepes, Maxime Trebitsch

17 pages, 15 figures, submitted to MNRAS, comments welcome !

Using the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework, we explore the impact of environmental density and radiative feedback on the assembly of galaxies and their host halos during the Epoch of Reionization. The ASTRAEUS framework allows us to study the evolution of galaxies with masses ($\rm 10^{8.2}M_\odot < M_{\rm h} < 10^{13}M_\odot$) in wide variety of environment ($-0.5 < {\rm log}(1+\delta) < 1.3$ averaged over $(2~{\rm cMpc})^3$). We find that : (i) there exists a mass- and redshift- dependent "characteristic" environment (${\rm log} (1+\delta_a(M_{\rm h}, z)) = 0.021\times M_{\rm h}^{0.16} + 0.07 z -1.12$, up to $z\sim 10$) at which galaxies are most efficient at accreting dark matter, e.g at a rate of $0.2\%$ of their mass every Myr at $z=5$; (ii) the number of minor and major mergers and their contributions to the dark matter assembly increases with halo mass at all redshifts and is mostly independent of the environment; (iii) at $z=5$ minor mergers contribute slightly more (by up to $\sim 10\%$) to the dark matter assembly while for the stellar assembly, major mergers dominate the contribution from minor mergers for $M_{\rm h}\lesssim 10^{11.5}M_\odot$ galaxies; (iv) radiative feedback quenches star formation more in low-mass galaxies ($M_{\rm h} \lesssim 10^{9.5}M_\odot$) in over-dense environments (${\rm log}(1+\delta) > 0.5$); dominated by their major branch, this yields star formation histories biased towards older ages with a slower redshift evolution.

Quasars can be used to measure baryon acoustic oscillations at high redshift, which are considered as direct tracers of the most distant large-scale structures in the Universe. It is fundamental to select quasars from observations before implementing the above research. This work focuses on creating a catalog of quasar candidates based on photometric data to provide primary priors for further object classification with spectroscopic data in the future, such as The Dark Energy Spectroscopic Instrument (DESI) Survey. We adopt a machine learning algorithm (Random Forest, RF) for quasar identification. The training set includes $651,073$ positives and $1,227,172$ negatives, in which the photometric information are from DESI Legacy Imaging Surveys (DESI-LIS) \& Wide-field Infrared Survey Explore (WISE), and the labels are from a database of spectroscopically confirmed quasars based on Sloan Digital Sky Survey (SDSS) and the Set of Identifications \& Measurements and Bibliography for Astronomical Data (SIMBAD). The trained RF model is applied to point-like sources in DESI-LIS Data Release 9. To quantify the classifier's performance, we also inject a testing set into the to-be-applied data. Eventually, we obtained $1,953,932$ Grade-A quasar candidates and $22,486, 884$ Grade-B quasar candidates out of $425,540,269$ sources ($\sim 5.7\%$). The catalog covers $\sim 99\%$ of quasars in the to-be-applied data by evaluating the completeness of the classification on the testing set. The statistical properties of the candidates agree with that given by the method of color-cut selection. Our catalog can intensely decrease the workload for confirming quasars with the upcoming DESI data by eliminating enormous non-quasars but remaining high completeness. All data in this paper is publicly available online.

Vincenzo Sapienza, Marco Miceli, Aya Bamba, Satoru Katsuda, Tsutomu Nagayoshi, Yukikatsu Terada, Fabrizio Bocchino, Salvatore Orlando, Giovanni Peres

16 pages, 7 figures, 5 tables, Accepted for publication in ApJ

Synchrotron X-ray emission in young supernova remnants (SNRs) is a powerful diagnostic tool to study the population of high energy electrons accelerated at the shock front and the acceleration process. We performed a spatially resolved spectral analysis of NuSTAR and XMM-Newton observations of the young Kepler's SNR, aiming to study in detail its non-thermal emission in hard X-rays. We selected a set of regions all around the rim of the shell and extracted the corresponding spectra. The spectra were analyzed by adopting a model of synchrotron radiation in the loss-limited regime, to constrain the dependence of the cutoff energy of the synchrotron radiation on the shock velocity. We identify two different regimes of particle acceleration, characterized by different Bohm factors. In the north, where the shock interacts with a dense circumstellar medium (CSM), we found a more efficient acceleration than in the south, where the shock velocity is higher and there are no signs of shock interaction with dense CSM. Our results suggest an enhanced efficiency of the acceleration process in regions where the shock-CSM interaction generates an amplified and turbulent magnetic field. By combining hard X-ray spectra with radio and $\gamma-$ray observations of Kepler's SNR, we modelled the spectral energy distribution. In the light of our results we propose that the observed $\gamma-$ray emission is mainly hadronic, and originates in the northern part of the shell.

N.-E. Nèmec, A. I. Shapiro, E. Işık, K. Sowmya, S. K. Solanki, N. A. Krivova, R. H. Cameron, L. Gizon

10 pages, 6 Figures, accepted for publication in ApJL

Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres.

Maciek Wielgus, Nicola Marchili, Ivan Marti-Vidal, Garrett K. Keating, Venkatessh Ramakrishnan, Paul Tiede, Ed Fomalont, Sara Issaoun, Joey Neilsen, Michael A. Nowak, Lindy Blackburn, Charles F. Gammie, Ciriaco Goddi, Daryl Haggard, Daeyoung Lee, Monika Moscibrodzka, Alexandra J. Tetarenko, Geoffrey C. Bower, Chi-Kwan Chan, Koushik Chatterjee, Paul M. Chesler, Jason Dexter, Sheperd S. Doeleman, Boris Georgiev, Mark Gurwell, Michael D. Johnson, Daniel P. Marrone, Alejandro Mus, Dimitrios Psaltis, Bart Ripperda, Gunther Witzel + the Event Horizon Telescope Collaboration

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

Salvador Bará, Carmen Bao-Varela, Raul C. Lima

26 pages, 4 figures. Preprint

Low Earth orbit radiometers allow monitoring nighttime anthropogenic light emissions in wide areas of the planet. In this work we describe a simple model for assessing significant outdoor lighting changes at the municipality level using on-orbit measurements complemented with ground-truth information. We apply it to evaluate the transformation effected in the municipality of Ribeira (42{\deg} 33 23 N, 8{\deg} 59 32 W) in Galicia, which in 2015 reduced the amount of installed lumen in its publicly-owned outdoor lighting system from 93.2 to 28.7 Mlm. This significant cutback, with the help of additional controls, allowed to reduce from 0.768 to 0.208 Mlm/km2 the lumen emission density averaged across the territory. In combination with the VIIRS-DNB annual composite readings these data allow to estimate that the relative weight of the emissions of the public streetlight system with respect to the total emissions of light in the municipality changed from an initial value of 74.86% to 44.68% after the transformation. The effects of the sources spectral shift and the photon calibration factor on the radiance reported by the VIIRS-DNB are also evaluated.

Maria I. Carnerero, Claudia M. Raiteri, Lorenzo Rimoldini, Deborah Busonero, Enrico Licata, Nami Mowlavi, Isabelle Lecoeur-Taïbi, Marc Audard, Berry Holl, Panagiotis Gavras, Krzysztof Nienartowicz, Grégory Jevardat de Fombelle, Ruth Carballo, Gisella Clementini, Ludovic Delchambre, Sergei Klioner, Mario G. Lattanzi, Laurent Eyer

19 pages, 31 figures, 2 table. This paper is part of Gaia Data Release 3 (DR3). In press for A&A

One of the novelties of the Gaia-DR3 with respect to the previous data releases is the publication of the multiband light curves of about 1 million AGN. The goal of this work was the creation of a catalogue of variable AGN, whose selection was based on Gaia data only. We first present the implementation of the methods to estimate the variability parameters into a specific object study module for AGN. Then we describe the selection procedure that led to the definition of the high-purity variable AGN sample and analyse the properties of the selected sources. We started from a sample of millions of sources, which were identified as AGN candidates by 11 different classifiers based on variability processing. Because the focus was on the variability properties, we first defined some pre-requisites in terms of number of data points and mandatory variability parameters. Then a series of filters was applied using only Gaia data and the Gaia Celestial Reference Frame 3 (Gaia-CRF3) sample as a reference.The resulting Gaia AGN variable sample, named GLEAN, contains about 872000 objects, more than 21000 of which are new identifications. We checked the presence of contaminants by cross-matching the selected sources with a variety of galaxies and stellar catalogues. The completeness of GLEAN with respect to the variable AGN in the last Sloan Digital Sky Survey quasar catalogue is about 47%, while that based on the variable AGN of the Gaia-CRF3 sample is around 51%. From both a comparison with other AGN catalogues and an investigation of possible contaminants, we conclude that purity can be expected to be above 95%. Multiwavelength properties of these sources are investigated. In particular, we estimate that about 4% of them are radio-loud. We finally explore the possibility to evaluate the time lags between the flux variations of the multiple images of strongly lensed quasars, and show one case.

Christophe Pixius, Safak Celik, Matthias Bartelmann

26 pages, 4 figures

We present recent improvements in the perturbative treatment of particle interactions in Kinetic Field Theory (KFT) for inertial Zel'dovich trajectories. KFT has been developed for the systematic analytical calculation of non-linear cosmic structure formation on the basis of microscopic phase-space dynamics. We improve upon the existing treatment of the interaction operator by deriving a more rigorous treatment of phase-space trajectories of particles in an expanding universe. We then show how these results can be applied to KFT perturbation theory by calculating corrections to the late-time dark matter power spectrum at second order in the interaction operator. We find that the modified treatment of interactions w.r.t. inertial Zel'dovich trajectories improves the agreement of KFT with simulation results on intermediate scales compared to earlier results. Additionally, we illustrate that including particle interactions up to second order leads to a systematic improvement of the non-linear power spectrum compared to the first-order result.

Giorgos Leloudas, Mattia Bulla, Aleksandar Cikota, Lixin Dai, Lars L. Thomsen, Justyn R. Maund, Panos Charalampopoulos, Nathaniel Roth, Iair Arcavi, Katie Auchettl, Daniele B. Malesani, Matt Nicholl, Enrico Ramirez-Ruiz

Author's version of paper to appear in Nature Astronomy. In the journal version the detailed discussion on the ISP determination will be moved from the Methods section to a Supplementary Information section. 58 pages in double spacing format, including 5 Figures, 10 Extended Data Figures and 2 Tables

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events (TDEs) have now been identified in the optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three TDEs. The continuum polarization is independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with optical photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code we generate predictions for the degree of polarization as a function of disk mass and viewing angle, and we show that the predicted levels are compatible with the observations, for extended reprocessing envelopes of $\sim$1000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for TDE models, and opens an important new line of exploration in the study of TDEs.

Gábor Galgóczi, Jean-Paul Breuer, Valentina Fioretti, Jakub Zlámal, Norbert Werner, Vojtěch Čalkovský, Nathalie Boudin, Ivo Ferreira, Matteo Guainazzi, Andreas von Kienlin, Simone Lotti, Teresa Mineo, Silvano Molendi, Emanuele Perinati

SPIE conference proceeding

X-ray telescopes opened up a new window into the high-energy universe. However, the last generation of these telescopes encountered an unexpected problem: their optics focused not only X-rays but low-energy (so called soft) protons as well. These protons are very hard to model and can not be distinguished from X-rays. For example, 40\% of XMM-Newton observations is significantly contaminated by soft proton induced background flares. In order to minimize the background from such low-energy protons the Advanced Telescope for High ENergy Astrophysics (ATHENA) satellite introduced a novel concept, the so called Charged Particle Diverter (CPD). It is an array of magnets in a Hallbach design, which deflects protons below 76 keV before they would hit the Wide Field Imager (WFI) detector. In this work, we investigate the effect of scattering of the deflected protons with the CPD walls and the inner surfaces of the WFI detector assembly. Such scattered protons can loose energy, change direction and still hit the WFI. In order to adopt the most realistic instrument model, we imported the CAD model of both the CPD and the WFI focal plane assembly. Soft protons corresponding to $\approx$2.5 hours of exposure to the L1 solar wind are simulated in this work. The inhomogeneous magnetic field of the CPD is included in the simulation. We present a preliminary estimate of the WFI residual background induced by soft proton secondary scattering, in the case of the optical blocking filter present in the field of view. A first investigation of the volumes responsible for scattering the protons back into the field of view is reported.

G. Del Zanna, B. Mondal, Y.K.Rao, N. P. S. Mithun, S. V. Vadawale, K. K. Reeves, H. E. Mason, A. Sarkar, P. Janardhan, A. Bhardwaj

accepted for publication

We have reviewed the first year of observations of the Solar X-ray Monitor (XSM) onboard Chandrayaan-2, and the available multi-wavelength observations to complement the XSM data, focusing on Solar Dynamics Observatory AIA and Hinode XRT, EIS observations. XSM has provided disk-integrated solar spectra in the 1--15 keV energy range, observing a large number of microflares. We present an analysis of multi-wavelength observations of AR 12759 during its disk crossing. We use a new radiometric calibration of EIS to find that the quiescent AR core emission during its disk crossing has a distribution of temperatures and chemical abundances that does not change significantly over time. An analysis of the XSM spectra confirms the EIS results, and shows that the low First Ionization Potential (FIP) elements are enhanced, compared to their photospheric values. The frequent microflares produced by the AR did not affect the abundances of the quiescent AR core. We also present an analysis of one of the flares it produced, SOL2020-04-09T09:32. The XSM analysis indicates isothermal temperatures reaching 6 MK. The lack of very high-T emission is confirmed by AIA. We find excellent agreement between the observed XSM spectrum and the one predicted using an AIA DEM analysis. In contrast, the XRT Al-Poly / Be-thin filter ratio gives lower temperatures for the quiescent and flaring phases. We show that this is due to the sensitivity of this ratio to low temperatures, as the XRT filter ratios predicted with a DEM analysis based on EIS and AIA gives values in good agreement with the observed ones.

P. Christian Schneider, H. Moritz Günther, Sabina Ustamujic

Invited chapter for the "Handbook of X-ray and Gamma-ray Astrophysics" (Eds. C. Bambi and A. Santangelo, Springer Nature, 2022), accepted (34 pages, 11 figures)

Low-mass pre-main sequence (PMS) stars are strong X-ray sources, because they possess hot corona like their older main-sequence counterparts. Unique to young stars, however, are X-rays from accretion and outflows, and both processes are of pivotal importance for star and planet formation. We describe how X-ray data provide important insight into the physics of accretion and outflows. First, mass accreted from a circumstellar disk onto the stellar surface reaches velocities up to a few hundred km/s, fast enough to generate soft X-rays in the post-shock region of the accretion shock. X-ray observations together with laboratory experiments and numerical simulations show that the accretion geometry is complex in young stars. Specifically, the center of the accretion column is likely surrounded by material shielding the inner flow from view but itself also hot enough to emit X-rays. Second, X-rays are observed in two locations of protostellar jets: an inner stationary emission component probably related to outflow collimation and outer components, which evolve withing years and are likely related to working surfaces where the shock travels through the jet. Jet-powered X-rays appear to trace the fastest jet component and provide novel information on jet launching in young stars. We conclude that X-ray data will continue to be highly important for understanding star and planet formation, because they directly probe the origin of many emission features studied in other wavelength regimes. In addition, future X-ray missions will improve sensitivity and spectral resolution to probe key model parameters (e.g. velocities) in large samples of PMS stars.

Francesco Iacovelli, Michele Mancarella, Stefano Foffa, Michele Maggiore

18 pages, 6 figures, $\texttt{GWFAST}$ available at $\href{ this https URL }{\rm this\ link}$, $\texttt{WF4Py}$ available at $\href{ this https URL }{\rm this\ link}$

This is a supplement to "Forecasting the detection capabilities of third-generation gravitational-wave detectors using $\texttt{GWFAST}$", where the detection capabilities of the second and third generation of ground-based gravitational-wave detectors are studied. The software used to produce these results is $\texttt{GWFAST}$ (https://github.com/CosmoStatGW/gwfast), a Fisher information $\texttt{Python}$ code that allows us to easily and efficiently estimate signal-to-noise ratios and parameter measurement errors for large catalogs of resolved sources observed by networks of gravitational-wave detectors. In particular, $\texttt{GWFAST}$ includes the effects of the Earth's motion during the evolution of the signal, supports parallel computation, and relies on automatic differentiation rather than on finite differences techniques, which allows the computation of derivatives with accuracy close to machine precision. We also release the library $\texttt{WF4Py}$ (https://github.com/CosmoStatGW/WF4Py) implementing state-of-the-art gravitational-wave waveforms in $\texttt{Python}$. In this supplement we provide a documentation of $\texttt{GWFAST}$ and $\texttt{WF4Py}$ with practical examples and tests of performance and reliability.

We investigated the incompressible and compressible magnetohydrodynamic (MHD) energy cascade rates in the solar wind at different heliocentric distances. We used in situ magnetic field and plasma observations provided by the Parker Solar Probe (PSP) mission and exact relations in fully developed turbulence. To estimate the compressible cascade rate, we applied two recent exact relations for compressible isothermal and polytropic MHD turbulence, respectively. Our observational results show a clear increase of the compressible and incompressible cascade rates as we get closer to the Sun. Moreover, we obtained an increase in both isothermal and polytropic cascade rates with respect to the incompressible case as compressibility increases in the plasma. Further discussion about the relation between the compressibility and the heliocentric distance is carried out. Finally, we compared both exact relations as compressibility increases in the solar wind and although we note a slightly trend to observe larger cascades using a polytropic closure, we obtained essentially the same cascade rate in the range of compressibility observed.

Abhijit Kayal, Veeresh Singh, C.H. Ishwara Chandra, Yogesh Wadadekar, Sushant Dutta

18 pages, 10 figures, 4 tables, accepted for publication in the Journal of Astrophysics and Astronomy (to appear in the special issue of Indian participation in the SKA)

Radio observations being insensitive to the dust-obscuration, have been exploited to unveil the population of Active Galactic Nuclei (AGN) residing in galaxies with large dust content. In this paper, we investigate the radio characteristics of 321 dust-obscured galaxies (DOGs; S$_{\rm 24 {\mu}m}$/S$_{\rm r band}$ $\geq$1000) by using mainly deep band-3 (250$-$550 MHz) observations from the upgraded Giant Metrewave Radio Telescope (uGMRT) and 1.5 GHz Jansky Very Large Array (JVLA) observations. We find that, for our sample of DOGs, deep (median noise-rms $=$ 30 $\mu$Jy beam$^{-1}$) 400 MHz band-3 uGMRT observations yield the highest detection rate (28 per cent) among those obtained with the JVLA, and the LOw Frequency ARray (LOFAR) radio observations and the ${\it XMM-N}$ X-ray observations. The radio characteristics of our sample sources, ${\it i.e.,}$ linear extent ($<$40 kpc at $z$ $<$1.2), bimodal spectral index (${\alpha}_{\rm 400 MHz}^{\rm 1.5 GHz}$) distribution and the radio luminosities (L$_{\rm 1.5 GHz}$ $>$5.0 $\times$ 10$^{23}$ W Hz$^{-1}$), suggest them to be mainly consist of Compact$-$Steep$-$Spectrum (CSS) or Peaked-Spectrum (PS) sources representing an early phase of the AGN-jet activity in dust-obscured environments. With stacking we find the existence of faint radio emission (S$_{\rm 400 MHz}$ = 72.9 $\mu$Jy beam$^{-1}$ and S$_{\rm 1.5 GHz}$ = 29 $\mu$Jy beam$^{-1}$ with signal-to-noise ratio $\sim$ 20) in otherwise radio-undetected DOGs. Our study revealing the faint emission at a few tens of $\mu$Jy level in high$-z$ DOGs can be used as a test-bed for the deeper radio continuum surveys planned with the Square-Kilometer Array (SKA) and its pathfinders.

Yue Hu, Siyao Xu, James M. Stone, Alex Lazarian

17 pages, 15 figures, submitted to ApJ

Magnetic fields on the order of 100 $\mu$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven by the strongest preshock density contrast and cascades to smaller scales. The resulting turbulence amplifies the postshock magnetic field. The magnetic fields' time evolution agrees with the prediction of the nonlinear dynamo theory in Xu & Lazarian (2016). When the initial weak magnetic field is perpendicular to shock normal, the total amplification of the field's strength achieves a factor of $\approx200$, which is twice larger than the one in the parallel shock case. However, the strongest magnetic field has a low volume filling factor and is limited by the turbulent energy due to the reconnection diffusion taking place in a turbulent and magnetized fluid. The magnetic field strength averaged along the shock surface is reduced by a factor $\gtrsim10$. We decompose the turbulent velocity and magnetic field into solenoidal and compressive modes. The solenoidal mode is dominant and follows the Kolmogorov scaling, even though the preshock density distribution has a shallow spectrum. When the preshock density distribution has a Kolmogorov spectrum, the fraction of the compressive component increases. We find that the perpendicular shock exhibits a smaller turbulent Alfv\'en Mach number in the vicinity of the shock front than the parallel shock.

Dingshan Deng, Yang Sun, Tianding Wang, Yuxi Wang, Biwei Jiang

Accepted for publication in the ApJ (July 13, 2022). 17 pages, 9 figures

The infrared excess from OB stars are commonly considered as contributions from ionized stellar wind or circumstellar dust. With the newly published LAMOST-OB catalog and GOSSS data, this work steps further on understanding the infrared excess of OB stars. Based on a forward modeling approach comparing the spectral slope of observational Spectral Energy Distributions (SED) and photospheric models, 1147 stars are found to have infrared excess from 7818 stars with good-quality photometric data. After removing the objects in the sightline of dark clouds, 532 ($\sim7\%$) B-type stars and 118 ($\sim23\%$) O-type stars are identified to be true OB stars with circumstellar infrared excess emission. The ionized stellar wind model and the circumstellar dust model are adopted to explain the infrared excess, and Bayes Factors are computed to quantitatively compare the two. It is shown that the infrared excess can be accounted for by the stellar wind for about 65\% cases in which 33\% by free-free emission and 32\% by synchrotron radiation. Other 30\% sources could have and 4\% should have a dust component or other mechanisms to explain the sharply increase flux at $\lambda > 10\mu$m. The parameters of dust model indicate a large-scale circumstellar halo structure which implies the origin of the dust from the birthplace of the OB stars. A statistical study suggests that the proportion with infrared excess in OB stars increases with stellar effective temperature and luminosity, and that there is no systematic change of the mechanism for infrared emission with stellar parameters.

M. Rybak, J. A. Hodge, T. R. Greve, D. Riechers, I. Lamperti, J. van Marrewijk, F. Walter, J. Wagg, P. P. van der Werf

Submitted to A&A

Dusty star-forming galaxies (DSFGs) at redshift z$\geq$1 are among the most vigorously star-forming galaxies in the Universe. However, their dense ($\geq$10$^5$ cm$^{-3}$ ) gas phase - typically traced by HCN(1-0) - remains almost entirely unexplored: only two DSFGs have been detected in HCN(1-0) to date. We present results of a JVLA survey of the J=1-0 transition of HCN, HCO+, and HNC(1-0) in six strongly lensed DSFGs at z = 2.5 - 3.3, effectively doubling the number of DSFGs with deep observations of these lines. We detect HCN(1-0) emission in one source (J1202+5354, 4.4$\sigma$), with a tentative HCO+ (1-0) detection in another (J1609+6045, 3.3$\sigma$). Spectral stacking yields strict upper limits on the HCN/FIR ($\leq$3.6$\times$10$^{-4}$) and HCN/CO(1-0) ratios ($\leq$0.045). The inferred HCN/FIR ratios (a proxy for the star-formation efficiency) are consistent with those in z$\sim$0 FIR-luminous starbursts. However, the HCN/CO ratios - a proxy for the dense-gas fraction - are a factor of a few lower than suggested by the two previous DSFG detections. Our results imply that most DSFGs have low dense-gas fractions. A comparison with Krumholz & Thompson (2007) models of star-forming galaxies indicates that the bulk of gas in DSFGs is at lower densities ($\approx$10$^2$ cm$^{-3}$ ), similar to "normal" star-forming galaxies, rather than ultraluminous starbursts.

Thomas Navarro, Timothy M. Merlis, Nicolas B. Cowan, Natalya Gomez

Temperate terrestrial planets orbiting low-mass stars are subject to strong tidal forces. The effects of gravitational tides on the solid planet and that of atmospheric thermal tides have been studied, but the direct impact of gravitational tides on the atmosphere itself has so far been ignored. We first develop a simplified analytic theory of tides acting on the atmosphere of a planet. We then implement gravitational tides into a general circulation model of a static-ocean planet in a short-period orbit around a low-mass star -- the results agree with our analytic theory. Because atmospheric tides and solid-body tides share a scaling with the semi-major axis, we show that there is a maximum amplitude of the atmospheric tide that a terrestrial planet can experience while still having a solid surface; Proxima Centauri b is the poster child for a planet that could be geophysically Earth-like but with atmospheric tides more than 500$\times$ stronger than Earth's. In this most extreme scenario, we show that atmospheric tides significantly impact the planet's meteorology -- but not its climate. Two possible modest climate impacts are enhanced longitudinal heat transport and cooling of the lowest atmospheric layers. The strong radiative forcing of such planets dominates over gravitational tides, unlike moons of cold giant planets, such as Titan. We speculate that atmospheric tides could be climatologically important on planets where the altitude of maximal tidal forcing coincides with the altitude of cloud formation and that the effect could be detectable for non-Earth-like planets subject to even greater tides.

Bahar Bidaran, Francesco La Barbera, Anna Pasquali, Reynier Peletier, Glenn van de Ven, Eva K. Grebel, Jesus Falcón-Barroso, Agnieszka Sybilska, Dimitri A. Gadotti, Lodovico Coccato

17 pages, 10 figures, Accepted for publication in MNRAS

Using MUSE spectra, we investigate how pre-processing and accretion onto a galaxy cluster affect the integrated stellar population properties of dwarf early-type galaxies (dEs). We analyze a sample of nine dEs with stellar masses of $\rm \sim 10^9 \, M_\odot$, which were accreted ($\sim$ 2-3 Gyr ago) onto the Virgo cluster as members of a massive galaxy group. We derive their stellar population properties, namely age, metallicity ([M/H]), and the abundance ratio of $\alpha$ elements ([$\alpha$/Fe]), by fitting observed spectral indices with a robust, iterative procedure, and infer their star formation history (SFH) by means of full spectral fitting. We find that these nine dEs are more metal-poor (at the 2-3$\sigma$ level) and significantly more $\alpha$-enhanced than dEs in the Virgo and Coma clusters with similar stellar mass, cluster-centric distance, and infall time. Moreover, for six dEs, we find evidence for a recent episode of star formation during or right after the time of accretion onto Virgo. We interpret the high [$\alpha$/Fe] of our sample of dEs as the result of the previous exposure of these galaxies to an environment hostile to star formation, and/or the putative short burst of star formation they underwent after infall into Virgo. Our results suggest that the stellar population properties of low-mass galaxies may be the result of the combined effect of pre-processing in galaxy groups and environmental processes (such as ram-pressure triggering star formation) acting during the early phases of accretion onto a cluster.

Harrison F. Agrusa, Fabio Ferrari, Yun Zhang, Derek C. Richardson, Patrick Michel

21 pages, 19 figures, published in PSJ

Previous efforts have modeled the Didymos system as two irregularly shaped rigid bodies, although it is likely that one or both components are in fact rubble piles. Here, we relax the rigid-body assumption to quantify how this affects the spin and orbital dynamics of the system following the DART impact. Given known fundamental differences between our simulation codes, we find that faster rigid-body simulations produce nearly the same result as rubble-pile models in scenarios with a moderate value for the momentum enhancement factor, $\beta$ ($\beta{\sim}3$) and an ellipsoidal secondary. This indicates that the rigid-body approach is likely adequate for propagating the post-impact dynamics necessary to meet DART Mission requirements. Although, if Dimorphos has a highly-irregular shape or structure, or if $\beta$ is unexpectedly large, then rubble-pile effects may become important. If Dimorphos's orbit and spin state are sufficiently excited, then surface particle motion is also possible. However, these simulations are limited in their resolution and range of material parameters, so they serve as a demonstration of principle, and Future work is required to fully understand the likelihood and magnitude of surface motion.

We present a new method for the detection of double-lined spectroscopic binaries (SB2) $v\sin{i}$ values from spectral fits. The method is tested on synthetic and real spectra from LAMOST-MRS. It can reliably detect SB2 candidates for double-lined binaries with $v\sin{i}_1+v\sin{i}_2<300\,km\,s^{-1}$ if the radial velocity separation is large enough. Using this method, we detect 2460 SB2 candidates, 1410 of which are new discoveries. We confirm the correlation between the radial velocity separation estimated by the binary model and $v\sin{i}_0$ estimated by the single star model using the selected sample. Additionally, our method finds one new SB2 candidate in open cluster M~11.

Emily O. Garvin, Sandor Kruk, Claude Cornen, Rachana Bhatawdekar, Raoul Cañameras, Bruno Merín

24 page, 14 figures, 5 tables, accepted for publication in A&A June 28 2022

The Hubble Space Telescope (HST) archives constitute a rich dataset of high resolution images to mine for strong gravitational lenses. While many HST programs specifically target strong lenses, they can also be present by coincidence in other HST observations. We aim to identify non-targeted strong gravitational lenses in almost two decades of images from the ESA it Hubble Space Telescope archive (eHST), without any prior selection on the lens properties. We used crowdsourcing on the Hubble Asteroid Hunter (HAH) citizen science project to identify strong lenses, alongside asteroid trails, in publicly available large field-of-view HST images. We visually inspected 2354 objects tagged by citizen scientists as strong lenses to clean the sample and identify the genuine lenses. We report the detection of 252 strong gravitational lens candidates, which were not the primary targets of the HST observations. 198 of them are new, not previously reported by other studies, consisting of 45 A grades, 74 B grades and 79 C grades. The majority are galaxy-galaxy configurations. The newly detected lenses are, on average, 1.3 magnitudes fainter than previous HST searches. This sample of strong lenses with high resolution HST imaging is ideal to follow-up with spectroscopy, for lens modelling and scientific analyses. This paper presents an unbiased search of lenses, which enabled us to find a high variety of lens configurations, including exotic lenses. We demonstrate the power of crowdsourcing in visually identifying strong lenses and the benefits of exploring large archival datasets. This study shows the potential of using crowdsourcing in combination with artificial intelligence for the detection and validation of strong lenses in future large-scale surveys such as ESA's future mission Euclid or in JWST archival images.

Derek C. Richardson, Harrison F. Agrusa, Brent Barbee, William F. Bottke, Andrew F. Cheng, Siegfried Eggl, Fabio Ferrari, Masatoshi Hirabayashi, Özgür Karatekin, Jay McMahon, Stephen R. Schwartz, Ronald-Louis Ballouz, Adriano Campo Bagatin, Elisabetta Dotto, Eugene G. Fahnestock, Oscar Fuentes-Muñoz, Ioannis Gkolias, Douglas P. Hamilton, Seth A. Jacobson, Martin Jutzi, Josh Lyzhoft, Rahil Makadia, Alex J. Meyer, Patrick Michel, Ryota Nakano, Guillaume Noiset, Sabina D. Raducan, Nicolas Rambaux, Alessandro Rossi, Paul Sánchez, Daniel J. Scheeres, Stefania Soldini, Angela M. Stickle, Paolo Tanga, Kleomenis Tsiganis, Yun Zhang

23 pages, 13 figures, published in PSJ

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $\beta$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $\beta$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous.

Yikun Zhang, Rafael S. de Souza, Yen-Chi Chen

20 pages, 9 figures, 2 tables

The latticework structure known as the cosmic web provides a valuable insight into the assembly history of large-scale structures. Despite the variety of methods to identify the cosmic web structures, they mostly rely on the assumption that galaxies are embedded in a Euclidean geometric space. Here we present a novel cosmic web identifier called SCONCE (Spherical and CONic Cosmic wEb finder) that inherently considers the 2D (RA,DEC) spherical or the 3D (RA,DEC,$z$) conic geometry. The proposed algorithms in SCONCE generalize the well-known subspace constrained mean shift (SCMS) method and primarily address the predominant filament detection problem. They are intrinsic to the spherical/conic geometry and invariant to data rotations. We further test the efficacy of our method with an artificial cross-shaped filament example and apply it to the SDSS galaxy catalogue, revealing that the 2D spherical version of our algorithms is robust even in regions of high declination. Finally, using N-body simulations from Illustris, we show that the 3D conic version of our algorithms is more robust in detecting filaments than the standard SCMS method under the redshift distortions caused by the peculiar velocities of halos. Our cosmic web finder is packaged in python as SCONCE-SCMS and has been made publicly available.

We analyze 152 large confined flares (GOES class $\geq$M1.0 and $\leq$$45^{\circ}$ from disk center) during 2010$-$2019, and classify them into two types according to the criterion taken from the work of Li et al. (2019). "Type I" flares are characterized by slipping motions of flare loops and ribbons and a stable filament underlying the flare loops. "Type II" flares are associated with the failed eruptions of the filaments, which can be explained by the classical 2D flare model. A total of 59 flares are "Type I" flares (about 40\%) and 93 events are "Type II" flares (about 60\%). There are significant differences in distributions of the total unsigned magnetic flux ($\Phi$$_\mathrm{AR}$) of active regions (ARs) producing the two types of confined flares, with "Type I" confined flares from ARs with a larger $\Phi$$_{AR}$ than "Type II". We calculate the mean shear angle $\Psi$$_\mathrm{HFED}$ within the core of an AR prior to the flare onset, and find that it is slightly smaller for "Type I" flares than that for "Type II" events. The relative non-potentiality parameter $\Psi$$_\mathrm{HFED}$/$\Phi$$_\mathrm{AR}$ has the best performance in distinguishing the two types of flares. About 73\% of "Type I" confined flares have $\Psi$$_\mathrm{HFED}$/$\Phi$$_\mathrm{AR}$$<$1.0$\times$$10^{-21}$ degree Mx$^{-1}$, and about 66\% of "Type II" confined events have $\Psi$$_\mathrm{HFED}$/$\Phi$$_\mathrm{AR}$$\geq$1.0$\times$$10^{-21}$ degree Mx$^{-1}$. We suggest that "Type I" confined flares cannot be explained by the standard flare model in 2D/3D, and the occurrence of multiple slipping magnetic reconnections within the complex magnetic systems probably leads to the observed flare.

C. de la Fuente Marcos, R. de la Fuente Marcos

37 pages, 14 figures, 10 tables, 3 appendices. Accepted for publication in Celestial Mechanics and Dynamical Astronomy

The region where the main asteroid belt is now located may have started empty, to become populated early in the history of the Solar system with material scattered outward by the terrestrial planets and inward by the giant planets. These dynamical pathways toward the main belt may still be active today. Here, we present results from a data mining experiment aimed at singling out present-day members of the main asteroid belt that may have reached the belt during the last few hundred years. Probable newcomers include 2003 BM1, 2007 RS62, 457175 (2008 GO98), 2010 BG18, 2010 JC58, 2010 JV52, 2010 KS6, 2010 LD74, 2010 OX38, 2011 QQ99, 2013 HT149, 2015 BH103, 2015 BU525, 2015 RO127, 2015 RS139, 2016 PC41, 2016 UU231, 2020 SA75, 2020 UO43, and 2021 UJ5, all of them in the outer belt. Some of these candidates may have been inserted in their current orbits after experiencing relatively recent close encounters with Jupiter. We also investigated the likely source regions of such new arrivals. Asteroid 2020 UO43, if real, has a non-negligible probability of having an origin in the Oort cloud or even interstellar space. Asteroid 2003 BM1 may have come from the neighborhood of Uranus. However, most newcomers -- including 457175, 2011 QQ99, and 2021 UJ5 -- might have had an origin in Centaur orbital space. The reliability of these findings is assessed within the context of the uncertainties of the available orbit determinations.

N. Walters, J. Farihi, T. R. Marsh, E. Breedt, P. W. Cauley, T. von Hippel, J. J. Hermes

14 pages, 8 figures, 3 tables. Submitted to MNRAS

This study presents observations of two possible substellar survivors of post-main sequence engulfment, currently orbiting white dwarf stars. Infrared and optical spectroscopy of GD 1400 reveal a 9.98 h orbital period, where the benchmark brown dwarf has $M_2=65\pm3$ M$_{\rm Jup}$, $T_{\rm eff}\approx1900$ K, and an age close to 1 Gyr. The $0.558\pm0.005$ M$_{\odot}$ mass of GD 1400A suggests the substellar companion escaped contact with the RGB host, but was later enveloped during the AGB. The known infrared excess towards PG 0010+281 is consistent with a substellar companion, yet no radial velocity or photometric variability is found despite multiple instrument searches. Three independent mass determinations for PG 0010+281 all suggest enhanced mass loss associated with binary evolution, where the youngest total age for an isolated star is $7.5\pm2.5$ Gyr. A possible solution to this conundrum is the cannibalization of one or more giant planets, which enhanced mass loss post-main sequence, but were ultimately destroyed. Thus, PG 0010+281 is likely orbited by a debris disk that is comfortably exterior to the Roche limit, adding to the growing number of non-canonical disks orbiting white dwarfs. At present, only L-type (brown) dwarfs are known to survive direct engulfment during the post-main sequence, whereas T- and Y-type substellar companions persist at wide separations. These demographics suggest that roughly 50 M$_{\rm Jup}$ is required to robustly avoid post-main sequence annihilation, suggesting all closely-orbiting giant planets are consumed, which may contribute to mass loss and magnetic field generation in white dwarfs and their immediate progenitors.

G. Del Zanna, V. Polito, J. Dudík, P. Testa, H.E. Mason, E. Dzifčáková

Accepted for publication in ApJ

We present joint Hinode/EIS and IRIS observations of Fe XII lines in active regions, both on-disk and off-limb. We use an improved calibration for the EIS data, and find that the 192.4 A / 1349 A observed ratio is consistent with the values predicted by CHIANTI and the coronal approximation in quiescent areas, but not in all active region observations, where the ratio is often lower than expected by up to a factor of about two. We investigate a number of physical mechanisms that could affect this ratio, such as opacity and absorption from cooler material. We find significant opacity in the EIS Fe XII 193 and 195 A lines, but not in the 192.4 A line, in agreement with previous findings. As we cannot rule out possible EUV absorption by H, He and He II in the on-disk observations, we focus on an off-limb observation where such absorption is minimal. After considering these, as well as possible non-equilibrium effects, we suggest that the most likely explanation for the observed low Fe XII 192.4 A / 1349 A ratio is the presence of non-Maxwellian electron distributions in the active regions. This is in agreement with previous findings based on EIS and IRIS observations independently.

Grace E. Chesmore, Kathleen Harrington, Carlos E. Sierra, Patricio A. Gallardo, Shreya Sutariya, Tommy Alford, Alexandre E. Adler, Tanay Bhandarkar, Gabriele Coppi, Nadia Dachlythra, Joseph Golec, Jon Gudmundsson, Saianeesh K. Haridas, Bradley R. Johnson, Anna M. Kofman, Jeffrey Iuliano, Jeff McMahon, Michael D. Niemack, John Orlowski-Scherer, Karen Perez Sarmiento, Roberto Puddu, Max Silva-Feaver, Sara M. Simon, Julia Robe, Edward J. Wollack, Zhilei Xu

in process

We present near-field radio holography measurements of the Simons Observatory Large Aperture Telescope Receiver optics. These measurements demonstrate that radio holography of complex millimeter-wave optical systems comprising cryogenic lenses, filters, and feed horns can provide detailed characterization of wave propagation before deployment. We used the measured amplitude and phase, at 4K, of the receiver near-field beam pattern to predict two key performance parameters: 1) the amount of scattered light that will spill past the telescope to 300K and 2) the beam pattern expected from the receiver when fielded on the telescope. These cryogenic measurements informed the removal of a filter, which led to improved optical efficiency and reduced side-lobes at the exit of the receiver. Holography measurements of this system suggest that the spilled power past the telescope mirrors will be less than 1% and the main beam with its near side-lobes are consistent with the nominal telescope design. This is the first time such parameters have been confirmed in the lab prior to deployment of a new receiver. This approach is broadly applicable to millimeter and sub-millimeter instruments.

R. Aafreen, R. Abhishek, B. Ajithkumar, Arunkumar M. Vaidyanathan, Indrajit. V. Barve, Sahana Bhattramakki, Shashank Bhat, B. S. Girish, Atul Ghalame, Y. Gupta, Harshal G. Hayatnagarkar, P. A. Kamini, A. Karastergiou, L. Levin, S. Madhavi, M. Mekhala, M. Mickaliger, V. Mugundhan, Arun Naidu, J. Oppermann, B. Arul Pandian, N. Patra, A. Raghunathan, Jayanta Roy, Shiv Sethi, Benjamin Shaw, K. Sherwin, O. Sinnen, S. K. Sinha, K. S. Srivani, B. Stappers, C. R.Subrahmanya, Thiagaraj Prabu, C. Vinutha, Y. G. Wadadekar, Haomiao Wang, C. Williams

Accepted for JoAA, SKA Special issue on SKA (2022)

This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwidths, forming several 100s of beams, and processing the data in the SKA real-time pulsar search pipelines. A brief account of the different architectures of the accelerators, primarily the new generation Field Programmable Gate Array-based accelerators, showing their critical roles to achieve high-performance computing and in handling the enormous data volume problems of the SKA is presented here. It also presents the power-performance efficiency of this emerging technology and presents potential future scenarios.

Jin Beniyama, Shigeyuki Sako, Ryou Ohsawa, Satoshi Takita, Naoto Kobayashi, Shin-ichiro Okumura, Seitaro Urakawa, Makoto Yoshikawa, Fumihiko Usui, Fumi Yoshida, Mamoru Doi, Yuu Niino, Toshikazu Shigeyama, Masaomi Tanaka, Nozomu Tominaga, Tsutomu Aoki, Noriaki Arima, Ko Arimatsu, Toshihiro Kasuga, Sohei Kondo, Yuki Mori, Hidenori Takahashi, Jun-ichi Watanabe

This article is published in PASJ as open access, published by OUP ( this https URL ). 27 pages, 16 figures

We report the results of video observations of tiny (diameter less than 100 m) near-Earth objects (NEOs) with Tomo-e Gozen on the Kiso 105 cm Schmidt telescope. A rotational period of a tiny asteroid reflects its dynamical history and physical properties since smaller objects are sensitive to the YORP effect. We carried out video observations of 60 tiny NEOs at 2 fps from 2018 to 2021 and successfully derived the rotational periods and axial ratios of 32 NEOs including 13 fast rotators with rotational periods less than 60 s. The fastest rotator found during our survey is 2020 HS7 with a rotational period of 2.99 s. We statistically confirmed that there is a certain number of tiny fast rotators in the NEO population, which have been missed with any previous surveys. We have discovered that the distribution of the tiny NEOs in a diameter and rotational period (D-P) diagram is truncated around a period of 10 s. The truncation with a flat-top shape is not explained well either by a realistic tensile strength of NEOs or suppression of YORP by meteoroid impacts. We propose that the dependence of the tangential YORP effect on the rotational period potentially explains the observed pattern in the D-P diagram.

John Orlowski-Scherer, Saianeesh K. Haridas, Luca Di Mascolo, Karen Perez Sarmiento, Charles E. Romero, Simon Dicker, Tony Mroczkowski, Tanay Bhandarkar, Eugene Churazov, Tracy E Clarke, Mark Devlin, Massimo Gaspari, Ian Lowe, Brian Mason, Craig L Sarazin, Jonathon Sievers, Rashid Sunyaev

8 pages, 3 figures

Mechanical feedback from active galactic nuclei (AGN) is thought to be the dominant feedback mechanism quenching cooling flows and star formation in galaxy cluster cores. However, the mechanisms by which AGN couple to the intracluster medium (ICM) are not well understood. The nature of pressure supporting the cavities is not known. Using the MUSTANG-2 instrument on the Green Bank Telescope (GBT), we aimed to measure thermal Sunyaev-Zeldovich (SZ) effect signals associated with the X-ray cavities in MS0735.6+7421, a moderate mass cluster hosting one of the most energetic AGN outbursts known. We use these measurements to infer the level of non-thermal sources of pressure, such as magnetic fields and turbulence, as well as relativistic and cosmic ray components, supporting the cavities. We used preconditioned gradient descent to fit a model for the cluster, cavities, and central point source directly to the time ordered data of the MUSTANG-2 signal. We use this model to probe the thermodynamic state of the cavities. We have shown that the SZ signal associated with the cavities is suppressed compared to the expectations for a thermal plasma with the temperature $\sim$few tens keV. The smallest value of the suppression factor $f$ that is consistent with the data is $\sim$0.4, lower than inferred in earlier work. Larger values of $f$ are possible once the contribution of the cocoon shock surrounding the bubbles is taken into account. The baseline model with this particular geometrical setup yields best-fitting value f~0.5, which at face value implies a mix of thermal and non-thermal pressure support. Larger values of $f$ (up to 1, i.e. no tSZ signal from the bubbles) are still possible when allowing for variations in the line-of-sight geometry.

Guillaume Mahler, Mathilde Jauzac, Johan Richard, Benjamin Beauchesne, Harald Ebeling, David Lagattuta, Priyamvada Natarajan, Keren Sharon, Hakim Atek, Adélaïde Claeyssens, Benjamin Clément, Dominique Eckert, Alastair Edge, Brenda Frye, Jean-Paul Kneib, Anna Niemiec

link to the lensing outputs - this https URL

Exploiting the fundamentally achromatic nature of gravitational lensing, we present a lens model for the massive galaxy cluster SMACS\,J0723.3$-$7323 (SMACS\,J0723, $z=0.388$) that significantly improve upon earlier work. Building on strong-lensing constraints identified in prior Hubble Space Telescope (HST) observations, the mass model utilizes 21 multiple-image systems, 16 of which were newly discovered in Early Release Observational (ERO) data from the James Webb Space Telescope (JWST). The resulting lens model maps the cluster mass distribution to an RMS spatial precision of \rms and is publicly available https://www.dropbox.com/sh/3iatmz5k4hafzqf/AAAh0JvLgpBVoLp6qsxYZkFGa?dl=0. Consistent with previous analysis, our study shows SMACS\,J0723.3$-$7323 to be well described by a single cluster-scale component centered on the location of the brightest cluster galaxy, but deviates by adding two more diffuse components West of the cluster. A comparison of the galaxy distribution, the mass distribution, and gas distribution in the core of SMACS\,J0723 based on \HST - \JWST, and {\it Chandra} data finds a fairly concentrated regular elliptical profile but also signs of recent merger activity, possibly close to our line of sight. The exquisite sensitivity of NIRCAM \JWST\ reveals in spectacular fashion both the extended intra-cluster-light distribution and numerous star-forming clumps in magnified background galaxies. The high-precision lens model derived here for SMACS\,0723 demonstrated impressively the power of combining \HST\ and \JWST\ data for studies of structure formation and evolution in the distant Universe.

Massimo Pascale, Brenda Frye, Jose Diego, Lukas Furtak, Adi Zitrin, Tom Broadhurst, Chirstopher Conselice, Liang Dai, Leonardo Ferreira, Nathan Adams, Jeremy Lim, Asheesh Meena, Rogier Windhorst, Stephen Wilkins

The first deep field images from the \textit{James Webb Space Telescope} (JWST) of the galaxy cluster SMACS~J0723.3-7327 reveal a wealth of new lensed images at uncharted infrared wavelengths, with unprecedented depth and resolution. Here we securely identify 13 new sets of multiply lensed galaxies totalling 39 images, adding to the five sets of bright and multiply-imaged galaxies already known from Hubble data. We find examples of arcs crossing critical curves with magnification factors of at least 150, allowing detailed community follow-up, including JWST spectroscopy for precise redshift determinations, chemical abundances and detailed internal gas dynamics of very distant, young galaxies. We also detect an Einstein cross candidate only visible thanks to JWST's superb resolution. Our \emph{parametric} lens model is available through at https://www.dropbox.com/sh/gwup2lvks0jsqe5/AAC2RRSKce0aX-lIFCc9vhBXa?dl=0 , and will be regularly updated using additional spectroscopic redshifts. The model reproduces very well the multiple images, and allows for accurate magnification estimates of high-redshift galaxies. This work represents a first taste of the enhanced power JWST will have for lensing-related science.

We show that Dense Neural Networks can be used to accurately model the cooling of high-energy particles in the early universe, in the context of the public code package DarkHistory. DarkHistory self-consistently computes the temperature and ionization history of the early universe in the presence of exotic energy injections, such as might arise from the annihilation or decay of dark matter. The original version of DarkHistory uses large pre-computed transfer function tables to evolve photon and electron spectra in redshift steps, which require a significant amount of memory and storage space. We present a light version of DarkHistory that makes use of simple Dense Neural Networks to store and interpolate the transfer functions, which performs well on small computers without heavy memory or storage usage. This method anticipates future expansion with additional parametric dependence in the transfer functions without requiring exponentially larger data tables.

Raimon Luna, Gabriele Bozzola, Vitor Cardoso, Vasileios Paschalidis, Miguel Zilhão

12 pages, 10 figures

We compute the emission of linear momentum (kicks) by both gravitational and electromagnetic radiation in fully general-relativistic numerical evolutions of quasi-circular charged black hole binaries. We derive analytical expressions for slowly moving bodies and explore numerically a variety of mass ratios and charge-to-mass ratios. We find that for the equal mass case our analytical expression is in excellent agreement with the observed values and, contrarily to what happens in the vacuum case, we find that in presence of electromagnetic fields there is emission of momentum by gravitational waves. We also find that the strong gravitational kicks of binaries with unequal masses affect the electromagnetic kicks, causing them to strongly deviate from Keplerian predictions. For the values of charge-to-mass ratio considered in this work, we observe that magnitudes of the electromagnetic kicks are always smaller than the gravitational ones.

Yuki Fujimoto, Kenji Fukushima, Larry D. McLerran, Michal Praszalowicz

6 pages, 3 figures

We discuss an interpretation that a peak in the sound velocity in neutron star matter, as suggested by the observational data, signifies strongly-coupled conformal matter. The normalized trace anomaly is a dimensionless measure of conformality leading to the derivative and the non-derivative contributions to the sound velocity. We find that the peak in the sound velocity is attributed to the derivative contribution from the trace anomaly that steeply approaches the conformal limit. Smooth continuity to the behavior of high-density QCD implies that the matter part of the trace anomaly may be positive definite. We discuss a possible implication of the positivity condition of the trace anomaly on the $M$-$R$ relation of the neutron stars.

Serge Koutchmy

40 pages, 14 Figures. For the special issue on "Refraction atmospherique et astronomie" of the Compte-Rendus Physique Acad. Sc. Paris, Ed. L. Dettwiller and P. L\'ena

The main observations of 1761 by M. Lomonossov and those that followed are recalled by extending the discussion to other remarkable visual observations of the passages, then with more and more powerful imagers producing images in profusion. The modern treatment of parasitic effects is briefly recalled by focusing on the expert observation of 1761 which has recently been widely commented on and criticized. It included a spurious effect called the "black drop effect". The shell or aureole or atmospheric ring of Venus observed outside the solar disk is considered with reference to the today parameters of the Venus atmosphere. The contacts during the transit are discussed taking into account effects of scattering, absorption and the dominant effects of the refraction at the small angular distances found to be comparable to a fraction of the angular dimension of the planet. Modern observations of the 2004 and the 2012 transit are tentatively discussed to elucidate what is the arc of Lomonossov?

Brian Batell, Nikita Blinov, Christopher Hearty, Robert McGehee

21 pages, 8 figures; contribution to Snowmass 2021

A broad program of searches at high intensity experiments during the coming decade and beyond will sensitively probe new light mediator particles interacting through the minimal renormalizable vector, Higgs, and neutrino portals as well as higher-dimension axion-like particle portals. These portals may link the visible and dark sectors and play a critical role in many proposed solutions to some of the big open questions in particle physics and cosmology. In this whitepaper, we survey the theoretical and experimental progress, status, and prospects in the study of minimal dark sector portals.

Raghuveer Garani, Michel H.G. Tytgat, Jérôme Vandecasteele

10+19 pages, 6+8 figures

We explore the possible phases of a condensed dark matter (DM) candidate taken to be in the form of a fermion with a Yukawa coupling to a scalar particle, at zero temperature but at finite density. This theory essentially depends on only four parameters, the Yukawa coupling, the fermion mass, the scalar mediator mass, and the DM density. At low fermion densities we delimit the Bardeen-Cooper-Schrieffer (BCS), Bose-Einstein Condensate (BEC) and crossover phases as a function of model parameters using the notion of scattering length. We further study the BCS phase by consistently including emergent effects such as the scalar density condensate and superfluid gaps. Within the mean field approximation, we derive the consistent set of gap equations, retaining their momentum dependence, and valid in both the non-relativistic and relativistic regimes. We present numerical solutions to the set of gap equations, in particular when the mediator mass is smaller and larger than the DM mass. Finally, we discuss the equation of state (EoS) and possible astrophysical implications for asymmetric DM.

Recent studies of quantum field theory in FLRW spacetime suggest that the cause of the speeding up of the universe is the quantum vacuum, no need to introduce ad hoc scalar fields. Appropriate renormalization of the energy-momentum tensor shows that the vacuum energy density is a smooth function of the Hubble rate and its derivatives: $\rho_{\rm vac}=\rho_{\rm vac}(H, \dot{H},\ddot{H},...)$. This is because in QFT the quantum scaling of $\rho_{\rm vac}$ with the renormalization point turns into cosmic evolution with $H$. As a result, any two nearby points of the cosmic expansion during the standard FLRW epoch are smoothly related through $\delta\rho_{\rm vac}\sim {\cal O}(H^2)$. In this scenario, no fine tuning is needed at all. What we call the `cosmological constant' $\Lambda$ is just the nearly sustained value of $\rho_{\rm vac}(H)$ around (any) given epoch. In the very early universe, higher (even) powers $\rho_{\rm vac}\sim{\cal O}(H^N)$ ($N=4,6,..$) triggered fast inflation during a short period in which $H=$const. In it, the equation of state (EoS) of the vacuum is very close to $w_{\rm vac}=-1$, but this ceases to be true during the FLRW era. Amazingly, the quantum vacuum acts as a formidable cosmic chameleon: it subsequently adopts the EoS of matter during the relativistic ($w_{\rm vac}=1/3$) and non-relativistic ($w_{\rm vac}=0$) epochs, and in the late universe it mimics quintessence, $w_{\rm vac}\gtrsim-1$, only to tend again to $-1$ in the remote future. In the transit, the quantum vacuum helps to solve the $H_0$ and $\sigma_8$ tensions.