Abstract visibility
Abstract text size

Papers for Friday, Apr 23 2021

Papers with local authors

Long Wang, David M. Hernandez

8 pages, 6 figures, submitted to MNRAS

0 votes
Paper 9 — arXiv:2104.10843
0 votes
Paper 9 — arXiv:2104.10843

It is well known that numerical errors grow exponentially in $N$-body simulations of gravitational bound stellar systems, but it is not well understood how the accuracy parameters of algorithms affect the physical evolution in simulations. By using the hybrid $N$-body code, PeTar, we investigate how escapers and the structure evolution of collisional stellar systems (e.g., star clusters) depend on the accuracy of long-range and short-range interactions. We study a group of simulations of ideal low-mass star clusters in which the accuracy parameters are varied. We find that although the number of escapers is different in individual simulations, its distribution from all simulations can be described by Poisson statistics. The density profile also has a similar feature. By using a self-consistent set-up of the accuracy parameters for long- and short-range interactions, such that orbits are resolved well enough, the physical evolution of the models is identical. But when the short-range accuracy is too low, a nonphysical dynamical evolution can easily occur; this is not the case for long-range interactions. This strengthens the need to include a proper algorithm (e.g. regularization methods) in the realistic modelling of collisional stellar systems. We also demonstrate that energy-conservation is not a good indicator to monitor the quality of the simulations. The energy error of the system is controlled by the hardest binary, and thus, it may not reflect the ensemble error of the global system.

D. S. Doelman, F. Snik, E. H. Por, S. P. Bos, G.P.P.L. Otten, M. Kenworthy, S. Y. Haffert, M. Wilby, A. J. Bohn, B. J. Sutlieff, K. Miller, M. Ouellet, J. de Boer, C. U. Keller, M. J. Escuti, S. Shi, N.Z. Warriner, K. J. Hornburg, J. L. BirkbyJ. Males, K. M. Morzinski, L.M. Close, J. Codona, J. Long, L. Schatz, J. Lumbres, A. Rodack, K. Van Gorkom, A Hedglen, O. Guyon, J. Lozi, T. Groff, J. Chilcote, N. Jovanovic, S. Thibault, C. de Jonge, G. Allain, C. Vallée, D. Patel, O. Côté, C. Marois, P. Hinz, J. Stone, A. Skemer, Z. Briesemeister, A. Boehle, A. M. Glauser, W. Taylor, P. Baudoz, E. Huby, O. Absil, B. Carlomagno, C. Delacroix

38 pages, 17 figures, accepted for publication in Applied Optics

0 votes
Paper 39 — arXiv:2104.11211
0 votes
Paper 39 — arXiv:2104.11211

Over the last decade, the vector-apodizing phase plate (vAPP) coronagraph has been developed from concept to on-sky application in many high-contrast imaging systems on 8-m class telescopes. The vAPP is an geometric-phase patterned coronagraph that is inherently broadband, and its manufacturing is enabled only by direct-write technology for liquid-crystal patterns. The vAPP generates two coronagraphic PSFs that cancel starlight on opposite sides of the point spread function (PSF) and have opposite circular polarization states. The efficiency, that is the amount of light in these PSFs, depends on the retardance offset from half-wave of the liquid-crystal retarder. Using different liquid-crystal recipes to tune the retardance, different vAPPs operate with high efficiencies ($>96\%$) in the visible and thermal infrared (0.55 $\mu$m to 5 $\mu$m). Since 2015, seven vAPPs have been installed in a total of six different instruments, including Magellan/MagAO, Magellan/MagAO-X, Subaru/SCExAO, and LBT/LMIRcam. Using two integral field spectrographs installed on the latter two instruments, these vAPPs can provide low-resolution spectra (R$\sim$30) between 1 $\mu$m and 5 $\mu$m. We review the design process, development, commissioning, on-sky performance, and first scientific results of all commissioned vAPPs. We report on the lessons learned and conclude with perspectives for future developments and applications.

Papers reserved for later discussion

These are papers reserved by people for discussion at a later date. All reservations are kept for 2 days after the date of the reservation.

Hajime Fukushima, Hidenobu Yajima

26 pages, 23 figures, submitted to MNRAS

0 votes
04/23/2021: arXiv:2104.10892
0 votes
04/23/2021: arXiv:2104.10892

By performing three-dimensional radiation hydrodynamics simulations, we study the formation of young massive star clusters (YMCs, $M_{*}>10^4~M_{\odot}$) in clouds with the surface density ranging from $\Sigma_{\rm cl} = 80$ to $3200~M_{\odot}\;{\rm pc^{-2}}$. We find that photoionization feedback suppresses star formation significantly in clouds with low surface density. Once the initial surface density exceeds $\sim 100~M_{\odot}\;{\rm pc^{-2}}$ for clouds with $M_{\rm cl}=10^{6}~M_{\odot}$ and $Z= Z_{\odot}$, most of the gas is converted into stars because the photoionization feedback is inefficient in deep gravitational potential. In this case, the star clusters are massive and gravitationally bounded as YMCs. The transition surface density increases as metallicity decreases, and it is $\sim 350~M_{\odot}\;{\rm pc^{-2}}$ for $Z=10^{-2}~Z_{\odot}$. We show that more than 10 percent of star-formation efficiency (SFE) is needed to keep a star cluster gravitationally bounded even after the disruption of a cloud. Also, we develop a semi-analytical model reproducing the SFEs obtained in our simulations. We find that the SFEs are fit with a power-law function with the dependency $\propto \Sigma_{\rm cl}^{1/2}$ for low-surface density and rapidly increases at the transition surface densities. The conditions of the surface density and the metallicity match recent observations of giant molecular clouds forming YMCs in nearby galaxies.

Sigurd Naess, Simone Aiola, Nick Battaglia, Richard J. Bond, Erminia Calabrese, Steve K. Choi, Nicholas F. Cothard, Mark Halpern, J. Colin Hill, Brian J. Koopman, Mark Devlin, Jeff McMahon, Simon Dicker, Adriaan J. Duivenvoorden, Jo Dunkley, Alexander Van Engelen, Valentina Fanfani, Simone Ferraro, Patricio A. Gallardo, Yilun Guan, Dongwon Han, Matthew Hasselfield, Adam D. Hincks, Kevin Huffenberger, Arthur B. Kosowsky, Thibaut Louis, Amanda Macinnis, Mathew S. Madhavacheril, Federico Nati, Michael D. Niemack, Lyman Page, Maria Salatino, Emmanuel Schaan, John Orlowski-Scherer, Alessandro Schillaci, Benjamin Schmitt, Neelima Sehgal, Cristóbal Sifón, Suzanne Staggs, Edward J. Wollack

23 pages, 10 figures, 5 tables

0 votes
04/22/2021: arXiv:2104.10264
0 votes
04/22/2021: arXiv:2104.10264

We use Atacama Cosmology Telescope (ACT) observations at 98 GHz (2015--2019), 150 GHz (2013--2019) and 229 GHz (2017--2019) to perform a blind shift-and-stack search for Planet 9. The search explores distances from 300 AU to 2000 AU and velocities up to 6.3 arcmin per year, depending on the distance. For a 5 Earth-mass Planet 9 the detection limit varies from 325 AU to 625 AU, depending on the sky location. For a 10 Earth-mass planet the corresponding range is 425 AU to 775 AU. The search covers the whole 18,000 square degrees of the ACT survey, though a slightly deeper search is performed for the parts of the sky consistent with Planet 9's expected orbital inclination. No significant detections are found, which is used to place limits on the mm-wave flux density of Planet 9 over much of its orbit. Overall we eliminate roughly 17% and 9% of the parameter space for a 5 and 10 Earth-mass Planet 9 respectively. We also provide a list of the 10 strongest candidates from the search for possible follow-up.

All other papers

Alexander D. Rathcke, Ryan J. MacDonald, Joanna K. Barstow, Jayesh M. Goyal, Mercedes Lopez-Morales, João M. Mendonça, Jorge Sanz-Forcada, Gregory W. Henry, David K. Sing, Munazza K. Alam, Nikole K. Lewis, Katy L. Chubb, Jake Taylor, Nikolay Nikolov, Lars A. Buchhave

27 pages, 10 figures. Resubmitted to AJ after referee report

We present a new optical transmission spectrum of the hot Jupiter WASP-79b. We observed three transits with the STIS instrument mounted on HST, spanning 0.3 - 1.0 um. Combining these transits with previous observations, we construct a complete 0.3 - 5.0 um transmission spectrum of WASP-79b. Both HST and ground-based observations show decreasing transit depths towards blue wavelengths, contrary to expectations from Rayleigh scattering or hazes. We infer atmospheric and stellar properties from the full near-UV to infrared transmission spectrum of WASP-79b using three independent retrieval codes, all of which yield consistent results. Our retrievals confirm previous detections of H$_{2}$O (at 4.0$\sigma$ confidence), while providing moderate evidence of H$^{-}$ bound-free opacity (3.3$\sigma$) and strong evidence of stellar contamination from unocculted faculae (4.7$\sigma$). The retrieved H$_{2}$O abundance ($\sim$ 1$\%$) suggests a super-stellar atmospheric metallicity, though stellar or sub-stellar abundances remain consistent with present observations (O/H = 0.3 - 34$\times$ stellar). All three retrieval codes obtain a precise H$^{-}$ abundance constraint: log(X$_{\rm{H^{-}}}$) $\approx$ -8.0 $\pm$ 0.7. The potential presence of H$^{-}$ suggests that JWST observations may be sensitive to ionic chemistry in the atmosphere of WASP-79b. The inferred faculae are $\sim$ 500 K hotter than the stellar photosphere, covering $\sim$ 15$\%$ of the stellar surface. Our analysis underscores the importance of observing UV - optical transmission spectra in order to disentangle the influence of unocculted stellar heterogeneities from planetary transmission spectra.

Elizabeth J. Gonzalez, Facundo Rodriguez, Manuel Merchán, Diego García Lambas, Martín Makler, Martín Chalela, Maria E. S. Pereira, Bruno Moraes, HuanYuan Shan

13 pages, 8 figures, plus appendices. Accepted in MNRAS

Galaxy group masses are important to relate these systems with the dark matter halo hosts. However, deriving accurate mass estimates is particularly challenging for low-mass galaxy groups. Moreover, calibration of bservational mass-proxies using weak-lensing estimates have been mainly focused on massive clusters. We present here a study of halo masses for a sample of galaxy groups identified according to a spectroscopic catalogue, spanning a wide mass range. The main motivation of our analysis is to assess mass estimates provided by the galaxy group catalogue derived through an abundance matching luminosity technique. We derive total halo mass estimates according to a stacking weak-lensing analysis. Our study allows to test the accuracy of mass estimates based on this technique as a proxy for the halo masses of large group samples. Lensing profiles are computed combining the groups in different bins of abundance matching mass, richness and redshift. Fitted lensing masses correlate with the masses obtained from abundance matching. However, when considering groups in the low- and intermediate-mass ranges, masses computed according to the characteristic group luminosity tend to predict higher values than the determined by the weak-lensing analysis. The agreement improves for the low-mass range if the groups selected have a central early-type galaxy. Presented results validate the use of mass estimates based on abundance matching techniques which provide good proxies to the halo host mass in a wide mass range.

Pablo Villanueva-Domingo, Kiyotomo Ichiki

18 pages, 8 figures

Primordial black holes (PBHs) as part of the Dark Matter (DM) would modify the evolution of large-scale structures and the thermal history of the universe. Future 21 cm forest observations, sensitive to small scales and the thermal state of the Inter Galactic Medium (IGM), could probe the existence of such PBHs. In this article, we show that the shot noise isocurvature mode on small scales induced by the presence of PBHs can enhance the amount of low mass halos, or minihalos, and thus, the number of 21 cm absorption lines. However, if the mass of PBHs is as large as $M_{\rm PBH}\gtrsim 10 \, M_\odot$, with an abundant enough fraction of PBHs as DM, $f_{\rm PBH}$, the IGM heating due to accretion onto the PBHs counteracts the enhancement due to the isocurvature mode, reducing the number of absorption lines instead. The concurrence of both effects imprints distinctive signatures in the number of absorbers, allowing to bound the abundance of PBHs. We compute the prospects for constraining PBHs with future 21 cm forest observations, finding achievable competitive upper limits on the abundance as low as $f_{\rm PBH} \sim 10^{-3}$ at $M_{\rm PBH}= 100 \, M_\odot$, or even lower at larger masses, in unexplored regions of the parameter space by current probes. The impact of astrophysical X-ray sources on the IGM temperature is also studied, which could potentially weaken the bounds.

Turner Woody, Kevin C. Schlaufman

20 pages, 4 figures, and 5 tables in aastex63 format; accepted for publication in AJ

Globular clusters can form inside their host galaxies at high redshift when gas densities were higher and gas-rich mergers were common. They can also form inside lower-mass galaxies that have since been accreted and tidally disrupted, leaving their globular cluster complement bound to higher-mass halos. We argue that the age-metallicity-specific orbital energy relation in a galaxy's globular cluster system can be used to identify its origin. Gas-rich mergers should produce tightly bound systems in which metal-rich clusters are younger than metal-poor clusters. Globular clusters formed in massive disks and then scattered into a halo should have no relationship between age and specific orbital energy. Accreted globular clusters should produce weakly bound systems in which age and metallicity are correlated with each other but inversely correlated with specific orbital energy. We use precise relative ages, self-consistent metallicities, and space-based proper motion-informed orbits to show that the Milky Way's metal-poor globular cluster system lies in a plane in age-metallicity-specific orbital energy space. We find that relatively young or metal-poor globular clusters are weakly bound to the Milky Way, while relatively old or metal-rich globular clusters are tightly bound to the Galaxy. While metal-rich globular clusters may be formed either in situ or ex situ, our results suggest that metal-poor clusters formed outside of the Milky Way in now-disrupted dwarf galaxies. We predict that this relationship between age, metallicity, and specific orbital energy in a $L^{*}$ galaxy's globular cluster system is a natural outcome of galaxy formation in a $\Lambda$CDM universe.

Wei Yan, Ryan C. Hickox, Chien-Ting J. Chen, Claudio Ricci, Alberto Masini, Franz E. Bauer, David M. Alexander

8 pages, 5 figures, accepted by ApJ

In active galactic nuclei (AGN), fluorescent Fe K$\alpha$ (iron) line emission is generally interpreted as originating from obscuring material around a supermassive black hole (SMBH) on the scale of a few parsecs (pc). However, recent Chandra studies indicate the existence of iron line emission extending to kpc scales in the host galaxy. The connection between iron line emission and large-scale material can be spatially resolved directly only in nearby galaxies, but could be inferred in more distant AGNs by a connection between line emission and star-forming gas and dust that is more extended than the pc-scale torus. Here we present the results from a stacking analysis and X-ray spectral fitting performed on sources in the Chandra Deep Field South (CDFS) 7 Ms observations. From the deep stacked spectra, we select sources with stellar mass $\log(M_*/M_\odot)>10$ at $0.5<z<2$, obtaining 25 sources with high infrared luminosity ($ {\rm SFR}_{\rm FIR} \geq 17\;M_{\odot}\;{\rm yr}^{-1}$) and 32 sources below this threshold. We find that the equivalent width of the iron line EW(Fe) is a factor of three higher with 3$\sigma$ significance for high infrared luminosity measured from Herschel observations, indicating a connection between iron line emission and star-forming material on galaxy scales. We show that there is no significant dependence in EW(Fe) on $M_*$ or X-ray luminosity, suggesting the reflection of AGN X-ray emission over large scales in their host galaxies may be widespread.

John Chambers

Accepted for publication in The Astrophysical Journal

The formation of gas-giant planets within the lifetime of a protoplanetary disk is challenging especially far from a star. A promising model for the rapid formation of giant-planet cores is pebble accretion in which gas drag during encounters leads to high accretion rates. Most models of pebble accretion consider disks with a monotonic, radial pressure profile. This causes a continuous inward flux of pebbles and inefficient growth. Here we examine planet formation in a disk with multiple, intrinsic pressure bumps. In the outer disk, pebbles become trapped near these bumps allowing rapid growth under suitable conditions. In the inner disk, pebble traps may not exist because the inward gas advection velocity is too high. Pebbles here are rapidly removed. In the outer disk, growth is very sensitive to the initial planet mass and the strength of turbulence. This is because turbulent density fluctuations raise planetary eccentricities, increasing the planet-pebble relative velocity. Planetary seeds above a distance-dependent critical mass grow to a Jupiter mass in 0.5--3 million years out to at least 60 AU in a 0.03 solar-mass disk. Smaller bodies remain near their initial mass, leading to a sharp dichotomy in growth outcomes. For turbulent alpha = 1e-4, the critical masses are 1e-4 and 1e-3 Earth masses at 9 and 75 AU, respectively. Pressure bumps in disks may explain the large mass difference between the giant planets and Kuiper belt objects, and also the existence of wide-orbit planets in some systems.

The topic of the present study is combining a dynamic model of a protoplanetary disk with the computations of radiation transfer for obtaining synthetic spectra and disk images suitable for immediate comparison of the model with observations. Evolution of the disk was computed using the FEOSAD hydrodynamic model, which includes a self-consistent calculation of the dynamics of dust and gas in the 2D thin disk approximation. Radiation transfer was simulated by the open code RADMC-3D. Three phases of disk evolution were considered: a young gravitationally unstable disk, a disk during an accretion luminosity burst, and an evolved disk. For these stages, the influence of various processes upon the disk's thermal structure was analyzed, as well as the differences between the temperatures obtained in the initial dynamic model and in the model with a detailed calculation of the radiation transfer. It is shown that viscous heating in the inner regions and adiabatic heating in the disk spirals can be important sources of heating. On the basis of the calculated spectral energy distributions, using SED-fitter software package used for the observations, physical parameters of the model disks were reconstructed. A significant spread between reconstructed parameters and initial characteristics of the disk indicates verification necessity of the models within the framework of spatially resolved observations of disks in the different spectral ranges

A new $\gamma$-ray supernova remnant (SNR) G206.9+2.3 is first reported in this study, analyzing 12.4 years of observation data of the Fermi Large Area Telescope (Fermi-LAT). Its $\gamma$-ray spatial distribution did not appear extended feature. Its photon flux was (3.68$\pm$1.21) 10$^{-10}$ cm$^{-2}$ s$^{-1}$, and its power-law spectral index was 3.37$\pm$0.57 in the 0.8-500 GeV energy band. By subtracting the background around the SNR G206.9+2.3 region, we found that its TS values of the global fit from the four different energy bands were greater than 9. We identify that this is a real $\gamma$-ray signal. Furthermore, we found that its GeV spatial location was in good agreement with that of its radio band. Its 12.4 years of the light curve (LC) did not exhibit significant variability. We suggest that the new $\gamma$-ray source is a likely counterpart to SNR G206.9+2.3.

Qianwei Qu, Sergei N. Yurchenko, Jonathan Tennyson

13 pages, 14 figures, 2 tables

An accurate line list, called XABC, is computed for nitric oxide which covers its pure rotational, vibrational and rovibronic spectra. A mixture of empirical and theoretical electronic transition dipole moments are used for the final calculation of $^{14}\mathrm{N}^{16}\mathrm{O}$ rovibronic $\mathrm{A}\,^2\Sigma^+$ -- $\mathrm{X}\,^2\Pi$, $\mathrm{B}\,^2\Pi$ -- $\mathrm{X}^2\Pi$ and $\mathrm{C}\,^2\Pi$ -- $\mathrm{X}\,^2\Pi$ which correspond to the $\gamma$, $\beta$ and $\delta$ band systems, respectively, as well as minor improvements to transitions within the $\mathrm{X}\,^2\Pi$ ground state. The work is a major update of the ExoMol NOname line list. It provides a high-accuracy NO ultraviolet line list covering the complicated regions where the $\mathrm{B}\,^2\Pi$-$\mathrm{C}\,^2\Pi$ states interact. XABC provides comprehensive data for the lowest four doublet states of NO in the region of $\lambda > 160 ~ \mathrm{nm}$ ($\tilde{\nu} < 63~000~\mathrm{cm}^{-1}$) for the analysis of atmospheric NO on Earth, Venus or Mars, other astronomical observations and applications. The data is available via www.exomol.com.

Alex J. Kemp, Amanda I. Karakas, Andrew R. Casey, Robert G. Izzard, Ashley J. Ruiter, Floor S. Broekgaarden, Karel D. Temmink

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code \binaryc. This treatment of novae represents a significant improvement on the `averaging' treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 \pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14\pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year.

Louise Harra, Vincenzo Andretta, Thierry Appourchaux, Frédéric Baudin, Luis Bellot-Rubio, Aaron C. Birch, Patrick Boumier, Robert H. Cameron, Matts Carlsson, Thierry Corbard, Jackie Davies, Andrew Fazakerley, Silvano Fineschi, Wolfgang Finsterle, Laurent Gizon, Richard Harrison, Donald M. Hassler, John Leibacher, Paulett Liewer, Malcolm MacDonald, Milan Maksimovic, Neil Murphy, Giampiero Naletto, Giuseppina Nigro, Christopher Owen, Valentín Martínez-Pillet, Pierre Rochus, Marco Romoli, Takashi Sekii, Daniele Spadaro, Astrid Veronig, Werner Schmutz

29 pages, 9 figures, White paper for ESA Voyage 2050

A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by any other mission. One of the major mysteries of the Sun is the solar cycle. The activity cycle of the Sun drives the structure and behaviour of the heliosphere and is, of course, the driver of space weather. In addition, solar activity and variability provides fluctuating input into the Earth climate models, and these same physical processes are applicable to stellar systems hosting exoplanets. One of the main obstructions to understanding the solar cycle, and hence all solar activity, is our current lack of understanding of the polar regions. In this White Paper, submitted to the European Space Agency in response to the Voyage 2050 call, we describe a mission concept that aims to address this fundamental issue. In parallel, we recognise that viewing the Sun from above the polar regions enables further scientific advantages, beyond those related to the solar cycle, such as unique and powerful studies of coronal mass ejection processes, from a global perspective, and studies of coronal structure and activity in polar regions. Not only will these provide important scientific advances for fundamental stellar physics research, they will feed into our understanding of impacts on the Earth and other planets' space environment.

T. Mishenina, N. Basak, V. Adibekyan, C. Soubiran, V. Kovtyukh

25 pages, 11 figures, Accepted in MNRAS

Stellar parameters of 25 planet-hosting stars and abundances of Li, C, O, Na, Mg, Al, S, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Ni, Zn, Y, Zr, Ba, Ce, Pr, Nd, Sm and Eu, were studied based on homogeneous high resolution spectra and uniform techniques. The iron abundance [Fe/H] and key elements (Li, C, O, Mg, Si) indicative of the planet formation, as well as the dependencies of [El/Fe] on $T_{cond}$, were analyzed. The iron abundances determined in our sample stars with detected massive planets range within -0.3<[Fe/H]<0.4. The behaviour of [C/Fe], [O/Fe], [Mg/Fe] and [Si/Fe] relative to [Fe/H] is consistent with the Galactic Chemical Evolution trends. The mean values of C/O and [C/O] are <C/O>= 0.48 +/-0.07 and <[C/O]>=-0.07 +/-0.07, which are slightly lower than solar ones. The Mg/Si ratios range from 0.83 to 0.95 for four stars in our sample and from 1.0 to 1.86 for the remaining 21 stars. Various slopes of [El/Fe] vs. Tcond were found. The dependencies of the planetary mass on metallicity, the lithium abundance, the C/O and Mg/Si ratios, and also on the [El/Fe]-Tcond slopes were considered.

Chen Zhengyi, Gu Qiu-Sheng, Rubén García-Benito, Zhang Zhi-Yu, Ge Xue, Xiao Mengyuan, Yu Xiaoling

16 pages, 13 figures, accepted for publication in ApJ

Lenticular galaxies (S0s) were considered mainly as passive evolved spirals due to environmental effects for a long time; however, most S0s in the field cannot fit into this common scenario. In this work, we study one special case, SDSS J120237.07+642235.3 (PGC 38025), a star-forming field S0 galaxy with an off-nuclear blue core. We present optical integral field spectroscopic (IFS) observation with the 3.5 meter telescope at Calar Alto (CAHA) Observatory, and high-resolution millimeter observation with the NOrthern Extended Millimeter Array (NOEMA). We estimated the star formation rate (SFR = 0.446 $M_\odot yr^{-1}$) and gaseous metallicity (12 + log(O/H) = 8.42) for PGC 38025, which follows the star formation main sequence and stellar mass - metallicity relation. We found that the ionized gas and cold molecular gas in PGC 38025 show the same spatial distribution and kinematics, whilst rotating misaligned with stellar component. The off-nuclear blue core is locating at the same redshift as PGC 38025 and its optical spectrum suggest it is \rm H\,{\sc ii} region. We suggest that the star formation in PGC 38025 is triggered by a gas-rich minor merger, and the off-nuclear blue core might be a local star-formation happened during the accretion/merger process.

Oliver J. Hall, Guy R. Davies, Jennifer van Saders, Martin B. Nielsen, Mikkel N. Lund, William J. Chaplin, Rafael A. García, Louis Amard, Angela A. Breimann, Saniya Khan, Victor See, Jamie Tayar

48 pages, 10 figures, 1 table. Published in Nature Astronomy

Studies using asteroseismic ages and rotation rates from star-spot rotation have indicated that standard age-rotation relations may break down roughly half-way through the main sequence lifetime, a phenomenon referred to as weakened magnetic braking. While rotation rates from spots can be difficult to determine for older, less active stars, rotational splitting of asteroseismic oscillation frequencies can provide rotation rates for both active and quiescent stars, and so can confirm whether this effect really takes place on the main sequence. We obtained asteroseismic rotation rates of 91 main sequence stars showing high signal-to-noise modes of oscillation. Using these new rotation rates, along with effective temperatures, metallicities and seismic masses and ages, we built a hierarchical Bayesian mixture model to determine whether the ensemble more closely agreed with a standard rotational evolution scenario, or one where weakened magnetic braking takes place. The weakened magnetic braking scenario was found to be 98.4% more likely for our stellar ensemble, adding to the growing body of evidence for this stage of stellar rotational evolution. This work represents the largest catalogue of seismic rotation on the main sequence to date, opening up possibilities for more detailed ensemble analysis of rotational evolution with Kepler.

Aurélien Stcherbinine, Mathieu Vincendon, Franck Montmessin, Pierre Beck

Submitted to Icarus

Mars northern polar latitudes are known to harbor an enhanced 3 ${\mu}$m spectral signature when observed from orbit. This may indicate a greater amount of surface adsorbed or bound water, although it has not yet been possible to easily reconcile orbital observations with ground measurements by Phoenix. Here we reprocessed OMEGA/Mars Express observations acquired during the Northern summer to further characterize this 3 ${\mu}$m absorption band increase. We identify the presence of a new specific spectral signature composed of an additional narrow absorption feature centered at 3.03 ${\mu}$m coupled with an absorption at ${\lambda}$ ${\geq}$ 3.8 ${\mu}$m. This signature is homogeneously distributed over a bright albedo open ring surrounding the circumpolar low-albedo terrains between ~ 68{\deg}N and 76{\deg}N and ~ 0{\deg}E and 270{\deg}E. This location includes the Phoenix landing site. This feature shows no time variability and can be confidently attributed to a seasonally stable surface component. All together, the stability, spectral shape and absence of significant correlation with other signatures in the 1 $-$ 2.5 ${\mu}$m range discard interpretations relying on water ice or easily exchangeable adsorbed water. The exact full spectral shape cannot be easily reproduced by pure minerals samples, although sulfates, notably lowly hydrated Ca-sulfates, provide interesting comparisons. A modification of the chemical or physical properties of the soil, potentially involving additional sulfates contaminants, or modification of the hydration state of sulfates, and/or modification of their grains size, seems a plausible explanation to this observation, which may then indicate geologically recent water alteration at high northern latitudes.

Elimboto Yohana, Yin-Zhe Ma, Di Li, Xuelei Chen, Wei-Ming Dai

14 pages, 11 figures, accepted in MNRAS

The 21-cm intensity mapping (IM) of neutral hydrogen (HI) is a promising tool to probe the large-scale structures. Sky maps of 21-cm intensities can be highly contaminated by different foregrounds, such as Galactic synchrotron, free-free emission, extragalactic point sources, and atmospheric noise. We here present a model of foreground components and a method of removal, especially to quantify the potential of Five-hundred-meter Aperture Spherical radio Telescope (FAST) for measuring HI IM. We consider 1-year observational time with the survey area of $20,000\,{\rm deg}^{2}$ to capture significant variations of the foregrounds across both the sky position and angular scales relative to the HI signal. We first simulate the observational sky and then employ the Principal Component Analysis (PCA) foreground separation technique. We show that by including different foregrounds, thermal and $1/f$ noises, the value of the standard deviation between reconstructed 21-cm IM map and the input pure 21-cm signal is $\Delta T = 0.034\,{\rm mK}$, which is well under control. The eigenmode-based analysis shows that the underlying HI eigenmode is just less than $1$ per cent level of the total sky components. By subtracting the PCA cleaned foreground+noise map from the total map, we show that PCA method can recover HI power spectra for FAST with high accuracy.

Yajie Chen, Damien Przybylski, Hardi Peter, Hui Tian, F. Auchère, D. Berghmans

Accepted by A&A (Solar Orbiter special issue)

Context. Recent observations by the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter have characterized prevalent small-scale transient brightenings in the corona above the quiet Sun termed campfires. Aims. In this study we search for comparable brightenings in a numerical model and then investigate their relation to the magnetic field and the processes that drive these events. Methods. We use the MURaM code to solve the 3D radiation MHD equations in a box that stretches from the upper convection zone to the corona. The model self-consistently produces a supergranular network of the magnetic field and a hot corona above this quiet Sun. For the comparison with the model we synthesize the coronal emission as seen by EUI in its 174 {\AA} channel, isolate the seven strongest transient brightenings, and investigate (the changes of) the magnetic field in and around these in detail. Results. The transients we isolate have a lifetime of about 2 minutes and are elongated loop-like features with lengths around 1Mm to 4 Mm. They tend to occur at heights of about 2Mm to 5Mm above the photosphere a bit offset from magnetic concentrations that mark the bright chromospheric network and they reach temperatures of above 1 MK. With this they very much resemble the (larger) campfires found in observations. In our model most events are energised by component reconnection between (bundles of) field lines that interact at coronal heights. In one case we find that untwisting of a highly twisted flux rope initiates the heating. Conclusions. Based on our study we propose that the majority of campfire events found by EUI are driven by component reconnection and our model suggests that this process contributes significantly to the heating of the corona above the quiet Sun.

Bhuwan Joshi (USO/PRL, India), Prabir K. Mitra (USO/PRL, India, Dept. of Physics, Gujarat Univ., India), R. Bhattacharyya (USO/PRL, India), Kushagra Upadhyay (USO/PRL, India), Divya Oberoi (NCRA/TIFR, India), K. Sasikumar Raja (IIA, India), Christian Monstein (IRSOL, Switzerland)

30 pages, 14 figures, 1 table; accepted for publication in Solar Physics

We present a multi-wavelength investigation of a C-class flaring activity that occurred in the active region NOAA 12734 on 8 March 2019. The investigation utilises data from AIA and HMI on board the SDO and the Udaipur-CALLISTO solar radio spectrograph of the Physical Research Laboratory. This low intensity C1.3 event is characterised by typical features of a long duration event (LDE), viz. extended flare arcade, large-scale two-ribbon structures and twin coronal dimmings. The eruptive event occurred in a coronal sigmoid and displayed two distinct stages of energy release, manifested in terms of temporal and spatial evolution. The formation of twin dimming regions are consistent with the eruption of a large flux rope with footpoints lying in the western and eastern edges of the coronal sigmoid. The metric radio observations obtained from Udaipur-CALLISTO reveals a broad-band ($\approx$50-180 MHz), stationary plasma emission for $\approx$7 min during the second stage of the flaring activity that resemble a type IV radio burst. A type III decametre-hectometre radio bursts with starting frequency of $\approx$2.5 MHz precedes the stationary type IV burst observed by Udaipur-CALLISTO by $\approx$5 min. The synthesis of multi-wavelength observations and Non-Linear Force Free Field (NLFFF) coronal modelling together with magnetic decay index analysis suggests that the sigmoid flux rope underwent a zipping-like uprooting from its western to eastern footpoints in response to the overlying asymmetric magnetic field confinement. The asymmetrical eruption of the flux rope also accounts for the observed large-scale structures viz. apparent eastward shift of flare ribbons and post flare loops along the polarity inversion line (PIL), and provides an evidence for lateral progression of magnetic reconnection site as the eruption proceeds.

Elisa Frattin, I. Bertini, S.L.Ivanovski, F. Marzari, M. Fulle, V.V. Zakharov, F. Moreno, G. Naletto, M. Lazzarin, P. Cambianica, G. Cremonese, S. Ferrari, F. Ferri, C. Guettler, F. La Forgia, A. Lucchetti, M. Pajola, L. Penasa, A. Rotundi, H. Sierks, C. Tubiana

In this work we aim to characterise the dust motion in the inner coma of comet 67P/Churyumov- Gerasimenko to provide constraints for theoretical 3D coma models. The OSIRIS camera onboard the Rosetta mission was able for the first time to acquire images of single dust particles from inside the cometary coma, very close to the nucleus. We analyse a large number of particles, performing a significant statistic of their behaviour during the post perihelion period, when the spacecraft covered distances from the nucleus ranging between 80 and 400 km. We describe the particle trajectories, investigating their orientation and finding highly radial motion with respect to the nucleus. Then, from the particle brightness profiles, we derive a particle rotational frequency of v < 3.6 Hz, revealing that they are slow rotators and do not undergo fragmentation. We use scattering models to compare the observed spectral radiance of the particles with the simulated ones in order to estimate their size, finding values that range from millimetres up to centimetres. The statistics performed in this paper provide useful parameters to constrain the cometary coma dynamical models.

Kornél Kapás, Tamás Bozóki, Gergely Dálya, János Takátsy, László Mészáros, András Pál

Due to the advancement of nano-satellite technology, CubeSats and fleets of CubeSats can form an alternative to high-cost large-size satellite missions with the advantage of extended spatial coverage. One of these initiatives is the Cubesats Applied for MEasuring and LOcalising Transients (CAMELOT) mission concept, aimed at detecting and localizing gamma-ray bursts with an efficiency and accuracy comparable to large gamma-ray space observatories. While precise attitude control is not necessary for such a mission, attitude determination is an important issue in the interpretation of scintillator detector data as well as optimizing downlink telemetry. The employment of star trackers is not always a viable option for such small satellites, hence another alternative is necessary. A new method is proposed in this series of papers, utilizing thermal imaging sensors to provide simultaneous measurement of the attitude of the Sun and the horizon by employing a homogeneous array of such detectors. The combination with Sun and horizon detection w.r.t. the spacecraft would allow the full 3-DoF recovery of its attitude. In this paper we determine the spherical projection function of the MLX90640 infrasensors planned to be used for this purpose. We apply a polynomial transformation with radial corrections to map the spatial coordinates to the sensor plane. With the determined projection function the location of an infrared point source can be determined with an accuracy of ~40', well below the design goals of a nano-satellite designed for gamma-ray detection.

Enrico Peretti, Giovanni Morlino, Pasquale Blasi, Pierre Cristofari

11 pages, 5 figures, 2 tables, submitted to MNRAS, comments are welcome

The enhanced star forming activity, typical of starburst galaxies, powers strong galactic winds expanding on kpc scales and characterized by bubble structures. Here we discuss the possibility that particle acceleration may take place at the termination shock of such winds. We calculate the spectrum of such particles and their maximum energy, that turns out to be in the range $\sim 10-10^2$ PeV for typical values of the parameters. Cosmic rays accelerated at the termination shock can be advected towards the edge of the wind bubble and eventually escape into extragalactic space. We also calculate the flux of gamma rays and neutrinos produced by hadronic interactions in the bubble as well as the diffuse flux resulting from the superposition of the contribution of starburst galaxies on cosmological scales. Finally we compute the diffuse flux of cosmic rays from starburst bubbles and compare it with existing data.

Álvaro Pastor-Gutiérrez, Harm Schoorlemmer, Robert Daniel Parsons, Michael Schmelling

In the sub-TeV regime, the most widely used hadronic interaction models disagree significant lying their predictions for post-first interaction and ground-level particle spectra from cosmic ray induced air showers. These differences generate an important source of systematic uncertainty in their experimental use. We investigate the nature and impact of model uncertainties through a simultaneous analysis of ground level particles and first interaction scenarios. We focus on air shower primaries with energies close to the transition between high and low energy hadronic interaction models, where the dissimilarities have been shown to be the largest and well within the range of accelerator measurements. Interaction models are shown to diverge as several shower scenarios are compared, reflecting intrinsic differences in the model theoretical frameworks. Finally, we discuss the importance of interactions in the energy regime where the switching between models occurs (<1 TeV) and the effect of the choice of model on the number of hadronic interactions within cosmic ray induced air showers of higher energies.

A. Borghese, N. Rea, R. Turolla, M. Rigoselli, J.A.J. Alford, E.V. Gotthelf, M. Burgay, A. Possenti, S. Zane, F. Coti Zelati, R. Perna, P. Esposito, S. Mereghetti, D. Viganó, A. Tiengo, D. Götz, A. Ibrahim, G.L. Israel, J. Pons, R. Sathyaprakash

14 pages, 12 figures, 4 tables. Accepted for publications in MNRAS

After 15 years, in late 2018, the magnetar XTE J1810-197 underwent a second recorded X-ray outburst event and reactivated as a radio pulsar. We initiated an X-ray monitoring campaign to follow the timing and spectral evolution of the magnetar as its flux decays using Swift, XMM-Newton, NuSTAR, and NICER observations. During the year-long campaign, the magnetar reproduced similar behaviour to that found for the first outburst, with a factor of two change in its spin-down rate from $\sim7.2\times10^{-12}$ s s$^{-1}$ to $\sim1.5\times10^{-11}$ s s$^{-1}$ after two months. Unique to this outburst, we confirm the peculiar energy-dependent phase shift of the pulse profile. Following the initial outburst, the spectrum of XTE J1810-197 is well-modelled by multiple blackbody components corresponding to a pair of non-concentric, hot thermal caps surrounded by a cooler one, superposed to the colder star surface. We model the energy-dependent pulse profile to constrain the viewing and surface emission geometry and find that the overall geometry of XTE J1810-197 has likely evolved relative to that found for the 2003 event.

Haritma Gaur, Prashanth Mohan, Ashwani Pandey

13 pages, 1 figure, 4 tables; accepted for publication in ApJ

The flat spectrum radio quasars (FSRQs) are a sub-class of blazars characterised by prominent optical emission lines and a collimated large-scale jet along the observer line of sight. An X-ray spectral flattening has been reported in FSRQs (at relatively high redshifts), attributable to either absorption from gas along the line of sight or intrinsic jet based radiative processes. We study a sample of 16 high redshift FSRQs (z of 1.1 -- 4.7; rest frame energy upto 50 keV) observed with XMM-Newton and Swift satellites spanning 29 epochs. The X-ray spectra are fit with a power law including free excess absorption and one multiplied by an exponential roll off to account for the intrinsic jet based processes. A statistical analysis is used to distinguish between these models to understand the origin of the spectral flattening. The model selection is unable to distinguish between them in ten of the sixteen FSRQs. Intrinsic jet based radiative processes are indicated in four FSRQs where we infer energy breaks consistent with their expectation from the external Compton scattering of low energy ambient photons. Two of the FSRQs indicate mixed results, supportive of either scenario, illustrating the difficulty in identifying X-ray absorption signatures. A clear detection can be employed to disentangle the relative contributions from the inter-galactic medium and the intra-cluster medium, the methodology of which is outlined and applied to the latter two sources.

V. L. Oknyansky, M. S. Brotherton, S. S. Tsygankov, A. V. Dodin, D.-W. Bao, B.-X. Zhao, P. Du, M. A. Burlak, N. P. Ikonnikova, A. M. Tatarnikov, A. A. Belinski, A. A. Fedoteva, N. I. Shatsky, E. O. Mishin, S. G. Zheltouhov, S. A. Potanin, J.-M. Wang, J. N. McLane, H. A. Kobulnicky, D. A. Dale, T. E. Zastrocky, J. Maithil, K. A. Olson, C. Adelman, Z. Carter, A. M. Murphree, M. Oeur, S. Schonsberg, T. Roth

18 pages, 18 figures, accepted by the MNRAS

We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC 3516 carried out in 2018 to 2020 covering the wavelength range from the X-ray to the optical. The facilities included the telescopes of the CMO SAI MSU, the 2.3-m WIRO telescope, and the XRT and UVOT of Swift. We found that NGC 3516 brightened to a high state and could be classified as Sy1.5 during the late spring of 2020. We have measured time delays in the responses of the Balmer and He II 4686 lines to continuum variations. In the case of the best-characterized broad H-beta line, the delay to continuum variability is about 17 days in the blue wing and is clearly shorter, 9 days, in the red, which is suggestive of inflow. As the broad lines strengthened, the blue side came to dominate the Balmer lines, resulting in very asymmetric profiles with blueshifted peaks during this high state. During the outburst the X-ray flux reached its maximum on 1 April 2020 and it was the highest value ever observed for NGC 3516 by the Swift observatory. The X-ray hard photon index became softer, about 1.8 in the maximum on 21 Apr 2020 compared to the mean about 0.7 during earlier epochs before 2020. We have found that the UV and optical variations correlated well (with a small time delay of 1-2 days) with the X-ray until the beginning of April 2020, but later, until the end of Jun. 2020, these variations were not correlated. We suggest that this fact may be a consequence of partial obscuration by Compton-thick clouds crossing the line of sight.

Gavin P Lamb, D. Alexander Kann, Joseph John Fernández, Ilya Mandel, Andrew J. Levan, Nial R. Tanvir

11 pages, 5 figures, submitted to MNRAS

We investigate the shape of the jet break in within-beam gamma-ray burst (GRB) optical afterglows for various lateral jet structure profiles. We consider cases with and without lateral spreading and a range of inclinations within the jet core half-opening angle, $\theta_c$. We fit model and observed afterglow lightcurves with a smoothly-broken power-law function with a free-parameter $\kappa$ that describes the sharpness of the break. We find that the jet break is sharper ($\kappa$ is greater) when lateral spreading is included than in the absence of lateral spreading. For profiles with a sharp-edged core, the sharpness parameter has a broad range of $0.1\lesssim\kappa\lesssim4.6$, whereas profiles with a smooth-edged core have a narrower range of $0.1\lesssim\kappa\lesssim2.2$ when models both with and without lateral spreading are included. For sharp-edged jets, the jet break sharpness depends strongly on the inclination of the system within $\theta_c$, whereas for smooth-edged jets, $\kappa$ is more strongly dependent on the size of $\theta_c$. Using a sample of 20 GRBs we find nine candidate smooth-edged jet structures and eight candidate sharp-edged jet structures, while the remaining three are consistent with either. The shape of the jet break, as measured by the sharpness parameter $\kappa$, can be used as an initial check for the presence of lateral structure in within-beam GRBs where the afterglow is well-sampled at and around the jet-break time.

Thomas J. Maccarone (Texas Tech University)

9 pages, accepted to Journal of the American Association of Variable Star Observers

As new classes of transients and variable stars are discovered, and theoretical models are established to work or not to work for a few members of the class, it is often the case that some researchers will make arguments on the basis of Occam's razor that all members of the class must be produced by whichever mechanism first successfully explained one of the objects. It is also frequent that this assumption will be more more implicitly. Retrospective analysis shows rather clearly that this argument fails a large fraction of the time, and in many cases, this search for false consistency has led to more fundamental astrophysical errors, a few of which are quite prominent in the history of astronomy. A corollary of this is that on numerous occasions, theoretical models to explain transients have turned out to be models that describe real (but often not yet discovered) phenomenona other than the ones to which they have first been applied, albeit with minor errors that caused the model to appear to fit to a known phenomenon it did not describe. A set of examples of such events is presented here (some of which will be quite familiar to most astronomers), along with a discussion of why this phenomenon occurs, how it may be manifesting itself at the present time. Some discussion will also be made of why and when survey designs have led to immediate separation of various transient mechanisms, generally by being overpowered in some way relative to what is needed to {\it detect} a new class of objects.

Brad Koplitz, Jared Johnson, Benjamin F. Williams, Knox S. Long, William P. Blair, Jeremiah W. Murphy, Andrew Dolphin, Tristan Hillis

28 pages, 13 figures, 4 tables

We constrained the progenitor masses for 169 supernova remnants, 8 historically observed supernovae, and the black hole formation candidate in NGC 6946, finding that they are consistent with originating from a standard initial mass function. Additionally, there were 16 remnants that showed no sign of nearby star formation consistent with a core-collapse supernova, making them good Type Ia candidates. Using $Hubble$ $Space$ $Telescope$ broadband imaging, we measured stellar photometry of ACS/WFC fields in F435W, F555W, F606W, and F814W filters as well as WFC3/UVIS fields in F438W, F606W, and F814W. We then fitted this photometry with stellar evolutionary models to determine the ages of the young populations present at the positions of the SNRs and SNe. We then infer a progenitor mass probability distribution from the fitted age distribution. For 37 SNRs we tested how different filter combinations affected the inferred masses. We find that filters sensitive to H$\alpha$, [N II], and [S II] gas emission can bias mass estimates for remnants that rely on our technique. Using a KS-test analysis on our most reliable measurements, we find the progenitor mass distribution is well-matched by a power-law index of $-2.6^{+0.5}_{-0.6}$, which is consistent with a standard initial mass function.

Cyril Pitrou, Alain Coc, Jean-Philippe Uzan, Elisabeth Vangioni

2 pages

Nuclear physics experiments give reaction rates that, via modelling and comparison with primordial abundances, constrain cosmological parameters. The error bars of a key reaction, \dpg, were tightened in 2020, bringing to light discrepancies between different analyses and calling for more accurate measurements of other reactions.

W. N. Alston, C. Pinto, D. Barret, A. D'Ai, M. Del Santo, H. Earnshaw, A. C. Fabian, F. Fuerst, E. Kara, P. Kosec, M. J. Middleton, M. L. Parker, F. Pintore, A. Robba, T. Roberts, R. Sathyaprakash, D. Walton, E. Ambrosi

8 pages, 6 figures, submitted to MNRAS

Most ultraluminous X-ray sources (ULXs) are believed to be stellar mass black holes or neutron stars accreting beyond the Eddington limit. Determining the nature of the compact object and the accretion mode from broadband spectroscopy is currently a challenge, but the observed timing properties provide insight into the compact object and details of the geometry and accretion processes. Here we report a timing analysis for an 800 ks XMM-Newton campaign on the supersoft ultraluminous X-ray source, NGC 247 ULX-1. Deep and frequent dips occur in the X-ray light curve, with the amplitude increasing with increasing energy band. Power spectra and coherence analysis reveals the dipping preferentially occurs on $\sim 5$ ks and $\sim 10$ ks timescales, indicating the mechanism is located at $\sim 10^4$ $R_g$ from the compact object. The dips can be caused by either the occultation of the central X-ray source by an optically thick structure, such as warping of the accretion disc, or from obscuration by a wind launched from the accretion disc, or both. This behaviour supports the idea that supersoft ULXs are viewed close to edge-on to the accretion disc.

C. Pinto, R. Soria, D. Walton, A. D'Ai, F. Pintore, P. Kosec, W. N. Alston, F. Fuerst, M. J. Middleton, T. P. Roberts, M. Del Santo, D. Barret, E. Ambrosi, A. Robba, H. Earnshaw, A. Fabian

20 pages, 17 figures, submitted to MNRAS

Most ULXs are believed to be powered by super-Eddington accreting neutron stars and, perhaps, black holes. Above the Eddington rate the disc is expected to thicken and to launch powerful winds through radiation pressure. Winds have been recently discovered in several ULXs. However, it is yet unclear whether the thickening of the disc or the wind variability causes the switch between the classical soft and supersoft states observed in some ULXs. In order to understand such phenomenology and the overall super-Eddington mechanism, we undertook a large (800 ks) observing campaign with XMM-Newton to study NGC 247 ULX-1, which shifts between a supersoft and classical soft ULX state. The new observations show unambiguous evidence of a wind in the form of emission and absorption lines from highly-ionised ionic species, with the latter indicating a mildly-relativistic outflow (-0.17c) in line with the detections in other ULXs. Strong dipping activity is observed in the lightcurve and primarily during the brightest observations, which is typical among soft ULXs, and indicates a close relationship between the accretion rate and the appearance of the dips. The latter is likely due to a thickening of the disc scale-height and the wind as shown by a progressively increasing blueshift in the spectral lines.

Dongwoo T. Chung, Patrick C. Breysse, Håvard Tveit Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Kieran A. Cleary, Hans Kristian Eriksen, Marie Kristine Foss, Joshua Ott Gundersen, Laura C. Keating, Jonas Gahr Sturtzel Lunde, Nils-Ole Stutzer, Marco P. Viero, Duncan J. Watts, Ingunn Kathrine Wehus

24 pages + appendix and bibliography (29 pages total), 13 figures, 2 tables; submitted to ApJ

Line-intensity mapping observations will find fluctuations of integrated line emission are attenuated by varying degrees at small scales due to the width of the line emission profiles. This attenuation may significantly impact estimates of astrophysical or cosmological quantities derived from measurements. We consider a theoretical treatment of the effect of line broadening on both the clustering and shot-noise components of the power spectrum of a generic line-intensity power spectrum using a halo model. We then consider possible simplifications to allow easier application in analysis, particularly in the context of inferences that require numerous, repeated, fast computations of model line-intensity signals across a large parameter space. For the CO Mapping Array Project (COMAP) and the CO(1-0) line-intensity field at $z\sim3$ serving as our primary case study, we expect a $\sim10\%$ attenuation of the spherically averaged power spectrum on average at relevant scales of $k\approx0.2$-$0.3$ Mpc$^{-1}$, compared to $\sim25\%$ for the interferometric Millimetre-wave Intensity Mapping Experiment (mmIME) targeting shot noise from CO lines at $z\sim1$-$5$ at scales of $k\gtrsim1$ Mpc$^{-1}$. We also consider the nature and amplitude of errors introduced by simplified treatments of line broadening, and find that while an approximation using a single effective velocity scale is sufficient for spherically-averaged power spectra, a more careful treatment is necessary when considering other statistics such as higher multipoles of the anisotropic power spectrum or the voxel intensity distribution.

Eva Sextl, Rolf-Peter Kudritzki, Jochen Weller, Miguel A. Urbaneja, Achim Weiss

To be published in the Astrophysical Journal

We calculate models of stellar evolution for very massive stars and include the effects of modified gravity to investigate the influence on the physical properties of blue supergiant stars and their use as extragalactic distance indicators. With shielding and fifth force parameters in a similar range as in previous studies of Cepheid and tip of the red giant branch (TRGB) stars we find clear effects on stellar luminosity and flux-weighted gravity. The relationship between flux weighted gravity, g_F = g/Teff^4, and bolometric magnitude M_bol (FGLR), which has been used successfully for accurate distance determinations, is systematically affected. While the stellar evolution FGLRs show a systematic offset from the observed relation, we can use the differential shifts between models with Newtonian and modified gravity to estimate the influence on FGLR distance determinations. Modified gravity leads to a distance increase of 0.05 to 0.15 magnitudes in distance modulus. These change are comparable to the ones found for Cepheid stars. We compare observed FGLR and TRGB distances of nine galaxies to constrain the free parameters of modified gravity. Not accounting for systematic differences between TRGB and FGLR distances shielding parameters of 5*10^-7 and 10^-6 and fifth force parameters of 1/3 and 1 can be ruled out with about 90% confidence. Allowing for potential systematic offsets between TRGB and FGLR distances no determination is possible for a shielding parameter of 10^-6. For 5*10$^-7 a fifth force parameter of 1 can be ruled out to 92% but 1/3 is unlikely only to 60%.

Maxime Devogèle, Marin Ferrais, Emmanuel Jehin, Nicholas Moskovitz, Brian A. Skiff, Stephen E. Levine, Annika Gustafsson, Davide Farnocchia, Marco Micheli, Colin Snodgrass, Galin Borisov, Jean Manfroid, Youssef Moulane, Zouhair Benkhaldoun, Artem Burdanov, Francisco J. Pozuelos, Michael Gillon, Julien de Wit, Simon F. Green, Philippe Bendjoya, Jean-Pierre Rivet, Luy Abe, David Vernet, Colin Orion Chandler, Chadwick A. Trujillo

15 pages, 11 figures; Accepted in MNRAS

In December 2018, the main-belt asteroid (6478)~Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observations over two apparitions to monitor its activity, rotation period, composition, and possible non-gravitational orbital evolution. We find that Gault has a rotation period of $P = 2.4929 \pm 0.0003$ hours with a lightcurve amplitude of $0.06$ magnitude. This short rotation period close to the spin barrier limit is consistent with Gault having a density no smaller than $\rho = 1.85$~g/cm$^3$ and its activity being triggered by the YORP spin-up mechanism. Analysis of the Gault phase curve over phase angles ranging from $0.4^{\circ}$ to $23.6^{\circ}$ provides an absolute magnitude of $H = 14.81 \pm 0.04$, $G1=0.25 \pm 0.07$, and $G2= 0.38 \pm 0.04$. Model fits to the phase curve find the surface regolith grain size constrained between 100-500 $\rm{\mu}$m. Using relations between the phase curve and albedo we determine that the geometrical albedo of Gault is $p_{\rm v} = 0.26 \pm 0.05$ corresponding to an equivalent diameter of $D = 2.8^{+0.4}_{-0.2}$ km. Our spectroscopic observations are all consistent with an ordinary chondrite-like composition (S, or Q-type in the Bus-DeMeo taxonomic classification). A search through archival photographic plate surveys found previously unidentified detections of Gault dating back to 1957 and 1958. Only the latter had been digitized, which we measured to nearly double the observation arc of Gault. Finally, we did not find any signal of activity during the 2020 apparition or non-gravitational effects on its orbit.

R. Modzelewska, G.A. Bazilevskaya, M. Boezio, S.V. Koldashov, M.B. Krainev, N. Marcelli, A.G. Mayorov, M.A. Mayorova, R. Munini, I.K. Troitskaya, R.F. Yulbarisov, X. Luo, M.S. Potgieter, O.P.M. Aslam

13 pages, 9 figures

Using measurements from the PAMELA and ARINA spectrometers onboard the RESURS DK-1 satellite, we have examined the 27-day intensity variations in galactic cosmic ray (GCR) proton fluxes in 2007-2008. The PAMELA and ARINA data allow for the first time a study of time profiles and the rigidity dependence of the 27-day variations observed directly in space in a wide rigidity range from ~300 MV to several GV. We find that the rigidity dependence of the amplitude of the 27-day GCR variations cannot be described by the same power-law at both low and high energies. A flat interval occurs at rigidity R = <0.6-1.0> GV with a power-law index gamma = - 0.13+/-0.44 for PAMELA, whereas for R >= 1 GV the power-law dependence is evident with index gamma = - 0.51+/-0.11. We describe the rigidity dependence of the 27-day GCR variations for PAMELA and ARINA data in the framework of the modulation potential concept using the force-field approximation for GCR transport. For a physical interpretation, we have considered the relationship between the 27-day GCR variations and solar wind plasma and other heliospheric parameters. Moreover, we have discussed possible implications of MHD modeling of the solar wind plasma together with a stochastic GCR transport model concerning the effects of corotating interaction regions.

David S. Doelman, Joost P. Wardenier, Peter Tuthill, Michael P. Fitzgerald, Jim Lyke, Steph Sallum, Barnaby Norris, N. Zane Warriner, Christoph Keller, Michael J. Escuti, Frans Snik

17 pages, 21 figures, 7 tables, accepted for publication in A&A

We report on the design, construction, and commissioning of a prototype aperture masking technology implemented at the Keck OSIRIS Imager: the holographic aperture mask. Holographic aperture masking (HAM) aims at (i) increasing the throughput of sparse aperture masking (SAM) by selectively combining all subapertures across a telescope pupil in multiple interferograms using a phase mask, and (ii) adding low-resolution spectroscopic capabilities. Using liquid-crystal geometric phase patterns, we manufacture a HAM mask that uses an 11-hole SAM design as the central component and a holographic component comprising 19 different subapertures. Thanks to a multilayer liquid-crystal implementation, the mask has a diffraction efficiency higher than 96% from 1.1 to 2.5 micron. We create a pipeline that extracts monochromatic closure phases from the central component as well as multiwavelength closure phases from the holographic component. We test the performance of the HAM mask in the laboratory and on-sky. The holographic component yields 26 closure phases with spectral resolutions between R$\sim$6.5 and R$\sim$15. On April 19, 2019, we observed the binary star HDS 1507 in the Hbb filter ($\lambda_0 = 1638$ nm and $\Delta \lambda = 330$ nm) and retrieved a constant separation of 120.9 $\pm 0.5$ mas for the independent wavelength bins, which is in excellent agreement with literature values. For both the laboratory measurements and the observations of unresolved reference stars, we recorded nonzero closure phases -- a potential source of systematic error that we traced to polarization leakage of the HAM optic. We propose a future upgrade that improves the performance, reducing this effect to an acceptable level. Holographic aperture masking is a simple upgrade of SAM with increased throughput and a new capability of simultaneous low-resolution spectroscopy that provides new differential observables.

William Lake, Smadar Naoz, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Kyle Kremer

17 pages, 8 figures

Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where these relative velocities are coherent. However, it is challenging to study SIGOs using simulations on large scales due to the varying relative velocities at scales larger than a few Mpc. Here, we study SIGO abundances semi-analytically: using perturbation theory, we predict the number density of SIGOs analytically, and compare these results to small-box numerical simulations. We use the agreement between the numerical and analytic calculations to extrapolate the large-scale variation of SIGO abundances over different stream velocities. As a result, we predict similar large-scale variations of objects with high gas densities before reionization that could possibly be observed by JWST. If indeed SIGOs are progenitors of globular clusters, then we expect a similar variation of globular cluster abundances over large scales. Significantly, we find that the expected number density of SIGOs is consistent with observed globular cluster number densities. As a proof-of-concept, and because globular clusters were proposed to be natural formation sites for gravitational wave sources from binary black hole (BBH) mergers, we show that SIGOs should imprint an anisotropy on the gravitational wave signal on the sky, consistent with SIGOs' distribution.

A. N. Lasenby, W. J. Handley, D. J. Bartlett, C.S. Negreanu

28 pages, 32 figures, prepared for submission to PRD

The concordance model of cosmology predicts a universe which finishes in a finite amount of conformal time at a future conformal boundary. We show that for particular cases we study, the background variables and perturbations may be analytically continued beyond this boundary and that the "end of the universe" is not necessarily the end of their physical development. Remarkably, these theoretical considerations of the end of the universe might have observable consequences today: perturbation modes consistent with these boundary conditions have a quantised power spectrum which may be relevant to features seen in the large scale cosmic microwave background. Mathematically these cosmological models may either be interpreted as a palindromic universe mirrored in time, a reflecting boundary condition, or a double cover, but are identical with respect to their observational predictions and stand in contrast to the predictions of conformal cyclic cosmologies.

Modified gravity (MG) theories predict, in general, that the ratio of gravitational wave (GW) to electromagnetic (EM) luminosity distances, $\Xi$, differs from its general relativity (GR) value of unity at cosmological scales, thus providing another perturbative probe to MG. In this paper, we introduce new phenomenological parametrizations for both the Friedmann-Lema\^itre-Robertson-Walker (FLRW) background evolution of $f(R)$ models, via the dark energy equation of state parameter, $w_{\text{DE}}$, and for $\Xi$ in this class of theories. We simulate a mock dataset for the Einstein Telescope (ET) of 1000 GW signals from binary neutron star (BNS) mergers and redshift information from their EM counterpart, exploring the consequent constraints on the relevant gravitational, cosmological and phenomenological parameters. As a model of particular interest, we take $\gamma$-gravity theory and investigate whether it could be distinguished from GR's $\Lambda$CDM model. We then combine our results with actual data from type Ia supernovae (SNIa) and combined baryon acoustic oscillations (BAO) and cosmic microwave background (CMB) observations. Additionally, we also investigate the potential bounds to $f_{R0}$ for any viable $f(R)$ whose background evolution is indistinguishable from the standard model of cosmology above a certain redshift, showing that, for a $\Lambda$CDM fiducial model, ET data would provide $|f_{R0}|<10^{-2}$ at a 95\% level. We conclude altogether that probing the redshift evolution of the GW luminosity distance from detections of the ET in its first running decade will not substantially help constraining $f(R)$ theories of gravity.

In environments with high dense neutrino gases, such as in core-collapse supernovae (CCSNe), the neutrinos can experience collective neutrino oscillation due to their self-interactions. In particular, fast flavor conversion driven by the crossings in the neutrino angular distribution can affect explosion mechanism, nucleosynthesis, and neutrino observation. We perform the numerical computation of nonlinear flavor evolution on the neutrino angular distribution with tiny crossings expected to be generated in the preshock region. We demonstrate that the fast instability is triggered and a cascade develops under a realistic three-flavor model considering muon production and weak magnetism in the SN dynamics. The tiny crossing excites specific spatial modes, and then the flavor instability propagates into other modes which otherwise remain stable due to the nonlinear effects. Our results indicate that fast flavor conversion can rise in the preshock region and have a sufficient impact on the flavor contents.

Ibrahima Bah, Iosif Bena, Pierre Heidmann, Yixuan Li, Daniel R. Mayerson

26 pages + Appendices, 8 figures

We construct a family of non-supersymmetric extremal black holes and their horizonless microstate geometries in four dimensions. The black holes can have finite angular momentum and an arbitrary charge-to-mass ratio, unlike their supersymmetric cousins. These features make them and their microstate geometries astrophysically relevant. Thus, they provide interesting prototypes to study deviations from Kerr solutions caused by new horizon-scale physics. In this paper, we compute the gravitational multipole structure of these solutions and compare them to Kerr black holes. The multipoles of the black hole differ significantly from Kerr as they depend non-trivially on the charge-to-mass ratio. The horizonless microstate geometries have the same multipoles as their corresponding black hole, with small deviations set by the scale of their microstructure.

Magnetic reconnection is a primary driver of particle acceleration processes in space and astrophysical plasmas. Understanding how particles are accelerated and the resulting particle energy spectra is among the central topics in reconnection studies. We review recent advances in addressing this problem in nonrelativistic reconnection that is relevant to space and solar plasmas and beyond. We focus on particle acceleration mechanisms, particle transport due to 3D reconnection physics, and their roles in forming power-law particle energy spectra. We conclude by pointing out the challenges in studying particle acceleration and transport in a large-scale reconnection layer and the relevant issues to be addressed in the future.

Salvatore Capozziello, Shin'ichi Nojiri, Sergei D. Odintsov

20 pages

Scalar modes are considered for the propagation of gravitational waves in cosmological background taking into account of thermal radiation and near cosmological singularities. In particular, we point out that the contribution of thermal radiation can heavily affect the dynamics of gravitational waves giving enhancement or dissipation effects both at quantum and classical level. These effects are considered both in General Relativity and in modified theories like $F(R)$ gravity. The possible detection and the disentanglement of the various gravitational modes on the stochastic background are also discussed.

Sergey Yu. Vernov, Ekaterina O. Pozdeeva

13 pages, 2 figures

De Sitter solutions play an important role in cosmology because the knowledge of unstable de Sitter solutions can be useful to describe inflation, whereas stable de Sitter solutions are often used in models of late-time acceleration of the Universe. The Einstein-Gauss-Bonnet gravity cosmological models are actively used both as inflationary models and as dark energy models. To modify the Einstein equations one can add a nonlinear function of the Gauss-Bonnet term or a function of the scalar field multiplied on the Gauss-Bonnet term. The effective potential method essentially simplifies the search and stability analysis of de Sitter solutions, because the stable de Sitter solutions correspond to minima of the effective potential.

The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient $\beta_{21}$ with a favorable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($\beta^u_{21}$) and magnetic ($\beta^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $\beta^b_{21}$, the kinetic SCE $\beta^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal{O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $\beta^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $\beta^u_{21}$ can dominate even at energy equipartition $u_\text{rms}\simeq b_\text{rms}$, resulting in a positive total $\beta_{21}$ even though $\beta^b_{21}<0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation (SOCA) versus the spectral-$\tau$ closure (STC), for which opposite signs for $\beta^u_{21}$ have been reported, with the same sign for $\beta^b_{21}<0$. The results also offer an explanation for the simulations that find $\beta^u_{21}>0$ and an inconclusive overall sign of $\beta_{21}$ for $\mathcal{O}(10)$ Reynolds numbers. The transient behavior of $\beta^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work.