Abstract visibility
Abstract text size

Papers for Monday, Jul 05 2021

Papers with local authors

Jiachen Ding, Lifan Wang, Peter Brown, Ping Yang

21 pages, 20 figures, accepted to ApJ

0 votes
Paper 31 — arXiv:2107.01026
0 votes
Paper 31 — arXiv:2107.01026

We use a Monte Carlo radiative transfer model (MCRTM) to simulate the UBVRI light curves, images and linear polarization of a light echo from supernova SN$~$1987A in the Large Magellanic Cloud (LMC) using various dust cloud shapes, sizes, and optical properties. We compare the theoretical simulations to the observations of AT2019xis, a light echo detected at a large angular distance (4.05$^{'}$) from SN$~$1987A. We estimate the size and optical thickness of the dust cloud based on the simulation results and the observations of Optical Gravitational Lensing Experiment (OGLE-IV) Transient Detection System (OTDS) I-band light curve. The mass of the dust cloud is calculated using the estimated size, optical thickness and extinction coefficient. If the dust cloud is assumed to correspond to a gas-to-dust ratio of 300, the total mass of the dust cloud is approximately 7.8-9.3 $M_{\odot}$. Based on these theoretical models, we show that the morphological shapes of the light echoes in the wavelength range in or shorter than the U-band to be very different from those in the longer wavelength bands, and the difference carries important information on the early UV radiation of SN$~$1987A.

M. Amenomori, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He, K. Hibino, N. Hotta, Haibing Hu, H. B. Hu, J. Huang, H. Y. Jia, L. Jiang, P. Jiang, H. B. Jin, K. Kasahara, Y. Katayose, C. Kato, S. Kato, K. Kawata, M. Kozai, D. Kurashige, Labaciren, G. M. Le, A. F. Li, H. J. Li, W. J. Li, Y. Li, Y. H. Lin, B. Liu, C. Liu, J. S. Liu, L. Y. Liu, M. Y. Liu, W. Liu, X. L. Liu, Y.-Q. Lou, H. Lu, X. R. Meng, K. Munakata, H. Nakada, Y. Nakamura, Y. Nakazawa, H. Nanjo, C. C. Ning, M. Nishizawa, M. Ohnishi, T. Ohura, S. Okukawa, S. Ozawa, L. Qian, X. Qian, X. L. Qian, X. B. Qu, T. Saito, M. Sakata, T. Sako, et al.

Accepted for publication in the Physical Review Letters

0 votes
Paper 36 — arXiv:2107.01064
0 votes
Paper 36 — arXiv:2107.01064

We report observations of gamma-ray emissions with energies in the 100 TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1 with the pulsar moving away from its original birthplace situated around the centroid of the observed gamma-ray emission. This work would stimulate further studies of particle acceleration mechanisms at these gamma-ray sources.

All other papers

Andrija Kostić, Jens Jasche, Doogesh Kodi Ramanah, Guilhem Lavaux

5 + 3 pages, 4 figures. Submitted to A&A Letters. Comments welcome

We present maps revealing the expected information content of cosmic large-scale structures concerning cosmological physics. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. This approach has become feasible through the recent development of causal inference machinery that is informed on the physics of cosmic structure formation. Specifically, we measure the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contributions to the Fisher information. Our physical forward modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We showcase the potential of our approach by studying the information content of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative. The results presented in this work are the first of their kind and elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward to perform efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework.

Michał J. Michałowski, Przemek Mróz

ApJ Letters, accepted; 9 pages, 4 figures, 1 table

Gravitational lensing is an important prediction of general relativity, providing both its test and a tool to detect faint but amplified sources and to measure masses of lenses. For some applications (e.g. testing the theory), a point source lensed by a point-like lens would be more advantageous. However, until now only one gravitationally lensed star has been resolved. The future telescopes will resolve very small lensing signatures for stars orbiting the supermassive black hole (SMBH) in the center of the Milky Way. The lensing signatures should however be easier to detect for background stars. We predict that ELT, TMT, and GMT will resolve the lensed images of around 100 (60) stars in the disk and 30 (20) stars in the bulge in the background of the SMBH, down to 28 (27) mag (Vega) limits at K-band, requiring 5 (1) hr of integration. In order to detect several such stars one needs the limit of at least 24 mag. With a decade-long monitoring one can also detect the rotation of the lensed images. The detection of elongated images will not be possible, because this would require a nearly perfect source-lens alignment. JWST will likely be limited by the confusion caused by stars near the Galactic center. The detection of such lensed images will provide a very clean test of general relativity, when combined with the SMBH mass measurement from orbital motions of stars, and accurate measurements of the SMBH properties, because both the source and the lens can be considered point-like.

M. Llerena, R. Amorín, F. Cullen, L. Pentericci, A. Calabrò, R. McLure, A. Carnall, E. Pérez-Montero, F. Marchi, A. Bongiorno, M. Castellano, A. Fontana, D. J. McLeod, M. Talia, N. P. Hathi, P. Hibon, F. Mannucci, A. Saxena, D. Schaerer, G. Zamorani

Submitted to A&A. Comments are welcome

We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at $2<z<4$, selected from a parent sample of $\sim$750 main-sequence star-forming galaxies in the VANDELS survey. These CIII] emitters have a broad range of UV luminosities, thus allowing a detailed stacking analysis to characterize their stellar mass, star formation rate (SFR) and stellar metallicity, as a function of the UV emission line ratios, EWs, and the carbon-to-oxygen (C/O) abundance ratio. Reliable CIII] detections represent $\sim$30% of the parent sample. Extreme CIII] emitters (EW(CIII])$\gtrsim$8\r{A}) are exceedingly rare ($\sim$3%) in VANDELS. The UV line ratios of the sample suggest no ionization source other than massive stars. Stacks with larger EW(CIII]) show larger EW(Ly$\alpha$) and lower metallicity, but not all CIII] emitters are Ly$\alpha$ emitters. The stellar metallicities of CIII] emitters are not significantly different from that of the parent sample, increasing from $\sim$10% to $\sim$40% solar for stellar masses $\log$(M$_{\star}$/M$_{\odot})\sim$9-10.5. The stellar mass-metallicity relation of the CIII] emitters is consistent with previous works showing strong evolution from $z=0$ to $z\sim3$. The C/O abundances of the sample range 35%-150% solar, with a noticeable increase with FUV luminosity and a smooth decrease with the CIII] and Ly$\alpha$ EWs. We discuss the CIII] emitters in the C/O-Fe/H and the C/O-O/H planes and find they follow stellar and nebular abundance trends consistent with those of Milky Way halo and thick disc stars and local HII galaxies, respectively. A qualitative agreement is also found with chemical evolution models, which suggests that CIII] emitters at $z\sim$3 are experiencing an active phase of chemical enrichment.

D. M. Rowan, K. Z. Stanek, Z. Way, C. S. Kochanek, T. Jayasinghe, Todd A. Thompson, H. Barker, F.-J. Hamsch, T. Bohlsen, Stella Kafka, B. J. Shappee, T. W.-S. Holoien, J. L. Prieto

4 pages, 1 figure, published in RNAAS

We use ASAS V-band and ASAS-SN g-band observations to model the long-period detached eclipsing binary ASASSN-21co. ASAS observations show an eclipse of depth V ~ 0.6 mag in April of 2009. ASAS-SN g-band observations from March of 2021 show an eclipse of similar duration and depth, suggesting an orbital period of 11.9 years. We combine the g-band observations with additional BVRI photometry taken during the eclipse to model the eclipse using PHOEBE. We find that the system is best described by two M giants with a ratio of secondary radius to primary radius of ~0.61. Optical spectra taken during the eclipse are consistent with at least one component of the binary being an M giant, and we find no temporal changes in the spectral features. The eclipse itself is asymmetric, showing an increase in brightness near mid-eclipse, likely due to rotational variability that is too low amplitude to be observed out-of-eclipse.

Mateo Prgomet, Martin P. Rey, Eric P. Andersson, Alvaro Segovia Otero, Oscar Agertz, Florent Renaud, Andrew Pontzen, Justin I. Read

Main text 7 pages. Submitted to MNRAS, comments welcome

We study how an observationally-motivated, metallicity-dependent initial mass function (IMF) affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions. We complement the EDGE galaxy formation model from Agertz et al. (2020) with a new prescription for IMF variations according to Geha et al. (2013). At the low metallicities typical of faint dwarf galaxies, the IMF becomes top-heavy, increasing the efficiency of supernova and photo-ionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass-metallicity relation.

Nicholas M. Law, Hank Corbett, Nathan W. Galliher, Ramses Gonzalez, Alan Vasquez, Glenn Walters, Lawrence Machia, Jeff Ratzloff, Kendall Ackley, Chris Bizon, Christopher Clemens, Steven Cox, Steven Eikenberry, Ward S. Howard, Amy Glazier, Andrew W. Mann, Robert Quimby, Daniel Reichart, David Trilling

17 pages, 5 figures, 2 tables

New mass-produced, wide-field, small-aperture telescopes have the potential to revolutionize ground-based astronomy by greatly reducing the cost of collecting area. In this paper, we introduce a new class of large telescope based on these advances: an all-sky, arcsecond-resolution, 1000-telescope array which builds a simultaneously high-cadence and deep survey by observing the entire sky all night. As a concrete example, we describe the Argus Array, a 5m-class telescope with an all-sky field of view and the ability to reach extremely high cadences using low-noise CMOS detectors. Each 55 GPix Argus exposure covers 20% of the entire sky to g=19.6 each minute and g=21.9 each hour; a high-speed mode will allow sub-second survey cadences for short times. Deep coadds will reach g=23.6 every five nights over 47% of the sky; a larger-aperture array telescope, with an \'etendue close to the Rubin Observatory, could reach g=24.3 in five nights. These arrays can build two-color, million-epoch movies of the sky, enabling sensitive and rapid searches for high-speed transients, fast-radio-burst counterparts, gravitational-wave counterparts, exoplanet microlensing events, occultations by distant solar system bodies, and myriad other phenomena. An array of O(1,000) telescopes, however, would be one of the most complex astronomical instruments yet built. Standard arrays with hundreds of tracking mounts entail thousands of moving parts and exposed optics, and maintenance costs would rapidly outpace the mass-produced-hardware cost savings compared to a monolithic large telescope. We discuss how to greatly reduce operations costs by placing all optics in a thermally controlled, sealed dome with a single moving part. Coupled with careful software scope control and use of existing pipelines, we show that the Argus Array could become the deepest and fastest Northern sky survey, with total costs below $20M.

Tom Crossland, Pontus Stenetorp, Daisuke Kawata, Sebastian Riedel, Thomas D. Kitching, Anurag Deshpande, Tom Kimpson, Choong Ling Liew-Cain, Christian Pedersen, Davide Piras, Monu Sharma

23 pages, 14 figures. Submitted to Monthly Notices of the Royal Astronomical Society. Astronomical measurement database available at this http URL

We develop a new model for automatic extraction of reported measurement values from the astrophysical literature, utilising modern Natural Language Processing techniques. We use this model to extract measurements present in the abstracts of the approximately 248,000 astrophysics articles from the arXiv repository, yielding a database containing over 231,000 astrophysical numerical measurements. Furthermore, we present an online interface (Numerical Atlas) to allow users to query and explore this database, based on parameter names and symbolic representations, and download the resulting datasets for their own research uses. To illustrate potential use cases we then collect values for nine different cosmological parameters using this tool. From these results we can clearly observe the historical trends in the reported values of these quantities over the past two decades, and see the impacts of landmark publications on our understanding of cosmology.

Alaina Henry, Marc Rafelski, Ben Sunnquist, Norbert Pirzkal, Camilla Pacifici, Hakim Atek, Micaela Bagley, Ivano Baronchelli, Guillermo Barro, Andrew J. Bunker, James Colbert, Y. Sophia Dai, Bruce G. Elmegreen, Debra Meloy Elmegreen, Steven Finkelstein, Dale Kocevski, Anton Koekemoer, Matthew Malkan, Crystal L. Martin, Vihang Mehta, Anthony Pahl, Casey Papovich, Michael Rutkowski, Jorge Sanchez Almeida, Claudia Scarlata, Gregory Snyder, Harry Teplitz

Accepted for publication in ApJ; 41 pages, 20 figures

We present a new measurement of the gas-phase mass-metallicity relation (MZR), and its dependence on star formation rates (SFRs) at 1.3 < z < 2.3. Our sample comprises 1056 galaxies with a mean redshift of z = 1.9, identified from the Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopy in the Cosmic Assembly Near-Infrared Deep Extragalactic Survey (CANDELS) and the WFC3 Infrared Spectroscopic Parallel Survey (WISP). This sample is four times larger than previous metallicity surveys at z ~ 2, and reaches an order of magnitude lower in stellar mass (10^8 M_sun). Using stacked spectra, we find that the MZR evolves by 0.3 dex relative to z ~ 0.1. Additionally, we identify a subset of 49 galaxies with high signal-to-noise (SNR) spectra and redshifts between 1.3 < z < 1.5, where H-alpha emission is observed along with [OIII] and [OII]. With accurate measurements of SFR in these objects, we confirm the existence of a mass-metallicity-SFR (M-Z-SFR) relation at high redshifts. These galaxies show systematic differences from the local M-Z-SFR relation, which vary depending on the adopted measurement of the local relation. However, it remains difficult to ascertain whether these differences could be due to redshift evolution, as the local M-Z-SFR relation is poorly constrained at the masses and SFRs of our sample. Lastly, we reproduced our sample selection in the IllustrisTNG hydrodynamical simulation, demonstrating that our line flux limit lowers the normalization of the simulated MZR by 0.2 dex. We show that the M-Z-SFR relation in IllustrisTNG has an SFR dependence that is too steep by a factor of around three.

Compact astrophysical objects have been considered in the literature as dark matter (DM) probes, via the observational effects of annihilating captured DM. In this respect, Population III (Pop III) stars are particularly interesting candidates, since they form at high redshifts, in very high DM density environments. It is customary to assume such a star would form roughly at the center of a DM halo, and, as such, have no rotational velocity. In this paper, we break from this assumption and explore the effects we can expect to observe if a Pop III star forms at some distance away from the center of the halo and thus has a non-zero rotational velocity. The capture rate of DM in such a star is suppressed by a predictable amount. We develop and validate an analytical expression for the capture rate suppression factor and re-evaluate the bounds placed on the DM-nucleon cross section as a result of DM capture. We find that our previous results, obtained under the assumption of star formed within the central 10 AU of the DM mini-halo are essentially unchanged, even when considering the possible rotational velocities for those central stars.

M. Witt, A. Psaltis, H. Yasin, C. Horn, M. Reichert, T. Kuroda, M. Obergaulinger, S. M. Couch, A. Arcones

16 pages, 15 figures, submitted to ApJ

We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_\odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.

B. Ansarinejad (1, 2), T. Shanks (1), R. M. Bielby (1, 3), N. Metcalfe (1), L. Infante (4, 5, 6), D. N. A. Murphy (7), D. J. Rosario (1), S. M. Stach (1) ((1) Durham University, (2) The University of Melbourne, (3) Department for Education, (4) Las Campanas Observatory, (5) Universidad Diego Portales, (6) Pontificia Universidad Catolica de Chile, (7) University of Cambridge)

14 pages, 7 figures, submitted to MNRAS

We investigate the dust heating mechanisms of sub-mm galaxies (SMGs) to consider the contribution of Active Galactic Nuclei (AGN) compared to star-formation. We have used ALMA at $0.''1$ resolution to image a complete sample of seven sub-mm sources previously shown to have spectral energy distributions (SEDs) that were as well-fitted by obscured AGN as star-forming galaxy templates. Indeed, two sub-mm sources were known to be quasars from their absorbed X-ray emission. We now find that the sub-mm sizes of all SMGs are small ($\approx1-2$ kpc) and generally $>3\times$ smaller than any host detected in the near-Infra-Red. In all cases, the five SMGs are comparable in sub-mm size to the two known quasars and four $z\approx6$ quasars, also observed with ALMA. We detect no evidence of diffuse spiral arms in this complete sample. We then convert the Far-Infra-Red (FIR) luminosities to star-formation rate (SFR) surface densities and find that the SMGs occupy the same range as the known quasars in our sample. We conclude that in terms of sub-mm size, extent relative to host and SFR density as well as luminosity and mid-IR colour, there is little distinction between the SMGs and sub-mm bright quasars. In light of these results, we continue to suggest that dust+gas absorbed quasars may simultaneously dominate the FIR and hard X-ray backgrounds.

S. K. Leggett, Pascal Tremblin, Mark W. Phillips, Trent J. Dupuy, Mark Marley, Caroline Morley, Adam Schneider, Dan Caselden, Colin Guillaume, Sarah E. Logsdon

40 pages which includes 16 Figures and 10 Tables. The Journal publication will include data behind the Figures for Figures 5, 8 and 9, and a machine readable version of Table 10

Cold, low-mass, field brown dwarfs are important for constraining the terminus of the stellar mass function, and also for optimizing atmospheric studies of exoplanets. In 2020 new model grids for such objects were made available: Sonora-Bobcat and ATMO 2020. Also, new candidate cold brown dwarfs were announced, and new spectroscopic observations at lambda ~4.8 um were published. In this paper we present new infrared photometry for some of the coldest brown dwarfs, and put the new data and models together to explore the properties of these objects. We reconfirm the importance of mixing in these atmospheres, which leads to CO and NH_3 abundances that differ by orders of magnitude from chemical equilibrium values. We also demonstrate that the new models retain the known factor >~3 discrepancy with observations at 2 <~ lambda um <~ 4, for brown dwarfs cooler than 600 K. We show that the entire 1 <~ lambda um <~ 20 energy distribution of six brown dwarfs with 260 <= T_eff K <= 475 can be well reproduced, for the first time, by model atmospheres which include dis-equilibrium chemistry as well as a photospheric temperature gradient which deviates from the standard radiative/convective equilibrium value. This change to the pressure-temperature profile is not unexpected for rotating and turbulent atmospheres which are subject to diabatic processes. A limited grid of modified-adiabat model colors is generated, and used to estimate temperatures and metallicities for the currently known Y dwarfs. A compilation of the photometric data used here is given in the Appendix.

Supermassive black holes (with $\mathrm{M_{BH} \sim 10^9 M_{\odot}}$) are observed in the first Gyr of the Universe, and their host galaxies are found to contain unexpectedly large amounts of dust and metals. In light of the two empirical facts, we explore the possibility of supercritical accretion and early black hole growth occurring in dusty environments. We generalise the concept of photon trapping to the case of dusty gas and analyse the physical conditions leading to dust photon trapping. Considering the parameter space dependence, we obtain that the dust photon trapping regime can be more easily realised for larger black hole masses, higher ambient gas densities, and lower gas temperatures. The trapping of photons within the accretion flow implies obscured active galactic nuclei (AGNs), while it may allow a rapid black hole mass build-up at early times. We discuss the potential role of such dust photon trapping in the supercritical growth of massive black holes in the early Universe.

Brian Punsly, Sandor Frey, Cormac Reynolds, Paola Marziani, Alexander Pushkarev, Sina Chen, Shang Li, Preeti Kharb

To appear in ApJ

PKS\,1351$-$018 at a redshift of $z=3.71$ is one of the most luminous, steady synchrotron sources with a luminosity $> 10^{47}$\,erg~s$^{-1}$. The synchrotron luminosity does not seem to vary by more than $\sim 25\%$ over 35 years. In order to appreciate this remarkable behavior, if it were at $z=0.5$, it would have a flux density at 15 GHz in a range of $110 - 137$\,Jy over 11 yrs. In spite of this steady behavior, two strong $\gamma$-ray flares $\lesssim 10^{49}$\,erg~s$^{-1}$ were detected in 2011 and 2016. There is blazar-like behavior coexisting with the steady behavior. This study is aimed at elucidating the dual nature of this source. We find that the radio source is extremely compact with a bright core and a steep spectrum secondary component, 12\,mas away, that appears to be constant in position and flux density in six epochs from 1995 to 2018. We estimate that a jet with a time averaged power of $(5.2 \pm 3.2) \times 10^{45}$\,erg~s$^{-1}$ terminates in this lobe that is advancing $\gtrsim 0.9 c$ at a deprojected distance of $1-3$\,kpc from the central engine. This is the rare case of a young ($\sim 6000$\,yr), very powerful radio source that is viewed a few degrees from the jet axis. We find evidence of a high velocity (4000\,km~s$^{-1}$), high ionization wind emanating form a luminous quasar. The young radio jet appears to experience modest bending as it navigates through the intense quasar environment.

Shivan Khullar, Christoph Federrath, Mark R. Krumholz, Christopher D. Matzner

18 pages, 10 figures, Accepted in MNRAS

We conduct numerical experiments to determine the density probability distribution function (PDF) produced in supersonic, isothermal, self-gravitating turbulence of the sort that is ubiquitous in star-forming molecular clouds. Our experiments cover a wide range of turbulent Mach number and virial parameter, allowing us for the first time to determine how the PDF responds as these parameters vary, and we introduce a new diagnostic, the dimensionless star formation efficiency versus density ($\epsilon_{\rm ff}(s)$) curve, which provides a sensitive diagnostic of the PDF shape and dynamics. We show that the PDF follows a universal functional form consisting of a log-normal at low density with two distinct power law tails at higher density; the first of these represents the onset of self-gravitation, and the second reflects the onset of rotational support. Once the star formation efficiency reaches a few percent, the PDF becomes statistically steady, with no evidence for secular time-evolution at star formation efficiencies from about five to 20 percent. We show that both the Mach number and the virial parameter influence the characteristic densities at which the log-normal gives way to the first power-law, and the first to the second, and we extend (for the former) and develop (for the latter) simple theoretical models for the relationship between these density thresholds and the global properties of the turbulent medium.

Fabio Cafardo, Rodrigo Nemmen

Accepted for publication in The Astrophysical Journal. 21 pages, 9 figures, 4 tables

Sagittarius A* (Sgr A*) -- the supermassive black hole (SMBH) in the center of our galaxy -- has been observed in most of the electromagnetic spectrum, from radio to X-rays. Diffuse gamma-ray emission has been observed around Sgr A* and a gamma-ray point source has been detected coinciding with the SMBH's position, although there is no definitive association between the two. In this work, we have used ~11 years of Fermi Large Area Telescope (LAT) observations of the point source 4FGL J1745.6-2859 and performed a detailed imaging analysis across four energy bands. Our goal is to elucidate the nature of the gamma-ray emission at the Galactic Center (GC) and whether it is associated with the SMBH. We find that the centroid of the emission approaches Sgr A*'s location as the energy increases. Assuming that the gamma-ray point source is located at the GC, we estimate a luminosity of 2.61E36 erg/s in the 100 MeV to 500 GeV energy range. This is consistent with Sgr A*'s bolometric luminosity. Based on the point source properties, we ruled out several potential candidates for its nature and favor a cosmic ray origin either from protons, electrons or both, accelerated by -- or in the vicinity of -- the SMBH. Our results indicate that the point source at the GC is indeed the gamma-ray counterpart of Sgr A* in the GeV range.

P. Kornecki, E. Peretti, S. del Palacio, P. Benaglia, L. J. Pellizza

Submitted to Astronomy & Astrophysics (on 12/05/2021)

Star-forming galaxies (SFGs) emit non-thermal radiation from radio to gamma-rays. We aim to investigate the main mechanisms of global CR transport and cooling in SFGs. The way they contribute in shaping the relations between non-thermal luminosities and SFR could shed light onto their nature. We develop a model to compute the CR populations of SFGs, taking into account their production, transport, and cooling. The model is parameterised only through global galaxy properties, and describes the non-thermal emission in both radio and gamma-rays. We focus on the role of diffusive and advective transport by galactic winds, either driven by turbulent or thermal instabilities. We compare model predictions to observations, for which we compile a homogeneous set of luminosities in these radio bands, and update those available in gamma-rays. Our model reproduces reasonably well the observed relations between the gamma-ray or 1.4 GHz radio luminosities and the SFR, assuming a single power-law scaling of the magnetic field with the latter with index beta=0.3, and winds blowing either at Alfvenic speeds or typical starburst wind velocities. Escape of CR is negligible for > 30 Mo/yr. A constant ionisation fraction of the interstellar medium fails to reproduce the 150 MHz radio luminosity throughout the whole SFR range. Our results reinforce the idea that galaxies with high SFR are CR calorimeters, and that the main mechanism driving proton escape is diffusion, whereas electron escape also proceeds via wind advection. They also suggest that these winds should be CR or thermally-driven at low and intermediate SFR, respectively. Our results globally support that magnetohydrodynamic turbulence is responsible for the dependence of the magnetic field strength on the SFR and that the ionisation fraction is strongly disfavoured to be constant throughout the whole SFR range.

Presence of large-scale surface magnetic field in early-type stars leads to several unique electromagnetic phenomena producing radiation over X-ray to radio bands. Among them, the rarest type of emission is electron cyclotron maser emission (ECME) observed as periodic, circularly polarized radio pulses. The phenomenon was first discovered in the hot magnetic star CU Vir. Past observations of this star led to the consensus that the star produces only right circularly polarized ECME, suggesting that only one magnetic hemisphere takes part in the phenomenon. Here we present the first ultra-wideband (0.4$-$4 GHz) study of this star using the upgraded Giant Metrewave Radio telescope and the Karl G. Jansky Very Large Array, which led to the surprising discovery of ECME of both circular polarizations up to around 1.5 GHz. The GHz observations also allowed us to infer that the upper ECME cut-off frequency is at $\gtrsim 5\,\mathrm{GHz}$. The sub-GHz observation led to the unexpected observation of more than two pairs of ECME pulses per rotation cycle. In addition, we report the discovery of a `giant pulse', and transient enhancements, which are potentially the first observational evidence of `centrifugal breakout' of plasma from the innermost part of the stellar magnetosphere. The stark contrast between the star's behavior at GHz and sub-GHz frequencies could either be due to propagation effects, a manifestation of varying magnetic field topology as a function of height, or a signature of an additional `ECME engine'.

Leonid Kitchatinov, Anna Khlystova

12 pages, 9 figures, accepted in ApJ

Observations reveal a relatively small but statistically significant North-South (NS) asymmetry in sunspot activity varying on a time scale of several solar cycles. This paper proposes a dynamo model for the phenomenon of long-term NS asymmetry. The model separates dynamo equations for magnetic fields of dipolar and quadrupolar equatorial parity. The NS asymmetry results from the superposition of dipolar and quadrupolar fields. Model computations confirm the formerly proposed excitation of the quadrupolar dynamo mode by a dominant dipolar mode mediated by the equator-symmetric fluctuations in the $\alpha$-effect as a mechanism for the long-term NS asymmetry. An analytically solvable example of oscillations excited by short-term random forcing is given to justify the numerical result of NS asymmetry coherent on a time scale of several (about 6 in the present model) solar cycles resulting from random variations in the $\alpha$-effect on a time scale of one solar rotation. The model computations show the phase locking phenomenon of dipolar and quadrupolar fields oscillating predominantly in phase (northern type asymmetry) or in antiphase (southern type asymmetry) with relatively short irregular transitions between these two states. Large asymmetry in the simulated Grand minima is found and explained by weak magnetic quenching of the $\alpha$-effect during the minima. The possibility of polar field asymmetry in activity minima as a precursor of sunspot asymmetry in the following activity cycles is discussed based on the dynamo model and observations.

Molecular clouds have complex density structures produced by processes including turbulence and gravity. We propose a triangulation-based method to dissect the density structure of a molecular cloud and study the interactions between dense cores and their environments. In our {approach}, a Delaunay triangulation is constructed, which consists of edges connecting these cores. Starting from this construction, we study the physical connections between neighboring dense cores and the ambient environment in a systematic fashion. We apply our method to the Cygnus-X massive GMC complex and find that the core separation is related to the mean surface density by $\Sigma_{\rm edge} \propto l_{\rm core }^{-0.28 }$, which can be explained by {fragmentation controlled by a scale-dependent turbulent pressure (where the pressure is a function of scale, e.g. $p\sim l^{2/3}$)}. We also find that the masses of low-mass cores ($M_{\rm core} < 10\, M_{\odot}$) are determined by fragmentation, whereas massive cores ($M_{\rm core} > 10\, M_{\odot}$) grow mostly through accretion. The transition from fragmentation to accretion coincides with the transition from a log-normal core mass function (CMF) to a power-law CMF. By constructing surface density profiles measured along edges that connect neighboring cores, we find evidence that the massive cores have accreted a significant fraction of gas from their surroundings and thus depleted the gas reservoir. Our analysis reveals a picture where cores form through fragmentation controlled by scale-dependent turbulent pressure support, followed by accretion onto the massive cores, {and the method can be applied to different regions to achieve deeper understandings in the future.

We compute the cosmological boost factor at high redshifts of $z=10$-$100$ by integrating the non-linear matter power spectrum measured from high-resolution cosmological $N$-body simulations. An accurate boost factor is required to estimate the energy injection from dark matter annihilation, which may affect the cosmological re-ionization process. We combined various box-size simulations (side lengths of $1 \, {\rm kpc}$-$10 \, {\rm Mpc}$) to cover a wide range of scales, i.e. $k=1$-$10^7 \, {\rm Mpc}^{-1}$. The boost factor is consistent with the linear theory prediction at $z \gtrsim 50$ but strongly enhanced at $z \lesssim 50$ as a result of non-linear matter clustering. Although dark matter free-streaming damping was imposed at $k_{\rm fs}=10^6 \, {\rm Mpc}^{-1}$ in the initial power spectrum, the damping disappears at later times of $z\lesssim40$ as a result of the power transfer from large to small scales. Because the simulations do not explore very small-scale clustering at $k>10^7 \, {\rm Mpc}^{-1}$, our result is a lower bound on the boost factor at $z \lesssim 40$. A simple fitting function of the boost factor is also presented.

Krisztina É. Gabányi, Sándor Frey, Krisztina Perger

7 pages, 3 figures, accepted to MNRAS

The coevolution of galaxies and their central supermassive black holes is a subject of intense research. A class of objects, the dust-obscured galaxies (DOGs) are particularly interesting in this respect as they are thought to represent a short evolutionary phase when violent star formation activity in the host galaxy may coexist with matter accretion onto the black hole powering the active nucleus. Here we investigate different types of DOGs classified by their mid-infrared spectral energy distributions to reveal whether they can be distinguished by their arcsec-scale radio properties. Radio emission is unaffected by dust obscuration and may originate from both star formation and an active nucleus. We analyse a large sample of 661 DOGs complied from the literature and find that only a small fraction of them ($\sim 2$ per cent) are detected with flux densities exceeding $\sim 1$ mJy in the Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. These radio-detected objects are almost exclusively `power-law' DOGs. Stacking analysis of the FIRST image cutouts centred on the positions of individually radio-undetected sources suggests that weak radio emission is present in `power-law' DOGs. On the other hand, radio emission from `bump' DOGs is only marginally detected in the median-stacked FIRST image.

M. Fabrizio, V.F. Braga, J. Crestani, G. Bono, I. Ferraro, G. Fiorentino, G. Iannicola, G.W. Preston, C. Sneden, F. Thévenin, G. Altavilla, B. Chaboyer, M. Dall'Ora, R. da Silva, E. K. Grebel, C.K. Gilligan, H. Lala, B. Lemasle, D. Magurno, M. Marengo, S. Marinoni, P.M. Marrese, C. E. Martìnez-Vàzquez, N. Matsunaga, M. Monelli, J. P. Mullen, J. Neeley, M. Nonino, Z. Prudil, M Salaris, P. B. Stetson, E. Valenti, M. Zoccali

22 pages, 13 figures, 3 tables. Accepted for publication in ApJ

We discuss the largest and most homogeneous spectroscopic dataset of field RR Lyrae variables (RRLs) available to date. We estimated abundances using both high-resolution and low-resolution ({\Delta S} method) spectra for fundamental (RRab) and first overtone (RRc) RRLs. The iron abundances for 7,941 RRLs were supplemented with similar literature estimates available, ending up with 9,015 RRLs (6,150 RRab, 2,865 RRc). The metallicity distribution shows a mean value of <[Fe/H]> = -1.51\pm0.01, and {\sigma}(standard deviation)= 0.41 dex with a long metal-poor tail approaching [Fe/H] = -3 and a sharp metal-rich tail approaching solar iron abundance. The RRab variables are more metal-rich (<[Fe/H]>ab = -1.48\pm0.01, {\sigma} = 0.41 dex) than RRc variables (<[Fe/H]>c = -1.58\pm0.01, {\sigma} = 0.40 dex). The relative fraction of RRab variables in the Bailey diagram (visual amplitude vs period) located along the short-period (more metal-rich) and the long-period (more metal-poor) sequences are 80% and 20\%, while RRc variables display an opposite trend, namely 30\% and 70\%. We found that the pulsation period of both RRab and RRc variables steadily decreases when moving from the metal-poor to the metal-rich regime. The visual amplitude shows the same trend, but RRc amplitudes are almost two times more sensitive than RRab amplitudes to metallicity. We also investigated the dependence of the population ratio (Nc/Ntot) of field RRLs on the metallicity and we found that the distribution is more complex than in globular clusters. The population ratio steadily increases from ~0.25 to ~0.36 in the metal-poor regime, it decreases from ~0.36 to ~0.18 for -1.8 < [Fe/H] < -0.9 and it increases to a value of ~0.3 approaching solar iron abundance.

V.F. Braga (1, 2), J. Crestani (1, 3, 4), M. Fabrizio (1, 2), G. Bono (1, 4), G.W. Preston (5), C. Sneden (6), J. Storm (7), S. Kamann (8), M. Latour (9), H. Lala (10), B. Lemasle (10), Z. Prudil (11), G. Altavilla (1, 2), B. Chaboyer (12), M. Dall'Ora (13), I. Ferraro (1), C.K. Gilligan (12), G. Fiorentino (1), G. Iannicola (1), L. Inno (14), S. Kwak (4), M. Marengo (15), S. Marinoni (1, 2), P.M. Marrese (1, 2), C.E. Martínez-Vázquez (16), M. Monelli (17), J.P. Mullen (15), N. Matsunaga (18), J. Neeley (19), P.B. Stetson (20), E. Valenti (21), M. Zoccali (22, 23) ((1) INAF-Osservatorio Astronomico di Roma, Monte Porzio Catone, Italy, (2) Space Science Data Center, Roma, Italy, (3) Departamento de Astronomia, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil, (4) Dipartimento di Fisica, Università di Roma Tor Vergata, Roma, Italy, (5) The Observatories of the Carnegie Institution for Science, Pasadena, CA, USA, (6) Department of Astronomy and McDonald Observatory, The University of Texas, Austin, TX, USA, (7) Leibniz-Institut für Astrophysik Potsdam, Potsdam, Germany, (8) Astrophysics Research Institute, Liverpool John Moores University, Liverpool, UK, (9) Institute for Astrophysics, Georg-August-University Göttingen, Germany, (10) Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Heidelberg, Germany, (11) Université de Nice Sophia-antipolis, CNRS, Observatoire de la Côte d'Azur, Laboratoire Lagrange, Nice, France, (12) Department of Physics and Astronomy, Dartmouth College, Hanover, USA, (13) INAF-Osservatorio Astronomico di Capodimonte, Napoli, Italy, (14) Università degli Studi di Napoli "Parthenope'', Napoli, Italy, (15) Department of Physics and Astronomy, Iowa State University, Ames, IA, USA, (16) Cerro Tololo Inter-American Observatory, NSF's National Optical-Infrared Astronomy Research Laboratory, La Serena, Chile, (17) Instituto de Astrofísica de Canarias, La Laguna, Tenerife, Spain, (18) Department of Astronomy, The University of Tokyo, Tokyo, Japan, (19) Department of Physics, Florida Atlantic University, Boca Raton, FL, USA, (20) Herzberg Astronomy and Astrophysics, National Research Council, Victoria, BC, Canada, (21) European Southern Observatory, Garching bei Munchen, Germany, (22) Instituto Milenio de Astrofísica, Santiago, Chile, (23) Pontificia Universidad Catolica de Chile, Instituto de Astrofisica, Santiago, Chile)

We collected the largest spectroscopic catalog of RR Lyrae (RRLs) including $\approx$20,000 high-, medium- and low-resolution spectra for $\approx$10,000 RRLs. We provide the analytical forms of radial velocity curve (RVC) templates. These were built using 36 RRLs (31 fundamental -- split into three period bins -- and 5 first overtone pulsators) with well-sampled RVCs based on three groups of metallic lines (Fe, Mg, Na) and four Balmer lines (H$_\alpha$, H$_\beta$, H$_\gamma$, H$_\delta$). We tackled the long-standing problem of the reference epoch to anchor light curve and RVC templates. For the $V$-band, we found that the residuals of the templates anchored to the phase of the mean magnitude along the rising branch are $\sim$35\% to $\sim$45\% smaller than those anchored to the phase of maximum light. For the RVC, we used two independent reference epochs for metallic and Balmer lines and we verified that the residuals of the RVC templates anchored to the phase of mean RV are from 30\% (metallic lines) up to 45\% (Balmer lines) smaller than those anchored to the phase of minimum RV. We validated our RVC templates by using both the single- and the three-phase points approach. We found that barycentric velocities based on our RVC templates are two-three times more accurate than those available in the literature. We applied the current RVC templates to Balmer lines RVs of RRLs in the globular NGC~3201 collected with MUSE at VLT. We found the cluster barycentric RV of $V_{\gamma}$=496.89$\pm$8.37(error)$\pm$3.43 (standard deviation) km/s, which agrees well with literature estimates.

M. V. Popov, N. Bartel, M. S. Burgin, V. A. Soglasnov

6 pages, 4 figures. Accepted for publication in MNRAS

We report on slow phase variations of the response of the space-ground radio interferometer RadioAstron during observations of pulsar B0329+54. The phase variations are due to the ionosphere and clearly distinguishable from effects of interstellar scintillation. Observations were made in a frequency range of 316-332~MHz with the 110-m Green Bank Telescope and the 10-m RadioAstron telescope in 1-hour sessions on 2012 November 26, 27, 28, and 29 with progressively increasing baseline projections of about 60, 90, 180, and 240 thousand kilometres. Quasi-periodic phase variations of interferometric scintles were detected in two observing sessions with characteristic time-scales of 12 and 10 minutes and amplitudes of up to 6.9~radians. We attribute the variations to the influence of medium-scale Travelling Ionospheric Disturbances. The measured amplitude corresponds to variations in vertical total electron content in ionosphere of about $0.1\times10^{16}\, \mathrm{m}^{-2}$. Such variations would noticeably constrain the coherent integration time in VLBI studies of compact radio sources at low frequencies.

Alison Wong, Benjamin Pope, Louis Desdoigts, Peter Tuthill, Barnaby Norris, Chris Betters

The principal limitation in many areas of astronomy, especially for directly imaging exoplanets, arises from instability in the point spread function (PSF) delivered by the telescope and instrument. To understand the transfer function, it is often necessary to infer a set of optical aberrations given only the intensity distribution on the sensor - the problem of phase retrieval. This can be important for post-processing of existing data, or for the design of optical phase masks to engineer PSFs optimized to achieve high contrast, angular resolution, or astrometric stability. By exploiting newly efficient and flexible technology for automatic differentiation, which in recent years has undergone rapid development driven by machine learning, we can perform both phase retrieval and design in a way that is systematic, user-friendly, fast, and effective. By using modern gradient descent techniques, this converges efficiently and is easily extended to incorporate constraints and regularization. We illustrate the wide-ranging potential for this approach using our new package, Morphine. Challenging applications performed with this code include precise phase retrieval for both discrete and continuous phase distributions, even where information has been censored such as heavily-saturated sensor data. We also show that the same algorithms can optimize continuous or binary phase masks that are competitive with existing best solutions for two example problems: an Apodizing Phase Plate (APP) coronagraph for exoplanet direct imaging, and a diffractive pupil for narrow-angle astrometry. The Morphine source code and examples are available open-source, with a similar interface to the popular physical optics package Poppy.

Paolo Bonfini, Andreas Zezas, Matthew L. N. Ashby, Steven P. Willner, Alexandros Maragkoudakis, Konstantinos Kouroumpatzakis, Paul. H. Sell, Konstantinos Kovlakas

33 pages, 21 figures, published in MNRAS, 504, 3831-3861 (2021)

We constrain the mass distribution in nearby, star-forming galaxies with the Star Formation Reference Survey (SFRS), a galaxy sample constructed to be representative of all known combinations of star formation rate (SFR), dust temperature, and specific star formation rate (sSFR) that exist in the Local Universe. An innovative two-dimensional bulge/disk decomposition of the 2MASS/$K_{s}$-band images of the SFRS galaxies yields global luminosity and stellar mass functions, along with separate mass functions for their bulges and disks. These accurate mass functions cover the full range from dwarf galaxies to large spirals, and are representative of star-forming galaxies selected based on their infra-red luminosity, unbiased by AGN content and environment. We measure an integrated luminosity density $j$ = 1.72 $\pm$ 0.93 $\times$ 10$^{9}$ L$_{\odot}$ $h^{-1}$ Mpc$^{-3}$ and a total stellar mass density $\rho_{M}$ = 4.61 $\pm$ 2.40 $\times$ 10$^{8}$ M$_{\odot}$ $h^{-1}$ Mpc$^{-3}$. While the stellar mass of the \emph{average} star-forming galaxy is equally distributed between its sub-components, disks globally dominate the mass density budget by a ratio 4:1 with respect to bulges. In particular, our functions suggest that recent star formation happened primarily in massive systems, where they have yielded a disk stellar mass density larger than that of bulges by more than 1 dex. Our results constitute a reference benchmark for models addressing the assembly of stellar mass on the bulges and disks of local ($z = 0$) star-forming galaxies.

Konstantinos Kouroumpatzakis, Andreas Zezas, Alexandros Maragkoudakis, Steven P. Willner, Paolo Bonfini, Matthew L. N. Ashby, Paul H. Sell, Thomas H. Jarrett

20 pages, 21 Figures, accepted for publication in MNRAS

We present new H$\alpha$ photometry for the Star-Formation Reference Survey (SFRS), a representative sample of star-forming galaxies in the local Universe. Combining these data with the panchromatic coverage of the SFRS, we provide calibrations of H$\alpha$-based star-formation rates (SFRs) with and without correction for the contribution of [$\rm N_{^{II}}$] emission. We consider the effect of extinction corrections based on the Balmer decrement, infrared excess (IRX), and spectral energy distribution (SED) fits. We compare the SFR estimates derived from SED fits, polycyclic aromatic hydrocarbons, hybrid indicators such as 24 $\mu$m + H$\alpha$, 8 $\mu$m + H$\alpha$, FIR + FUV, and H$\alpha$ emission for a sample of purely star-forming galaxies. We provide a new calibration for 1.4 GHz-based SFRs by comparing to the H$\alpha$ emission, and we measure a dependence of the radio-to-H$\alpha$ emission ratio based on galaxy stellar mass. Active galactic nuclei introduce biases in the calibrations of different SFR indicators but have only a minimal effect on the inferred SFR densities from galaxy surveys. Finally, we quantify the correlation between galaxy metallicity and extinction.

Sonia Fornasier, Jules Bourdelle de Micas, Pedro H. Hasselmann, Van H. Hoang, Maria Antonietta Barucci, Holger Sierks

14 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

We investigated Wosret, a region located on the small lobe of the 67P/Churyumov-Gerasimenko comet subject to strong heating during the perihelion passage. This region includes Abydos, the final landing site of the Philae lander. We analyzed high-resolution images of the Wosret region acquired between 2015 and 2016 by the OSIRIS instrument on board the Rosetta spacecraft. We observed a few morphological changes in Wosret, related to local dust coating removal with an estimated depth of $\sim$ 1 m, along with the formation of a cavity measuring 30 m in length and 6.5 m in depth, for a total estimated mass loss of 1.2 $\times$ 10$^6$ kg. The spectrophotometry of the region is typical of medium-red regions of comet 67P, with spectral slope values of 15-16 \%/(100 nm) in pre-perihelion data acquired at phase angle 60$^o$. Wosret has a spectral phase reddening of 0.0546 $\times 10^{-4}$ nm$^{-1} deg^{-1}$, which is about a factor of 2 lower than what was determined for the nucleus northern hemisphere regions, possibly indicating a reduced surface micro-roughness due to the lack of widespread dust coating. A few tiny bright spots are observed. Morphological features such as "goosebumps" or clods are widely present and larger in size than similar features located in the big lobe. Compared to Anhur and Khonsu, two southern hemisphere regions in the big lobe which are also exposed to high insolation during perihelion, Wosret exhibits fewer exposed volatiles and less morphological variations due to activity events. Our analysis indicates that the small lobe has different physical and mechanical properties than the big one and a lower volatile content, at least in its uppermost layers. These results support the hypothesis that comet 67P originated from the merging of two distinct bodies in the early Solar System.

Supriya Hebbur Dayananda, Javier Trujillo Bueno, Ángel de Vicente, Tanausú del Pino Alemán

20 pages, 14 figures; accepted for publication in The Astrophysical Journal

The near-Earth space weather is driven by the quick release of magnetic free energy in the solar corona. Probing this extremely hot and rarified region of the extended solar atmosphere requires modeling the polarization of forbidden and permitted coronal lines. To this end, it is important to develop efficient codes to calculate the Stokes profiles that emerge from given three-dimensional (3D) coronal models, and this should be done taking into account the symmetry breaking produced by the presence of magnetic fields and non-radial solar wind velocities. We have developed such a tool with the aim of theoretically predicting and interpreting spectropolarimetric observations of the solar corona in permitted and forbidden lines. In this paper we show the results of a theoretical investigation of the linear polarization signals produced by scattering processes in the H I Ly-${\alpha}$ line at 1216 A and in the He II Ly-${\alpha}$ line at 304 A using 3D coronal models by Predictive Science Inc. These spectral lines have very different critical magnetic fields for the onset of the Hanle effect (53 G and 850 G, respectively), as well as different sensitivities to the Doppler effect caused by the solar wind velocities. We study under which circumstances simultaneous observations of the scattering polarization in these Ly-${\alpha}$ lines can facilitate the determination of magnetic fields and macroscopic velocities in the solar corona.

Fernanda Duplancic, Diego G. Lambas, Sol Alonso, Georgina V. Coldwell

11 pages, 10 figures, accepted for publication in MNRAS

We perform a comparative study of AGNs in pairs, triplets and groups. To this end we use the Duplancic et al. catalogue of small galaxy systems and consider BPT and WHAN diagnostic diagrams to select optical AGNs. Also we identify mid-IR AGNs by using WISE data. We performed a comparison between the different AGN classification methodologies and study the AGN fraction in pairs, triplets, and groups with four to six members. We also analyse the main properties of Optical and mid-IR AGN hosts and the influence of environment on the active nuclei phenomena in these small galaxy systems. Our results show that, regardless the specifically adopted classification scheme, the fraction of AGN in pairs and triplets is always higher than the corresponding fraction in groups. Moreover, the fraction of powerful AGNs in pair and triplets is about twice the fraction of regular AGNs. We also find a remarkable difference between Optical and mid-IR AGNs in groups, where host galaxies of WISE AGNs are less massive and concentrated, with young stellar populations and blue colours. Also all WISE AGNs in groups have a very close companions and reside in an intermediate global density environment. Galaxy triplets show a larger AGN fraction for galaxies with a close nearest neighbours, while pairs present a nearly constant AGN fraction regardless the distance to the nearest companion. Our studies highlight the important role of interactions, besides the global environment dependence, in the activation of the AGN phenomenon in small galaxy systems.

Hayato Uehara, Kazuhito Dobashi, Shingo Nishiura, Tomomi Shimoikura, Takahiro Naoi

Accepted by The Astrophysical Journal, 21 pages, 11 figures, 6 tables

We have studied optical properties of interstellar dust around the Orion A molecular cloud to investigate the size distribution and the composition of dust grains. Orion A is one of the most studied molecular clouds in the solar vicinity ($d \simeq 400\ \rm{pc}$). In this paper, we used optical and near-infrared photometric data. The optical data were obtained by $BVRI$ bands imaging observations. The near-infrared data consisting of $JHK_{S}$ bands were taken from 2MASS point source catalog. We produced some color excess maps around Orion A, and measured their ratios such as $E(R-I)/E(B-V)$. In order to investigate dust properties, we compared the observed ratios with results of simulation performed by Naoi T. et al. (2021) who calculated the extinction in the optical to near-infrared wavelengths based on a standard dust model; they assumed a power-law grain-size distribution with an upper cutoff radius and assumed $graphite$ and $silicate$ as dominant components. As a result, we found that the upper cutoff radius around Orion A is $\simeq 0.3$ $\rm{\mu}$m, and $silicate$ predominates compared with $graphite$ (with the fraction of $silicate$ grater than $93$ %). In addition, we further derived the total-to-selective extinction ratio $R_{V}$ from the observed extinction of $A_{V}$ and the color excess $E(B-V)$, and compared it with the model calculations. Dust properties (i.e., the upper cutoff radius and the ratio of $graphite/silicate$) derived from $R_{V}$ is almost consistent with those derived from the color excess ratios.

Rafaela Gsponer, Johannes Noller

17 pages, 9 figures

A number of stability criteria exist for dark energy theories, associated with requiring the absence of ghost, gradient and tachyonic instabilities. Tachyonic instabilities are the least well explored of these in the dark energy context and we here discuss and derive criteria for their presence and size in detail. Our findings suggest that, while the absence of ghost and gradient instabilities is indeed essential for physically viable models and so priors associated with the absence of such instabilities significantly increase the efficiency of parameter estimations without introducing unphysical biases, this is not the case for tachyonic instabilities. Even strong such instabilities can be present without spoiling the cosmological validity of the underlying models. Therefore, we caution against using exclusion priors based on requiring the absence of (strong) tachyonic instabilities in deriving cosmological parameter constraints. We illustrate this by explicitly computing such constraints within the context of Horndeski theories, while quantifying the size and effect of related tachyonic instabilities.

Tania M. Barone, Francesco D'Eugenio, Nicholas Scott, Matthew Colless, Sam P. Vaughan, Arjen van der Wel, Amelia Fraser-McKelvie, Anna de Graaff, Jesse van de Sande, Po-Feng Wu, Rachel Bezanson, Sarah Brough, Eric Bell, Scott M. Croom, Luca Cortese, Simon Driver, Anna R. Gallazzi, Adam Muzzin, David Sobral, Joss Bland-Hawthorn, Julia J. Bryant, Michael Goodwin, Jon S. Lawrence, Nuria P. F. Lorente, Matt S. Owers

15 pages, 8 figures, submitted to MNRAS

We investigate the change in mean stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift ($0.60\leq z\leq0.76$) using the LEGA-C Survey, to low redshift ($0.014\leq z\leq0.10$) using the SAMI Galaxy Survey. We find that, similarly to their low-redshift counterparts, the stellar metallicity of quiescent galaxies at $0.60\leq z\leq 0.76$ closely correlates with $M_*/R_\mathrm{e}$ (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals, by determining the escape velocity required by metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. On the other hand, we find no correlation between stellar age and $M_*/R_\mathrm{e}^2$ (stellar mass surface density $\Sigma$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$\Sigma$ relation in the context of the redshift evolution of the star-forming and quiescent populations in the mass--size plane, and find our results can be explained as a consequence of galaxies forming more compactly at higher redshifts, and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The $z\sim 0$ age--$\Sigma$ relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and so a range of surface densities.

J. Malfait, W. Homan, S. Maes, J. Bolte, L. Siess, F. De Ceuster, L. Decin

23 pages, 22 figures

The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularly, there is a lack of studies investigating the structure formation in eccentric systems. We aim to improve our understanding of the wind morphology of eccentric AGB binary systems by investigating the mechanism responsible for the different small-scale structures and global morphologies that arise in a polytropic wind with different velocities. Using the smoothed particle hydrodynamics (SPH) code Phantom, we generate nine different high-resolution, 3D simulations of an AGB star with a solar-mass companion with various wind velocity and eccentricity combinations. The models assume a polytropic gas, with no additional cooling. We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, that is flattened. In the case of a lower wind velocity, the stronger interaction results in the formation of a high-energy region and bow-shock structure that can shape the wind into an irregular morphology if instabilities arise. High-eccentricity models show a complex, phase-dependent interaction leading to wind structures that are irregular in three dimensions. However, the significant interaction with the companion compresses matter into an equatorial density enhancement, irrespective of eccentricity.

Valeri V. Makarov, Norbert Zacharias, Charles T. Finch

published in RNAAS

The USNO ground-based astrometric program URAT-Bright in combination with the Hipparcos mission epoch astrometry provides precise proper motions of a thousand bright stars in the southern hemisphere on a time basis of about 25 years. Small but statistically significant differences between these proper motions and Gaia EDR3 data can reveal long-period exoplanets similar to Jupiter in the nearest star systems. The presence of such a planet orbiting the magnetically active dwarf $\epsilon$ Eri is confirmed from both URAT--Hipparcos--EDR3 data and Hipparcos--EDR3 data with a corresponding projected velocity of $(+5,+8)$ and $(+6,+13)$ m s$^{-1}$, respectively. These signals are formally significant at a 0.989 and 1.0 confidence. We conclude that the newest astrometric results confirm the existence of a long-period exoplanet orbiting $\epsilon$ Eri, which was marginally detected from precision radial velocity measurements some 20 years ago.

R. Abuter, A. Amorim, M. Bauböck, F. Baganoff, J.P. Berge, H. Boyce, H. Bonnet, W. Brandner, Y. Clénet, R. Davies, P.T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, A. Eckart, F. Eisenhauer, G.G. Fazio, N.M. Förster Schreiber, K. Foster, C. Gammie, P. Garcia, F. Gao, E. Gendron, R. Genzel, G. Ghisellini, S. Gillessen, M.A. Gurwell, M. Habibi, D. Haggard, C. Hailey, F. A. Harrison, X. Haubois, G. Heißel, T. Henning, S. Hippler, J.L. Hora, M. Horrobin, A. Jiménez-Rosales, L. Jochum, L. Jocou, A. Kaufer, P. Kervella, S. Lacour, V. Lapeyrère, J.-B. Le Bouquin, P. Léna, P.J. Lowrance, D. Lutz, S. Markoff, K. Mori, M.R. Morris, J. Neilsen, M. Nowak, T. Ott, T. Paumard, K. Perraut, G. Perrin, G. Ponti, O. Pfuhl, S. Rabien, G. Rodríguez-Coira, J. Shangguan, et al. (22 additional authors not shown)

accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $\nu L_\nu\propto \nu^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $\nu L_\nu \propto \nu^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $\gamma_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $\gamma_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $\gamma_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $\gamma_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare.

Marcin Hajduk, Marijke Haverkorn, Timothy Shimwell, Mateusz Olech, Joseph R. Callingham, Harish K. Vedantham, Glenn J. White, Marco Iacobelli, Alexander Drabent

20 pages, 10 figures, ApJ, accepted

We present observations of planetary nebulae with the LOw Frequency ARray (LOFAR) between 120 and 168 MHz. The images show thermal free-free emission from the nebular shells. We have determined the electron temperatures for spatially resolved, optically thick nebulae. These temperatures are 20 to 60% lower than those estimated from collisionally excited optical emission lines. This strongly supports the existence of a cold plasma component, which co-exists with hot plasma in planetary nebulae. This cold plasma does not contribute to the collisionally excited lines, but does contribute to recombination lines and radio flux. Neither of the plasma components are spatially resolved in our images, although we infer that the cold plasma extends to the outer radii of planetary nebulae. However, more cold plasma appears to exist at smaller radii. The presence of cold plasma should be taken into account in modeling of radio emission of planetary nebulae. Modelling of radio emission usually uses electron temperatures calculated from collisionally excited optical and/or infrared lines. This may lead to an underestimate of the ionized mass and an overestimate of the extinction correction from planetary nebulae when derived from the radio flux alone. The correction improves the consistency of extinction derived from the radio fluxes when compared to estimates from the Balmer decrement flux ratios.

Á. Kereszturi, V. Barta, I. Bondár, Cs. Czanik, A. Igaz, P. Mónus, D. Rezes, L. Szabados, B. D. Pál

Accepted for publication in MNRAS

Joint evaluation of different meteor observation types support the better understanding of both the meteor phenomenon and the terrestrial atmosphere. Two types of examples are presented in this work, linking ionospheric effects to specific meteors, where almost one third of the meteors emerged at high altitudes were simultaneously recorded with an optical camera. Very few such observations have been realized yet. With a daytime fireball, the recorded infrasound effect and the atmospheric blast produced shock wave related small earthquakes were identified by a network of ground stations. The overview of these observational types highlights specific topics where substantial improvements and discoveries are expected in the near future.

David Abarca, Kyle Parfrey, Włodek Kluźniak

9 pages, 4 figures, Submitted to ApJ Letters

We perform a global 2.5D general-relativistic radiation magnetohydrodynamic simulation of super-critical accretion onto a neutron star with a $2\times 10^{10}$ G dipolar magnetic field, as a model of a neutron-star-powered ultraluminous X-ray source (ULX). We compute a lower limit on the total luminosity of $\sim 2.5\,L_\mathrm{Edd}$, and find the radiation to be highly beamed by the accretion disk outflows. The apparent isotropic luminosity, which is a function of the viewing angle, reaches a maximum above $100\,L_\mathrm{Edd}$, consistent with the luminosities observed in ULXs.

Hideki Umehata, Ian Smail, Charles C. Steidel, Matthew Hayes, Douglas Scott, A. M. Swinbank, R.J. Ivison, Toru Nagao, Mariko Kubo, Kouichiro Nakanishi, Yuichi Matsuda, Soh Ikarashi, Yoichi Tamura, J. E. Geach

23 pages, 11 figures, accepted for publication in ApJ

We present observations of a giant Lyman-alpha blob in the SSA22 proto-cluster at z=3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Dust continuum, along with [C II]158um, and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z=3.0987-3.1016 in the surroundings are identified in [C II] and dust continuum emission, with two of them potential companions or tidal structures associated with the most massive galaxies. Spatially resolved [C II] and infrared luminosity ratios for the widely distributed media (L[C II]/LIR~0.01-0.001) suggest that the observed extended interstellar media are likely to have originated from star-formation activity and the contribution from shocked gas is probably not dominant. LAB1 is found to harbour a total molecular gas mass Mmol=(8.7+/-2.0)e+10 Msun, concentrated in the core region of the Ly-alpha-emitting area. While (primarily obscured) star-formation activity in the LAB1 core is one of the most plausible power sources for the Ly-alpha emission, multiple major-mergers found in the core may also play a role in making LAB1 exceptionally bright and extended in Ly-alpha as a result of cooling radiation induced by gravitational interactions.

We examine the evolution of intragroup gas X-ray scaling relations for group-sized halos ($M_{500}=10^{12.3-15}M_{\odot}$) in the Simba galaxy formation simulation. X-ray luminosity $L_X$ vs $M_{500}$ shows increasing deviation from self-similarity from $z=3\to 0$, with $M_{500}<10^{13.5} M_{\odot}$ halos exhibiting a large reduction in $L_X$ and slight increase in X-ray luminosity-weighted temperature $T_X$. These shifts are driven by a strong drop in $f_{\rm gas}$ with time for these halos, and coincides with the onset of black hole jet feedback in these systems at $z\sim 1.5$ in Simba. The connection with black hole feedback is corroborated by $f_{BH}\equiv M_{BH}/M_{500}$ in $M_{500}<10^{13.5} M_{\odot}$ halos being strongly anti-correlated with $L_X$ and $f_{\rm gas}$ at $z\la 1.5$. This is further reflected in the scatter of $L_X-T_X$: halos with small $f_{BH}$ lie near self-similarity, while those with the highest $f_{BH}$ lie furthest below. Turning off jet feedback results in mostly self-similar behaviour down to $z=0$. For the X-ray weighted metallicity $Z_X$, stellar feedback impacts the enrichment of halo gas. Finally, halo profiles show that jet feedback flattens the electron density and entropy profiles, and introduces a core in X-ray surface brightness particularly at $M_{500}<10^{13.5} M_{\odot}$. This argues that intragroup X-ray evolution is largely driven by jet feedback removing hot gas from the cores of massive groups, and expelling gas altogether in less massive groups.

Lisa Benato, Erik Buhmann, Martin Erdmann, Peter Fackeldey, Jonas Glombitza, Nikolai Hartmann, Gregor Kasieczka, William Korcari, Thomas Kuhr, Jan Steinheimer, Horst Stöcker, Tilman Plehn, Kai Zhou

13 pages, 5 figures, 5 tables

We introduce a collection of datasets from fundamental physics research -- including particle physics, astroparticle physics, and hadron- and nuclear physics -- for supervised machine learning studies. These datasets, containing hadronic top quarks, cosmic-ray induced air showers, phase transitions in hadronic matter, and generator-level histories, are made public to simplify future work on cross-disciplinary machine learning and transfer learning in fundamental physics. Based on these data, we present a simple yet flexible graph-based neural network architecture that can easily be applied to a wide range of supervised learning tasks in these domains. We show that our approach reaches performance close to state-of-the-art dedicated methods on all datasets. To simplify adaptation for various problems, we provide easy-to-follow instructions on how graph-based representations of data structures, relevant for fundamental physics, can be constructed and provide code implementations for several of them. Implementations are also provided for our proposed method and all reference algorithms.

Mohsen Khodadi, Gaetano Lambiase, David F. Mota

31 pages, 14 figures

Thanks to the release of the extraordinary EHT image of shadow attributed to the M87* supermassive black hole (SMBH), we have a novel window to assess the validity of fundamental physics in the strong-field regime. Motivated by this, we consider Johannsen \& Psaltis metric parameterized by mass, spin, and an additional dimensionless hair parameter $\epsilon$ which in the high rotation regimes is able to provide a suitable framework for the test of the no-hair theorem (NHT) using the EHT data. Incorporating the $\epsilon$ into the standard Kerr spacetime enrich it in the sense that, depending on setting the positive and negative values for that, we deal with alternative compact objects: deformed Kerr naked singularity and Kerr BH solutions, respectively. Shadows associated with these two possible solutions indicate that the deformation parameter $\epsilon$ affects the geometry shape of standard shadow such that it become more oblate and prolate with $\epsilon<0$ and $\epsilon>0$, respectively. By scanning the window associated with three shadow observables: oblateness, deviation from circularity, and shadow diameter, we perform a numerical analysis within the range $a_*=0.9\mp0.1$ of the dimensionless rotation parameter, to find the constraints on the hair parameter $\epsilon$ in both possible solutions. For both possible signs of $\epsilon$, we extract a variety of upper bounds that are in interplay with $a_*$. Our analysis suggests that as the rotation parameter approaches the extreme limit, although the allowable range of both hair parameters becomes narrower, the hairy Kerr BH solution is a more promising candidate to play a role of the alternative compact object instead of standard Kerr BH. The lack of tension between hairy Kerr BH with the current observation of the EHT shadow of the M87* SMBH carries this message that the possibility of NHT violation is not excluded.

Fermion dark matter particles can aggregate to form extended dark matter structures via a first-order phase transition in which the particles get trapped in the false vacuum. We study Fermi balls created in a phase transition induced by a generic quartic thermal effective potential. We show that for Fermi balls of mass, $3\times 10^{-12}M_\odot \lesssim M_{\rm FB} \lesssim 10^{-5}M_\odot$, correlated observations of gravitational waves produced during the phase transition (at SKA/THEIA), and gravitational microlensing caused by Fermi balls (at Subaru-HSC), can be made.

We study the behaviour of Yukawa and Newtonian gravitational forces in a cubic box with fully periodic boundaries commonly encountered in N-body simulations of the structure formation. Placing a single gravitating body at the origin of coordinates, we reveal the scales at which non-negligible deviation from the Yukawa law occurs when the Newtonian approximation is employed. We discuss the results in terms of the corresponding physical distances today as well as earlier, back at the matter-dominated stage. Revisiting the problem for free boundaries, we also compare the periodic and plain gravitational forces for Yukawa-type interactions.

Thomas Celora, Nils Andersson, Ian Hawke, Gregory L. Comer

24 pages, 1 figure, comments are welcome

Models of turbulent flows require the resolution of a vast range of scales, from large eddies to small-scale features directly associated with dissipation. As the required resolution is not within reach of large scale numerical simulations, standard strategies involve a smoothing of the fluid dynamics, either through time averaging or spatial filtering. These strategies raise formal issues in general relativity, where the split between space and time is observer dependent. To make progress, we develop a new covariant framework for filtering/averaging based on the fibration of spacetime associated with fluid elements and the use of Fermi coordinates to facilitate a meaningful local analysis. We derive the resolved equations of motion, demonstrating how "effective" dissipative terms arise because of the coarse-graining, and pay particular attention to the thermodynamical interpretation of the resolved quantities. Finally, as the smoothing of the fluid-dynamics inevitably leads to a closure problem, we propose a new closure scheme inspired by recent progress in the modelling of dissipative relativistic fluids, and crucially, demonstrate the linear stability of the proposed model.

D Alesini, D Babusci, C Barone, B Buonomo, M M Beretta, L Bianchini, G Castellano, F Chiarello, D Di Gioacchino, P Falferi, G Felici, G Filatrella, L G Foggetta, A Gallo, C Gatti, F Giazotto, G Lamanna, F Ligabue, N Ligato, C Ligi, G Maccarrone, B Margesin, F Mattioli, E Monticone, L Oberto, S Pagano, F Paolucci, M Rajteri, A Rettaroli, L Rolandi, P Spagnolo, A Toncelli, G Torrioli

7 pages, 9 figures

Josephson junctions, in appropriate configurations, can be excellent candidates for detection of single photons in the microwave frequency band. Such possibility has been recently addressed in the framework of galactic axion detection. Here are reported recent developments in the modelling and simulation of dynamic behaviour of a Josephson junction single microwave photon detector. For a Josephson junction to be enough sensitive, small critical currents and operating temperatures of the order of ten of mK are necessary. Thermal and quantum tunnelling out of the zero-voltage state can also mask the detection process. Axion detection would require dark count rates in the order of 0.001 Hz. It is, therefore, is of paramount importance to identify proper device fabrication parameters and junction operation point.

Motivated by electromagnetic-field confinement due to plasma near accreting black holes, we continue our exploration of the linear dynamics of an electromagnetic field propagating in curved spacetime in the presence of plasma by including three effects that were neglected in our previous analysis: collisions in the plasma, thermal corrections, and the angular momentum of the background black-hole spacetime. We show that: (i) the plasma-driven long-lived modes survive in a collisional plasma except when the collision timescale is unrealistically small; (ii) thermal effects, which might be relevant for accretion disks around black holes, do not affect the axial long-lived modes; (iii) in the case of a spinning black hole the plasma-driven modes become superradiantly unstable at the linear level; (iv) the polar sector in the small-frequency regime admits a reflection point due to the resonant properties of the plasma. Dissipative effects such as absorption, formation of plasma waves, and nonlinear dynamics play a crucial role in the vicinity of this resonant point.

Yu Seon Jeong, Weidong Bai, Milind Diwan, Maria Vittoria Garzelli, Fnu Karan Kuma, Mary Hall Reno

8 pages, 6 figures. Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021)

Theoretical predictions of the prompt atmospheric neutrino flux have large uncertainties associated with charm hadron production, by far the dominant source of prompt neutrinos in the atmosphere. The flux of cosmic rays, with its steeply falling energy spectrum, weights the forward production of charm in the evaluation of the atmospheric neutrino flux at high energies. The current LHCb experiment at CERN constrains charm production in kinematic regions relevant to the prompt atmospheric neutrino flux. The proposed Forward Physics Facility has additional capabilities to detect neutrino fluxes from forward charm production at the LHC. We discuss the implications of the current and planned experiments on the development of theoretical predictions of the high energy atmospheric neutrino flux.