Abstract visibility
Abstract text size

Papers for Monday, Jun 27 2022

Papers with local authors

A. K. Saydjari, E. F. Schlafly, D. Lang, A. M. Meisner, G. M. Green, C. Zucker, I. Zelko, J. S. Speagle, T. Daylan, A. Lee, F. Valdes, D. Schlegel, D. P. Finkbeiner

40 pages, 33 figures, submitted to ApJ

0 votes
Paper 6 — arXiv:2206.11909
0 votes
Paper 6 — arXiv:2206.11909

Deep optical and near-infrared imaging of the entire Galactic plane is essential for understanding our Galaxy's stars, gas, and dust. The second data release of the DECam Plane Survey (DECaPS2) extends the five-band optical and near-infrared survey of the southern Galactic plane to cover $6.5\%$ of the sky, |b| < 10{\deg} and 6{\deg} > l > -124{\deg}, complementary to coverage by Pan-STARRS1. Typical single-exposure effective depths, including crowding effects and other complications, are 23.5, 22.6, 22.1, 21.6, and 20.8 mag in $g$, $r$, $i$, $z$, and $Y$ bands, respectively, with around 1 arcsecond seeing. The survey comprises 3.32 billion objects built from 34 billion detections in 21.4 thousand exposures, totaling 260 hours open shutter time on the Dark Energy Camera (DECam) at Cerro Tololo. The data reduction pipeline features several improvements, including the addition of synthetic source injection tests to validate photometric solutions across the entire survey footprint. A convenient functional form for the detection bias in the faint limit was derived and leveraged to characterize the photometric pipeline performance. A new post-processing technique was applied to every detection to de-bias and improve uncertainty estimates of the flux in the presence of structured backgrounds, specifically targeting nebulosity. The images and source catalogs are publicly available at this http URL

Wenting Wang, Ling Zhu, Zhaozhou Li, Yang Chen, Jiaxin Han, Feihong He, Xiaohu Yang, Yipeng Jing, Carlos Frenk, Jialu Nie, Hao Tian, Chao Liu, Yanan Cao, Xiaoqing Qiu, John Helly, Robert J. J. Grand, Facundo A. Gomez

submitted to ApJ

0 votes
Paper 25 — arXiv:2206.12121
0 votes
Paper 25 — arXiv:2206.12121

Mock member stars for 28 dwarf galaxies are constructed from the cosmological Auriga simulation, which reflect the true dynamical status of realistic tracers. The axis-symmetric Jeans Anisotropic Multi-Gaussian Expansion (JAM) modeling is applied to 6,000 star particles for each, to recover the underlying matter distribution. The stellar or dark matter component individually is poorly recovered, but the total profile is constrained reasonably. The mass within the half-mass radius of tracers is recovered the tightest, and the mass between 200 and 300~pc, $M(200-300\mathrm{pc})$, is constrained ensemble unbiasedly, with a scatter of 0.169~dex. Using a smaller sample of 2,000 particles and only line-of-sight velocities with errors, the scatter in $M(200-300\mathrm{pc})$ is increased by $\sim$40%. Quiescent Sagittarius dSph-like systems and star-forming systems with outflows show distinct features, with $M(200-300\mathrm{pc})$ mostly under-estimated for the former, and more likely over-estimated for the latter. The biases correlate with the dynamical status, which is a result of contractions due to tidal effects in quiescent systems or prominent galactic winds in star-forming systems, driving the systems out of equilibrium. After including Gaia DR3 proper motion errors, we find proper motions are as useful as line-of-sight velocities for nearby systems at $<\sim$60~kpc. The expected proper motion errors of the China Space Station Telescope is useful at $<\sim$20~kpc. Extrapolating the density profiles and the dynamical constraints down to scales below the resolution, we find the mass within 150~pc can be constrained ensemble unbiasedly, with a scatter of $\sim$0.3~dex. In the end, we show that the contraction of member stars in nearby systems is detectable with Gaia DR3 proper motion errors.

All other papers

Arindam Mal, Sarbani Palit, Sisir Roy

8 pages,2 figure

Recent observational evidence in extra galactic astronomy, the interpretation of the nature of quasar redshift continues to be research interest. Spectrum observation of high redshift quasar is young in nature. Observational evidence discuss on physical interpretation of redshift periodicity with statistical confirmation. Karlsson observed redshift periodicity at integer multiples of 0.089 in log scale and Burbidge observed redshift periodicity integer multiple of 0.061 in linear scale .Data analysis is important in order to form correct interpretations of the observed phenomena. Since Singular value decomposition (SVD) based periodicity estimation is known to be superior for noisy data sets, especially when the data contains multiple harmonics and overtones, mainly irregular in nature, we have chosen it to be our primary tool for analysis of the quasar-galaxy pair redshift data. Kernel density estimation has been performed for estimating the bin width as proper computation of this quantity is crucial for the correctness of the analysis and prevention of over smoothing of the data.We observed fundamental periodicity to be an integer multiple of 0.063 and 0.0604 using method1 and method2 in the transformed quasar redshift data with 95% confidence interval in linear scale. Our results clearly establish that redshift is quantized for quasar-galaxy pair data and its histogram exhibits periodic peak(s). At last briefly discussed on physical interpretation of quantized redshift for quasar and galaxy.Hoyle Narlikar theory of gravity explain the Mystery in recent observation.

R. Baier-Soto, R. Herrera-Camus, N. M. Förster Schreiber, A. Contursi, R. Genzel, D. Lutz, L. Tacconi

Accepted for publication in Astronomy & Astrophysics Letters, 8 pages, 6 figures

GN-108036 is a star-forming galaxy at $z=7.21$, and one of the most distant known sources in the Northern hemisphere. Based on observations from the NOrthern Extended Millimeter Array (NOEMA), here we report the tentative detection of the [CII] line at $\approx4\sigma$ significance. The integrated [CII] line emission is spatially offset about $\sim4$ kpc from the rest-frame ultraviolet (UV) emission. The total [CII] luminosity ($L_{\rm [CII]}=2.7\times10^8~L_{\odot}$) is consistent with the relation between [CII] luminosity and star formation rate (SFR) observed in nearby and high-$z$ star forming galaxies. More interestingly, the [CII] line is blueshifted with respect to the Ly$\alpha$ line by $980\pm10$ km s$^{-1}$. If confirmed, this corresponds to the largest velocity offset reported to date between the Ly$\alpha$ line and a non-resonant line at $z\gtrsim6$. According to trends observed in other high redshift galaxies, the large Ly$\alpha$ velocity offset in GN-108036 is consistent with its low Ly$\alpha$ equivalent width and high UV absolute magnitude. Based on Ly$\alpha$ radiative transfer models of expanding shells, the large Ly$\alpha$ velocity offset in GN-108036 could be interpreted as the presence of a large column density of hydrogen gas, and/or an outflow with a velocity of $v_{\rm out}\sim\Delta v_{\rm Ly \alpha}/2\sim500$ km s$^{-1}$. We also report the 3$\sigma$ detection of a potential galaxy companion located $\sim30$ kpc east of GN-108036, at a similar systemic velocity, and with no counterpart rest-frame UV emission.

Václav Pavlík, Enrico Vesperini

9 pages, 13 figures, accepted for publication in MNRAS

This paper is the third in a series investigating, by means of $N$-body simulations, the implications of an initial radially anisotropic velocity distribution on the dynamics of star clusters. Such a velocity distribution may be imprinted during a cluster's early evolutionary stages and several observational studies have found examples of old globular clusters in which radial anisotropy is still present in the current velocity distribution. Here we focus on its influence on mass segregation and the dynamics of primordial binary stars (disruptions, ejections, and component exchanges). The larger fraction of stars on radial/highly eccentric orbits in the outer regions of anisotropic clusters lead to an enhancement in the dynamical interactions between inner and outer stars that affects both the process of mass segregation and the evolution of primordial binaries. The results of our simulations show that the time scale of mass segregation of the initially anisotropic cluster is longer in the core and shorter in the outer regions, when compared to the initially isotropic system. The evolution of primordial binaries is also significantly affected by the initial velocity distribution and we find that the rate of disruptions, ejections and exchange events affecting the primordial binaries in the anisotropic clusters is higher than in the isotropic ones.

Akos Bogdan, Mark Vogelsberger

Accepted chapter of the Section on Galaxies for the Handbook of X-ray and Gamma-ray Astrophysics (Springer, eds. Banbi & Santangelo)

The presence of gaseous X-ray halos around massive galaxies is a basic prediction of all past and modern structure formation simulations. The importance of these X-ray halos is further emphasized by the fact that they retain signatures of the physical processes that shape the evolution of galaxies from the highest redshift to the present day. In this review, we overview our current observational and theoretical understanding of hot gaseous X-ray halos around nearby massive galaxies and we also describe the prospects of observing X-ray halos with future instruments.

Ananthan Karunakaran, Kristine Spekkens, Rhys Carroll, David J. Sand, Paul Bennet, Denija Crnojević, Michael G. Jones, Burçin Mutlu-Pakdil

12 pages, 5 Figures, 3 Tables; Submitted to MNRAS; Comments welcome!

We present neutral atomic hydrogen (HI) observations using the Robert C. Byrd Green Bank Telescope (GBT) along the lines of sight to 49 dwarf satellite galaxy candidates around eight Local Volume systems (M104, M51, NGC1023, NGC1156, NGC2903, NGC4258, NGC4565, NGC4631). We detect the HI reservoirs of two candidates (dw0934+2204 and dw1238$-$1122) and confirm them as background sources relative to their nearest foreground host systems. The remaining 47 satellite candidates are not detected in HI, and we place stringent $5\sigma$ upper limits on their HI mass. We note that some (15/47) of our non-detections stem from satellites being occluded by their putative host's HI emission. In addition to these new observations, we compile literature estimates on the HI mass for an additional 17 satellites. We compare the HI properties of these satellites to those within the Local Group, finding broad agreement between them. Crucially, these observations probe a ``transition'' region between $-10\gtrsim M_V \gtrsim -14$ where we see a mixture of gas-rich and gas-poor satellites and where quenching processes shift from longer timescales (i.e. via starvation) to shorter ones (i.e. via stripping). While there are many gas-poor satellites within this region, some are gas rich and this suggests that the transition towards predominantly gas-rich satellites occurs at $L_{V}\sim10^{7}L_{\odot}$, in line with simulations. The observations presented here are a key step toward characterizing the properties of dwarf satellite galaxies around Local Volume systems and future wide-field radio surveys with higher angular resolution (e.g.~WALLABY) will vastly improve upon the study of such systems.

Andrew J. Winter, Thomas J. Haworth

Community feedback requested for this solicited review, submitted to EPJ+ - 47 pages, 17 figures

Planet-forming disc evolution is not independent of the star formation and feedback process in giant molecular clouds. In particular, OB stars emit UV radiation that heats and disperses discs in a process called 'external photoevaporation'. This process is understood to be the dominant environmental influence acting on planet-forming discs in typical star forming regions. Our best studied discs are nearby, in sparse stellar groups where external photoevaporation is less effective. However the majority of discs are expected to reside in much stronger UV environments. Understanding external photoevaporation is therefore key to understanding how most discs evolve, and hence how most planets form. Here we review our theoretical and observational understanding of external photoevaporation. We also lay out key developments for the future to address existing unknowns and establish the full role of external photoevaporation in the disc evolution and planet formation process.

Dillon Dong (Caltech), Gregg Hallinan (Caltech)

31 pages, 11 figures; submitted to ApJ; comments very much welcome!

We report the discovery and follow-up observations of VT 1137-0337: an unusual radio transient found in our systematic search for extragalactic explosions in the VLA Sky Survey (VLASS). VT 1137-0337 is located in the brightest region of a dwarf starburst galaxy (stellar mass $\sim 10^{8.3} M_{\odot}$, star formation rate $\sim 0.5 M_{\odot}$ yr$^{-1}$) at a luminosity distance of 121.6 Mpc. Its 3 GHz luminosity of $\sim 2.5 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$ is comparable to luminous radio supernovae associated with dense circumstellar interaction and relativistic outflows. However, its broadband radio spectrum - a featureless power law $\propto \nu^{-0.35 \pm 0.02}$ over a range of $\gtrsim$10$\times$ in frequency and fading at a rate of $\sim$ 5% per year over 4 years - cannot be directly explained by the shock of a stellar explosion. Jets launched by various classes of accreting black holes also struggle to account for VT 1137-0337's combination of observational properties. Instead, we propose that VT 1137-0337 is a $\sim$decades old pulsar wind nebula that has recently emerged from within the free-free opacity of its surrounding supernova ejecta. If the nebula is powered by spindown, the central neutron star should be highly magnetized, with a surface dipole field of $\sim 10^{13} - 10^{14}$ G and a present-day spin period of $\sim 10 - 100$ ms. Alternatively, the nebula may be powered by the release of magnetic energy from a magnetar. Magnetar nebulae have been proposed to explain the persistent radio sources associated with the repeating fast radio bursts FRB 121102 and FRB 190520B. These FRB persistent sources have not previously been observed as transients, but do bear a striking resemblance to VT 1137-0337 in their radio luminosity, spectral index, and host galaxy properties.

Corentin Cadiou, Andrew Pontzen, Hiranya V. Peiris

8 pages, 3 figures. Submitted to MNRAS. Comments welcome

The angular momentum of galaxies controls the kinematics of their stars, which in turn drives observable quantities such as the apparent radius, the bulge fraction, and the alignment with other nearby structures. To show how angular momentum of galaxies is determined, we build high (${35}\,\mathrm{pc}$) resolution numerical experiments in which we increase or decrease the angular momentum of the Lagrangian patches in the early universe. We simulate three galaxies over their histories from $z=200$ to $z=2$, each with five different choices for the angular momentum (fifteen simulations in total). Our results show that altering early-universe angular momentum changes the timing and orbital parameters of mergers, which in turn changes the total stellar angular momentum within a galaxy's virial radius in a predictable manner. Of our three galaxies, one has no large satellite at $z=2$; in this case, the specific angular momentum is concentrated in the central galaxy. We modify its stellar angular momentum over $0.7\,\mathrm{dex}$ (from $61$ to $320\,\mathrm{kpc.km.s^{-1}}$) and show that this causes its effective radius to grow by $40\,\%$, its $v/\sigma$ parameter to grow by a factor $\times 2.6$ and its bulge fraction to decrease from $0.72$ to $0.57$. The ability to control angular momentum will allow future studies to probe the causal origin of scaling relations between galaxy mass, angular momentum and morphology, and to better understand the origin of galactic intrinsic alignments.

Moritz Reichert, Martin Obergaulinger, Miguel-Angel Aloy, Michael Gabler, Almudena Arcones, Friedrich-Karl Thielemann

Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: i) a prompt ejection of matter right after core bounce, ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, iii) small amount of mass ejected with high entropies in the center of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as $^{26}$Al, $^{44}$Ti, $^{56}$Ni, and $^{60}$Fe. The obtained $^{56}$Ni masses vary between $0.01 - 1\,\mathrm{M_\odot}$. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few $\sim 10^{43} \,\mathrm{erg\,s^{-1}}$. Under certain conditions, we find a significant impact of the $^{66}$Ni decay chain that can raise the peak luminosity up to $\sim 38\%$ compared to models including only the $^{56}$Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements.

Lodovico Coccato (1), Amelia Fraser-McKelvie (2, 3), Yara L. Jaffé (4), Evelyn J. Johnston (5), Arianna Cortesi (6, 7), Diego Pallero (4) ((1) European Southern Observatory, Karl-Schwarzchild-str., 2, 85748 Garching b. M unchen, Germany, (2) International Centre for Radio Astronomy Research, The University of Western Australia, 35 Stirling Hwy, 6009 Crawley, WA, Australia, (3) ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), (4) Instituto de Fisica y Astronomia, Universidad de Valparaiso, Avda. Gran Bretaña 1111 Valparaiso, Chile, (5) Núcleo de Astronomía de la Facultad de Ingeniería y Ciencias, Universidad Diego Portales, Av. Ejército Libertador 441, Santiago, Chile, (6) Observatório do Valongo, Ladeira do Pedro Antônio 43, CEP:20080-090, Rio de Janeiro, RJ, Brazil, (7) Centro Brasileiro de Pesquisas Físicas, Rua Dr. Xavier Sigaud 150, CEP 22290-180, Rio de Janeiro, RJ, Brazil.)

12 pages, 9 figures, accepted for publication in MNRAS

It is well established that there are at least two main channels to form lenticular (or S0) galaxies. The first, which we name "faded spiral" scenario, includes quenching events that led to consumption or removal of gas from a spiral progenitor. The second, which we call "merger" scenario, includes merger-like events and interactions between galaxies. Each scenario leaves characteristic signatures in the newly-formed lenticular galaxy. However, the conditions that trigger one mechanism over another are still unknown. This paper is the third of a series aimed at understanding the role of the environment in the formation of lenticular galaxies. In this study, we combine the kinematics, morphology, and properties of the stellar populations of 329 S0s from the SAMI and MaNGA surveys in order to highlight the role of the environment in the process.We divide the S0s into two classes (A and B) according to their global properties, that we can associate to the products of a faded spiral scenario (class A) or a merger scenario (class B). We then study how the various classes are distributed within different environments. Our study reveals that the "faded spiral" pathway is the most efficient channel to produce S0s, and it becomes more efficient as the mass of the group or cluster or local density of galaxies increase. The merger pathway is also a viable channel, and its efficiency becomes higher with decreasing local density or environment mass.

Trevor Z. Dorn-Wallenstein, Emily M. Levesque, James R. A. Davenport, Kathryn F. Neugent, Brett M. Morris, K. Azalee Bostroem

Accepted to the Astrophysical Journal. Comments welcome

Fast yellow pulsating supergiants (FYPS) are a recently-discovered class of evolved massive pulsator. As candidate post-red supergiant objects, and one of the few classes of pulsating evolved massive stars, these objects have incredible potential to change our understanding of the structure and evolution of massive stars. Here we examine the lightcurves of a sample of 126 cool supergiants in the Magellanic Clouds observed by the Transiting Exoplanet Survey Satellite (TESS) in order to identify pulsating stars. After making quality cuts and filtering out contaminant objects, we examine the distribution of pulsating stars in the Hertzprung-Russel (HR) diagram, and find that FYPS occupy a region above $\log L/L_\odot \gtrsim 5.0$. This luminosity boundary corresponds to stars with initial masses of $\sim$18-20 $M_\odot$, consistent with the most massive red supergiant progenitors of supernovae (SNe) II-P, as well as the observed properties of SNe IIb progenitors. This threshold is in agreement with the picture that FYPS are post-RSG stars. Finally, we characterize the behavior of FYPS pulsations as a function of their location in the HR diagram. We find low frequency pulsations at higher effective temperatures, higher frequency pulsations at lower temperatures, with a transition between the two behaviors at intermediate temperatures. The observed properties of FYPS make them fascinating objects for future theoretical study.

Alex Gough, Cora Uhlemann

Submitted to the Open Journal of Astrophysics, 18 pages + appendices, 15 figures, comments welcome

The classical dynamics of collisionless cold dark matter, commonly described by fluid variables or a phase-space distribution, can be captured in a single semiclassical wavefunction. We illustrate how classical multi-streaming creates wave interference in a toy model corresponding to the dynamics of the Zel'dovich approximation and link it to diffraction optics. Wave interference dresses the classical skeleton of cold dark matter with universal features akin to the physical imprints of wavelike (or fuzzy) dark matter. We untangle this wave interference to obtain single-stream wavefunctions corresponding to the classical fluid streams, by writing the wavefunction in an integral form. Our wave decomposition captures the full phase-space information and isolates the multi-stream phenomena related to vorticity and velocity dispersion. We link the wave interference features of our system to the standard forms of diffraction catastrophe integrals, which produce bright caustics in optical fields analogous to the cold dark matter density field. Our two complementary descriptions of dark matter wave-fields present rich universal features that can unlock new ways of modelling and probing wavelike dark matter on the scales of the cosmic web.

Lewis Prole, Paul Clark, Ralf Klessen, Simon Glover, Ruediger Pakmor

Submitted to MNRAS

Population III stars form in groups due to the fragmentation of primordial gas. While uniform magnetic fields have been shown to support against fragmentation in present day star formation, it is unclear whether realistic k^3/2 primordial fields can have the same effect. We bypass the issues associated with simulating the turbulent dynamo by introducing a saturated magnetic field at equipartition with the velocity field when the central densities reaches 10-13 g cm-3. We test a range of sink particle creation densities from 10-10-10-8 g cm-3. Within the range tested, the fields did not suppress fragmentation of the gas and hence could not prevent the degree of fragmentation from increasing with increased resolution. The number of sink particles formed and total mass in sink particles was unaffected by the magnetic field across all seed fields and resolutions. The magnetic pressure remained sub-dominant to the gas pressure except in the highest density regions of the simulation box, where it became equal to but never exceeded gas pressure. Our results suggest that the inclusion of magnetic fields in numerical simulations of Pop III star formation is largely unimportant.

Krishan Chand (ARIES), Gopal-Krishna (CEBS)

Accepted for publication in MNRAS Letter

Flat-spectrum radio quasars (FSRQs) whose brightness is dominated by a relativistically beamed core, are frequently found in 'blazar state' commonly inferred from a high optical polarization ( > 3%), and/or a large continuum variability. Here we use these two prime optical markers to investigate continuance of an FSRQ in blazar (or non-blazar) state over an exceptionally long time baseline spanning 4 decades. Our basic sample is a well-defined, unbiased set of 80 FSRQs whose blazar state stood confirmed during 1980s from optical polarimetry. Four decades later, blazar state of each FSRQ is ascertained here from variability of their optical light-curves of typical duration $\sim$ 3.5 years, a low noise (rms $\sim$ 2%) and good cadence ($\sim$ 3 days), obtained under the Zwicky Transient Facility project ongoing since 2018. For about 40% of these FSRQs, blazar state could be ascertained additionally from the opto-polarimetric survey RoboPol (2013-2017). From both these databases it is found that only $\sim$ 10% of the FSRQs have undergone a blazar $\leftrightarrow$ non-blazar state transition over the past 3 - 4 decades. This reinforces the case for a long-term stability of blazar state in individual FSRQs, despite their state fluctuating more commonly on year-like time scales.

Aishwarya R. Iyer, Michael R. Line, Philip S. Muirhead, Jonathan J. Fortney, Ehsan Gharib-Nezhad

Submitted to ApJ

About 70-80% of stars in our solar and galactic neighborhood are M dwarfs. They span a range of low masses and temperatures relative to solar-type stars, facilitating molecule formation throughout their atmospheres. Standard stellar atmosphere models primarily designed for FGK stars face challenges when characterizing broadband molecular features in spectra of cool stars. Here, we introduce SPHINX--a new 1-D self-consistent radiative-convective thermochemical equilibrium chemistry model grid of atmospheres and spectra for M dwarfs in low-resolution (R~250). We incorporate the latest pre-computed absorption cross-sections with pressure-broadening for key molecules dominant in late-K, early/main-sequence-M stars. We then validate our grid models by acquiring fundamental properties (Teff, log(g), [M/H], radius, and C/O) for 10 benchmark M+G binary stars with known host metallicities and 10 M dwarfs with interferometrically measured angular diameters. Incorporating a Gaussian-process inference tool Starfish, we account for correlated and systematic noise in low-resolution (spectral stitching of SpeX, SNIFS, and STIS) observations and derive robust estimates of fundamental M dwarf atmospheric parameters. Additionally, we assess the influence of photospheric heterogeneity on acquired [M/H] and find that it could explain some deviations from observations. We also probe whether the model-assumed convective mixing-length parameter influences inferred radii, effective temperature, and [M/H] and again find that may explain discrepancies between interferometry observations and model-derived stellar parameters for cooler M dwarfs. Mainly, we show the unique strength in leveraging broadband molecular absorption features occurring in low-resolution M dwarf spectra and demonstrate the ability to improve constraints on fundamental properties of exoplanet hosts and late brown dwarf companions.

X-ray binaries exhibit a soft spectral state comprising thermal blackbody emission at 1 keV and a power-law tail above 10 keV. Empirical models fit the high-energy power-law tail to radiation from a nonthermal electron distribution, but the physical location of the nonthermal electrons and the reason for their power-law index and high-energy cut-off are still largely unknown. Here, we propose that the nonthermal electrons originate from within the black hole's innermost stable circular orbit (the ''plunging region''). Using an analytic model for the plunging region dynamics and electron distribution function properties from particle-in-cell simulations, we outline a steady-state model that can reproduce the observed spectral features. In particular, our model reproduces photon indices of $\Gamma\gtrsim2$ and power-law luminosities on the order of a few percent of the disk luminosity for strong magnetic fields, consistent with observations of the soft state. Because the emission originates so close to the black hole, we predict that the power-law luminosity should strongly depend on the system inclination angle and black hole spin. This model could be extended to the power-law tails observed above 400 keV in the hard state of X-ray binaries.

Water-rich exoplanet is a type of terrestrial planet that is water-rich and its ocean depth can reach tens of to hundreds of kilo-meters with no exposed continents. Due to the lack of exposed continents, neither western boundary current nor coastal upwelling exists, and ocean overturning circulation becomes the most important way to return the nutrients deposited in deep ocean back to the thermocline and to the surface ocean.Here we investigate the depth of the thermocline in both wind-dominated and mixing-dominated systems on water-rich exoplanets using the global ocean model MITgcm. We find that the wind-driven circulation is dominated by overturning cells through Ekman pumping and subduction and by zonal (west--east) circum-longitudinal currents, similar to the Antarctic Circumpolar Current on Earth. The wind-influenced thermocline depth shows little dependence on the ocean depth, and under a large range of parameters, the thermocline is restricted within the upper layers of the ocean. The mixing-influenced thermocline is limited within the upper 10 km of the ocean and can not reach the bottom of the ocean even under extremely strong vertical mixing. The scaling theories for the thermocline depth on Earth are applicable for the thermocline depth on water-rich exoplanets. However, due to the lack of exposed continents, the zonal and meridional flow speeds are not in the same magnitude as that in the oceans of Earth, which results in the scaling relationships for water-rich exoplanets are a little different from that used on Earth.

F. H. Vincent, S. Gralla, A. Lupsasca, M. Wielgus

21 pages, 9 figure, submitted to A&A

High-frequency very-long-baseline interferometry (VLBI) observations can now resolve the horizon-scale emission from sources in the immediate vicinity of nearby supermassive black holes. Future space-VLBI observations will access highly lensed features of black hole images -- photon rings -- that will provide particularly sharp probes of strong-field gravity. Focusing on the particular case of the supermassive black hole M87*, our goal is to explore a wide variety of accretion flows onto a Kerr black hole and to understand their corresponding images and visibilities. We are particularly interested in the visibility on baselines to space, which encodes the photon ring shape and whose measurement could provide a stringent test of the Kerr hypothesis. We develop a fully analytical model of stationary, axisymmetric accretion flows with a variable disk thickness and a matter four-velocity that can smoothly interpolate between purely azimuthal rotation and purely radial infall. We then determine the observational appearance of such flows, taking care to include the effects of thermal synchrotron emission and absorption. Our images generically display a "wedding cake" structure composed of discrete, narrow photon rings (n=1,2,...) stacked on top of broader primary emission that surrounds a central brightness depression of model-dependent size. We find that the "black hole shadow" is a model-dependent phenomenon -- even for diffuse, optically thin sources -- and should not be regarded as a generic prediction of general relativity. At 230 GHz, the n=1 ring is always visible, but the n=2 ring is sometimes suppressed due to absorption. At 345 GHz, the medium is optically thinner and the n=2 ring displays clear signatures in both the image and visibility domains, identifying this frequency as more promising for future space-VLBI measurements of the photon ring shape.

Chengyuan Li, Zhenya Zheng, Xiaodong Li, Xiaoying Pang, Baitian Tang, Antonino P. Milone, Yue Wang, Haifeng Wang, Dengkai Jiang

21 pages, 16 figures, RAA accepted

Multiple stellar populations (MPs) in most star clusters older than 2 Gyr, as seen by lots of spectroscopic and photometric studies, have led to a significant challenge to the traditional view of star formation. In this field, space-based instruments, in particular the Hubble Space Telescope (HST), have made a breakthrough as they significantly improved the efficiency of detecting MPs in crowding stellar fields by images. The China Space Station Telescope (CSST) and the HST are sensitive to a similar wavelength interval, but it covers a field of view which is about 5-8 times wider than that of HST. One of its instruments, the Multi-Channel Imager (MCI), will have multiple filters covering a wide wavelength range from NUV to NIR, making the CSST a potentially powerful tool for studying MPs in clusters. In this work, we evaluate the efficiency of the designed filters for the MCI/CSST in revealing MPs in different color-magnitude diagrams (CMDs). We find that CMDs made with MCI/CSST photometry in appropriate UV filters are powerful tools to disentangle stellar populations with different abundances of He, C, N, O and Mg. On the contrary, the traditional CMDs are blind to multiple populations in globular clusters (GCs). We show that CSST has the potential of being the spearhead instrument for investigating MPs in GCs in the next decades.

M. Dixon (1), C. Lidman (2, 3), J. Mould (1), L. Kelsey (4, 5), D. Brout (6), A. Möller (1), P. Wiseman (5), M. Sullivan (5), L. Galbany (7, 8), T. M. Davis (9), M. Vincenzi (4, 5), D. Scolnic (10), G. F. Lewis (11), M. Smith (5), R. Kessler (12, 13), A. Duffy (1), E. Taylor (1), C. Flynn (1), T. M. C. Abbott (14), M. Aguena (15), S. Allam (16), F. Andrade-Oliveir (17), J. Annis (16), J. Asorey (18), E. Bertin (19, 20), S. Bocquet (21), D. Brooks (22), D. L. Burke (23, 24), A. Carnero Rosell (15, 25, 26), D. Carollo (27), M. Carrasco Kind (28, 29), J. Carretero (30), M. Costanzi (27, 31, 32), L. N. da Costa (15), M. E. S. Pereira (33), P. Doel (22), S. Everett (34), I. Ferrero (35), B. Flaugher (16), D. Friedel (28), J. Frieman (13, 16), J. García-Bellido (36), M. Gatti (37), D. W. Gerdes (17, 38), et al. (34 additional authors not shown)

15 pages, 10 figures. Submitting to MNRAS

We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features (< 1.3{\sigma}), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [OII] (4.4{\sigma}) and the Balmer emission lines (3{\sigma}). These trends are weaker than the well known trend between Hubble residuals and host galaxy stellar mass (7.2{\sigma}) that is derived from broad band photometry. After light curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3{\sigma}) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H\b{eta} and H{\gamma}. The trend is only present in the redder SNe Ia, suggesting that bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters.

Nayana A.J., Poonam Chandra, Anoop Krishna, G.C. Anupama

18 pages, 8 figures, accepted for publication in ApJ

We present extensive radio monitoring of a type IIb supernova (SN IIb), SN 2016gkg during $t \sim$ 8$-$1429 days post explosion at frequencies $\nu \sim$ 0.33$-$25 GHz. The detailed radio light curves and spectra are broadly consistent with self-absorbed synchrotron emission due to the interaction of the SN shock with the circumstellar medium. The model underpredicts the flux densities at $t \sim$ 299 days post-explosion by a factor of 2, possibly indicating a density enhancement in the CSM due to a non-uniform mass-loss from the progenitor. Assuming a wind velocity $v_{\rm w} \sim$ 200 km s$^{-1}$, we estimate the mass-loss rate to be $\dot{M} \sim$ (2.2, 3.6, 3.8, 12.6, 3.7, and 5.0) $\times$ 10$^{-6}$ $M_{\odot}$ yr$^{-1}$ during $\sim$ 8, 15, 25, 48, 87, and 115 years, respectively before the explosion. The shock wave from SN 2016gkg is expanding from $R \sim$ 0.5 $\times$ 10$^{16}$ to 7 $\times$ 10$^{16}$ cm during $t \sim$ 24$-$492 days post-explosion indicating a shock deceleration index, $m$ $\sim$ 0.8 ($R \propto t^m$), and mean shock velocity $v \sim$ 0.1c. The radio data being inconsistent with free-free absorption model and higher shock velocities are in support of a relatively compact progenitor for SN 2016gkg.

S.V. Logvinenko, S.A. Tyul'bashev, V.M. Malofeev

Published in Bulletin of the Lebedev Physics Institute and translated from russian version of paper by Yandex translator with correction of scientific lexis. 5 page, 2 figure

Two new rotating transients were detected in the 2020 observations carried out on the LPA LPI radio telescope. The dispersion measures of the found transients are DM=21 and 35 pc/cm^3, the pulse half-widths are We=18 and 35 ms for J1550+09 and J2047+13, respectively. The upper estimate of the period RRAT J2047+13 P=2.925s was obtained. The study shows the existence of rotating transients whose pulses appear less frequently than one pulse per 10 hours of observations.

O. S. Bayandina, C. L. Brogan, R. A. Burns, A. Caratti o Garatti, J. O. Chibueze, S. P. van den Heever, S. E. Kurtz, G. C. MacLeod, L. Moscadelli, A. M. Sobolev, K. Sugiyama, I. E. Val'tts, Y. Yonekura

The massive young stellar object (MYSO) G358.93-0.03-MM1 showed an extraordinary near-infrared- to (sub-)millimetre-dark and far-infrared-loud accretion burst, which is closely associated with flares of several class II methanol maser transitions, and, later, a 22 GHz water maser flare. Water maser flares provide an invaluable insight into ejection events associated with accretion bursts. Although the short timescale of the 22 GHz water maser flare made it impossible to carry out a very long baseline interferometry observation, we could track it with the Karl G. Jansky Very Large Array (VLA). The evolution of the spatial structure of the 22 GHz water masers and their association with the continuum sources in the region is studied with the VLA during two epochs, pre- and post-H2O maser flare. A drastic change in the distribution of the water masers is revealed: in contrast to the four maser groups detected during epoch I, only two newly formed clusters are detected during epoch II. The 22 GHz water masers associated with the bursting source MM1 changed in morphology and emission velocity extent. Clear evidence of the influence of the accretion burst on the ejection from G358.93-0.03-MM1 is presented. The accretion event has also potentially affected a region with a radius of ~2'' (~13 500 AU at 6.75 kpc), suppressing water masers associated with other point sources in this region.

Daria Teplykh, Valery Malofeev, Oleg Malov, Sergey Tyul'bashev

published in Open Astronomy, 6 pages, 8 figures

The analysis of radio emission of three new pulsars discovered at the Pushchino Radio Astronomy Observatory is presented. The detailed observations were carried out at a frequency of 111 MHz using the Large Phased Array (LPA) and the standard digital receiver with a total bandwidth is 2.245 MHz and time resolution is 2.46 or 5.12 ms. All pulsars exhibit features of their radiation, the subpulse drift is observed in J0220+3622, the flare activity is exhibited in J0303+2248 and the nulling phenomenon has been detected in J0810+3725.

M. Spinelli, G. Kyriakou, G. Bernardi, P. Bolli, L.J. Greenhill, A. Fialkov, H. Garsden

Accepted for publication in MNRAS

Cosmic Dawn, the onset of star formation in the early universe, can in principle be studied via the 21cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be {\it O}(10-100)\,mK. Detection requires separation of the 21cm signal from bright chromatic foreground emission due to Galactic structure, and the characterisation of how it couples to instrumental response. In this work, we present characterisation of antenna gain patterns for the Large-aperture Experiment to detect the Dark Ages (LEDA) via simulations, assessing the effects of the antenna ground-plane geometries used, and measured soil properties. We then investigate the impact of beam pattern uncertainties on the reconstruction of a Gaussian absorption feature. Assuming the pattern is known and correcting for the chromaticity of the instrument, the foregrounds can be modelled with a log-polynomial, and the 21cm signal identified with high accuracy. However, uncertainties on the soil properties lead to \textperthousand\ changes in the chromaticity that can bias the signal recovery. The bias can be up to a factor of two in amplitude and up to few \% in the frequency location. These effects do not appear to be mitigated by larger ground planes, conversely gain patterns with larger ground planes exhibit more complex frequency structure, significantly compromising the parameter reconstruction. Our results, consistent with findings from other antenna design studies, emphasise the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments.

M. Brož, P. Harmanec, P. Zasche, R. Catalan-Hurtado, B.N. Barlow, W. Frondorf, M. Wolf, H. Drechsel, R. Chini, A. Nasseri, J. Labadie-Bartz, G.W. Christie, W.S.G. Walker, M. Blackford, D. Blane, A.A. Henden, T. Bohlsen, H. Božić, J. Jonák

A&A, submitted

HD 93206 is early-type massive stellar system, composed of components resolved by direct imaging (Ab, Ad, B, C, D) as well as a compact sub-system (Aa1, Aa2, Ac1, Ac2). Its geometry was already determined on the basis of extensive photometric, spectroscopic and interferometric observations. However, the fundamental absolute parameters are still not known precisely enough. We use an advanced N-body model to account for all mutual gravitational perturbations among the four close components, and all observational data types, including: astrometry, radial velocities, eclipse timing variations, squared visibilities, closure phases, triple products, normalized spectra, and spectral-energy distribution (SED). The respective model has 38 free parameters, namely three sets of orbital elements, component masses, and their basic radiative properties ($T$, $\log g$, $v_{\rm rot}$). We revised the fundamental parameters of QZ Car as follows. For a model with the nominal extinction coefficient $R_V \equiv A_V/E(B-V) = 3.1$, the best-fit masses are $m_1 = 26.1\,M_{\rm S}$, $m_2 = 32.3\,M_{\rm S}$, $m_3 = 70.3\,M_{\rm S}$, $m_4 = 8.8\,M_{\rm S}$, with uncertainties of the order of $2\,M_{\rm S}$, and the system distance $d = (2800\pm 100)\,{\rm pc}$. In an alternative model, where we increased the weights of RV and TTV observations and relaxed the SED constraints, because extinction can be anomalous with $R_V \sim 3.4$, the distance is smaller, $d = (2450\pm 100)\,{\rm pc}$. This would correspond to that of Collinder 228 cluster. Independently, this is confirmed by dereddening of the SED, which is only then consistent with the early-type classification (O9.7Ib for Aa1, O8III for Ac1). Future modelling should also account for an accretion disk around Ac2 component.

X. Wang, H. Wang, H. J. Mo

15 pages, 12 figures, submitted to A&A

We use a large N-body simulation to study the characteristic scales in the density gradient profiles in and around halos with masses ranging from $10^{12}$ to $10^{15} h^{-1}{\rm M_\odot}$. We investigate the profiles separately along the major (T_1) and minor (T_3) axes of the local tidal tensor and how the characteristic scales depend on halo mass, formation time, and environment. We find two kinds of prominent characteristic features in the gradient profiles, a deep `valley' and a prominent `peak'. We use the Gaussian Process Regression to fit the gradient profiles and identify the local extrema to determine the scales associate with these features. Around the valley, we identify three types of distinct local minima, corresponding to caustics of particles orbiting around halos. The appearance and depth of the three caustics depend significantly on the direction defined by the local tidal field, formation time and environment of halos. The first caustic is located at a radius r>0.8R_{200}, corresponding to the splashback feature, and is dominated by particles at their first apocenter after infall. The second and third caustics, around 0.6R_{200} and 0.4R_{200} respectively, can be determined reliably only for old halos. The first caustic is always the most prominent feature along T_3, but may not be the case along T_1 or in azimuthally-averaged profiles, suggesting that caution must be taken when using averaged profiles to investigate the splashback radius. We find that the splashback feature is approximately isotropic when proper separations are made between the first and the other caustics. We also identify a peak feature located at $\sim$ 2.5R_{200} in the density gradient profile. This feature is the most prominent along T_1 and is produced by mass accumulations from the structure outside halos. We also discuss the origins of these features and their observational implications.

S. B. Shaulov, V. A. Ryabov, S. E. Pyatovsky, A. L. Shepetov, V. V. Zhukov

Submission to SciPost Phys. Proc

The conclusions of this work are based on the analysis of the characteristics of the EAS cores obtained using X-ray emulsion chambers. According to these data, a number of anomalous effects are observed in the knee region, such as scaling violation in the spectra of secondary hadrons, an excess of muons in EAS with gamma families and others. At the same energies equivalent to 1-100 PeV the laboratory system colliders show scaling behavior. So analysis of the data on the EAS cores suggests that the knee in their spectrum is formed by a component of cosmic rays of a non-nuclear nature, possibly consisting of stable (quasi-stable) particles of hypothetical strange quark matter.

Zhihong He, Kun Wang, Yangping Luo, Jing Li, Xiaochen Liu, Qingquan Jiang

16 figures, 2 tables (full version available at CDS), 886 figures in figure set, accepted for publication on ApJs

Although previous searches for star clusters have been very successful, many clusters are likely still omitted, especially at high Galactic latitude regions. In this work, based on the astrometry of Gaia EDR3, we searched nearby (parallax > 0.8 mas) all-sky regions, obtaining 886 star clusters, of which 270 candidates have not been cataloged before. At the same time, we have presented the physical parameters of the clusters by fitting theoretical isochrones to their optical magnitudes. More halo members and expanding structures in many star clusters were also found. Most of the new objects are young clusters that are less than 100 million years old. Our work greatly increased the sample size and physical parameters of star clusters in the solar neighborhood, in particular, 46 clusters are newly found with |b| > 20 deg, which represents an increase of nearly three fold of cluster numbers at high Galactic latitude regions. The cluster parameters and member stars are available at CDS via https://cdsarc.u-strasbg.fr/ftp/vizier.submit//hezh22b/, and the cluster figure sets are available via https://doi.org/10.12149/101133.

Gaia Collaboration, L. Galluccio, M. Delbo, F. De Angeli, T. Pauwels, P. Tanga, F. Mignard, A. Cellino, A. G. A. Brown, K. Muinonen, A. Penttila, S. Jordan, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren, X. Luri, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran, N. A. Walton, C. A. L. Bailer-Jones, U. Bastian, R. Drimmel, F. Jansen, D. Katz, M. G. Lattanzi, F. van Leeuwen, J. Bakker, C. Cacciari, J. Castaneda, C. Fabricius, M. Fouesneau, Y. Frémat, A. Guerrier, U. Heiter, E. Masana, R. Messineo, N. Mowlavi, C. Nicolas, K. Nienartowicz, F. Pailler, P. Panuzzo, F. Riclet, W. Roux, G. M. Seabroke, R. Sordo, F. Thévenin, et al. (385 additional authors not shown)

30 pages, 26 figures

The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was derived from measurements obtained by means of the Blue and Red photometers (BP/RP), which were binned in 16 discrete wavelength bands. We describe the processing of the Gaia spectral data of SSOs, explaining both the criteria used to select the subset of asteroid spectra published in Gaia DR3, and the different steps of our internal validation procedures. In order to further assess the quality of Gaia SSO reflectance spectra, we carried out external validation against SSO reflectance spectra obtained from ground-based and space-borne telescopes and available in the literature. For each selected SSO, an epoch reflectance was computed by dividing the calibrated spectrum observed by the BP/RP at each transit on the focal plane by the mean spectrum of a solar analogue. The latter was obtained by averaging the Gaia spectral measurements of a selected sample of stars known to have very similar spectra to that of the Sun. Finally, a mean of the epoch reflectance spectra was calculated in 16 spectral bands for each SSO. The agreement between Gaia mean reflectance spectra and those available in the literature is good for bright SSOs, regardless of their taxonomic spectral class. We identify an increase in the spectral slope of S-type SSOs with increasing phase angle. Moreover, we show that the spectral slope increases and the depth of the 1 um absorption band decreases for increasing ages of S-type asteroid families.

C. Sánchez Contreras, J. Alcolea, R. Rodríguez Cardoso, V. Bujarrabal, A. Castro-Carrizo, G. Quintana-Lacaci, L. Velilla-Prieto, M. Santander-García

30 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics (on 19/06/2022)

We present ALMA continuum and molecular line emission maps at $\sim$1 mm of OH 231.8, a well studied bipolar nebula around an AGB star. The excellent angular resolution of our maps ($\sim$20 mas) allows us to scrutinise the central nebular regions of OH 231.8, which hold the clues to unravel how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind, the discovery of a rotating circumbinary disk selectively traced by NaCl, KCl, and H$_2$O emission lines. This represents the first detection of KCl in an oxygen-rich AGB circumstellar envelope. The rotating disk, of radius $\sim$30 au, lies at the base of a young bipolar wind traced by SiO and SiS emission, which also presents signs of rotation at its base. The NaCl equatorial structure is characterised by a mean rotation velocity of $\sim$4 km s$^{-1}$ and extremely low expansion speeds, $\sim$3 km s$^{-1}$. The outflow has a predominantly expansive kinematics characterised by a constant radial velocity gradient of $\sim$65 km s$^{-1}$ arcsec$^{-1}$ at its base. Beyond $r$$\sim$350 au, the gas in the outflow continues radially flowing at a constant terminal speed of $\sim$16 km s$^{-1}$. Our continuum maps reveal a spatially resolved dust disk-like structure perpendicular to the outflow, with the NaCl, KCl and H$_2$O emission arising from the disk's surface layers. Within the disk, we also identify an unresolved point continuum source, which likely represents the central Mira-type star QX Pup enshrouded by a $\sim$3 $R_{\star}$ component of hot ($\sim$1400 K) freshly formed dust. The point source is slightly off-centered from the disk centroid, enabling us for the first time to place constraints to the orbital separation and period of the central binary system, $a$$\sim$20 au and $P_{\rm orb}$$\sim$55 yr, respectively. (abridged).

L. Ricci, B. Boccardi, E. Nokhrina, M. Perucho, N. MacDonald, G. Mattia, P. Grandi, E. Madika, T. P. Krichbaum, J. A. Zensus

12 pages plus appendix, 8 figures. Accepted for publication in A&A

Aims. Hot accretion flows are thought to be able to power the relativistic jets observed in Active Galactic Nuclei. They can present themselves as SANE (Standard And Normal Evolution) disks or MAD (Magnetically Arrested Disks), two states implying profound differences in the physical properties of the disks themselves and of the outflows they produce. Methods. In this paper we use a multi-frequency and multi-epoch data set to study the giant radio galaxy NGC 315, with the goal to explore the properties of its accretion disk and sub-parsec jet. We analyze the source maps with a pixel-based analysis and we use theoretical models to link the observational properties of the jet to the physical state of the accretion disk. Results. We propose that the bulk flow in NGC 315 accelerates on sub-pc scales, concurrently with the parabolic expansion. We show that this fast acceleration can be theoretically reconciled with a magnetically driven acceleration. Along the acceleration and collimation zone, we observe an unexpected spectral behavior, with very steep spectral index values $\alpha \sim -1.5$ ($S_\nu \propto \nu^\alpha$) between 22 GHz and 43 GHz. Based on the properties of this region, we predict the black hole of NGC 315 to be fast rotating and the magnetic flux threading the accretion disk to be in excellent agreement with that expected in the case of a MAD. Using a new formalism based on the core-shift effect, we model the magnetic field downstream a quasi-parabolic accelerating jet and we reconstruct it up to the event horizon radius. In the MAD scenario, we compare it with the expected magnetic saturation strengths in the disk, finding a good agreement.

Chloe Fisher, Kevin Heng

14 pages, 8 figures. Accepted for publication in ApJ

The construction and implementation of atmospheric model grids is a popular tool in exoplanet characterisation. These typically vary a number of parameters linearly, containing one model for every combination of parameter values. Here we investigate alternative methods of sampling parameters, including random sampling and Latin hypercube (LH) sampling, and how these compare to linearly sampled grids. We use a random forest to analyse the performance of these grids for two different models, as well as investigate the information content of the particular model grid from Goyal et al. 2019. We also use nested-sampling to implement mock atmospheric retrievals on simulated JWST transmission spectra by interpolating on linearly sampled model grids. Our results show that random or LH sampling out-performs linear sampling in parameter predictability for our higher dimensional models, requiring fewer models in the grid, and thus allowing for more computationally intensive forward models to be used. We also find that using a traditional retrieval with interpolation on a linear grid can produce biased posterior distributions, especially for parameters with non-linear effects on the spectrum. In particular, we advise caution when performing linear interpolation on the C/O ratio, cloud properties, and metallicity. Finally, we find that the information content analysis of the grid from Goyal et al. 2019 is able to highlight key areas of the spectra where the presence or absence of certain molecules can be detected, providing good indicators for parameters such as temperature and C/O ratio.

Jack White, Karel Adamek, Jayanta Roy, Sofia Dimoudi, Scott M. Ransom, Wesley Armour

The Fourier Domain Acceleration Search (FDAS) is an effective technique for detecting faint binary pulsars in large radio astronomy datasets. This paper quantifies the sensitivity impact of reducing numerical precision in the GPU accelerated FDAS pipeline of the AstroAccelerate software package. The prior implementation used IEEE-754 single-precision in the entire binary pulsar detection pipeline, spending a large fraction of the runtime computing GPU accelerated FFTs. AstroAccelerate has been modified to use bfloat16 (and IEEE754 double-precision to provide a "gold standard" comparison) within the Fourier domain convolution section of the FDAS routine. Approximately 20,000 synthetic pulsar filterbank files representing binary pulsars were generated using SIGPROC with a range of physical parameters. They have been processed using bfloat16, single and double-precision convolutions. All bfloat16 peaks are within 3% of the predicted signal-to-noise ratio of their corresponding single-precision peaks. Of 14,971 "bright" single-precision fundamental peaks above a power of 44.982 (our experimentally measured highest noise value), 14,602 (97.53%) have a peak in the same acceleration and frequency bin in the bfloat16 output plane, whilst in the remaining 369 the nearest peak is located in the adjacent acceleration bin. There is no bin drift measured between the single and double-precision results. The bfloat16 version of FDAS achieves a speedup of approximately 1.6x compared to single-precision. A comparison between AstroAccelerate and the PRESTO software package is presented using observations collected with the GMRT of PSR J1544+4937, a 2.16ms black widow pulsar in a 2.8 hour compact orbit.

Tomohiro Fujita, Kyohei Mukaida, Yuichiro Tada

15pages, 4 figures

We develop the stochastic formalism for $\mathrm{U}(1)$ gauge fields that has the Chern-Simons coupling to a rolling pseudo-scalar field during inflation. The Langevin equations for the physical electromagnetic fields are derived and the analytic solutions are studied. Using numerical simulation we demonstrate that the electromagnetic fields averaged over the Hubble scale continuously change their direction and their amplitudes fluctuate around the analytically obtained expectation values. Though the isotropy is spontaneously broken by picking up a particular local Hubble patch, each Hubble patch is understood independent and the isotropy is conserved globally by averaging all the Hubble patches.

H. Zhao, M. Schultheis, T. Zwitter, C.A.L. Bailer-Jones, P. Panuzzo, P. Sartoretti, G.M. Seabroke, A. Recio-Blanco, P. de Laverny, G. Kordopatis, O.L. Creevey, T.E. Dharmawardena, Y. Frémat, R. Sordo, R. Drimmel, D.J. Marshall, P.A. Palicio, G. Contursi, M.A. Álvarez, S. Baker, K. Benson, M. Cropper, C. Dolding, H.E. Huckle, M. Smith, O. Marchal, C. Ordenovic, F. Pailler, I. Slezak

7 pages, 4 figures, submitted to A&A

Studies on the correlation between different diffuse interstellar bands (DIBs) are important for exploring their origins. However, the Gaia-RVS spectral window between 846 and 870 nm has few DIBs, the strong DIB at 862 nm being the only convincingly confirmed one. Here we attempt to confirm the existence of a broad DIB around 864.8 nm and estimate its characteristics using the stacked Gaia-RVS spectra of a large number of stars. We study the correlations between the two DIBs at 862 nm and 864.8 nm, as well as the interstellar extinction. We obtain spectra of the interstellar medium (ISM) absorption by subtracting the stellar components using templates constructed from real spectra at high Galactic latitudes with low extinctions. We then stack the ISM spectra in Galactic coordinates, pixelized by the HEALPix scheme, to measure the DIBs. The stacked spectrum is modeled by the profiles of the two DIBs, Gaussian for $\lambda$862 and Lorentzian for $\lambda$864.8, and a linear continuum. We obtain 8458 stacked spectra in total, of which 1103 (13%) have reliable fitting results after applying numerous conservative filters. This work is the first one to fit and measure $\lambda$862 and $\lambda$864.8 simultaneously in cool-star spectra. We find that the EWs and CDs of the two DIBs are well-correlated with each other. The full width at half maximum (FWHM) of $\lambda$864.8 is estimated as $1.62 \pm 0.33$ nm which compares to $0.55 \pm 0.06$ nm for $\lambda$862. We also measure the vacuum rest-frame wavelength of $\lambda$864.8 to be $\lambda_0 = 864.53 \pm 0.14$ nm, smaller than previous estimates. We make a solid confirmation of the existence of the DIB around 864.8 nm by exploring its correlation with $\lambda$862 and estimating its FWHM. $\lambda$862 correlates better with E(BP-RP) than $\lambda$864.8.

Roland Diehl, Andreas Korn, Bruno Leibundgut, Maria Lugaro, Anton Wallner

44 pages, 55 figures, 2 tables. Invited review. Accepted for publication in 'Progress in Particle and Nuclear Physics'

The origins of the elements and isotopes of cosmic material is a critical aspect of understanding the evolution of the universe. Nucleosynthesis typically requires physical conditions of high temperatures and densities. These are found in the Big Bang, in the interiors of stars, and in explosions with their compressional shocks and high neutrino and neutron fluxes. Many different tools are available to disentangle the composition of cosmic matter, in material of extraterrestrial origins such as cosmic rays, meteorites, stardust grains, lunar and terrestrial sediments, and through astronomical observations across the electromagnetic spectrum. Understanding cosmic abundances and their evolution requires combining such measurements with approaches of astrophysical, nuclear theories and laboratory experiments, and exploiting additional cosmic messengers, such as neutrinos and gravitational waves. Recent years have seen significant progress in almost all these fields; they are presented in this review. Models are required to explore nuclear fusion of heavier elements. These have been confirmed by observations of nucleosynthesis products in the ejecta of stars and supernovae, as captured by stardust grains and by characteristic lines in spectra seen from these objects, and also by ejecta material captured by Earth over millions of years in sediments. All these help to piece together how cosmic materials are transported in interstellar space and re-cycled into and between generations of stars. Our description of cosmic compositional evolution needs observational support, as it rests on several assumptions that appear challenged. This overview presents the flow of cosmic matter and the various sites of nucleosynthesis, as understood from combining many techniques and observations, towards the current knowledge of how the universe is enriched with elements.

Mina Pak, Joon Hyeop Lee, Hyunjin Jeong, Woong-Seob Jeong

10 pages, 6 figures, Accepted for publication in AJ

We investigate the stellar and ionized gas kinematics, and stellar populations of NGC3182 galaxy using integral field spectrograph data from the Calar Alto Legacy Integral Field Area survey. We try to clarify the nature of the ring structure in NGC 3182. We find a negative stellar age gradient out to the ring, while [{\alpha}/Fe] considerably enhanced in the ring. The stellar metallicity shows a smooth negative gradient. From the line ratio diagnostic diagrams, we confirm that NGC 3182 is a Seyfert galaxy from emission line flux ratio, while the gas in the inner ring is ionized mostly by young stars. However, any obvious feature of outflows is not found in its gas kinematics. In the ring, star formation seems to have recently occurred and the gas metallicity is slightly enhanced compared to the center. From our results, we conclude that star formation has occurred in the circumnuclear region within a short period and this may result from a positive feedback by AGN radiation pressure.

D. Mesa, M. Bonavita, S. Benatti, R. Gratton, S. Marino, P. Kervella, V. D'Orazi, S.Desidera, T. Henning, M. Janson, M. Langlois, E. Rickman, A. Vigan, A. Zurlo, J.-L. Baudino, B. Biller, A. Boccaletti, M. Bonnefoy, W. Brandner, E. Buenzli, F. Cantalloube, D. Fantinel, C. Fontanive, R. Galicher, C. Ginski, J. Girard, J. Hagelberg, T. Kopytova, C. Lazzoni, H. Le Coroller, R. Ligi, M. Llored, A.-L. Maire, D. Mouillet, C. Perrot, S. Rochat, C. Romero, D. Rouan, M. Samland, T.O.B. Schmidt, E. Sissa, F. Wildi

11 pages, 7 figures, accepted for publications on A&A

Aims. This work aims at constraining the masses and separations of potential substellar companions to five accelerating stars (HIP 1481, HIP 88399, HIP 96334, HIP 30314 and HIP 116063) using multiple data sets acquired with different techniques. Methods. Our targets were originally observed as part of the SPHERE/SHINE survey, and radial velocity (RV) archive data were also available for four of the five objects. No companions were originally detected in any of these data sets, but the presence of significant proper motion anomalies (PMa) for all the stars strongly suggested the presence of a companion. Combining the information from the PMa with the limits derived from the RV and SPHERE data, we were able to put constraints on the characteristics of the unseen companions. Results. Our analysis led to relatively strong constraints for both HIP 1481 and HIP 88399, narrowing down the companion masses to 2-5 M_Jup and 3-5 M_Jup and separations within 2-15 au and 3-9 au, respectively. Because of the large age uncertainties for HIP 96334, the poor observing conditions for the SPHERE epochs of HIP 30314 and the lack of RV data for HIP 116063, the results for these targets were not as well defined, but we were still able to constrain the properties of the putative companions within a reasonable confidence level. Conclusions. For all five targets, our analysis has revealed that the companions responsible for the PMa signal would be well within reach for future instruments planned for the ELT (e.g., MICADO), which would easily achieve the required contrast and angular resolution. Our results therefore represent yet another confirmation of the power of multi-technique approaches for both the discovery and characterisation of planetary systems.

Silke Maes, Lionel Siess, Ward Homan, Jolien Malfait, Frederik De Ceuster, Thomas Ceulemans, Dion Donne, Mats Esseldeurs, Leen Decin

Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary interaction. In order to improve our understanding about these systems, theoretical studies are needed in the form of numerical simulations. Currently, a handful of simulations exist, albeit they mainly focus on the hydrodynamics of the outflow. Hence, we here present the pathway to more detailed and accurate modelling of companion-perturbed outflows with Phantom, by discussing the missing but crucial physical and chemical processes. With these state-of-the-art simulations we aim to make a direct comparison with observations to unveil the true identity on the embedded systems.

Kendall Sullivan (UT Austin), Adam L. Kraus (UT Austin), Andrew W. Mann (UNC Chapel Hill)

Accepted to ApJ. 28 pages, 15 figures, 4 tables (including 1 appendix). The code to reproduce the analysis is available on GitHub (see this https URL ) and the first author is happy to assist with any issues!

To fully leverage the statistical strength of the large number of planets found by projects such as the Kepler survey, the properties of planets and their host stars must be measured as accurately as possible. One key population for planet demographic studies is circumstellar planets in close binaries ($\rho < 50$au), where the complex dynamical environment of the binary inhibits most planet formation, but some planets nonetheless survive. Accurately characterizing the stars and planets in these complex systems is a key factor in better understanding the formation and survival of planets in binaries. Toward that goal, we have developed a new Markov Chain Monte Carlo fitting algorithm to retrieve the properties of binary systems using unresolved spectra, unresolved photometry, and resolved contrasts. We have analyzed 8 Kepler Objects of Interest in M star binary systems using literature data, and have found that the temperatures of the primary stars (and presumed planet hosts) are revised upward by an average of 200 K. The planetary radii should be revised upward by an average of 20% if the primary star is the host, and 80% if the secondary star is the planet host. The average contrast between stellar components in the Kepler band is 0.75 mag, which is small enough that neither star in any of the binaries can be conclusively ruled out as a potential planet host. Our results emphasize the importance of accounting for multiplicity when measuring stellar parameters, especially in the context of exoplanet characterization.

Zhe-Yu Daniel Lin, Zhi-Yun Li, Haifeng Yang, Olga Muñoz, Leslie Looney, Ian Stephens, Charles L. H. Hull, Manuel Fernández-López, Rachel Harrison

13 pages, 8 figures, submitted to MNRAS

The size of dust grains plays a fundamental role in determining the physical and chemical processes in circumstellar disks, but observational constraints of grain size, $a$, remain challenging. (Sub)millimeter continuum observations often show a percent-level polarization parallel to the disk minor axis, which is generally interpreted as coming from scattering by $\sim 100\mu$m-sized spherical grains (with a size parameter $x \equiv 2\pi a / \lambda < 1$, where $\lambda$ is the wavelength). Larger spherical grains (with $x$ greater than unity) would produce polarization parallel to the disk major axis, in contradiction to the observed pattern. The inferred grain size of $\sim 100\mu$m is in tension with the opacity index $\beta$ that point to larger mm/cm-sized grains. However, grains are expected to be non-spherical and irregular under realistic conditions. Here we investigate the scattering-produced polarization by large irregular grains with $x$ much larger than unity with optical properties obtained from laboratory experiments. Using the radiation transfer code, RADMC-3D, we simulate polarization images and find that large irregular grains still produce polarization parallel to the disk minor axis. Our results suggest that scattering grains at (sub)millimeter wavelengths in disks are not limited to $\sim 100 \mu$m, but can easily have millimeter sizes or even larger, thus alleviating the long-standing tension between the grain sizes inferred from scattering and other means. Additionally, if large irregular grains are not settled to the midplane, their strong forward scattering can produce asymmetries between the near and far side of an inclined disk, which can be used to infer their presence.

Manali Jeste, Yan Gong, Ka Tat Wong, Karl M. Menten, Tomasz Kamiński, Friedrich Wyrowski

36 pages, 26 figures, accepted for publication in A&A

Context. HCN is the most abundant molecule after H$_{2}$ and CO in the circumstellar envelopes (CSEs) of carbon-rich asymptotic giant branch (AGB) stars. Its rotational lines within vibrationally excited states are exceptional tracers of the innermost region of C-rich CSEs. Aims. We aim to constrain the physical conditions of CSEs of C-rich stars using thermal lines of the HCN molecule. Additionally, we also search for new HCN masers and probe the temporal variations for HCN masers, which should shed light on their pumping mechanisms. Methods. We observed 16 C-rich AGB stars in various HCN rotational transitions within the ground and 12 vibrationally excited states with the APEX 12 m submillimeter telescope. Results. We detect 68 vibrationally excited HCN lines from 13 C-rich stars, including 39 thermal transitions and 29 maser lines, which suggests that vibrationally excited HCN lines are ubiquitous in C-rich stars. Population diagrams constructed, for two objects from the sample, for thermal transitions from different vibrationally excited states give excitation temperature around 800-900 K, confirming that they arise from the hot innermost regions of CSEs (i.e., r<20$R_{*}$). Among the detected masers, 23 are newly detected, and the results expand the total number of known HCN masers lines toward C-rich stars by 47%. In particular, the J=2-1 (0, 3$^{1e}$, 0), J=3-2 (0, 2, 0), J=4-3 (0, 1$^{1f}$, 0) masers are detected in an astronomical source for the first time. Our observations confirm temporal variations of the 2-1 (0, 1$^{1e}$, 0) maser on a timescale of a few years. Our analysis of the data suggests that all detected HCN masers are unsaturated. A gas kinetic temperature of ${\gtrsim}$700K and an H$_{2}$ number density of >10$^{8}$cm$^{-3}$ are required to excite the HCN masers. In some ways, HCN masers in C-rich stars might be regarded as an analogy of SiO masers in O-rich stars.

Louis Perenon, Matteo Martinelli, Roy Maartens, Stefano Camera, Chris Clarkson

10 pages, 6 figures

We combine cosmic chronometer and growth of structure data to derive the redshift evolution of the dark energy equation of state $w$, using a novel agnostic approach. The background and perturbation equations lead to two expressions for $w$, one purely background-based and the other relying also on the growth rate of large-scale structure. We compare the features and performance of the growth-based $w$ to the background $w$, using Gaussian Processes for the reconstructions. We find that current data is not precise enough for robust reconstruction of the two forms of $w$. By using mock data expected from next-generation surveys, we show that the reconstructions will be robust enough and that the growth-based $w$ will out-perform the background $w$. Furthermore, any disagreement between the two forms of $w$ will provide a new test for deviations from the standard model of cosmology.

Florian Lalande, Alessandro Alberto Trani

10 pages, 6 figures, submitted to ApJ

Understanding the long-term evolution of hierarchical triple systems is challenging due to its inherent chaotic nature, and it requires computationally expensive simulations. Here we propose a convolutional neural network model to predict the stability of hierarchical triples by looking at their evolution during the first $5 \times 10^5$ inner binary orbits. We employ the regularized few-body code \textsc{tsunami} to simulate $5\times 10^6$ hierarchical triples, from which we generate a large training and test dataset. We develop twelve different network configurations that use different combinations of the triples' orbital elements and compare their performances. Our best model uses 6 time-series, namely, the semimajor axes ratio, the inner and outer eccentricities, the mutual inclination and the arguments of pericenter. This model achieves an area under the curve of over $95\%$ and informs of the relevant parameters to study triple systems stability. All trained models are made publicly available, allowing to predict the stability of hierarchical triple systems $200$ times faster than pure $N$-body methods.

In this study, we use a single-mode squeezed thermal vacuum state formalism and examine the nature of a massive homogeneous scalar field, minimally coupled to the gravity in the framework of semiclassical gravity in the Bianchi type-I universe. We have obtained an estimate leading solution to the semiclassical Einstein equation for the Bianchi type-I universe shows, each scale factor in its respective direction obeys t^(2/3) power-law expansion. The mechanism of the nonclassical thermal cosmological particle production is also analyzed in the Bianchi type-I universe.

Mike Walmsley, Inigo Val Slijepcevic, Micah Bowles, Anna M. M. Scaife

Accepted at the ICML 2022 Workshop on Machine Learning for Astrophysics. Data: www.github.com/mwalmsley/pytorch-galaxy-datasets. Please reach out to share your labelled data - all contributions will be credited in future work

New astronomical tasks are often related to earlier tasks for which labels have already been collected. We adapt the contrastive framework BYOL to leverage those labels as a pretraining task while also enforcing augmentation invariance. For large-scale pretraining, we introduce GZ-Evo v0.1, a set of 96.5M volunteer responses for 552k galaxy images plus a further 1.34M comparable unlabelled galaxies. Most of the 206 GZ-Evo answers are unknown for any given galaxy, and so our pretraining task uses a Dirichlet loss that naturally handles unknown answers. GZ-Evo pretraining, with or without hybrid learning, improves on direct training even with plentiful downstream labels (+4% accuracy with 44k labels). Our hybrid pretraining/contrastive method further improves downstream accuracy vs. pretraining or contrastive learning, especially in the low-label transfer regime (+6% accuracy with 750 labels).

P.S. Koliogiannis, G.A. Tsalis, C.P. Panos, Ch.C. Moustakidis

v1: 9 pages, 3 figures

We systematically explore the statement that the configurational entropy provides an alternative approach to studying gravitational stability of compact objects, carried out in the previous work of M. Gleiser and N. Jiang, Phys. Rev. D {\bf 92}, 044046 (2015). To guarantee the accuracy of the calculations, we first employ two analytical solutions of the Tolman-Oppenheimer-Volkoff equations, i.e. the Uniform and the Tolman VII. The predictions of the analytical solutions are compared to the corresponding numerical calculations, paying special attention on the use of the Fourier transform, which plays the major role for the connection of the configurational entropy to the bulk neutron star properties. Afterwards, we extend our study in order to include a large set of realistic equations of state. All of them are used extensively in the literature for compact star studies. In particular, we focus on neutron stars, quark stars, as well as on the the third family of compact stars (hybrid stars), where a possible phase transition may lead to the existence of twin stars (stars with equal mass but different radius). We found that a general rule to relate the stability region obtained from the traditional perturbation methods to the one obtained by the minimum of configurational entropy, does not hold. We found only one case which confirms this statement, that is the configuration corresponding to the free Fermi gas equation of state. We conclude that the suggested prediction of the stability by the minimization of the configurational entropy, concerning various compact objects, is rather a conjecture that does not have any basis on theoretical arguments and, even more, is not confirmed either empirically.

Photon initiated chemistry, \textit{i.e.} the interaction of light with chemical species, is a key factor in the evolution of the atmosphere of exoplanets. For planets orbiting stars in UV-rich environments, photodissociation induced by high energy photons dominates the atmosphere composition and dynamics. The rate of photodissociation can be highly dependent on atmospheric temperature, as increased temperature leads to increased population of vibrational excited states and the consequent lowering of the photodissociation threshold. This paper inaugurates a new series of papers presenting computed temperature-dependent photodissociation cross sections with rates generated for different stellar fields. Cross sections calculations are performed by solving the time-independent Schr\"{o}dinger equation for each electronic state involved in the process. Here photodissociation cross sections for hydrogen chloride and hydrogen fluoride are computed for a grid of 34 temperatures between 0 and 10~000 K. Use of different radiation fields shows that for the Sun and cooler stars the photodissociation rate can increase exponentially for molecular temperatures above 1000 K; conversely the photodissociation rates in UV rich fields instead are almost insensitive to the temperature of the molecule. Furthermore, these rates show extreme sensitivity to the radiation model used for cool stars, suggesting that further work on these may be required. The provision of an ExoMol database of cross sections is discussed.

Enrico Cannizzaro, Laura Sberna, Andrea Caputo, Paolo Pani

10+3 pages, 4 figures

Black-hole superradiance has been used to place very strong bounds on a variety of models of ultralight bosons such as axions, new light scalars, and dark photons. It is common lore to believe that superradiance bounds are broadly model independent and therefore pretty robust. In this work we show however that superradiance bounds on dark photons can be challenged by simple, compelling extensions of the minimal model. In particular, if the dark photon populates a larger dark sector and couples to dark fermions playing the role of dark matter, then superradiance bounds can easily be circumvented, depending on the mass and (dark) charge of the dark matter.

In this paper we investigate a non-minimal, space-time derivative dependent, coupling between the $k$-essence field and a relativistic fluid using a variational approach. The derivative coupling term couples the space-time derivative of the $k$-essence field with the fluid 4-velocity via an inner product. The inner product has a coefficient whose form specifies the various models of interaction. By introducing a coupling term at the Lagrangian level and using the variational technique we obtain the $k$-essence field equation and the Friedmann equations in the background of a spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric. Explicitly using the dynamical analysis approach we analyze the dynamics of this coupled scenario in the context of two kinds of interaction models. The models are distinguished by the form of the coefficient multiplying the derivative coupling term. In the simplest approach we work with an inverse square law potential of the $k$-essence field. Both of the models are not only capable of producing a stable accelerating solution, they can also explain different phases of the evolutionary universe.