Abstract visibility
Abstract text size

Papers for Wednesday, Jan 19 2022

Papers with local authors

Kareem El-Badry, Charlie Conroy, Eliot Quataert, Hans-Walter Rix, Jonathan Labadie-Bartz, Tharindu Jayasinghe, Todd Thompson, Phillip Cargile, Keivan G. Stassun, Ilya Ilyin

25 pages, 19 figures, submitted to MNRAS

0 votes
Paper 2 — arXiv:2201.05614
0 votes
Paper 2 — arXiv:2201.05614

Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6 \,M_{\odot}$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\,M_{\odot}$), subgiant companion on a 5.47-day orbit. The emission lines originate in an accretion disk caused by ongoing mass transfer, not from a decretion disk as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor's core: the surface helium fraction is $Y_{\rm He}\approx 0.6$, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system's properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\,M_{\odot}$ donor leaves the main sequence, with the secondary becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdOB star. Discovery of one object in this short-lived ($\sim$1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdOB binaries. We infer that $(28_{-16}^{+27})\,\%$ of Be stars have stripped companions, most of which are faint. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdOB binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars.

Xin Li, Jian Li, Zhihong Jeff Xia, Nikolaos Georgakarakos

12 pages, 9 figures, accepted for pubblication on Monthly Notices of the Royal Astronomical Society

0 votes
Paper 73 — arXiv:2201.06743
0 votes
Paper 73 — arXiv:2201.06743

Most recently, machine learning has been used to study the dynamics of integrable Hamiltonian systems and the chaotic 3-body problem. In this work, we consider an intermediate case of regular motion in a non-integrable system: the behaviour of objects in the 2:3 mean motion resonance with Neptune. We show that, given initial data from a short 6250 yr numerical integration, the best-trained artificial neural network (ANN) can predict the trajectories of the 2:3 resonators over the subsequent 18750 yr evolution, covering a full libration cycle over the combined time period. By comparing our ANN's prediction of the resonant angle to the outcome of numerical integrations, the former can predict the resonant angle with an accuracy as small as of a few degrees only, while it has the advantage of considerably saving computational time. More specifically, the trained ANN can effectively measure the resonant amplitudes of the 2:3 resonators, and thus provides a fast approach that can identify the resonant candidates. This may be helpful in classifying a huge population of KBOs to be discovered in future surveys.

Lin Wang, Ye-Zhao Yu, Feifei Kou, Kuo Liu, Xinxin Wang, Bo Peng
0 votes
Paper 81 — arXiv:2201.06815
0 votes
Paper 81 — arXiv:2201.06815

We present simultaneous broad-band radio observations on the abnormal emission mode from PSR B1859$+$07 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). This pulsar shows peculiar emission phenomena, which are occasional shifts of emission to an early rotational phase and mode change of emission at the normal phase. We confirm all these three emission modes with our datasets, including the B (burst) and Q (quiet) modes of the non-shift pulses and the emission shift mode with a quasi-periodicity of 155 pulses. We also identify a new type of emission shift event, which has emission at the normal phase during the event. We studied polarisation properties of these emission modes in details, and found that they all have similar polarisation angle (PA) curve, indicating the emission of all these three modes are from the same emission height.

Jie Zheng, Yun Chen, Tengpeng Xu, Zong-Hong Zhu

9 pages, 3 figures, 3 tables. arXiv admin note: text overlap with arXiv:2107.08916

0 votes
Paper 93 — arXiv:2201.07011
0 votes
Paper 93 — arXiv:2201.07011

We explore and compare the capabilities of the recent observations of standard cosmological probes and the future observations of gravitational-wave (GW) standard sirens on constraining cosmological parameters. It is carried out in the frameworks of two typical dynamical models of cosmology, i.e., the $\omega_0\omega_a$CDM model with $\omega(z) = \omega_0 +\omega_a*z/(1+z)$, and the $\xi$-index model with $\rho_X\propto\rho_ma^{\xi}$, where $\omega(z)$ is the dark energy equation of state, and $\rho_X$ and $\rho_m$ are the energy densities of dark energy and matter, respectively. In the cosmological analysis, the employed data sets include the recent observations of the standard cosmological probes, i.e., Type Ia supernovae (SNe Ia), baryon acoustic oscillation (BAO) and cosmic microwave background (CMB), and the mock GW standard siren sample with 1000 merging neutron star events anticipated from the third-generation detectors. In the scenarios of both $\omega_0\omega_a$CDM and $\xi$-index models, it turns out that the mock GW sample can reduce the uncertainty of the Hubble constant $H_0$ by about 50% relative to that from the joint SNe+BAO+CMB sample; nevertheless, the SNe+BAO+CMB sample demonstrates better performance on limiting other parameters. Furthermore, the Bayesian evidence is applied to compare the dynamical models with the $\Lambda$CDM model. The Bayesian evidences computed from the SNe+BAO+CMB sample reveal that the $\Lambda$CDM model is the most supported one; moreover, the $\omega_0\omega_a$CDM model is more competitive than the $\xi$-index model.

Papers with votes

Todd M. Tripp (University of Massachusetts - Amherst)

Accepted for publication in MNRAS

1 vote
Paper 3 — arXiv:2201.05615
1 vote
Paper 3 — arXiv:2201.05615

The high-velocity clouds (HVCs) in the outer Milky Way at $20^{\circ} < l < 190^{\circ}$ have similar spatial locations, metallicities, and kinematics. Moreover, their locations and kinematics are coincident with several extraplanar stellar streams. The HVC origins may be connected to the stellar streams, either stripped directly from them or precipitated by the aggregate dynamical roiling of the region by the stream progenitors. This paper suggests that these HVCs are "misty" precipitation in the stream wakes based on the following observations. New high-resolution (2.6 km/s) ultraviolet spectroscopy of the QSO H1821+643 resolves what appears to be a single HVC absorption cloud (at 7 km/s resolution) into five components with $T \lesssim 3\times 10^{4}$ K. Photoionization models can explain the low-ionization components but require some depletion of refractory elements by dust, and model degeneracies allow a large range of metallicity. High-ionization absorption lines (SiIV, CIV, and OVI) are kinematically aligned with the lower-ionization lines and cannot be easily explained with photoionization or equilibrium collisional ionization; these lines are best matched by non-equilibrium rapidly cooling models, i.e., condensing/precipitating gas, with high metallicity and a significant amount of HI. Both the low- and high-ionization phases have low ratios of cooling time to freefall time and cooling time to sound-crossing time, which enables fragmentation and precipitation. The H1821+643 results are corroborated by spectroscopy of six other nearby targets that likewise show kinematically correlated low- and high-ionization absorption lines with evidence of dust depletion and rapid cooling.

All other papers

Ying Qin, Xinwen Shu, Shuangxi Yi, Yuan-Zhu Wang

Comments welcome! Accepted for publication in Research in Astronomy and Astrophysics

Recent observations of AdLIGO and Virgo have shown that the spin measurements in binary black hole (BH) systems are typically small, which is consistent with the predictions by the classical isolated binary evolution channel. In this standard formation channel, the progenitor of the first-born BH is assumed to have efficient angular momentum transport. The BH spins in high-mass X-ray binaries (HMXBs), however, have been found consistently to be extremely high. In order to explain the high BH spins, the inefficient angular momentum transport inside the BH progenitor is required. This requirement, however, is incompatible with the current understanding of conventional efficient angular momentum transport mechanism. We find that this tension can be highly alleviated as long as the hypercritical accretion is allowed. We show that, for a case study of Cygnus X-1, the hypercritical accretion cannot only be a good solution for the inconsistent assumption upon the angular momentum transport within massive stars, but match its other properties reported recently.

Ruobing Dong, Hauyu Baobab Liu, Nicolas Cuello, Christophe Pinte, Peter Abraham, Eduard Vorobyov, Jun Hashimoto, Agnes Kospal, Eugene Chiang, Michihiro Takami, Lei Chen, Michael Dunham, Misato Fukagawa, Joel Green, Yasuhiro Hasegawa, Thomas Henning, Yaroslav Pavlyuchenkov, Tae-Soo Pyo, Motohide Tamura

Published in Nature Astronomy. Here is the authors' version with the Supplementary Information integrated into the Methods section

Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the million year old binary protostar Z CMa in dust continuum and molecular line emission. A point source ~4700 au from the binary has been discovered at both millimeter and centimeter wavelengths. It is located along the extension of a ~2000 au streamer structure previously found in scattered light imaging, whose counterpart in dust and gas emission is also newly identified. Comparison with simulations shows signposts of a rare flyby event in action. Z CMa is a "double burster", as both binary components undergo accretion outbursts, which may be facilitated by perturbations to the host disk by flybys.

Alessia Franchini, Alessandro Lupi, Alberto Sesana

8 pages, submitted to ApJL on December 6th 2021

The study of the interaction of a massive black hole binary with its gaseous environment is crucial in order to be able to predict merger rates and possible electromagnetic counterparts of gravitational wave signals. The evolution of the binary semi-major axis resulting from this interaction has been recently debated, and a clear consensus is still missing, also because of several numerical limitations, i.e. fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. We use the 3D meshless finite mass method of the code GIZMO with Lagrangian hyper-refinement, to resolve for the first time the dynamics inside the cavity, and in particular the mini-discs that form around the two components of a live binary surrounded by a locally isothermal gaseous circumbinary disc. We show that the binary orbit decays with time for very cold and very warm discs and that the result of the interaction in the intermediate regime is strongly influenced by the disc viscosity as this essentially regulates the fraction of mass contained in the mini-discs as well as the fraction that is accreted by the binary. We find the balance between these two quantities to determine whether the binary semi-major axis decreases with time.

Sayantan Auddy, Shantanu Basu, Takahiro Kudoh

6 Pages, 3 Figures, submitted to ApJL

We study the magnetic field to density ($B-\rho$) relation in turbulent molecular clouds with dynamically important magnetic fields using nonideal three-dimensional magnetohydrodynamic simulations. Our simulations show that there is a distinguishable break density $\rho_{\rm T}$ between the relatively flat low density regime and a power-law regime at higher densities. We present an analytic theory for $\rho_{\rm T}$ based on the interplay of the magnetic field, turbulence, and gravity. The break density $\rho_{\rm T}$ scales with the strength of the initial Alfv\'en Mach number $\mathcal{M}_{\rm A0}$ for sub-Alfv\'enic ( $\mathcal{M}_{\rm A0}<1$) and trans-Alfv\'enic ($\mathcal{M}_{\rm A0} \sim 1$) clouds. We fit the variation of $\rho_{\rm T}$ for model clouds as a function of $\mathcal{M}_{\rm A0}$, set by different values of initial sonic Mach number $\mathcal{M_{\rm 0}}$ and the initial ratio of gas pressure to magnetic pressure $\beta_{\rm 0}$. This implies that $\rho_{\rm T}$, which denotes the transition in mass-to-flux ratio from the subcritical to supercritical regime, is set by the initial turbulent compression of the molecular cloud.

Foteini Oikonomou

16 pages, 7 figures, invited talk at the ICRC 2021

Active galactic nuclei (AGN) with relativistic jets are the most powerful persistent astrophysical sources of electromagnetic radiation in the Universe. Blazars are the most extreme subclass of AGN with jets directed along the line of sight of the observer. Their high-energy photon emission dominates the extragalactic gamma-ray sky and reaches multi-TeV energies. This demonstrates that blazars accelerate particles to very high energies. It has long been suspected that blazars may also accelerate protons to very high energies and thus be cosmic neutrino sources. Being extremely rare objects in addition to being bright, blazars are among the most readily testable neutrino candidate source classes. Several multi-messenger monitoring campaigns have recently been triggered in response to high-energy neutrinos observed with the IceCube Neutrino Observatory from the direction of blazars. In this contribution, I summarise the theoretical interpretation of these observations and give an overview of the possible role of blazars as neutrino sources in light of the experimental results.

Michael Ryan, David Radice

11 pages, 5 figures

Several dark matter models allow for the intriguing possibility of exotic compact object formation. These objects might have unique characteristics that set them apart from their baryonic counterparts. Furthermore, gravitational wave observations of their mergers may provide the only direct window on a potentially entirely hidden sector. Here we discuss dark white dwarfs, starting with an overview of the microphysical model and analytic scaling relations of macroscopic properties derived from the non-relativistic limit. We use the full relativistic formalism to confirm these scaling relations and demonstrate that dark white dwarfs, if they exist, would have masses and tidal deformabilities that are very different from those of baryonic compact objects. Further, and most importantly, we demonstrate that dark white dwarf mergers would be detectable by current or planned gravitational observatories across several orders of magnitude in the particle-mass parameter space. Lastly, we find universal relations analogous to the compactness-Love and binary Love relations in neutron star literature. Using these results, we show that gravitational wave observations would constrain the properties of the dark matter particles constituting these objects.

Jordan Simpson, David Jones, Roger Wesson, Jorge García-Rojas

3 pages, 1 figure

A30 belongs to a class of planetary nebulae identified as "born-again", containing dense, hydrogen-poor ejecta with extreme abundance discrepancy factors (ADFs), likely associated with a central binary system. We present intermediate-dispersion spectroscopy of one such feature-the J4 equatorial knot. We confirm the apparent physical and chemical segregation of the polar and equatorial knots observed in previous studies, and place an upper limit on the ADF for O$^{2+}$ of 35, significantly lower than that of the polar knots. These findings further reinforce the theory that the equatorial and polar knots originate from different events.

Alberto Roper Pol, Chiara Caprini, Andrii Neronov, Dmitri Semikoz

25 pages, 9 figures, 1 table. Submitted to Phys. Rev. D

The NANOGrav, Parkes, and European pulsar timing array (PTA) collaborations have reported evidence for a common-spectrum process that can potentially correspond to a stochastic gravitational wave background (SGWB) in the 1--100 nHz frequency range. We consider the scenario in which this signal is produced by magnetohydrodynamic (MHD) turbulence in the early universe, induced by a non-helical primordial magnetic field at the energy scale corresponding to the quark confinement phase transition. We perform MHD simulations to study the dynamical evolution of the magnetic field and compute the resulting SGWB. We show that the SGWB output from the simulations can be very well approximated by assuming that the magnetic anisotropic stress is constant in time, over a time interval related to the eddy turnover time. The analytical spectrum that we derive under this assumption features a change of slope at a frequency corresponding to the GW source duration that we confirm with the numerical simulations. We compare the SGWB signal with the PTA data to constrain the temperature scale at which the SGWB is sourced, as well as the amplitude and characteristic scale of the initial magnetic field. We find that the generation temperature is constrained to be in the 2-200 MeV range, the magnetic field amplitude must be $>1$\% of the radiation energy density at that time, and the magnetic field characteristic scale is constrained to be $>10$\% of the horizon scale. We show that the turbulent decay of this magnetic field will lead to a field at recombination that can help to alleviate the Hubble tension and can be tested by measurements in the voids of the Large Scale Structure with gamma-ray telescopes like the Cherenkov Telescope Array.

F. Capel, J. M. Burgess, D. J. Mortlock, P. Padovani

11 pages, 10 figures, submitted to A&A

Several marginally significant associations between high-energy neutrinos and potential astrophysical sources have been recently reported, but a conclusive identification of these sources remains challenging. We explore the use of Monte Carlo simulations to gain deeper insight into the implications of, in particular, the IC170922A-TXS 0506+056 observation. Assuming a null model, we find a 7.6% chance to mistakenly identify coincidences between flaring blazars and neutrino alerts in 10-year surveys. We confirm that a blazar-neutrino connection based on the ${\gamma}$-ray flux is required to find a low chance coincidence probability and, therefore, a significant IC170922A-TXS 0506+056 association. We then assume this blazar-neutrino connection for the whole population and find that the ratio of neutrino to ${\gamma}$-ray fluxes must be $\lesssim 10^{-2}$ in order not to overproduce the total number of neutrino alerts seen by IceCube. For the IC170922A-TXS 0506+056 association to make sense, we must either accept this low flux ratio or suppose that only some rare sub-population of blazars is capable of high-energy neutrino production. For example, if we consider neutrino production only in blazar flares, we expect the flux ratio of between $10^{-3}$ and $10^{-1}$ to be consistent with a single coincident observation of a neutrino alert and flaring blazar. These conclusions are robust with respect to the uncertainties in our modelling assumptions.

Claire Marie Guimond, John Rudge, Oliver Shorttle

33 pages, 9 figures, submitted to The Planetary Science Journal

Topography on a wet rocky exoplanet could raise land above its sea level. Although land elevation is the product of many complex processes, the large-scale topographic features on any geodynamically-active planet are the expression of the convecting mantle beneath the surface. This so-called "dynamic topography" exists regardless of a planet's tectonic regime or volcanism; its amplitude, with a few assumptions, can be estimated via numerical simulations of convection as a function of the mantle Rayleigh number. We develop new scaling relationships for dynamic topography on stagnant lid planets using 2D convection models with temperature-dependent viscosity. These scalings are applied to 1D thermal history models to explore how dynamic topography varies with exoplanetary observables over a wide parameter space. Dynamic topography amplitudes are converted to an ocean basin capacity, the minimum water volume required to flood the entire surface. Basin capacity increases less steeply with planet mass than does the amount of water itself, assuming a water inventory that is a constant planetary mass fraction. We find that dynamically-supported topography alone could be sufficient to maintain subaerial land on Earth-size stagnant lid planets with surface water inventories of up to approximately $10^{-4}$ times their mass, in the most favourable thermal states. By considering only dynamic topography, which has ~1-km amplitudes on Earth, these results represent a lower limit to the true ocean basin capacity. Our work indicates that deterministic geophysical modelling could inform the variability of land propensity on low-mass planets.

Anna C. Childs, Jason H. Steffen

Accepted for publication in MNRAS. 12 pages, 7 figures, and 5 tables

We present a fragmentation module and a composition tracking code for the $n$-body code REBOUND. Our fragmentation code utilises previous semi-analytic models and follows an implementation method similar to fragmentation for the $n$-body code MERCURY. In our $n$-body simulations with fragmentation, we decrease the collision and planet formation timescales by inflating the particle radii by an expansion factor $f$ and experiment with various values of $f$ to understand how expansion factors affect the collision history and final planetary system. As the expansion factor increases, so do the rate of mergers which produces planetary systems with more planets and planets at larger orbits. Additionally, we present a composition tracking code which follows the compositional change of homogeneous bodies as a function of mass exchange and use it to study how fragmentation and the use of an expansion factor affects volatile delivery to the inner terrestrial disc. We find that fragmentation enhances radial mixing relative to perfect merging and that on average, as $f$ increases so does the average water mass fraction of the planets. Radial mixing decreases with increasing $f$ as collisions happen early on, before the bodies have time to grow to excited orbits and move away from their original location

D.S. Akerib, S. Alsum, H.M. Araújo, X. Bai, J. Balajthy, J. Bang, A. Baxter, E.P. Bernard, A. Bernstein, T.P. Biesiadzinski, E.M. Boulton, B. Boxer, P. Brás, S. Burdin, D. Byram, N. Carrara, M.C. Carmona-Benitez, C. Chan, J.E. Cutter, L. de Viveiros, E. Druszkiewicz, J. Ernst, A. Fan, S. Fiorucci, R.J. Gaitskell, C. Ghag, M.G.D. Gilchriese, C. Gwilliam, C.R. Hall, S.J. Haselschwardt, S.A. Hertel, D.P. Hogan, M. Horn, D.Q. Huang, C.M. Ignarra, R.G. Jacobsen, O. Jahangir, W. Ji, K. Kamdin, K. Kazkaz, D. Khaitan, E.V. Korolkova, S. Kravitz, V.A. Kudryavtsev, E. Leason, B.G. Lenardo, K.T. Lesko, J. Liao, J. Lin, A. Lindote, M.I. Lopes, A. Manalaysay, R.L. Mannino, N. Marangou, D.N. McKinsey, D.-M. Mei, J.A. Morad, A.St.J. Murphy, A. Naylor, C. Nehrkorn, H.N. Nelson, et al. (38 additional authors not shown)

We present the results from combining machine learning with the profile likelihood fit procedure, using data from the Large Underground Xenon (LUX) dark matter experiment. This approach demonstrates reduction in computation time by a factor of 30 when compared with the previous approach, without loss of performance on real data. We establish its flexibility to capture non-linear correlations between variables (such as smearing in light and charge signals due to position variation) by achieving equal performance using pulse areas with and without position-corrections applied. Its efficiency and scalability furthermore enables searching for dark matter using additional variables without significant computational burden. We demonstrate this by including a light signal pulse shape variable alongside more traditional inputs such as light and charge signal strengths. This technique can be exploited by future dark matter experiments to make use of additional information, reduce computational resources needed for signal searches and simulations, and make inclusion of physical nuisance parameters in fits tractable.

Marcos A. van Dam, Sylvain Oberti, Johann Kolb, Jim Lyke, Sylvain Cetre, Benoit Neichel

Infrared tip-tilt sensors (IR TTSs) have been deployed on three different laser guide star adaptive optics (AO) systems on three different telescopes. These IR TTS benefit from the high-order loop PSF sharpening in the near infrared, hence they provide a low tip-tilt residual and a good sky coverage. Nevertheless, these IR TTS are challenging and their use in AO is limited. In this paper, we outline existing IR TTS, provide on-sky performance results and describe our experience using IR TTS and along with plans for the near future. The second part of the paper deals with unresolved challenges for IR TTS. These include algorithms and loop stability in the low Strehl regime, using IR TTSs to measure higher-order modes and guiding on multiple guide stars with different magnitudes.

J. R. Dawson (1 and 2), George Hobbs (1), Yansong Gao (3), Seyit Camtepe (4), Josef Pieprzyk (4 and 5), Yi Feng (6 and 7), Luke Tranfa (1 and 2), Sarah Bradbury (1 and 8), Weiwei Zhu (6 and 7), Di Li (6 and 7 and 9). ((1) CSIRO Space and Astronomy, (2) Department of Physics and Astronomy and MQ Research Centre in Astronomy, Astrophysics and Astrophotonics, Macquarie University, Sydney (3) School of Computer Science and Engineering, Nanjing University of Science and Technology, (4) CSIRO Data 61, (5) Institute of Computer Science, Polish Academy of Sciences, (6) National Astronomical Observatories, Chinese Academy of Sciences, Beijing (7) University of Chinese Academy of Sciences, (8) School of Chemistry and Physics, Queensland University of Technology (QUT), (9) NAOC-UKZN Computational Astrophysics Centre, University of KwaZulu-Natal)

12 pages, 4 figures, accepted to Astronomy & Computing

We demonstrate how radio pulsars can be used as random number generators. Specifically, we focus on publicly verifiable randomness (PVR), in which the same sequence of trusted and verifiable random numbers is obtained by multiple parties. PVR is a critical building block for many processes and algorithms (including cryptography, scientific trials, electoral audits and international treaties). However, current approaches (based on number theory) may soon become vulnerable to quantum computers, motivating a growing demand for PVR based on natural physical phenomena. In this context, we explore pulsars as a potential physical PVR source. We first show that bit sequences extracted from the measured flux densities of a bright millisecond pulsar can pass standardised tests for randomness. We then quantify three illustrative methods of bit-extraction from pulsar flux density sequences, using simultaneous observations of a second pulsar carried out with the Parkes telescope in Australia and the Five-hundred-metre Aperture Spherical radio Telescope (FAST) in China, supported by numerical simulations. We demonstrate that the same bit sequence can indeed be obtained at both observatories, but the ubiquitous presence of radiometer noise needs to be accounted for when determining the expected bit error rate between two independent sequences. We discuss our results in the context of an imaginary use-case in which two mutually distrusting parties wish to obtain the same random bit sequence, exploring potential methods to mitigate against a malicious participant.

Margarita Safonova

5 pages, 2 figures, presented at First Belgo-Indian Network for Astronomy & Astrophysics (BINA) workshop, Nainital, India, 15-18 November 2016, this https URL

Gravitational lensing is brightening of background objects due to deflection of light by foreground sources. Rich clusters of galaxies are very effective lenses because they are centrally concentrated. Such natural Gravitational Telescopes provide us with strongly magnified galaxies at high redshifts otherwise too faint to be detected or analyzed. With a lensing boost, we can study galaxies shining at the end of the `Dark Ages'. We propose to exploit the opportunity provided by the large field of view and depth, to search for sources magnified by foreground clusters in the vicinity of the cluster critical curves, where enhancements can be of several tens in brightness. Another aspect is microlensing (ML), where we would like to continue our survey of a number of Galactic globular clusters over time-scales of weeks to years to search for ML events from planets to hypothesized central intermediate-mass black holes (IMBH).

Supriyo Ghosh (TIFR, Mumbai, India), J. P. Ninan (PSU, USA), D. K. Ojha (TIFR, Mumbai, India)

15 pages, 16 figures, 5 tables. Accepted for publication in MNRAS

We explored here the near-infrared $H$-band atmospheric window aiming to provide quantitative diagnostic tools for deriving stellar parameters, for instance, effective temperature ($T_{eff}$) and metallicity ([$Fe/H$]), of cool giants ($T_{eff}$ $<$ 5000 K) using low-resolution spectra. We obtained 177 cool giants from the X-shooter spectral library covering a wider metallicity range ($-$2.35 dex $<$ [$Fe/H$] $<$ 0.5 dex) than in earlier works. Degrading the spectral resolution to R$\sim$ 1200, we estimated equivalent widths of several important spectral features, and the behavior of spectral features with stellar parameters are studied. Also, the empirical relations for deriving $T_{eff}$ and [$Fe/H$] are established in the Bayesian framework. We found that $^{12}$CO at 1.56 $\mu$m and 1.62 $\mu$m, and $^{12}$CO+MgI at 1.71 $\mu$m are the best three $T_{eff}$ indicators with a typical accuracy of 153 K, 123 K and 107 K, respectively. The cubic Bayesian model provides the best metallicity estimator with a typical accuracy of 0.22 dex, 0.28 dex, and 0.24 dex for FeH at 1.62 $\mu$m, $^{12}$CO at 1.64 $\mu$m, and Fe I at 1.66 $\mu$m, respectively. We also showed a detailed quantitative metallicity dependence of $T_{eff}-$EWs correlations defining three metallicity groups, supersolar ([$Fe/H$] $>$ 0.0 dex), solar ($-$0.3 dex $<$ [$Fe/H$] $<$ 0.3 dex), and subsolar ([$Fe/H$] $<$ $-$0.3 dex), from Hierarchical Bayesian modelling. The difference between the solar and subsolar relationship is statistically significant, but such difference is not evident between the solar and supersolar groups.

H.E.S.S. Collaboration, F. Aharonian, F. Ait Benkhali, E.O. Anguner, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernloehr, M. Boettcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, T. Bulik, T. Bylund, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand, A. Chen, O. Chibueze, G. Cotter, P. Cristofari, J. Damascene Mbarubucyeye, J. Devin, A. Djannati-Atai, A. Dmytriiev, K. Egberts, S. Einecke, J.-P. Ernenwein, K. Feijen, A. Fiasson, G. Fichet de Clairfontaine, G. Fontaine, S. Funk, S. Gabici, Y.A. Gallant, S. Ghafourizadeh, G. Giavitto, L. Giunti, D. Glawion, J. F. Glicenstein, M.-H. Grondin, M. Hoerbe, W. Hofmann, T. L. Holch, M. Holler, D. Horns, Zhiqiu Huang, M. Jamrozy, I. Jung-Richardt, E. Kasai, K. Katarzynski, et al. (101 additional authors not shown)

submitted to Astronomy and Astrophysics

Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the gamma-ray luminosities of the remnants of Tycho's and Kepler's SNe are similar, then the latter is expected to be one of the faintest gamma-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6 sigma for a VHE signal from the remnant of Kepler's SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours. The measured integral flux above an energy of 226 GeV is ~0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a gamma-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho's SN, possibly indicating the same non-thermal emission processes acting in both these young remnants of thermonuclear SNe.

Sunspot groups observed in white-light appear as complex structures. Analysis of these structures is usually based on simple morphological descriptors which capture only generic properties and miss information about fine details. We present a machine learning approach to introduce a complete yet compact description of sunspot groups. The idea is to map sunspot group images into an appropriate lower-dimensional (latent) space. We apply a combination of Variational Autoencoder and Principal Component Analysis to obtain a set of 285 latent descriptors. We demonstrate that the standard descriptors are embedded into the latent ones. Thus, latent features can be considered as an extended description of sunspot groups and, in our opinion, can expand the possibilities for the research on sunspot groups. In particular, we demonstrate an application for estimation of the sunspot group complexity. The proposed parametrization model is generic and can be applied to investigation of other traces of solar activity observed in various spectrum lines. Key components of this work, which are the parametrization model, dataset of sunspot groups and latent vectors, are available in the public GitHub repository https://github.com/observethesun/sunspot groups and can be used to reproduce the results and for further research.

G. Alestas, L. Perivolaropoulos, K. Tanidis

9 pages, 8 figures. The numerical analysis files that lead to the production of the figures may be downloaded from this https URL

It has recently been pointed out that a gravitational transition taking place at a recent redshift $z_t$, reducing the effective gravitational constant $G_{\rm eff}$ by about $10\%$ for $z>z_t$, has the potential to lead to a resolution of the Hubble tension if $z_t\lesssim 0.01$. Since $H(z)^2\sim G_{\rm eff}$, such a transition would also lead to sharp change of the slope of the Hubble diagram at $z=z_t$ and a sharp decrease in the number of galaxies per redshift bin at $z_t$. Here we attempt to impose constraints on such a transition by using two robust low-z redshift survey datasets ($z<0.01$), taken from the Six-degree Field Galaxy Survey (6dFGS) as well as the 2MASS Redshift Survey (2MRS). In both surveys, we bin the data in redshift bins and focus on the number of galaxies in each bin ($\Delta N(z_i)$). We observe a peak in the distribution of galaxies near a distance of approximately 20 Mpc in both datasets. This feature could be attributed to galactic density fluctuations, to coherent peculiar velocities of galaxies or to an ultra late-time gravitational transition in the same era. In the context of the later scenario we show that this feature could have been induced by a sharp change of $G_{\rm eff}$ by $\Delta G_{\rm eff}/G_{\rm eff} \simeq 0.6$ at $z_t\simeq 0.005$. In a conservative approach we can then impose a bound on such gravitational transition as $\Delta G_{\rm eff}/G_{\rm eff} \lesssim 0.6$. We conclude that the required level of gravitational transition $\Delta G/G \simeq 0.1$ for the resolution of the Hubble tension can not be excluded by redshift survey data at $z<0.01$.

Vineet Ojha

9 pages, 2 figures, and a table. Accepted for publication in Journal of Astrophysics and Astronomy

SDSS J163401.94$+$480940.2 is a non-jetted radio-loud narrow-line Seyfert 1 (NLSy1) galaxy. Optical monitoring of this object was carried out in two intra-night sessions each $\geq$ 3 hrs with 3.6m DOT. Intra-night optical variability (INOV) characterization is presented for the first time for this source. We have detected an unexpected remarkable flare in one of two monitoring sessions of SDSS J163401.94$+$480940.2, whose rapid brightening phase implied a minute like doubling time of $\sim$ 22 minutes, thereby approaching to the extremely fast minute like variability, observed from FSRQ PKS 1222$+$21 at 400 GeV. The detection of a minute-like variability suggests the existence of relativistic jets with a small viewing angle. We briefly discuss the possible mechanisms for the non-detection of relativistic jets in its Very Long Baseline Array (VLBA) observations.

Stephan Rosswog

16 pages, 7 figures

Computational fluid dynamics is a crucial tool to theoretically explore the cosmos. In the last decade, we have seen a substantial methodological diversification with a number of cross-fertilizations between originally different methods. Here we focus on recent developments related to the Smoothed Particle Hydrodynamics (SPH) method. We briefly summarize recent technical improvements in the SPH-approach itself, including smoothing kernels, gradient calculations and dissipation steering. These elements have been implemented in the Newtonian high-accuracy SPH code MAGMA2 and we demonstrate its performance in a number of challenging benchmark tests. Taking it one step further, we have used these new ingredients also in the first particle-based, general-relativistic fluid dynamics code that solves the full set of Einstein equations, SPHINCS_BSSN. We present the basic ideas and equations and demonstrate the code performance at examples of relativistic neutron stars that are evolved self-consistently together with the spacetime.

Eckhart Spalding, Katie M. Morzinski, Phil Hinz, Jared Males, Michael Meyer, Sascha P. Quanz, Jarron Leisenring, Jennifer Power

Published in AJ

The Large Binocular Telescope (LBT) has two 8.4-m primary mirrors that produce beams that can be combined coherently in a "Fizeau" interferometric mode. In principle, the Fizeau PSF enables the probing of structure at a resolution up to three times better than that of the adaptive-optics-corrected PSF of a single 8.4-m telescope. In this work, we examined the nearby star Altair (5.13 pc, type A7V, $\sim$100s Myr to $\approx$1.4 Gyr) in the Fizeau mode with the LBT at Br-$\alpha$ (4.05 $\mu$m) and carried out angular differential imaging to search for companions. This work presents the first filled-aperture LBT Fizeau science dataset to benefit from a correcting mirror which provides active phase control. In the analysis of the $\lambda/D$ angular regime, the sensitivity of the dataset is down to $\approx$0.5 $M_{\odot}$ at 1" for a 1.0 Gyr system. This sensitivity remains limited by the small amount of integration time, which is in turn limited by the instability of the Fizeau PSF. However, in the Fizeau fringe regime we attain sensitivities of $\Delta m \approx 5$ at 0.2" and put constraints to companions of 1.3 $M_{\odot}$ down to an inner angle of $\approx$0.15", closer than any previously published direct imaging of Altair. This analysis is a pathfinder for future datasets of this type, and represents some of the first steps to unlocking the potential of the first ELT. Fizeau observations will be able to reach dimmer targets with upgrades to the instrument, in particular the phase detector.

The logotropic model [P.H. Chavanis, Eur. Phys. J. Plus {\bf 130}, 130 (2015)] may be an interesting alternative to the $\Lambda$CDM model. It is able to account for the present accelerating expansion of the universe while solving at the same time the core-cusp problem of the CDM model. In the logotropic model, there is a single dark fluid. Its rest-mass plays the role of dark matter and its internal energy plays the role of dark energy. We highlight two remarkable predictions of the logotropic model. It yields cored dark matter halos with a universal surface density equal to $\Sigma_0^{\rm th}=0.01955 c\sqrt{\Lambda}/G=133\, M_{\odot}/{\rm pc}^2$ without free parameter in very good agreement with the observational value $\Sigma_0^{\rm obs}=141_{-52}^{+83}\, M_{\odot}/{\rm pc}^2$. It also predicts the present ratio of dark energy and dark matter to be the pure number $\Omega_{\rm de,0}^{\rm th}/\Omega_{\rm dm,0}^{\rm th}=e=2.71828...$ in very good agreement with the observations giving $\Omega_{\rm de,0}^{\rm obs}/\Omega_{\rm dm,0}^{\rm obs}=2.669\pm 0.08$. Using the measured present proportion of baryonic matter $\Omega_{\rm b,0}^{\rm obs}=0.0486\pm 0.0010$, we find that the values of the present proportion of dark matter and dark energy are $\Omega_{\rm dm,0}^{\rm th}=\frac{1}{1+e}(1-\Omega_{\rm b,0})=0.2559$ and $\Omega_{\rm de,0}^{\rm th}=\frac{e}{1+e}(1-\Omega_{\rm b,0})=0.6955$ in very good agreement with the observational values $\Omega_{\rm dm,0}^{\rm obs}=0.2589\pm 0.0057$ and $\Omega_{\rm de,0}^{\rm obs}=0.6911\pm 0.0062$ within the error bars. These theoretical predictions are obtained by advocating a mysterious strong cosmic coincidence (dubbed "dark magic") implying that our epoch plays a particular role in the history of the universe.

Marcos A. van Dam, Yolanda Martín Hernando, Miguel Núñez Cagigal, Luzma M. Montoya

Multiconjugate adaptive optics (MCAO) systems have the potential to deliver diffraction-limited images over much larger fields of view than traditional single conjugate adaptive optics systems. In MCAO, the high altitude deformable mirrors (DMs) cause a distortion of the pupil plane and lead to a dynamic misregistration between the DM actuators and the wavefront sensors (WFSs). The problem is much more acute for solar astronomy than for night-time observations due to the higher spatial sampling of the WFSs and DMs, and the fact that the science observations are often made through stronger turbulence and at lower elevations. The dynamic misregistration limits the quality of the correction provided by solar MCAO systems. In this paper, we present PropAO, the first AO simulation tool (to our knowledge) to model the effect of pupil distortion. It takes advantage of the Python implementation of the optical propagation library PROPER. PropAO uses Fresnel propagation to propagate the amplitude and phase of an incoming wave through the atmosphere and the MCAO system. The resulting wavefront is analyzed by the WFSs and also used to evaluate the corrected image quality. We are able to reproduce the problem of pupil distortion and test novel non-linear reconstruction strategies that take the distortion into account. PropAO is shown to be an essential tool to study the behavior of the wavefront reconstruction and control for the European Solar Telescope.

H. A. Weaver, S. B. Porter, J. R. Spencer, The New Horizons Science Team

30 pages, 15 figures

Using the New Horizons LORRI camera, we searched for satellites near five Kuiper belt objects (KBOs): four cold classicals (CCs: 2011 JY31, 2014 OS393, 2014 PN70, 2011 HZ102) and one scattered disk object (SD: 2011 HK103). These objects were observed at distances of 0.092-0.290 au from the New Horizons spacecraft, achieving spatial resolutions of 136-430 km (resolution is ~2 camera pixels), much higher than possible from any other facilities. Here we report that CC 2011 JY31 is a binary system with roughly equal brightness components, CC 2014 OS393 is likely an equal brightness binary system, while the three other KBOs did not show any evidence of binarity. The 2011 JY31 binary has a semi-major axis of 198.6 +/- 2.9 km, an orbital inclination of 61.34 +/- 1.34 deg, and an orbital period of 1.940 +/- 0.002 d. The 2014 OS393 binary objects have an apparent separation of ~150 km, making 2011 JY31 and 2014 OS393 the tightest KBO binary systems ever resolved. Both 2011 HK103 and 2011 HZ102 were detected with SNR~10, and our observations rule out equal brightness binaries with separations larger than ~430 km and ~260 km, respectively. The spatial resolution for 2014 PN70 was ~200 km, but this object had SNR~2.5-3, which limited our ability to probe its binarity. The binary frequency for the CC binaries probed in our small survey (67%, not including 2014 PN70) is consistent with the high binary frequency suggested by larger surveys of CCs (Fraser et al. 2017, Noll et al. 2020) and recent planetesimal formation models (Nesvorny et al. 2021), but we extend the results to smaller orbit semi-major axes and smaller objects than previously possible.

We have collected from different surveys some significant spectroscopic data observed from star-forming galaxies in the local Universe. The objects showing a relatively rich spectrum in number of lines from different elements were selected in order to constrain the models. In particular, we looked at the relatively weak lines such as [OIII]4363, HeII4686, HeI4471 and HeI5876. We have modelled in detail the spectra by the coupled effect of photoionization from the stars and shocks. We have found that the abundances relative to H of most of the elements are lower than solar but not as low as those evaluated by the direct strong line and the Te methods. Sulfur which appears through the [SII]6717, 6731 and [SIII]6312 lines is not depleted, revealing a strong contribution from the ISM. We have added to the sample the optical-UV spectra of local low-metallicity dwarf galaxies containing the CIV/Hb and CIII]/Hb line ratios in order to determine with relative precision the C/H relative abundance. The results show He/H lower than solar in some objects and suggest that the geometrical thickness of the clouds constrains the HeII/Hb line ratios. We explain the low He/H by mixing of the wind from the star-forming region with ISM clouds.

De-Chang Dai, Glenn Starkman, Dejan Stojkovic

4 figures, 7 pages

A systematic downward trend in the weak gravity (large distances) segment of the radial acceleration in galactic rotation curves was found in the EAGLE simulation \cite{Dai:2017unr}, which was based on the standard $\Lambda$CDM model. A similar feature was discovered in \cite{Chae:2020omu} and interpreted as a MOND prediction. We analyze the existing rotation-curve data that extends to very large distances and low accelerations for the Milky Way and M31, and find that they are consistent with the $\Lambda$CDM EAGLE simulation, while their deviation from the generic MOND prediction would need to be ascribed to an external field effect, or possibly a {\it post facto} selected acceleration function $\mu(x)$.

Natalí S. M. de Santi, Natália V. N. Rodrigues, Antonio D. Montero-Dorta, L. Raul Abramo, Beatriz Tucci, M. Celeste Artale

Elucidating the connection between the properties of galaxies and the properties of their hosting haloes is a key element in theories of galaxy formation. When the spatial distribution of objects is also taken under consideration, investigating the halo-galaxy connection becomes very relevant for cosmological measurements. In this paper, we use machine learning (ML) techniques to analyze these intricate relations in the IllustrisTNG300 magnetohydrodynamical simulation. We employ four different algorithms: extremely randomized trees (ERT), K-nearest neighbors (kNN), light gradient boosting machine (LGBM), and neural networks (NN), along with a stacked model where we combine results from all four approaches. Overall, the different ML algorithms produce consistent results in terms of predicting galaxy properties from a set of input halo properties that include halo mass, concentration, spin, and halo overdensity. For stellar mass, the (predicted v. true) Pearson correlation coefficient is 0.98, dropping down to 0.7-0.8 for specific star formation rate (sSFR), colour, and size. In addition, we test an existing data augmentation technique, designed to alleviate the problem of unbalanced datasets, and show that it improves slightly the shape of the predicted distributions. We also demonstrate that our predictions are good enough to reproduce the power spectra of multiple galaxy populations, defined in terms of stellar mass, sSFR, colour, and size with high accuracy. Our results align with previous reports suggesting that certain galaxy properties cannot be reproduced using halo features alone.

Zhitong Li, Tianmeng Zhang, Xiaofeng Wang, Hanna Sai, Jujia Zhang, Juncheng Chen, Xulin Zhao, Shengyu Yan, Bo Wang, Mark M. Phillips, Eric Y. Hsiao, Nidia Morrell, Carlos Contreras, Christopher R. Burns, Christopher Ashall, Maximilian Stritzinger, Kevin Krisciunas, Jose Prieto, Hu Zou, Jiali Wang, Jun Ma, Jundan Nie, Suijian Xue, Xu Zhou, Zhimin Zhou, Danfeng Xiang, Gaobo Xi

In this paper, we present photometric and spectroscopic observations of a subluminous type Ia supernova (SN Ia) 2012ij, which has an absolute $B$-band peak magnitude $M_{B,\rm{max}}$ = $-$17.95 $\pm$ 0.15 mag. The $B$-band light curve exhibits a fast post-peak decline with $\Delta m_{15}(B)$ = 1.86 $\pm$ 0.05 mag. All the $R$ and $I$/$i$-band light curves show a weak secondary peak/shoulder feature at about 3 weeks after the peak, like some transitional subclass of SNe Ia, which could result from an incomplete merger of near-infrared (NIR) double peaks. The spectra are characterized by Ti~{\sc ii} and strong Si~{\sc ii} $\lambda$5972 absorption features that are usually seen in low-luminosity objects like SN 1999by. The NIR spectrum before maximum light reveals weak carbon absorption features, implying the existence of unburned materials. We compare the observed properties of SN 2012ij with those predicted by the sub-Chandrasekhar-mass and the Chandrasekhar-mass delayed-detonation models, and find that both optical and NIR spectral properties can be explained to some extent by these two models. By comparing the secondary maximum features in $I$ and $i$ bands, we suggest that SN 2012ij is a transitional object linking normal SNe Ia to typical 91bg-like ones. From the published sample of SNe Ia from the $Carnegie~Supernova~Project~II$ (CSP-II), we estimate that the fraction of SN 2012ij-like SNe Ia is not lower than $\sim$ 2%.

Kevin D. Crowley, Peter Dow, Jordan E. Shroyer, Bradley Dober, Jacob Spisak, Nicholas Galitzki, Tanay Bhandarkar, Mark J. Devlin, Simon Dicker, Patricio A. Gallardo, Kathleen Harrington, Bradley R. Johnson, Delwin Johnson, Anna M. Kofman, Akito Kusaka, Adrian Lee, Michele Limon, Jeffrey Iuliano, Federico Nati, John Orlowski-Scherer, Lyman Page, Michael Randall, Grant Teply, Tran Tsan, Edward J. Wollack, Zhilei Xu, Ningfeng Zhu

We present the design and measured performance of a new carbon fiber strut design that is used in a cryogenically cooled truss for the Simons Observatory Small Aperture Telescope (SAT). The truss consists of two aluminum 6061 rings separated by 24 struts. Each strut consists of a central carbon fiber tube fitted with two aluminum end caps. We tested the performance of the strut and truss by (i) cryogenically cycling and destructively pull-testing strut samples, (ii) non-destructively pull-testing the final truss, and (iii) measuring the thermal conductivity of the carbon fiber tubes. We found that the strut strength is limited by the mounting fasteners and the strut end caps, not the epoxy adhesive or the carbon fiber tube. This result is consistent with our numerical predictions. Our thermal measurements suggest that the conductive heat load through the struts (from 4 K to 1 K) will be less than 1 mW. This strut design may be a promising candidate for use in other cryogenic support structures.

Although first considered as too diluted for the formation of molecules in-situ and too harsh an environment for their survival, the interstellar medium has turned out to host a rich palette of molecular species: to date, 256 species have been identified. The last decade has seen an explosion of new detections, including those of a number of complex organic species, which may be dubbed as prebiotic. Organic molecules have been discovered not just in interstellar clouds from the Solar neighbourhood, but also throughout the Milky-Way, as well as in nearby galaxies, or some of the most distant quasars. These discoveries were made possible by the completion of large sub-millimetre and radio facilities. Equipped with new generation receivers, those instruments have provided the orders of magnitude leap in sensitivity required to detect the weak rotational lines that allowed the molecule identifications. Last two years, 30 prebiotic molecules have been detected in TMC-1, a dust-enshrouded gaseous cloud located at 400 light-years from the Sun in the Taurus constellation. Ten new molecular species, have been identified in the arm of a spiral galaxy 6 billion light-yr distant, and 12 molecular species observed in a quasar at 11 billion light-yr. We present the latest spectral observations of this outlying quasar and discuss the implications of those detections in these 3 archetypal sources. The basic ingredients involved in the Miller-Urey experiment and related experiments appeared early after the formation of the first galaxies and are widespread throughout the Universe. The chemical composition of the gas in distant galaxies seems not much different from that in the nearby interstellar clouds. It presumably comprises, like for TMC-1, aromatic rings and complex organic molecules putative precursors of the RNA nucleobases, except the lines of such species are too weak to be detected that far.

Liyi Gu, Chintan Shah, Ruitian Zhang

accepted for publication in Sensors, plasma diagnostics special issue

Relevant uncertainties on theoretical atomic data are vital to determine the accuracy of plasma diagnostics in a number of areas including in particular the astrophysical study. We present a new calculation of the uncertainties on the present theoretical ion-impact charge exchange atomic data and X-ray spectra based on a set of comparisons with the existing laboratory data obtained in historical merged-beam, cold-target recoil-ion momentum spectroscopy, and electron beam ion traps experiments. The average systematic uncertainties are found to be 35-88% on the total cross sections, and 57-75% on the characteristic line ratios. The model deviation increases as the collision energy decreases. The errors on total cross sections further induce a significant uncertainty to the calculation of ionization balance for low temperature collisional plasmas. Substantial improvements of the atomic database and dedicated laboratory measurements are needed to get the current models ready for the X-ray spectra from the next X-ray spectroscopic mission.

Shiying Lu, Qiusheng Gu, Xue Ge, Luis C. Ho, Yulong Gao, Zhengyi Chen, Ke Xu, Zhi-Yu Zhang, Yong Shi, Qirong Yuan, Min Bao

18pages, 13figures, 1table, accepted by ApJ

Comparing to the inactive and gas-poor normal lenticular galaxies (S0s) in the local universe, we study a barred star-forming S0 galaxy, PGC 34107, which has been observed by the Centro Astron\'{o}mico Hispano Alem\'{a}n (CAHA) 3.5-m telescope and the Northern Extended Millimeter Array (NOEMA). The spatially resolved ionized gas and molecular gas traced by $^{12}$CO(1-0), hereafter CO(1-0), show the similar distribution and kinematics to the stellar component with an off-center star-forming region, $\sim$380 pc away from the center. The resolved kinematics of molecular CO(1-0) emission reveals that there is a blueshifted (redshifted) velocity component on the receding (approaching) side of the galaxy along the stellar bar. This might provide a plausible evidence of non-circular motion, such as the bar-induced molecular gas inflow. The velocity of molecular gas inflow decreases with approaching towards the peak of the off-center star formation in the north, which might be associated with the inner Lindblad resonance (ILR). In addition to CO(1-0), we also detect the isotopic line of $^{13}$CO(1-0). Most $\rm H\alpha$, CO(1-0) and $^{13}$CO(1-0) emissions are concentrated on this northern star-forming region. We find that PGC 34107 follows the local stellar mass-metallicity relation, star-forming main sequence, and the Kennicutt-Schmidt law. The resolved and integrated molecular gas main sequence suggest that there is a higher gas fraction in the galaxy central region, which supports a scenario that the bar-induced gas reservoir provides the raw material, and subsequently triggers the central star formation.

Shu-Xu Yi, Fiorenzo Stoppa, Gijs Nelemans, Eric Cator

8 pages, 10 figures, submitted to A&A

We employ the Gravitational Waves-Universe toolbox to generate synthetic catalogues of detections of stellar mass binary black hole (BBH) mergers. Using the catalogues, we study how GW observations of BBHs can be used to constrain the merger rate as function of redshift and masses. We study advanced LIGO (aLIGO) and Einstein Telescope (ET) as two representatives for the 2nd and 3rd generation GW observatories. We also simulate the observations from a detector that is half as sensitive as the ET in design which represents an early phase of ET. Two methods are used to obtain the constraints on the source population properties from the catalogues: the first assumes a parametric differential merger rate model and applies a Bayesian inference; The other is non-parametric and uses weighted Kernel density estimators. The results show the overwhelming advantages of the 3rd generation detector over the 2nd generation for the study of BBH population properties, especially at a redshifts higher than $\sim2$, where the merger rate is believed to peak. With the LIGO catalogue, the parameteric Bayesian method can still give some constraints on the merger rate density and mass function beyond its detecting horizon, while the non-parametric method lose the constraining ability completely there. These two methods represent the two extreme situations of general population reconstruction. We also find that, despite the numbers of detection of the half-ET can be easily compatible with full ET after a longer observation duration, the catalogue from the full ET can still give much better constraints on the population properties, due to its smaller uncertainties on the physical parameters of the GW events.

Observations of black hole X-ray binaries and active galactic nuclei indicate that the accretion flows around black holes are composed of hot and cold gas, which have been theoretically described in terms of either a hot geometrically thick corona lying above and below a cold geometrically thin disk or an inner advection dominated accretion flow connected to an outer thin disk. This article reviews the accretion flows around black holes, with an emphasis on the physics that determines the configuration of hot and cold accreting gas, and how the configuration varies with the accretion rate and thereby produces various luminosity and spectra.

Bing-Qiang Qiao, Qing Luo, Qiang Yuan, Yi-Qing Guo

8 pages, 2 figures

The energy spectra and anisotropies are very important probes of the origin of cosmic rays. Recent measurements show that complicated but very interesting structures exist, at similar energies, in both the spectra and energy-dependent anisotropies, indicating a common origin of these structures. Particularly interesting phenomenon is that there is a reversal of the phase of the dipole anisotropies, which challenges a theoretical modeling. In this work, for the first time, we identify that there might be an additional phase reversal at $\sim 100$ GeV energies of the dipole anisotropies as indicated by a few underground muon detectors and the first direct measurement by the Fermi satellite, coincident with the hundreds of GV hardenings of the spectra. We propose that these two phase reversals, together with the energy-evolution of the amplitudes and spectra, can be naturally explained with a nearby source overlapping onto the diffuse background. As a consequence, the spectra and anisotropies can be understood as the scalar and vector components of this model, and the two reversals of the phases characterize just the competition of the cosmic ray streamings between the nearby source and the background. The alignment of the cosmic ray streamings along the local large-scale magnetic field may play an important but sub-dominant role in regulating the cosmic ray propagation. More precise measurements of the anisotropy evolution at both low energies by space detectors and high energies by air shower experiments for individual species will be essential to further test this scenario.

Julia C. Santos, Leonardo Bronfman, Edgar Mendoza, Jacques R. D. Lépine, Nicolas U. Duronea, Manuel Merello, Ricardo Finger

A spectral survey of methyl acetylene (CH3CCH) was conducted toward the hot molecular core/outflow G331.512-0.103. Our APEX observations allowed the detection of 41 uncontaminated rotational lines of CH3CCH in the frequency range between 172-356 GHz. Through an analysis under the local thermodynamic equilibrium assumption, by means of rotational diagrams, we determined Texc = 50 \pm 1 K, N(CH3CCH) = (7.5 \pm 0.4) x 10^{15} cm^{-2}, X[CH3CCH/H2] ~ (0.8-2.8) x 10^{-8} and X[CH3CCH/CH3OH] ~ 0.42 \pm 0.05 for an extended emitting region (~10 arcsec). The relative intensities of the K=2 and K=3 lines within a given K-ladder are strongly negatively correlated to the transitions' upper J quantum-number (r=-0.84). Pure rotational spectra of CH3CCH were simulated at different temperatures, in order to interpret this observation. The results indicate that the emission is characterized by a non-negligible temperature gradient with upper and lower limits of ~45 and ~60 K, respectively. Moreover, the line widths and peak velocities show an overall strong correlation with their rest frequencies, suggesting that the warmer gas is also associated with stronger turbulence effects. The K=0 transitions present a slightly different kinematic signature than the remaining lines, indicating that they might be tracing a different gas component. We speculate that this component is characterized by lower temperatures, and therefore larger sizes. Moreover, we predict and discuss the temporal evolution of the CH3CCH abundance using a two-stage zero-dimensional model of the source constructed with the three-phase Nautilus gas-grain code.

Berta Margalef-Bentabol, Christopher J. Conselice, Boris Haeussler, Kevin Casteel, Chris Lintott, Karen Masters, Brooke Simmons

Many aspects concerning the formation of spiral and disc galaxies remain unresolved, despite their discovery and detailed study over the past $150$ years. As such, we present the results of an observational search for proto-spiral galaxies and their earliest formation, including the discovery of a significant population of spiral-like and clumpy galaxies at $z>1$ in deep \textit{Hubble Space Telescope} CANDELS imaging. We carry out a detailed analysis of this population, characterising their number density evolution, masses, star formation rates and sizes. Overall, we find a surprisingly high overall number density of massive $M_{*} >10^{10}\mathrm{M}_{\odot}$ spiral-like galaxies (including clumpy spirals) at $z > 1$ of $0.18\,{\rm per}\, \mathrm{arcmin}^{-2}$. We measure and characterise the decline in the number of these systems at higher redshift using simulations to correct for redshift effects in identifications, finding that the true fraction of spiral-like galaxies grows at lower redshifts as $\sim$ $(1+z)^{-1.1}$. This is such that the absolute numbers of spirals increases by a factor of $\sim 10$ between $z = 2.5$ and $z = 0.5$. We also demonstrate that these spiral-like systems have large sizes at $z>2$, and high star formation rates, above the main-sequence, These galaxies represent a major mode of galaxy formation in the early universe, perhaps driven by the spiral structure itself. We finally discuss the origin of these systems, including their likely formation through gas accretion and minor mergers, but conclude that major mergers are an unlikely cause.

Hou-Jun Lü, Hao-Yu Yuan, Ting-Feng Yi, Xiang-Gao Wang, You-Dong Hu, Yong Yuan, Jian-Guo Wang, Jia-Xin Cao, De-Feng Kong, Emilio Fernandez-García, Alberto J.Castro-Tirado, Ji-Shun Lian, Wen-Pei Gan, Shan-Qin Wang, Li-Ping Xin, M.D. Caballero-García, Yu-Feng Fan, En-Wei Liang

23 pages, 7 figures, and 1 Table. Comments are welcome

Long-duration gamma-ray bursts (GRBs) that associated with supernova (SN) are believed to originate from massive star core-collapse events, whereas short-duration GRBs that are related to compact star mergers are expected to be accompanied by kilonova. GRB 211227A, lasting about 84 s, has a initial short/hard spike followed by a series of soft gamma-ray extended emission at redshift $z=$0.228. Moreover, lack of supernova associated, is down to very stringent limits at such low redshift, and is with a large physical offset ($20.47\pm14.47$ kpc) from the galaxy center. Those behaviors are similar to that of GRB 060614, and suggest that the progenitor of GRB 211227A is not favored with the death of massive stars. Hence, we propose that GRB 211227A has the same physical origin with that of GRB 060614 which is from compact star mergers. One calculate pseudo-kilonova emission of this case by adopting the typical parameters with that of GRB 060614 and GRB 170817A, and we find that this pseudo-kilonova is too faint to be detected. If this is the case, it can naturally interpret the characteristic of prompt emission, lack of SN and kilonova emission, and the large physical offset from the galaxy center.

Fabian Walter, Marcel Neeleman, Roberto Decarli, Bram Venemans, Romain Meyer, Axel Weiss, Eduardo Banados, Sarah E. I. Bosman, Chris Carilli, Xiaohui Fan, Dominik Riechers, Hans-Walter Rix, Todd A. Thompson

Accepted for publication in The Astrophysical Journal

We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc.

Min Li, Shichun Huang, Zhaohuan Zhu, Michail I. Petaev, Jason H. Steffen

9 pages, 11 figures

The maximum temperature and radial temperature profile in a protoplanetary disc are important for the condensation of different elements in the disc. We simulate the evolution of a set of protoplanetary discs from the collapse of their progenitor molecular cloud cores as well as the dust decoupling within the discs as they evolve. We show how the initial properties of the cloud cores affect the thermal history of the protoplanetary discs using a simple viscous disc model. Our results show that the maximum midplane temperature in the disc occurs within 0.5 AU. It increases with the initial cloud temperature and decreases with its angular velocity and the viscosity of the disc. From the observed properties of the molecular cloud cores we find the median value of the maximum temperature is around 1250 K, with roughly 90% of them being less than 1500 K - a value that is lower than the 50% condensation temperatures of most refractory elements. Therefore, only cloud cores with high initial temperatures or low angular velocities and/or low viscosities within the planet-forming discs will result in refractory-rich planetesimals. To reproduce the volatile depletion pattern of CM, CO, and CV chondrites and the terrestrial planets in Solar system, one must either have rare properties of the initial molecular cloud cores like high core temperature, or other sources of energy to heat the disc to sufficiently high temperatures. Alternatively, the volatile depletion observed in these chondrites may be inherited from the progenitor molecular cloud.

Fulvio Melia

8 pages, 1 figure. Accepted for publication in Astronomische Nachrichten / Astronomical Notes

The evidence for a Big Bang origin of the Universe is truly compelling, though its cause remains a complete mystery. As the cosmic spacetime is revealed to us with ever improving detail, however, we are beginning to refine the range of its possible initial conditions -- at least within the framework of current physical theories. The Universe, it seems, is spatially flat, and here we discuss in clear, straightforward terms why this trait implies a cosmos with zero kinetic plus gravitational energy, though apparently not zero total energy. Such an outcome has far reaching consequences because of the possibility that the Universe may have begun its existence as a quantum fluctuation. Was this from `nothing,' or perhaps a pre-existing vacuum? A non-zero total energy would seemingly preclude the former scenario, but not necessarily the latter, though this would then raise the question of how a fluctuation with non-zero energy could have lived long enough, or classicalized, for us to see it 13.5 billion years later. The high-precision measurement of the Universe's spatial curvature may thus constitute the first tangible piece of evidence impacting a possible quantum beginning.

A. Ranjan, R. Srianand, P. Petitjean, G. Shaw, Y.-K. Sheen, S. A. Balashev, N. Gupta, C. Ledoux, K. N. Telikova

16 pages, 12 figures, 16 page appendix, 14 figures in appendix

We present the results of a spectroscopic analysis of extremely strong damped Lyman-{\alpha} absorbers (ESDLAs, log N(Hi)>=21.7) observed with VLT-XShooter. ESDLAs probe gas from within the star-forming disk of the associated galaxies and thus ESDLAs provide a unique opportunity to study the interstellar medium of galaxies at high-redshift. We report column densities (N), equivalent widths (w), and the kinematic spread ({\Delta} v90) of species from neutral, singly ionised, and higher ionisation species. We find that, using the dust correction prescription, the measured metallicities are consistent for singly ionised gas species such as Pii, S ii, Si ii, Mnii and Crii, and Znii in all ESDLAs within 3-sigma uncertainty. We find that the distributions of N(Ari)/N(Hi) ratio in DLAs and ESDLAs are similar. We further report that ESDLAs do not show a strong deficiency of Ari relative to other {\alpha}-capture elements as is seen in DLAs. This supports the idea that the mentioned under-abundance of Ari in DLAs is possibly caused by the presence of background UV photons that penetrate the low N(Hi) clouds to ionise Ari, but they cannot penetrate deep enough in the high N(Hi) ESDLA environment. The w(Mgii lambda2796) distribution in ESDLAs is found to be similar to that of metal-rich Ci-selected absorbers, but the velocity spread of their Mgii profile is different. For higher ionisation species (such as C iv and Si iv), {\Delta} v90 is similar in the two populations, while the {\Delta} v90 of singly ionised species is smaller for ESDLAs. This suggests that the ESDLAs sample a different Hi region of their associated galaxy compared to the general DLA population. We further study the N(Cl i) distribution in high-redshift DLA and ESDLA sightlines, as Cl i is a good tracer of H2 gas. The N(Cl i)-N(H2) correlation is followed by all the clouds (ESDLAs and otherwise) having log N(H2)<22.

F. Marin, T. Barnouin, E. Lopez-Rodriguez

4 pages, 2 figures, to be published in the proceedings of the French SF2A

The Faint Object Camera (FOC) aboard the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet imaging polarimetry between 1990 and 2002. Tremendous progresses have been made thanks to those high spatial resolution, high signal-to-noise ratio observations, such as the identification of the location of hidden active nuclei and the three dimensional arrangement of polar material within the first hundred of parsecs around the central core. However, not all AGN observations have been reduced and analyzed, and none in a standardized framework. In this lecture note, we present our project of downloading, reducing and analyzing all the AGN HST/FOC observations that were achieved using a consistent, novel and open-access reduction pipeline. We briefly present our methodology and show the first, preliminary result from our reduction pipeline: NGC 1068.

In this comment of the article [arXiv:2002.01150] "Locating the source field lines of Jovian decametric radio emissions" by YuMing Wang et al., 2020, we discuss the assumptions used by the authors to compute the beaming angle of Jupiter s decametric emissions induced by the moon Io. Their method, relying on multi-point radio observations, was applied to a single event observed on 14th March 2014 by Wind and both STEREO A/B spacecraft from 5 to 16 MHz, and erroneously identified as a northern emission (Io-B type) instead of a southern one (Io-D type). We encourage the authors to update their results with the right hemisphere of origin and to test their method on a larger sample of Jupiter-Io emissions.

The primordial $^4$He abundance (Y$_p$) is one of the key characteristics of Primordial Nucleosynthesis processes that occurred in the first minutes after the Big Bang. Its value depends on the baryon/photon ratio $\eta\equiv n_b/n_{\gamma}$, and is also sensitive to the relativistic degrees of freedom which affect the expansion rate of the Universe at the radiation-dominated era. The most used method of the determination of Y$_p$ is the study of the metal deficient HII regions located in blue compact dwarf galaxies (BCDs). In this paper, we discuss in detail various methods of the determination of HII region metallicity in the context of Y$_p$ analyses. We show that some procedures used in the methods lead to biases in the metallicity estimates and underestimation of their uncertainties. We propose a modified method for the metallicity determination, as well as an additional criterion for selecting objects. We have selected 69 objects (26 objects with high quality spectra from the HeBCD+NIR database and 43 objects from the SDSS catalog), for which we estimate Y and O/H using the proposed method. We have estimated Y$_p=0.2470\pm0.0020$ which is one of the most accurate estimates obtained up to date. Its comparison with the value Y$_p=0.2470\pm0.0002$ obtained as a result of numerical modelling of Primordial Nucleosynthesis with the value of $\Omega_b$ taken from the analysis of the CMB anisotropy (Planck mission), is an important tool for studying the self-consistency of the Standard cosmological model (a possible discrepancy between these estimates could be an indicator of a new physics). The application of the proposed method allows one to more correctly estimate Y$_p$ and the slope $d$Y/$d$(O/H). Further analysis of the data from the SDSS catalog can significantly increase the statistics of objects for the regression analysis, which in turn can refine the Y$_p$ estimate.

Realizing the inflationary epoch driven by a pseudo-Nambu Goldstone boson (pNGB) could ensure the coveted flatness, and the sub-Planckian scales related to the dynamics of the paradigm. In this work, we have taken the most general form of such a scenario: Goldstone inflation proposed in \cite{croon}, and studied the model in Einstein-Gauss-Bonnet gravity. Natural inflation is a limiting case of this model, is also studied here. The specific form of the EGB coupling gives ample opportunity to study the rich phenomenology associated with inflation as well as the reheating epoch. Predicted values of the inflationary observables, tensor to scalar ratio ($r$), and spectral index ($n_s$) are in good agreement with the recent observations from $Planck'18$ \cite{Planck2018}. Thus, in the framework of EGB one can resurrect the model, which otherwise needs quite a bit of fine tuning or diversion from the canonical domain as studied in \cite{Bhattacharya:2018xlw}, to survive in the standard cold inflationary scenario. Finally, the era of reheating is studied for different choice of model parameters.

M. Giulietti, M. Massardi, A. Lapi, M. Bonato, A. F. M. Enia, M. Negrello, Q. D'Amato, M. Behiri, G. De Zotti

Accepted for publication in MNRAS. 12 pages, 6 figures

We investigate the radio-far infrared (FIR) correlation for a sample of $28$ bright high-redshift ($1 \lesssim z \lesssim 4$) star-forming galaxies selected in the FIR from the Herschel-ATLAS fields as candidates to be strongly gravitationally lensed. The radio information comes either from high sensitivity dedicated ATCA observations at $2.1$ GHz or from cross-matches with the FIRST survey at $1.4$ GHz. By taking advantage of source brightness possibly enhanced by lensing magnification, we identify a weak evolution with redshift out to $z\lesssim 4$ of the FIR-to-radio luminosity ratio $q_{\rm FIR}$. We also find that the $q_{\rm FIR}$ parameter as a function of the radio power $L_{1.4\,\rm GHz}$ displays a clear decreasing trend, similarly to what is observed for optically/radio selected lensed quasars found in literature, yet covering a complementary region in the $q_{\rm FIR}-L_{1.4\,\rm GHz}$ diagram. We interpret such a behavior in the framework of an in-situ galaxy formation scenario, as a result of the transition from an early dust-obscured star-forming phase (mainly pinpointed by our FIR selection) to a late radio-loud quasar phase (preferentially sampled by the optical/radio selection).

Claire L. Davies, Evan A. Rich, Tim J. Harries, John D. Monnier, Anna S. E. Laws, Sean M. Andrews, Jaehan Bae, David J. Wilner, Narsireddy Anugu, Jacob Ennis, Tyler Gardner, Stefan Kraus, Aaron Labdon, Jean-Baptiste le Bouquin, Cyprien Lanthermann, Gail H. Schaefer, Benjamin R. Setterholm, Theo ten Brummelaar

18 pages, 16 figures; Accepted for publication in MNRAS

We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim \lambda/2\pi$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au.

Irina A. Acharova, Margarita E. Sharina, Egor A. Kazakov

16 pages, 6 figures accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

A method is proposed for determining the properties of type Ia supernovae from short-lived precursors -- Prompt SNIa. This method is based on the assumption that this subtype of type Ia supernovae exploded into low-metallicity globular clusters (GCs), and is responsible for the enrichment of the high-metallicity subgroup of GCs and circumgalactic clouds (CGCs) with the iron peak elements. We justify that CGCs are the formation places of GCs of both subgroups. The accuracy of the method depends, first, on the number of GCs, the spectra of which have been studied in detail; second, on the number of chemical elements, the abundances of which have been worked out. Only those elements are of interest for this method that are produced in supernova explosions and are not produced at the previous stage of the stellar evolution. Our estimates of nucleosynthesis in low-metallicity supernova GCs are in the best agreement with the following Prompt SNIa model: Single Degenerate Pure Deflagration Models of white dwarfs (WDs) burning with masses in the range from 1.30 Msun to 1.31 Msun if carbon explodes in the centre of a WD with a low central density from 0.5*10^9 g/cm^3 to 10^9 g/cm^3.

C. C. Espaillat, G. J. Herczeg, T. Thanathibodee, C. Pittman, N. Calvet, N. Arulanantham, K. France, Javier Serna, J. Hernandez, A. Kospal, F.M. Walter, A. Frasca, W.J. Fischer, C.M. Johns-Krull, P.C. Schneider, C. Robinson, Suzan Edwards, P. Abraham, Min Fang, J. Erkal, C.F. Manara, J.M. Alcala, E. Alecian, R.D. Alexander, J. Alonso-Santiago, Simone Antoniucci, David R. Ardila, Andrea Banzatti, M. Benisty, Edwin A. Bergin, Katia Biazzo, Cesar Briceno, Justyn Campbell-White, L. Ilsedore Cleeves, Deirdre Coffey, Jochen Eisloffel, Stefano Facchini, D. Fedele, Eleonora Fiorellino, Dirk Froebrich, Manuele Gangi, Teresa Giannini, K. Grankin, Hans Moritz Gunther, Zhen Guo, Lee Hartmann, Lynne A. Hillenbrand, P.C. Hinton, Joel H. Kastner, Chris Koen, K. Mauco, I. Mendigutia, B. Nisini, Neelam Panwar, et al. (8 additional authors not shown)

accepted to ApJ

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

E. Cristo, N. C. Santos, O. Demangeon, J. H. C. Martins, P. Figueira, N. Casasayas-Barris, M. R. Zapatero Osorio, F. Borsa, S. G. Sousa, M. Oshagh, G. Micela, H. M. Tabernero, J.V. Seidel, S. Cristiani, F. Pepe, R. Rebolo, V. Adibekyan, R. Allart, Y. Alibert, V. Bourrier, A. Cabral, E. Esparza Borges, J. I. González Hernández, J. Lillo-Box, G. Lo Curto, C. Lovis, A. Manescau, P. Di Marcantonio, C. J.A.P. Martins, A. S. Mascareño, D. Mégevand, A. Mehner, N. J. Nunes, E. Palle, T. Silva, A. Sozzetti, S. Udry

In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO.

Kun Hu, Matthew G. Baring, Joseph A. Barchas, George Younes

26 pages, 9 figures, accepted for publication in The Astrophysical Journal

The surfaces of neutron stars are sources of strongly polarized soft X rays due to the presence of strong magnetic fields. Radiative transfer mediated by electron scattering and free-free absorption is central to defining local surface anisotropy and polarization signatures. Scattering transport is strongly influenced by the complicated interplay between linear and circular polarizations. This complexity has been captured in a sophisticated magnetic Thomson scattering simulation we recently developed to model the outer layers of fully-ionized atmospheres in such compact objects, heretofore focusing on case studies of localized surface regions. Yet, the interpretation of observed intensity pulse profiles and their efficacy in constraining key neutron star geometry parameters is critically dependent upon adding up emission from extended surface regions. In this paper, intensity, anisotropy and polarization characteristics from such extended atmospheres, spanning considerable ranges of magnetic colatitudes, are determined using our transport simulation. These constitute a convolution of varied properties of Stokes parameter information at disparate surface locales with different magnetic field strengths and directions relative to the local zenith. Our analysis includes full general relativistic propagation of light from the surface to an observer at infinity. The array of pulse profiles for intensity and polarization presented highlights how powerful probes of stellar geometry are possible. Significant phase-resolved polarization degrees in the range of 10-60% are realized when summing over a variety of surface field directions. These results provide an important background for observations to be acquired by NASA's new IXPE X-ray polarimetry mission.

N. Vasantharaju, P. Vemareddy, B. Ravindra, V. H. Doddamani

Accepted for publication in The Astrophysical Journal on 14 January, 2022. Article contains 16 pages, 7 figures and 1 Table

The abrupt and permanent changes of photospheric magnetic field in the localized regions of active regions during solar flares called magnetic imprints (MIs), have been observed for the past nearly three decades. The well known "coronal implosion" model is assumed to explain such flare associated changes but the complete physical understanding is still missing and debatable. In this study, we made a systematic analysis of flare-related changes of photospheric magnetic field during 21 flares (14 eruptive and 7 non-eruptive) using the high-cadence (\texttt{135s}) vector-magnetogram data obtained from Helioseismic and Magnetic Imager. The MI regions for eruptive flares are found to be strongly localised, whereas the majority of non-eruptive events ($>70~\%$) have scattered imprint regions. To quantify the strength of the MIs, we derived the integrated change of horizontal field and total change of Lorentz force over an area. These quantities correlate well with the flare strength, irrespective of whether flares being eruptive or not, short or long duration. Further, the free-energy (FE), determined from virial-theorem estimates, exhibits statistically significant downward trend which starts around the flare time is observed in majority of flares. The change of FE during flares do not depend on eruptivity but have a strong positive correlation ($\approx 0.8$) with the Lorentz force change, indicating that the part of FE released would penetrate into the photosphere. While these results strongly favor the idea of significant feedback from corona on the photospheric magnetic field, the characteristics of MIs are quite indistinguishable for flares being eruptive or not.

M. Sánchez-García, S. García-Burillo, M. Pereira-Santaella, L. Colina, A. Usero, M. Querejeta, A. Alonso-Herrero, A. Fuente

30 pages, 34 figures, 4 tables. Accepted for publication in Astronomy & Astrophysics

We analyse the influence of the dynamical environment on the star formation (SF) relations of the dense molecular gas in the starburst (SB) ring of the Seyfert 2 galaxy NGC 1068. We used ALMA to image the emission of the 1-0 transitions of HCN and HCO+ with a resolution of 56 pc. We also used ancillary data of CO(1-0) at a resolution of ~100 pc, and CO(3-2) and its underlying continuum emission at ~40 pc. These observations allow us to probe a range of molecular gas densities (n(H2)~10$^{3-5}cm^{-3}$). The SF rate (SFR) is derived from Pa$\alpha$ line emission imaged by HST/NICMOS. We analysed how SF relations change depending on the choice of aperture sizes and molecular gas tracer. The scatter in the Kennicutt-Schmidt relation is about a factor of two to three lower for the HCN and HCO+ lines compared to CO(1-0) for a common aperture. Correlations lose statistical significance below a critical spatial scale $\approx$300-400 pc. The SF efficiency of the dense molecular gas (SFEdense) shows a scattered distribution as a function of the HCN luminosity (L'(HCN)) around a mean value of $\simeq0.01$Myr$^{-1}$. An alternative prescription for SF relations, linking the SFEdense and the boundedness of the gas measured by the parameter b$\equiv\Sigma$dense/$\sigma^2$, where $\Sigma$dense is the dense molecular gas surface density and $\sigma$ the velocity dispersion, resolves the degeneracy associated with the SFEdense-L'(HCN) plot. We identify two branches in the SFEdense-b plot that correspond to two dynamical environments in the SB ring, which are defined by their proximity to the bar-ring interface region. This region corresponds to the crossing of two density wave resonances, where an increased rate of cloud-cloud collisions would favour an enhanced compression of molecular gas. Our results suggest that galactic dynamics plays a major role in the efficiency of the gas conversion into stars.

Roberto Silvotti, Péter Németh, John H. Telting, Andrzej S. Baran, Roy H. Østensen, Jakub Ostrowski, Sumanta K. Sahoo, Saskia Prins

Accepted for publication in MNRAS Main Journal

TIC137608661/TYC4544-2658-1/FBS0938+788 is a new sdBV+dM reflection-effect binary discovered by the TESS space mission with an orbital period of 7.21 hours. In addition to the orbital frequency and its harmonics, the Fourier transform of TIC137608661 shows many g-mode pulsation frequencies from the sdB star. The amplitude spectrum is particularly simple to interpret as we immediately see several rotational triplets of equally spaced frequencies. The central frequencies of these triplets are equally spaced in period with a mean period spacing of 270.12 s, corresponding to consecutive l=1 modes. From the mean frequency spacing of 1.25 {\mu}Hz we derive a rotation period of 4.6 days in the deep layers of the sdB star, significantly longer than the orbital period. Among the handful of sdB+dM binaries for which the sdB rotation was measured through asteroseismology, TIC137608661 is the non-synchronized system with both the shortest orbital period and the shortest core rotation period. Only NY Vir has a shorter orbital period but it is synchronized. From a spectroscopic follow-up of TIC137608661 we measure the radial velocities of the sdB star, determine its atmospheric parameters, and estimate the rotation rate at the surface of the star. This measurement allows us to exclude synchronized rotation also in the outer layers and suggests a differential rotation, with the surface rotating faster than the core, as found in few other similar systems. Furthermore, an analysis of the spectral energy distribution of TIC137608661, together with a comparison between sdB pulsation properties and asteroseismic models, gives us further elements to constrain the system.

Edward Belbruno, James Green

11 pages, 2 figures

A resultant gravitational force due to the current estimates of the virial mass of the Milky Way galaxy, dominated by dark matter, is estimated near the Sun and is described in two different analytical models yielding consistent results. One is a two-step Hernquist model, the other is a Navarro-Frenk-White model. The effect of this force is estimated on trajectories for spacecraft sufficiently far from the Sun. The difficulty of detecting this force is studied. It is concluded that its effect should be considered for certain spacecraft missions. Its effect on the Pioneer and New Horizons spacecrafts is discussed. A future mission is discussed that may be able to detect this force. Implications of this force are discussed with its impact for problems in planetary astronomy and astrophysics.

Dominic Adams, Vihang Mehta, Hugh Dickinson, Claudia Scarlata, Lucy Fortson, Sandor Kruk, Brooke Simmons, Chris Lintott

23 pages, 13 figures, 4 tables, submitted to ApJ

Massive, star-forming clumps are a common feature of high-redshift star-forming galaxies. How they formed, and why they are so rare at low redshift, remains unclear. In this paper we identify the largest yet sample of clumpy galaxies (7,052) at low redshift using data from the citizen science project \textit{Galaxy Zoo: Clump Scout}, in which volunteers classified over 58,000 Sloan Digital Sky Survey (SDSS) galaxies spanning redshift $0.02 < z < 0.15$. We apply a robust completeness correction by comparing with simulated clumps identified by the same method. Requiring that the ratio of clump-to-galaxy flux in the SDSS $u$ band be greater than 8\% (similar to clump definitions used by other works), we estimate the fraction of local galaxies hosting at least one clump ($f_{clumpy}$) to be $2.68_{-0.30}^{+0.33}\%$. We also compute the same fraction with a less stringent cut of 3\% ($11.33_{-1.16}^{+0.89}\%$), as the higher number count and lower statistical noise of this fraction permits sharper comparison with future low-redshift clumpy galaxy studies. Our results reveal a sharp decline in $f_{clumpy}$ over $0 < z < 0.5$. The minor merger rate remains roughly constant over the same span, so we suggest that minor mergers are unlikely to be the primary driver of clump formation. Instead, the rate of galaxy turbulence is a better tracer for $f_{clumpy}$ over $0 < z < 1.5$ for galaxies of all masses, which supports the idea that clump formation is primarily driven by violent disk instability for all galaxy populations during this period.

František Dinnbier, Pavel Kroupa, Richard I. Anderson

16 pages, 9 figures, accepted to A&A

Some of the youngest stars (age $\lesssim 10$ Myr) are clustered, while many others are observed scattered throughout star forming regions or in complete isolation. It has been intensively debated whether the scattered or isolated stars originate in star clusters, or if they form truly isolated, which could help constrain the possibilities how massive stars are formed. We adopt the assumption that all stars form in gravitationally bound star clusters embedded in molecular cloud cores ($\Gamma$-$1\;$ model), which expel their natal gas, and compare the fraction of stars found in clusters with observational data. The star clusters are modelled by the code nbody6, which includes stellar and circumbinary evolution, gas expulsion, and the external gravitational field of their host galaxy. We find that small changes in the assumptions in the current theoretical model estimating the fraction, $\Gamma$, of stars forming in embedded clusters have a large influence on the results, and we present a counterexample as an illustration. This calls into question theoretical arguments about $\Gamma$ in embedded clusters, and it suggests that there is no firm theoretical ground for low $\Gamma$ in galaxies with lower star formation rates (SFRs). Instead, the assumption that all stars form in embedded clusters is in agreement with observational data for the youngest stars (age $\lesssim 10$ Myr). In the $\Gamma$-$1\;$ scenario, the observed fraction of the youngest stars in clusters increases with the SFR only weakly; the increase is caused by the presence of more massive clusters in galaxies with higher SFRs, which release fewer stars to the field in proportion to their mass. The $\Gamma$-$1\;$ model yields a higher fraction of stars in clusters for older stars (age between $10$ and $300$ Myr) than what is observed. This discrepancy can be caused by interactions with molecular clouds.

KwangHo Park, Gen Chiaki, John H. Wise

13 pages, 7 figures, 1 table, submitted to ApJ

We explore the effect of dust on the growth of seed black holes (BHs) in the early universe. Previous 1D radiation-hydrodynamic (RHD) simulations show that increased radiation pressure on dust further suppresses the accretion rate than the case for the chemically pristine gas. Using the Enzo+Moray code, we perform a suite of 3D RHD simulations of accreting BHs in a dusty interstellar medium (ISM). We use the modified Grackle cooling library to consider dust physics in its non-equilibrium chemistry. The BH goes through an early evolutionary phase, where ionizing BH radiation creates an oscillating HII region as it cycles between accretion and feedback. As the simulations proceed, dense cold gas accumulates outside the ionized region where inflow from the neutral medium meets the outflow driven by radiation pressure. In the late phase, high-density gas streams develop and break the quasi-spherical symmetry of the ionized region, rapidly boosting the accretion rate. The late phase is characterized by the coexistence of strong ionized outflows and fueling high-density gas inflows. The mean accretion rate increases with metallicity reaching a peak at Z$\sim$0.01-0.1$\,Z_\odot$, one order of magnitude higher than the one for pristine gas. However, as the metallicity approaches the solar abundance, the mean accretion rate drops as the radiation pressure becomes strong enough to drive out the high-density gas. Our results indicate that a dusty metal-poor ISM can accelerate the growth rate of BHs in the early universe, however, can stun its growth as the ISM is further enriched toward the solar abundance.

Hiroya Umeda, Masami Ouchi, Kimihiko Nakajima, Yuki Isobe, Shohei Aoyama, Yuichi Harikane, Yoshiaki Ono, Akinori Matsumoto

Strong high-ionization lines such as HeII of young galaxies are puzzling at high and low redshift. Although recent studies suggest the existence of non-thermal sources, whether their ionizing spectra can consistently explain multiple major emission lines remains a question. Here we derive the general shapes of the ionizing spectra for three local extremely metal-poor galaxies (EMPGs) that show strong HeII$\lambda$4686. We parameterize the ionizing spectra composed of a blackbody and power-law radiation mimicking various stellar and non-thermal sources. We use photoionization models for nebulae, and determine seven parameters of the ionizing spectra and nebulae by Markov Chain Monte Carlo methods, carefully avoiding systematics of abundance ratios. We obtain the general shapes of ionizing spectra explaining $\sim 10$ major emission lines within observational errors with smooth connections from observed X-ray and optical continua. We find that an ionizing spectrum of one EMPG has a blackbody-dominated shape, while the others have convex downward shapes at $>13.6$ eV, which indicate a diversity of the ionizing spectrum shapes. We confirm that the convex downward shapes are fundamentally different from ordinary stellar spectrum shapes, and that the spectrum shapes of these galaxies are generally explained by the combination of the stellar and ultra-luminous X-ray sources. Comparisons with stellar synthesis models suggest that the diversity of the spectrum shapes arises from differences in the stellar age. If galaxies at $z\gtrsim 6$ are similar to the EMPGs, high energy ($>54.4$ eV) photons of the non-stellar sources negligibly contribute to cosmic reionization due to relatively weak radiation.

Francesca Annibali, Monica Tosi (INAF-OAS Bologna)

This review has been published by Nature Astronomy on January 13$^{th}$ 2022, and is part of the Collection of short and long articles on Dwarf Galaxies that the Journal is publishing since December 13$^{th}$, 2021. This version of the article has been accepted for publication after peer review, but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections

Dwarf galaxies are the least massive, most abundant, and most widely distributed type of galaxies. Hence, they are key to testing theories of galaxy and Universe evolution. Dwarf galaxies sufficiently close to have their gas and stellar components studied in detail are of particular interest, because their properties and their evolution can be inferred with accuracy. This Review summarizes what is known of the stellar and chemical properties of star-forming dwarf galaxies closer than $\sim$20 Mpc. Given their low metallicity, high gas content and ongoing star formation, these objects are supposed to resemble the first galaxies that formed at the earliest epochs, and may thus represent a window on the distant, early Universe. We describe the major results obtained in the past decade on the star formation histories, chemical abundances, galaxy formation and evolution of star-forming dwarfs, and the uncertainties still affecting these results.

R. S. Nascimento, A. Rodríguez-Ardila, L. Dahmer-Hahn, M. A. Fonseca-Faria, R. Riffel, M. Marinello, T. Beuchert

17 pages, 9 figures, 5 tables, accepted for publication in MNRAS

In this work, we study the optical properties of compact radio sources selected from the literature in order to determine the impact of the radio-jet in their circumnuclear environment. Our sample includes 58 Compact Steep Spectrum (CSS) and GigaHertz Peaked Spectrum (GPS) and 14 Megahertz-Peaked spectrum (MPS) radio sources located at $z\leq 1$. The radio luminosity ($L_R$) of the sample varies between Log\,L$_R\sim$ 23.2 and 27.7 W\,Hz$^{-1}$. We obtained optical spectra for all sources from SDSS-DR12 and performed a stellar population synthesis using the {\sc starlight} code. We derived stellar masses (M$_\star$), ages $\langle t_\star \rangle$, star formation rates (SFR), metallicities $\langle Z_\star \rangle$ and internal reddening A$_V$ for all young AGNs of our sample. A visual inspection of the SDSS images was made to assign a morphological class for each source. Our results indicate that the sample is dominated by intermediate to old stellar populations and there is no strong correlation between optical and radio properties of these sources. Also, we found that young AGNs can be hosted by elliptical, spiral and interacting galaxies, confirming recent findings. When comparing the optical properties of CSS/GPS and MPS sources, we do not find any significant difference. Finally, the Mid-Infrared WISE colours analysis suggest that compact radio sources defined as powerful AGNs are, in general, gas-rich systems.

Suma Murthy, Raffaella Morganti, Nissim Kanekar, Tom Oosterloo

15 pages, 6 figures. Accepted for publication in Astronomy & Astrophysics

We present a search for associated HI 21 cm absorption in a sample of 29 radio-loud active galactic nuclei (AGNs) at $0.7 < z < 1$, carried out with the upgraded Giant Metrewave Radio Telescope. We detect HI 21 cm absorption against none of our target AGNs, obtaining $3\sigma$ upper limits to the optical depth of $\lesssim$ 1% per 50 km s$^{-1}$ channel. The radio luminosity of our sources is lower than that of most AGNs searched for HI 21 cm absorption at similar redshifts in the literature, and, for all targets except two, the UV luminosity is below the threshold $10^{23}$ W Hz$^{-1}$, above which the HI in the AGN environment has been suggested to be completely ionised. We stacked the HI spectra to obtain a more stringent limit of $\approx 0.17$% per 50 km s$^{-1}$ channel on the average HI 21 cm optical depth of the sample. The sample is dominated by extended radio sources, 24 of which are extended on scales of tens of kiloparsecs. Including similar extended sources at $0.7 < z < 1.0$ from the literature, and comparing with a low-$z$ sample of extended radio sources, we find statistically significant ($\approx 3\sigma$) evidence that the strength of HI 21 cm absorption towards extended radio sources is weaker at $0.7<z<1.0$ than at $z < 0.25$, with a lower detection rate of HI 21 cm absorption at $0.7 < z < 1.0$. Redshift evolution in the physical conditions of HI is the likely cause of the weaker associated HI 21 cm absorption at high redshifts, due to either a low HI column density or a high spin temperature in high-$z$ AGN environments.

Alejandra Rocha-Solache, Iván Rodríguez-Montoya, David Sánchez-Argüelles, Itziar Aretxaga

26 pages, 12 figures, submitted to The Astrophysical Journal Supplement

The complex physics involved in atmospheric turbulence makes it very difficult for ground-based astronomy to build accurate scintillation models and develop efficient methodologies to remove this highly structured noise from valuable astronomical observations. We argue that a Deep Learning approach can bring a significant advance to treat this problem because of deep neural networks' inherent ability to abstract non-linear patterns over a broad scale range. We propose an architecture composed of long-short term memory cells and an incremental training strategy inspired by transfer and curriculum learning. We develop a scintillation model and employ an empirical method to generate a vast catalog of atmospheric noise realizations and train the network with representative data. We face two complexity axes: the signal-to-noise ratio (SNR) and the degree of structure in the noise. Hence, we train our recurrent network to recognize simulated astrophysical point-like sources embedded in three structured noise levels, with a raw-data SNR ranging from 3 to 0.1. We find that a slow and repetitive increase in complexity is crucial during training to obtain a robust and stable learning rate that can transfer information through different data contexts. We probe our recurrent model with synthetic observational data, designing alongside a calibration methodology for flux measurements. Furthermore, we implement a traditional matched filtering (MF) to compare its performance with our neural network, finding that our final trained network can successfully clean structured noise and significantly enhance the SNR compared to raw data and in a more robust way than traditional MF.

Although it is accepted that perfect-merging is not a realistic outcome of collisions, some researchers state that perfect-merging simulations can still be considered as quantitatively reliable representations of the final stage of terrestrial planet formation. Citing the work of Kokubo & Genda [ApJL, 714L, 21], they argue that the differences between the final planets in simulations with perfect-merging and those where collisions are resolved accurately are small, and it is, therefore, justified to use perfect-merging results as an acceptable approximation to realistic simulations. In this paper, we show that this argument does not stand. We demonstrate that when the mass lost during collisions is taken into account, the final masses of the planets will be so different from those obtained from perfect-merging that the latter cannot be used as a valid approximation. We carried out a large number of SPH simulations of embryo-embryo collisions and determined the amount of the mass and water lost in each impact. We applied the results to collisions in a typical perfect-merging simulation and showed that even when the mass-loss in each collision is as small as 10%, perfect-merging can, on average, overestimate the masses of the final planets by $\sim 35\%$ and their water-contents by more than 18%. Our analysis demonstrates that, while perfect-merging simulations are still a powerful tool in proving concepts, they cannot be used to make predictions, draw quantitative conclusions (especially about the past history of a planetary system) and serve as a valid approximation to, or in lieu of the simulations in which collisions are resolved accurately.

Raffaele D'Abrusco, David Zegeye, Giuseppina Fabbiano, Michele Cantiello, Maurizio Paolillo, Andreas Zezas

40 pages, 26 figures, accepted for publication in ApJ

We report the discovery of statistically significant spatial structures in the projected two-dimensional distributions of the Globular Cluster (GC) systems of 10 among the brightest galaxies in the Fornax Cluster. We use a catalog of GCs extracted from the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) Fornax Cluster Survey (ACSFCS) imaging data. We characterize the size, shape and location relative to the host galaxies of the GC structures and suggest a classification based on their morphology and location that is suggestive of different formation mechanisms. We also investigate the GC structures in the context of the positions of their host galaxies relative to the general spatial distributions of galaxies and intra-cluster globular clusters in the Fornax Cluster. We finally estimate the dynamical masses of the progenitors of some GC structures, under the assumption that they are the relics of past accretion events of satellite galaxies by their hosts.

Feiyu Yu, Xiangliang Kong, Fan Guo, Wenlong Liu, Zelong Jiang, Yao Chen, Joe Giacalone

Accepted for publication in The Astrophysical Journal Letters

Recent observations have shown that in many large solar energetic particle (SEP) events the event-integrated differential spectra resemble double power laws. We perform numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation, including protons and heavier ions. We find that for all ion species the energy spectra integrated over the simulation domain can be described by a double power law, and the break energy depends on the ion charge-to-mass ratio as $E_B \sim (Q/A)^\alpha$, with $\alpha$ varying from 0.16 to 1.2 by considering different turbulence spectral indices. We suggest that the double power law distribution may emerge as a result of the superposition of energetic particles from different source regions where the acceleration rates differ significantly due to particle diffusion. The diffusion and mixing of energetic particles could also provide an explanation for the increase of Fe/O at high energies as observed in some SEP events. Although further mixing processes may occur, our simulations indicate that either power-law break or rollover can occur near the Sun and predict that the spectral forms vary significantly along the shock front, which may be examined by upcoming near-Sun SEP measurements from Parker Solar Probe and Solar Orbiter.

Tomoya Kinugawa (ICRR), Takashi Nakamura (Kyoto U.), Hiroyuki Nakano (Ryukoku U.)

Comments welcome

Two neutron star (NS) black hole (BH) binaries, GW200105 and GW200115 found in the LIGO/Virgo O3b run have smaller BH mass of $6-9$$M_{\odot}$ which is consistent with Population I and II origin. Our population synthesis simulations using $10^6$ Population I and II binaries with appropriate initial parameters show consistent binary mass, event rate, and no detection of radio pulsar (PSR) and BH binaries in our Galaxy so far. Especially, we found possible progenitors of GW200105 and GW200115 which were formed at redshift $z=0.15$ and $z=1.6$ with binary mass of $(34M_{\odot},\, 9.2M_{\odot})$ and $(23.7M_{\odot},\, 10.6M_{\odot})$, respectively. The final masses of these binaries are $(6.85M_{\odot},\,2.14M_{\odot})$ and $(6.04M_{\odot},\,1.31M_{\odot})$ which look like $(9.0_{-1.7}^{+1.7}M_{\odot},\, 1.91_{-0.24}^{+0.33}M_{\odot})$ of GW200105 and $(5.9_{-2.5}^{+2.0}M_{\odot},\,1.44_{-0.29}^{+0.85}M_{\odot})$ of GW200115, respectively. We also estimate that $4-20$ PSR-BH binaries in our galaxy will be observed by SKA. The existence of NS-BH binaries in our galaxy can be confirmed in future SKA era.

Mahendran Sithamparam, Nirmell Satthiyasilan, Chen Chen, Tony Z Jia, Kuhan Chandru

21 pages, 2 figures

The Panspermia hypothesis posits that either life's building blocks (molecular Panspermia) or life itself (organism-based Panspermia) may have been interplanetary transferred to facilitate the Origins of Life (OoL) on a given planet, complementing several current OoL frameworks. Although many spaceflight experiments were performed in the past to test for potential terrestrial organisms as Panspermia seeds, it is uncertain whether such organisms will likely "seed" a new planet even if they are able to survive spaceflight. Therefore, rather than using organisms, using abiotic chemicals as seeds has been proposed as part of the molecular Panspermia hypothesis. Here, as an extension of this hypothesis, we introduce and review the plausibility of a polymeric material-based Panspermia seed (M-BPS) theoretical concept, where the type of polymeric material that can function as a M-BPS must be able to: 1) survive spaceflight, and 2) "function", i.e., contingently drive chemical evolution towards some form of abiogenesis once arriving on a foreign planet. We use polymeric gels as a model example of a potential M-BPS. Polymeric gels that can be prebiotically synthesized on one planet (such as polyester gels) could be transferred to another planet via meteoritic transfer, where upon landing on a liquid bearing planet, can assemble into structures containing cellular-like characteristics and functionalities. Such features presupposed that these gels can assemble into compartments through phase separation to accomplish relevant functions such as encapsulation of primitive metabolic, genetic and catalytic materials, exchange of these materials, motion, coalescence, and evolution. All of these functions can result in the gels' capability to alter local geochemical niches on other planets, thereby allowing chemical evolution to lead to OoL events.

Evaristus U. Iyida, Christian I. Eze, Finbarr C. Odo

11 Pages, 14 figures, 4 tables. Accepted for publication in Astrophysics and Space Science

The concept of the evolutionary sequence of jetted active galactic nuclei (AGNs) has been challenged in the last few decades since AGN subclasses are considered to be different due to their viewing angle. In this paper, we collected a sample of 1108 blazars (472 flat-spectrum radio quasars, FSRQs and 636 BL Lacertae objects, BL Lacs) and 120 Seyfert galaxies (SGs) with available redshifts and spectral properties in radio, optical, X-ray, and gamma-ray bands in order to compare the properties of SGs and the blazar subclasses and also explore their possible unification through evolution. It is found that the ratio of the relative difference of SGs and the blazar subclasses are approximately the same implying that they have an evolutionary relationship. We discovered from the results of a two-dimensional Kolmogorov-Smirnov (K.S.) test, that the probabilities (p) for the composite spectral indices: optical-X-ray radio-X-ray radio-optical and X-ray-gamma-ray to come from the same parent population is p, less than 0.05, implying that the null hypothesis cannot be rejected, thus, supports SGs-BL Lacs-FSRQs evolutionary sequence................. These results agree with the prediction that SGs have an evolutionary relationship with the blazar subclasses.

Planet 9 is currently a hypothetical planet the orbital parameters of which are based on anomalous orbits of Kuiper Belt objects. The orbital parameters are such that, if Planet 9 exists, the theory of solar barycentric dynamics would be profoundly altered. We show that, with Planet 9 included in the solar system, barycentric theory is a much more effective predictor of solar activity on the decadal, centennial and millennium time scales. In particular the most elementary quantity of barycentric theory, the Sun to planet barycentre distance, is more coherent with decadal solar activity cycles and grand solar activity minima than barycentric distance without Planet 9 included. Further, barycentric theory including Planet 9 contains strong components at periods corresponding to the Hallstaatt and Gleissberg cycles whereas barycentric theory without Planet 9 shows no evidence of these cycles. A challenge that emerged during this study was the absence of the strongest component in barycentric theory, the Jose component at ~ 178 year period, from the spectra of millennium scale solar activity. This conundrum was solved by demonstrating that, during the transformation from solar motion to solar activity, the Jose component in solar motion was, in the spectrum of solar activity, split into multiple sidebands due to phase modulation by lower frequency cycles. The excellent fit to solar activity at multiple time scales by barycentric theory with Planet 9 included is itself supporting evidence for the existence of Planet 9, specifically by providing an estimate of the current heliographic longitudinal location and orbital period.

Dacheng Lin, Jimmy A. Irwin, Edo Berger, Ronny Nguyen

13 pages, 9 figures, submitted version after addressing the referee's comments

It was proposed that a remnant stable magnetar could be formed in a binary neutron-star merger, leading to a fast X-ray transient (FXT) that can last for thousands of seconds. Recently, Xue et al. suggested that CDF-S XT2 was exactly such a kind of source. If confirmed, such emission can be used to search for electromagnetic counterparts to gravitational wave events from binary neutron-star mergers that have short gamma-ray bursts and the corresponding afterglows seen off-axis and thus too weak to be detected. Here we report the discovery of three new FXTs, XRT 170901, XRT 030511, and XRT 110919, from preliminary search over Chandra archival data. Similar to CDF-S XT2, these new FXTs had a very fast rise (less than a few ten seconds) and a plateau of X-ray flux $\sim$$1.0\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ lasting for 1-2 ks, followed by a steep decay. Their optical/IR counterparts, if present, are very weak, arguing against a stellar flare origin for these FXTs. For XRT 170901, we identified a faint host galaxy with the source at the outskirts, very similar to CDF-S XT2. Therefore, our newly discovered FXTs are also strong candidates for magnetar-powered X-ray transients resulting from binary neutron star mergers.

Minghao Yue, Xiaohui Fan, Jinyi Yang, Feige Wang

18 pages, 12 figures, accepted by AJ, comments welcomed

We present a mock catalog of gravitationally lensed quasars at $z_\text{qso}<7.5$ with simulated images for the Rubin Observatory Legacy Survey of Space and Time (LSST). We adopt recent measurements of quasar luminosity functions to model the quasar population, and use the CosmoDC2 mock galaxy catalog to model the deflector galaxies, which successfully reproduces the observed galaxy velocity dispersion functions up to $z_d\sim1.5$. The mock catalog is highly complete for lensed quasars with Einstein radius $\theta_E>0\farcs07$ and quasar absolute magnitude $M_{i}<-20$. We estimate that there are $\sim10^3$ lensed quasars discoverable in current imaging surveys, and LSST will increase this number to $\sim2.4\times10^3$. Most of the lensed quasars have image separation $\Delta \theta>0\farcs5$, which will at least be marginally resolved in LSST images with seeing of $\sim0\farcs7$. There will be $\sim200$ quadruply-lensed quasars discoverable in the LSST. The fraction of quad lenses among all discoverable lensed quasars is about $\sim10\%-15\%$, and this fraction decreases with survey depth. This mock catalog shows a large diversity in the observational features of lensed quasars, in terms of lensing separation and quasar-to-deflector flux ratio. We discuss possible strategies for a complete search of lensed quasars in the LSST era.

Very high energy {\gamma}-rays are one of the most important messengers of the non-thermal Universe. The major motivation of very high energy {\gamma}-ray astronomy is to find sources of high energy cosmic rays. Several astrophysical sources are known to accelerate cosmic rays to very high energies under extreme conditions. Very high energy {\gamma}-rays are produced at these astrophysical sites or near through interactions of cosmic rays in the surrounding medium close to the sources. Gamma-rays, being neutral, travel in a straight line and thus give us valuable information about the cosmic ray sources and their surroundings. Additionally, very high energy {\gamma}-ray astronomy can probe many fundamental physics questions. Ground-based {\gamma}-ray astronomy began its journey in 1989 when Whipple telescope detected TeV {\gamma}-rays from the Crab, a pulsar wind nebula in the Milky Way. In the last two decades, technological improvements have facilitated the development of the latest generation of very high energy detectors and telescopes which have delivered exciting new results. Until now over two hundred very high energy {\gamma}-ray sources, both galactic and extra-galactic has been detected. These observations have provided a deeper insight into a large number of important questions in high energy astrophysics and astroparticle physics. This review article is an attempt to enumerate the most important results in the exciting and rapidly developing field of very high energy {\gamma}-ray astronomy.

Isaac Legred, Katerina Chatziioannou, Reed Essick, Philippe Landry

20 pages, 12 figures

The rapid increase in the number and precision of astrophysical probes of neutron stars in recent years allows for the inference of their equation of state. Observations target different macroscopic properties of neutron stars which vary from star to star, such as mass and radius, but the equation of state allows for a common description of all neutron stars. To connect these observations and infer the properties of dense matter and neutron stars simultaneously, models for the equation of state are introduced. Parametric models rely on carefully engineered functional forms that reproduce a large array of realistic equations of state. Such models benefit from their simplicity but are limited because any finite-parameter model cannot accurately approximate all possible equations of state. Nonparametric models overcome this by increasing model freedom at the cost of increased complexity. In this study, we compare common parametric and nonparametric models, quantify the limitations of the former, and study the impact of modeling on our current understanding of high-density physics. We show that parametric models impose strongly model-dependent, and sometimes opaque, correlations between density scales. Such inter-density correlations result in tighter constraints that are unsupported by data and can lead to biased inference of the equation of state and of individual neutron star properties.

Moun Meenakshi, Dipanjan Mukherjee, Alexander Y. Wagner, Nicole P. H. Nesvadba, Raffaella Morganti, Reinier M. J. Janssen, Geoffrey V. Bicknell

Accepted for publication in MNRAS

We use the results of relativistic hydrodynamic simulations of jet-ISM interactions in a galaxy with a radio-loud AGN to quantify the extent of ionization in the central few kpcs of the gaseous galactic disc. We perform post-process radiative transfer of AGN radiation through the simulated gaseous jet-perturbed disc to estimate the extent of photo-ionization by the AGN with an incident luminosity of $10^{45}~\mathrm{erg\,s^{-1}}$. We also map the gas that is collisionally ionized due to shocks driven by the jet. The analysis was carried out for simulations with similar jet power ($10^{45}~\mathrm{erg\,s^{-1}}$) but different jet orientations with respect to the gas disc. We find that the shocks from the jets can ionize a significant fraction (up to 33$\%$) of dense gas ($n>100\,\mathrm{cm^{-3}}$) in the disc, and that the jets clear out the central regions of gas for AGN radiation to penetrate to larger distances in the disc. Jets inclined towards the disc plane couple more strongly with the ISM and ionize a larger fraction of gas in the disc as compared to the vertical jet. However, similar to previous studies, we find that the AGN radiation is quickly absorbed by the outer layers of dense clouds in the disc, and is not able to substantially ionize the disc on a global scale. Thus, compared to jet-ISM interactions, we expect that photo-ionization by the AGN radiation only weakly affects the star-formation activity in the central regions of the galactic disc ($\lesssim 1$ kpc), although the jet-induced shocks can spread farther out.

Rui Wang, Ying D. Liu, Shangbin Yang, Huidong Hu

Accepted for publication in ApJ. 31 pages, 12 figures

Homologous coronal mass ejections (CMEs) are an interesting phenomenon, and it is possible to investigate the formation of CMEs by comparing multi-CMEs under a homologous physical condition. AR 11283 had been present on the solar surface for several days when a bipole emerged on 2011 September 4. Its positive polarity collided with the pre-existing negative polarity belonging to a different bipole, producing recurrent solar activities along the polarity inversion line (PIL) between the colliding polarities, namely the so-called collisional PIL (cPIL). Our results show that a large amount of energy and helicity were built up in the form of magnetic flux ropes (MFRs), with recurrent release and accumulation processes. These MFRs were built up along the cPIL. A flux deficit method is adopted and shows that magnetic cancellation happens along the cPIL due to the collisional shearing scenario proposed by Chintzoglou et al. The total amount of canceled flux was $\sim$0.7$\times$10$^{21}$ Mx with an uncertainty of $\sim$13.2$\%$ within the confidence region of the 30$^\circ$ sun-center distance. The canceled flux amounts to 24$\%$ of the total unsigned flux of the bipolar magnetic region. The results show that the magnetic fields beside the cPIL are very sheared, and the average shear angle is above 70$^\circ$ after the collision. The fast expansion of the twist kernels of the MFRs and the continuous eruptive activities are both driven by the collisional shearing process. These results are important for better understanding the buildup process of the MFRs associated with homologous solar eruptions.

Fataneh Salahedin, Reza Pazhouhesh, Mohammad Malekjani

15 pages,Eur. Phys. J. Plus

We use different combinations of data samples to investigate the new generalized Chaplygin gas (NGCG) model in the context of dark energy (DE) cosmology. Using the available cosmological data, we put constraints on the the free parameters of NGCG model based on the statistical Markov chain Monte Carlo method. We then find the best fit values of cosmological parameters and those confidence regions in NGCG cosmology. Our result for the matter density parameter calculated in NGCG model is in excellent agreement with that of the standard CDM cosmology. We also find that the equation of state of DE of the model slightly favors the phantom regime. We show that the big tension between the low- and high-redshift observations appearing in CDM universe to predict the Hubble constant H 0 can be alleviated in NGCG model. However, from the statistical point of view, our results show that the standard CDM model fits the observations better than the NGCG cosmology.

Alexia M. Lopez, Roger G. Clowes, Gerard M. Williger

Submitted to MNRAS

We present the serendipitous discovery of a 'Giant Arc on the Sky' at $z \sim 0.8$. The Giant Arc (GA) spans $\sim 1$ Gpc (proper size, present epoch), and appears to be almost symmetrical on the sky. It was discovered via intervening MgII absorbers in the spectra of background quasars, using the catalogues of Zhu & M\'enard. The use of MgII absorbers represents a new approach to the investigation of large-scale structures (LSSs) at redshifts $0.45 \lesssim z \lesssim 2.25$. We present the observational properties of the GA, and we assess it statistically using methods based on: (i) single-linkage hierarchical clustering ($\sim 4.5\sigma$); (ii) the Cuzick-Edwards test ($\sim 3.0\sigma$); and (iii) power spectrum analysis ($\sim 4.8\sigma$). Each of these methods has distinctive attributes and powers, and we advise considering the evidence from the ensemble. The overdensity of the GA is $\delta \rho / \rho \sim 1.3 \pm 0.3$. The GA is the newest and one of the largest of a steadily accumulating set of very large LSSs that may (cautiously) challenge the Cosmological Principle, upon which the 'standard model' of cosmology is founded. Conceivably, the GA is the precursor of a structure like the Sloan Great Wall (but the GA is about twice the size), seen when the Universe was about half its present age.

V.G.Gurzadyan, N.N.Fimin, V.M.Chechetkin

8 pages, 1 fig., Eur.Phys.J.Plus (in press)

The emergence of one and two-dimensional configurations -- Zeldovich pancakes -- progenitors of the observed filaments and clusters and groups of galaxies, is predicted by means of a developed kinetic approach in analyzing the evolution of initial density perturbations. The self--consistent gravitational interaction described by Vlasov-Poisson set of equations with branching conditions is shown to predict two--dimensional structures as of layers of increased density and voids between them, i.e. the cellular macro-structure of the Universe. The modified potential of weak-field General Relativity is involved, which enables one to explain the Hubble tension, revealing the conceptual discrepancy in the local galactic flows and the cosmological expansion. This demonstrates the possible essential role of self-consistent gravity in the formation of the cosmic web.

Maryna V. Ishchenko, Margarita O. Sobolenko, Mukhagali T. Kalambay, Bekdaulet T. Shukirgaliyev, Peter P. Berczik

9 pages, 9 figures, accepted for publication in Reports of National academy of sciences of the Republic of Kazakhstan

Using the data from Gaia (ESA) Data Release 2 we performed the orbital calculations of globular clusters (GCs) of the Milky Way. To explore possible collisions between the GCs, using our developed highorder {\phi}-GRAPE code, we integrated (backwards and forward) the orbits of 119 objects with reliable positions and proper motions. In calculations, we adopted a realistic axisymmetric Galactic potential (bulge + disk + halo). Using different impact conditions, we found five pairs of the GCs that likely experienced collisions: Terzan 3 - NGC 6553, Terzan 3 - NGC 6218, Liller 1 - NGC 6522, Djorg 2 - NGC 6552 and NGC 6355 - NGC 6637. We analyzed the GCs interaction rates with the central supermassive black hole. Assuming the maximum 100 pc distance criteria for separation between them we estimated 11 close encounter events. From our numerical simulations, we estimate the close interaction rate as at least one event per Gyr with the impact parameter less than 30 pc; and one event per Myr with the impact parameter less than 60 pc. Our calculations show one very close encounter of NGC 6121 with the central SMBH near 5.5 pc (practically direct collision). Based on the extended literature search for the possible progenitor of our selected 11 GCs, we found that most of them have a Milky Way main bulge origin.

P. Vemareddy, P. Démoulin, K. Sasikumar Raja, J. Zhang, N. Gopalswamy, N. Vasantharaju

12 pages, 6 figures, Accepted to publish in The Astrophysical Journal

Using the observations from Solar Dynamics Observatory, we study an eruption of a hot-channel flux rope (FR) near the solar-limb on February 9, 2015. The pre-eruptive structure is visible mainly in EUV 131 $\mathring{\mathrm{A}}$ images with two highly-sheared loop structures. They undergo slow rise motion and then reconnect to form an eruptive hot-channel as in the tether-cutting reconnection model. The J-shaped flare-ribbons trace the footpoint of the FR which is identified as the hot-channel. Initially, the hot channel is observed to rise slowly at 40 km s$^{-1}$, followed by an exponential rise from 22:55 UT at a coronal height of 87$\pm$2 Mm. Following the onset of the eruption at 23:00 UT, the flare-reconnection adds to the acceleration process of the CME within 3 R$_\odot$. Later on, the CME continues to accelerate at 8 m s$^{-2}$ during its propagation period. Further, the eruption launched type-II followed by III, IVm radio bursts. The start and end times of type-IVm correspond to the CME core height of 1.5 and 6.1 R$_\odot$, respectively. Also the spectral index is negative suggesting the non-thermal electrons trapped in the closed loop structure. Accompanied with type-IVm, this event is unique in the sense that the flare ribbons are very clearly observed along with the erupting hot channel, which strongly supports that the hooked-part of J-shaped flare ribbons outlines the boundary of the erupting FR.

Animikh Roy (University of Sussex, UK), Andor Kesselman (Pathr.ai, USA)

Graph Theory and Topological Data Analytics, while powerful, have many drawbacks related to their sensitivity and consistency with TDA & Graph Network Analytics. In this paper, we aim to propose a novel approach for encoding vectorized associations between data points for the purpose of enabling smooth transitions between Graph and Topological Data Analytics. We conclusively reveal effective ways of converting such vectorized associations to simplicial complexes representing micro-states in a Phase-Space, resulting in filter specific, homotopic self-expressive, event-driven unique topological signatures which we have referred as Roy-Kesselman Diagrams or R-K Diagrams with persistent homology, which emerge from filter-based encodings of R-K Models. The validity and impact of this approach were tested specifically on high-dimensional raw and derived measures of Gravitational Wave Data from the latest LIGO datasets published by the LIGO Open Science Centre along with testing a generalized approach for a non-scientific use-case, which has been demonstrated using the Tableau Superstore Sales dataset. We believe the findings of our work will lay the foundation for many future scientific and engineering applications of stable, high-dimensional data analysis with the combined effectiveness of Topological Graph Theory transformations.

Davide Gualdi, Licia Verde

17 pages plus appendix, 11 figures, prepared for submission to JCAP, comments are welcome!

We present the first detection of the integrated trispectrum ($\mathit{i}$-trispectrum) monopole and quadrupoles signal from BOSS CMASS NGC DR12. Extending the FKP estimators formalism to the Fourier transform of the four-point correlation function, we test shot-noise subtraction, Gaussianity of the i-trispectrum data-vector, significance of the detection and similarity between the signal from the data and from the galaxy mock catalogues used to numerically estimate the covariance matrix. Using scales corresponding to modes from minimum $k_\mathrm{min}=0.03\,h/\mathrm{Mpc}$ to maximum $k_\mathrm{min}=0.15\,h/\mathrm{Mpc}$, we find a detection in terms of distance from the null hypothesis of $(10.4,5.2,8.3,1.1,3.1)$ $\sigma$-intervals for the i-trispectrum monopole and quadrupoles respectively. This quantifies the presence of the physical signal of the four-points statistics on BOSS data. For completeness the same analysis is also performed for power spectrum and bispectrum, both monopoles and quadrupoles.

Mathieu Barthelemy, Elisa Robert, Vladimir Kalegaev, Vincent Grennerat, Thierry Sequies, Guillaume Bourdarot, Etienne Le Coarer, Jean-Jacques Correia, Patrick Rabou

Submitted to IEEE-JMASS. Under review

AMICal sat, a dedicated 2U cubesat, has been developed, in order to monitor the auroral emissions, with a dedicated imager. It aims to help to reconstruct the low energy electrons fluxes up to 30 keV in Earth auroral regions. It includes an imager entirely designed in Grenoble University Space Center. The imager uses a 1.3 Mpixels sparse RGB CMOS detector and a wide field objective (f=22.5 mm). The satellite platform has been built by the polish company Satrevolution. Launched September, 3rd, 2020 from Kuru (French Guyana) on board the Vega flight 16, it produces its first images in October 2020. The aim of this paper is to describe the design of the payload especially the optics and the proximity electronics, to describe the use of the payload for space weather purpose. A preliminary analysis of a first image showing the relevance of such an instrument for auroral monitoring is performed. This analysis allowed to reconstruct from one of the first images the local electron input flux at the top of the atmosphere during the exposure time.

K. Medler, P. A. Mazzali, J. Teffs, C. Ashall, J.P. Anderson, I. Arcavi, S. Benetti, K. A. Bostroem, J. Burke, Y.-Z. Cai, P. Charalampopoulos, N. Elias-Rosa, M. Ergon, L. Galbany, M. Gromadzki, D. Hiramatsu, D. A. Howell, C. Inserra, P. Lundqvist, C. McCully, T. Müller-Bravo, M. Newsome, M. Nicholl, E. Padilla Gonzalez, E. Paraskeva, A. Pastorello, C. Pellegrino, P. J. Pessi, A. Requitti, T. M. Reynolds, R. Roy, G. Terreran, L. Tomasella, D. R. Young

18 pages, 13 figures

The Ultra-Violet (UV) and Near Infrared (NIR) photometric and optical spectroscopic observations of SN 2020acat covering $\sim \! \! 250$ days after explosion are presented here. Using the fast rising photometric observations, spanning from the UV to NIR wavelengths, a pseudo-bolometric light curve was constructed and compared to several other well-observed Type IIb supernovae (SNe IIb). SN 2020acat displayed a very short rise time reaching a peak luminosity of $\mathrm{Log_{10}}(L) = 42.49 \pm 0.15 \, \mathrm{erg \, s^{-1}}$ in only $\sim \! \! 14.6 \pm 0.3$ days. From modelling of the pseudo-bolometric light curve, we estimated a total mass of $^{56} \mathrm{Ni}$ synthesised by SN 2020acat of $0.13 \pm 0.02 \, \mathrm{M_{\odot}}$, with an ejecta mass of $2.3 \pm 0.3 \, \mathrm{M_{\odot}}$ and a kinetic energy of $1.2 \pm 0.2 \times 10^{51}$ erg. The optical spectra of SN 2020acat display hydrogen signatures well into the transitional period ($\gtrsim 100$ days), between the photospheric and the nebular phases. The spectra also display a strong feature around $4900 \, \r{A}$ that cannot be solely accounted for by the presence of the $\mathrm{Fe_{II}}$ $5018$ line. We suggest that the $\mathrm{Fe_{II}}$ feature was augmented by $\mathrm{He_{I}}$ $5016$ and possibly by the presence of $\mathrm{N_{II}}$ $5005$. From both photometric and spectroscopic analysis, we inferred that the progenitor of SN\,2020acat was an intermediate mass compact star with a $M_\mathrm{ZAMS}$ of $18 - 22 \, \mathrm{M_{\odot}}$.

Tidal heating is often used to interpret "radius anomaly" of hot Jupiters (i.e. radii of a large fraction of hot Jupiters are in excess of 1.2 Jupiter radius which cannot be interpreted by the standard theory of planetary evolution). In this paper we find that tidal heating induces another phenomenon "runaway inflation" (i.e. planet inflation becomes unstable and out of control when tidal heating rate is above its critical value). With sufficiently strong tidal heating, luminosity initially increases with inflation, but across its peak it decreases with inflation such that heating is stronger than cooling and runaway inflation occurs. In this mechanism, the opacity near radiative-convective boundary (RCB) scales approximately as temperature to the fourth power and heat cannot efficiently radiate away from planet interior, which induces runaway inflation (similar to a tight lid on a boiling pot). Based on this mechanism, we find that radii of hot Jupiters cannot exceed $2.2R_J$, which is in good agreement with the observations. We also give an upper limit for orbital eccentricity of hot Jupiters. Moreover, by comparison to the observations we infer that tidal heating locates near RCB.

Jaime Ruiz-Zapatero, Carlos García-García, David Alonso, Pedro G. Ferreira, Richard D.P. Grumitt

Prepared for submission to MNRAS

We constrain the expansion history of the Universe and the cosmological matter density fraction in a model-independent way by exclusively making use of the relationship between background and perturbations under a minimal set of assumptions. We do so by employing a Gaussian process to model the expansion history of the Universe from present time to the recombination era. The expansion history and the cosmological matter density are then constrained using recent measurements from cosmic chronometers, Type-Ia supernovae, baryon acoustic oscillations, and redshift-space distortion data. Our results show that the evolution in the reconstructed expansion history is compatible with the \textit{Planck} 2018 prediction at all redshifts. The current data considered in this study can constrain a Gaussian process on $H(z)$ to an average $9.4 \%$ precision across redshift. We find $\Omega_m = 0.224 \pm 0.066$, lower but statistically compatible with the \textit{Planck} 2018 cosmology. Finally, the combination of future DESI measurements with the CMB measurement considered in this work holds the promise of $8 \%$ average constraints on a model-independent expansion history as well as a five-fold tighter $\Omega_m$ constraint using the methodology developed in this work.

Wojciech Szewczuk, Przemysław Walczak, Jadwiga Daszyńska-Daszkiewicz, Dawid Moździerski

Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

KIC8264293 is a fast rotating B-type pulsator observed by Kepler satellite. Its photometric variability is mainly due to pulsations in high-order g modes. Besides, we detected a weak H$\alpha$ emission. Thus, the second source of variability is the fluctuation in a disk around the star. The pulsational spectrum of KIC8264293 reveals a frequency grouping and period spacing pattern. Here, we present the thorough seismic analysis of the star based on these features. Taking into account the position of the star in the HR diagram and fitting the 14 frequencies that form the period spacing we constrain the internal structure of the star. We conclude that the star barely left the ZAMS and the best seismic model has $M = 3.54\,\mathrm{M}_\odot$, $V_\mathrm{rot}=248\,\mathrm{km\,s}^{-1}$ and $Z = 0.0112$. We found the upper limit on the mixing at the edge of the convective core, with the overshooting parameter up to $f_\mathrm{ov} = 0.03$. On the other hand, we were not able to constrain the envelope mixing for the star. To excite the modes in the observed frequency range, we had to modify the opacity data. Our best seismic model with an opacity increase by 100% at the "nickel" bump $\log T=5.46$ explains the whole instability. KIC8264293 is the unique, very young star pulsating in high-order g modes with the Be feature. However, it is not obvious whether the source of this circumstellar matter is the ejection of mass from the underlying star or whether the star has retained its protostellar disk.

Theories of stellar convective core overshoot can be examined through analysis of pulsating stars. Better accuracy can be achieved by obtaining external constraints such as those provided by observing pulsating stars in eclipsing binary systems, but this requires that the binary parameters be identified so photometric variations of the pulsating component may be isolated from the binary periodicity. This study aims to uncover the physical parameters of three binaries observed by the Kepler spacecraft. We also seek to evaluate the feasibility of accurately constraining binaries using only readily available time-series photometry and distance estimates. Binary models were constructed using the Physics of Eclipsing Binaries (PHOEBE) software package. Markov Chain Monte Carlo methods were used to sample the parameter space of these models and provide estimates of the posterior distributions for these systems. An initial run using binned light curve data was performed to identify general parameter trends and provide initializing distributions for a subsequent analysis incorporating the full data set. We present theoretical models for all three binaries, along with posterior distributions from our MCMC analyses. Models for KIC 8314879 and KIC 10727668 produced a good match to the observed data, while the model of KIC 5957123 failed to generate an appropriate synthetic light curve. For the two successful models, we interpret the posterior distributions and discuss confidence in our parameter estimates and uncertainties. We also evaluate the feasibility of this procedure in various contexts, and propose several modifications to improve the success of future studies.

We report the first detection in the ${\it J}_{{\it K}}$ = 1$_{0}$ - 0$_{0}$ rotational transition line of CH$_{3}$D towards three Class 0/I proto-brown dwarfs (proto-BDs) from IRAM 30 m observations. Assuming a rotational temperature of 25 K, the CH$_{3}$D abundances (relative to H$_{2}$) are in the range of (2.3--14.5) $\times$ 10$^{-7}$. The CH$_{4}$ abundances derived from the CH$_{3}$D abundances and assuming the DCO$^+$/HCO$^+$ ratios are in the range of (0.05--4.8) $\times $10$^{-5}$. The gas-phase formation of CH$_{3}$D via CH$_2$D$^+$ is enhanced at high densities of 10$^{8}$--10$^{10}$ cm$^{-3}$ and our observations are likely probing the innermost dense and warm regions in proto-BDs. Thermal and/or non-thermal desorption can return the CH$_{3}$D and CH$_{4}$ molecules formed at an early stage on grain surfaces to the gas-phase. The gas phase abundances are indicative of warm carbon-chain chemistry in proto-BDs where carbon-chain molecules are synthesized in a lukewarm ($\sim$20-30 K) region close to the central source.

In the studies on pre-recombination early dark energy (EDE), the evolution of Universe after recombination is usually regarded as ${\Lambda}CDM$-like, which corresponds that the equation of state of dark energy responsible for current accelerated expansion is $w=-1$. However, in realistic models, $w$ might be evolving. We consider the parametrizations of $w$ with respect to the redshift $z$ in Axion-like EDE and AdS-EDE models, respectively. We performed the Monte Carlo Markov chain analysis with recent cosmological data, and found that the bestfit $w(z)$ is compatible with $w_0=-1,w_a=0$ (the cosmological constant) and the evolution of $w$ is only marginally favored, which so has little effect on lifting the bestfit value of ${H_0}$.

Ameek Sidhu (1 and 2), Josh Bazely (1), Els Peeters (1, 2 and 3), Jan Cami (1, 2, and 3) ((1) Department of Physics & Astronomy, University of Western Ontario, London, ON, Canada, (2) Institute for Earth and Space Exploration, University of Western Ontario, London, ON, Canada, (3) SETI Institute, Mountain View, CA, USA)

15 pages, 12 figures, accepted for publication in MNRAS

We carried out a principal component analysis (PCA) of the fluxes of five polycyclic aromatic hydrocarbon (PAH) bands at 6.2, 7.7, 8.6, 11.0, and 11.2 $\mu$m in the reflection nebula NGC 7023 comprising of the photodissociation region (PDR) and a cavity. We find that only two principal components (PCs) are required to explain the majority of the observed variance in PAH fluxes (98 %). The first PC ($PC_{1}$), which is the primary driver of the variance, represents the total PAH emission. The second PC ($PC_{2}$) is related to the ionization state of PAHs across the nebula. This is consistent with the results of a similar analysis of the PAH emission in NGC 2023. The biplots and the correlations of PCs with the various PAH ratios show that there are two subsets of ionic bands with the 6.2 and 7.7 $\mu$m bands forming one subset and the 8.6 and 11.0 $\mu$m bands the other. However, the distinction between these subsets is only present in the PDR. We have also carried out a separate PCA analysis of the PAH fluxes, this time only considering variations in the cavity. This shows that in the cavity, $PC_{2}$ is not related to the charge state of PAHs, but possibly to structural molecular changes.

L. Pizzuti, I. D. Saltas, A. Biviano, G. Mamon, L. Amendola

23 pages, two figures. Techinical manual to be submitted to JOSS

The \textsc{MG-MAMPOSSt} code is a license-free \textsc{Fortran95} code to perform tests of General Relativity (GR) through the analyses of kinematical data of galaxy clusters based on the Jeans' equation. The code is based on the \textsc{MAMPOSSt} method, and extends the original code through new parametrisations of the gravitational potential for general families of gravity theories beyond GR aimed to explain dark energy. \textsc{MG-MAMPOSSt} is further supplemented with a new capability to produce weak lensing forecasts for joint kinematic+lensing analysis. This document provides a technical description of the code's new features, functionality with respect to the original version, and instructions on its installation and use. Finally, we explain how the code could be further modified to include a wider family of gravity models and/or density profiles, that could allow its application in broader theoretical frameworks as well as other physical systems such as stellar clusters. A detailed forecast analysis for the modified gravity models currently implemented in the code can be found in the paper of Pizzuti et al., 2021.

Imaging a supermassive black hole and extracting physical information requires good knowledge of both the gravitational and the astrophysical conditions near the black hole. When the geometrical properties of the black hole are well understood, extracting information on the emission properties is possible. Similarly, when the emission properties are well understood, extracting information on the black-hole geometry is possible. At present however, uncertainties are present both in the geometry and in the emission, and this inevitably leads to degeneracies in the interpretation of the observations. We explore here the impact of varying geometry and emission coefficient when modelling the imaging of a spherically-accreting black hole. Adopting the Rezzolla-Zhidenko parametric metric to model arbitrary static black-holes, we first demonstrate how shadow-size measurements leave degeneracies in the multidimensional space of metric-deviation parameters, even in the limit of infinite-precision measurements. Then, at finite precision, we show that these degenerate regions can be constrained when multiple pieces of information, such as the shadow-size and the peak image intensity contrast, are combined. Such degeneracies can potentially be eliminated with measurements at increased angular-resolution and flux-sensitivity. While our approach is restricted to spherical symmetry and hence idealised, we expect our results to hold also when more complex geometries and emission processes are considered.

Isabell Stucke, Deborah Morgenstern, Achim Zeileis, Georg J. Mayr, Thorsten Simon, Gerhard Diendorfer, Wolfgang Schulz, Hannes Pichler

Despite its scarcity, upward lightning initiated from tall structures causes more damage than common downward lightning. One particular subtype with a continuous current only is not detectable by conventional lightning location systems (LLS) causing a significantly reduced detection efficiency. Upward lightning has become a major concern due to the recent push in the field of renewable wind energy generation . The growing number of tall wind turbines increased lightning related damages. Upward lightning may be initiated by the tall structure triggering the flash itself (self-triggered) or by a flash striking close by (other-triggered). The major objective of this study is to find the driving atmospheric conditions influencing whether an upward flash is self-triggered or other-triggered and whether it is of the undetectable subtype. We explore upward flashes directly measured at the Gaisberg Tower in Salzburg (Austria) between 2000 and 2015. These upward flashes are combined with atmospheric reanalysis data stratified into five main meteorological groups: cloud physics, mass field, moisture field, surface exchange and wind field. We use classification methods based on tree-structured ensembles in form of conditional random forests. From these random forests we assess the meteorological influence and find the most important atmospheric drivers for one event or the other, respectively.

We have built and operated a crystalline/vapor xenon TPC, with the goal of improving searches for dark matter. The motivation for this instrument is the fact that beta decays from the radon decay chain to the ground state presently limit the state-of-the-art liquid/vapor xenon experiments. In contrast, a crystalline xenon target has the potential to tag and reject radon-chain backgrounds, due to the time and energy signature of their decays. The present article is the first demonstration of a crystalline/vapor xenon TPC with electroluminescence (gas gain) for the electron signal readout. It also shows that the scintillation yield in crystalline xenon appears to be identical to that in liquid xenon, in contrast to previous results.

The renormalization of the vacuum energy in quantum field theory (QFT) is usually plagued with theoretical conundrums related not only with the renormalization procedure itself, but also with the fact that the final result leads usually to very large (finite) contributions incompatible with the measured value of $\Lambda$ in cosmology. Herein, we compute the zero-point energy (ZPE) for a nonminimally coupled (massive) scalar field in FLRW spacetime using the off-shell adiabatic renormalization technique employed in previous work. The on-shell renormalized result first appears at sixth adiabatic order, so the calculation is rather cumbersome. The general off-shell result yields a smooth function $\rho_{\rm vac}(H)$ made out of powers of the Hubble rate and/or of its time derivatives involving different (even) adiabatic orders $\sim H^n$ ($n=2,4,6,...)$, i.e. it leads to the running vacuum model (RVM) structure. We have verified the same result from the effective action formalism and used it to find the $\beta$-function of the running quantum vacuum. No undesired contributions $\sim m^4$ from particle masses appear and hence no fine-tuning of the parameters is needed in $\rho_{\rm vac}(H)$. Furthermore, we find that the higher power $\sim H^6$ could naturally drive RVM-inflation in the early universe. Our calculation also elucidates in detail the equation of state of the quantum vacuum: it proves to be not exactly $-1$. The form of $\rho_{\rm vac}(H)$ at low energies is also characteristic of the RVM and consists of an additive constant (essentially the cosmological term) together with a small dynamical component $\sim \nu H^2$ ($|\nu|\ll1$). Finally, we predict a slow ($\sim\ln H$) running of Newton's gravitational coupling $G(H)$. The physical outcome of our semiclassical QFT calculation is revealing: today's cosmic vacuum and the gravitational strength should be both mildly dynamical.

Shou-shan Bao, Shaoqi Hou, Hong Zhang

9 pages, 7 figures. Comments are welcome

Wormholes bridging distant places of the universe are well-known solutions of general relativity. In particular, traversable wormholes which allow interstellar traveling are also popular in science fiction. However, no hint of their existence has been found yet. In this work, we propose using the gravitational wave (GW) scattering off spherical wormholes to search for their existence. We carefully calculate the reflected and transmitted waveforms with time-independent scattering theory. Our results quantitatively show the echo signatures in the two universes on both sides of the wormhole. In a certain wormhole mass range, the transmitted wave has a unique isolated chirp without an inspiral waveform, and the reflected wave has the anti-chirp behavior, i.e., the missing of the chirping signal. We also calculate the searching range of the current and projected GW telescopes. Our method can be adapted to efficiently calculate the templates to search for wormholes.

Sebastián Nájera, Aram Aguilar, Geovanny A. Rave-Franco, Celia Escamilla-Rivera, Roberto A. Sussman

17 pages, accepted for publication in IJGMMP

In this paper we explore the possibility to find exact solutions for Teleparallel Gravity (TG) of the type of spherically symmetric Lema\^\i tre-Tolman-Bondi (LTB) dust models. We apply to the LTB metric the formalism of Teleparallel Gravity in its extension to $f(T,B)$ models, which can be seen it as the analagous from the Schwarzschild solution in General Relativity. An exact LTB solution is obtained which is compatible with a specific $f(T,B)$ model whose observational constraints are cosmological viable in a standard spatially flat Robertson-Walker geometry.

We study a modification of the Higgs inflation scenario where we introduce an extra scalar $\phi$, with mass $m$, coupled to the Ricci scalar as $g\phi^2 R$, and mixed with the Higgs field $h$ by means of the Lagrangian term $\mu \phi h^2$. Both fields participate in the inflation process in a unitary theory that predicts values of the cosmological observables in agreement with the results from the Planck/BICEP/Keck collaborations. In addition, by means of a $\mathcal{CP}$ odd effective operator that couples $\phi$ to the Chern-Simons term of the hypercharge gauge group as $f_\phi^{-1}\phi \,Y_{\mu\nu}\tilde Y^{\mu\nu}$, maximally helical magnetic fields are produced during the last e-folds of inflation. We found a window in the coupling $ f_\phi$ where these fields survive all constraints until the electroweak phase transition, and source the baryon asymmetry of the Universe through the Standard Model chiral anomaly. From a phenomenological perspective, the model can solve the Standard Model instability problem at the scale $\mathcal Q_I\simeq 10^{11}$ GeV, provided that $\mu\lesssim m \lesssim \mathcal Q_I$, and for $m\lesssim \mathcal{O}$(few) TeV, the $\phi$-$h$ mixing becomes sizable while the theory turns natural. The latter thus predicts modifications of the trilinear and quartic couplings that could be explored at the HE-LHC, as well as at future colliders, and allows for direct $\phi$ production at the LHC followed by decay into $hh$. Present results from ATLAS and CMS already put (mild) bounds on the mass of the heavy scalar as $m\gtrsim 0.55$ TeV at 95% CL.

Nicolas Wijsen, Angels Aran, Camilla Scolini, David Lario, Alexandr Afanasiev, Rami Vainio, Blai Sanahuja, Jens Pomoell, Stefaan Poedts

14 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

We model the energetic storm particle (ESP) event of 14 July 2012 using the energetic particle acceleration and transport model named PARADISE, together with the solar wind and coronal mass ejection (CME) model named EUHFORIA. The simulation results illustrate both the capabilities and limitations of the utilised models. We show that the models capture some essential structural features of the ESP event; however, for some aspects the simulations and observations diverge. We describe and, to some extent, assess the sources of errors in the modelling chain of EUHFORIA and PARADISE and discuss how they may be mitigated in the future. The PARADISE model evolves energetic particle distributions in a background solar wind generated by the ideal MHD module of EUHFORIA. The CME generating the ESP event is simulated by using the spheromak model of EUHFORIA, which approximates the CME's flux rope as a linear force-free spheroidal magnetic field. In addition, a tool was developed to trace CME-driven shock waves in the EUHFORIA simulation domain. This tool is used in PARADISE to (i) inject 50 keV protons continuously at the CME-driven shock and (ii) include a foreshock and a sheath region, in which the energetic particle parallel mean free path, $\lambda_\parallel$, decreases towards the shock wave. The value of $\lambda_\parallel$ at the shock wave is estimated from in situ observations of the ESP event. For energies below 1 MeV, the simulation results agree well with both the upstream and downstream components of the ESP event observed by the Advanced Composition Explorer (ACE). This suggests that these low-energy protons are mainly the result of interplanetary particle acceleration. In the downstream region, the sharp drop in the energetic particle intensities is reproduced at the entry into the following magnetic cloud, illustrating the importance of a magnetised CME model.

Kalle Ala-Mattinen, Matti Heikinheimo, Kimmo Kainulainen, Kimmo Tuominen

14 pages, 11 figures, 1 table

We present a framework for solving a coupled set of momentum-dependent Boltzmann equations for the phase space distribution of cosmic relic particles, without resorting to approximations of assuming kinetic equilibrium or neglecting back scattering or elastic interactions. Our framework is amendable to precision numerical computations. To test it, we consider two benchmark models where the momentum-dependence of dark matter distribution function is potentially important: a real singlet scalar extension near the Higgs resonance, and a sterile neutrino dark matter model with a singlet scalar mediator. The singlet scalar example shows that even near sharp resonances the kinetic equilibrium holds well enough to justify the use of integrated momentum-independent Boltzmann equation in preliminary parameter scans. However, the integrated method may underestimate the relic density by up to 40% in extreme cases. In the sterile neutrino dark matter model we studied how the inclusion of previously ignored elastic interactions and processes with initial state sterile neutrinos could affect the non-thermal nature of their resulting distributions. Here the effects turned out to be negligible, proving the robustness of the earlier predictions.

Eduardo Guendelman, Ramon Herrera

19 pages and 2 figures. arXiv admin note: text overlap with arXiv:1408.5344

A two scalar field model that incorporates the Two Measures Theory (TMT) is introduced, in order to unify the early and present universe. We define the scale invariant couplings of the scalar fields to the different measures through exponential potentials. Spontaneous breaking of scale invariance takes place when integrating the fields that define the measures. When going to the Einstein frame we obtain: (i) An effective potential for the scalar fields with three flat regions which allows for a unified description of both early universe inflation (in the higher energy density flat region) as well as of present dark energy epoch which can be realized with a double phase, i.e., in two flat regions. (ii) In the slow roll inflation, only one field combination the ``dilaton", which transforms under scale transformations, has non trivial dynamics, the orthogonal one, which is scale invariant remains constant. The corresponding perturbations of the dilaton are calculated. (iii) For a reasonable choice of the parameters the present model perturbations conforms to the Planck Collaboration data. (iv) In the late universe we define scale invariant couplings of Dark Matter to the dilaton.(v) We calculate the evolution of the late universe under these conditions with the realization of two different possible realizations of $\Lambda$CDM type scenarios depending of the flat region in the late universe. These two phases could appear at different times in the history of the universe.(vi) From the Planck data, we find the constraints on the parameters during the inflationary epoch and these values are used to obtain constraints relevant to the present epoch.

We consider the propagation of gravitational waves in our late time Universe with the presence of structure. Gravitational waves emitted from distant sources have to traverse through regions that are far from smooth and homogeneous, before detection. We investigate the effect of inhomogeneities on the observables associated with the gravitational waves sources. In particular, we evaluate the impact of inhomogeneities on gravitational wave propagation employing the Buchert's framework of averaging. In context of a toy model within the above framework, it is first shown how the redshift versus distance relation, as well as the redshift drift get affected through the averaging process. We then study the variation of the redshift dependent part of the observed gravitational wave amplitude for different combination of our model parameters. We show that the variation of the gravitational wave amplitude with respect to redshift can deviate significantly in comparison with that in the Lambda-CDM model. Our result signifies the importance of local inhomogeneities on precision measurements of parameters of gravitational wave sources.

Tomoyo Namigata, C. J. Horowitz, R. Widmer-Schnidrig

six pages, 6 figures

Primordial black holes, with masses comparable to asteroids, are an attractive possibility for dark matter. In addition, other forms of dark matter could form compact dark objects (CDO). We search for small tidal accelerations from low mass black holes or CDOs orbiting near the Earth, and find none. Using about 10 years of data from the superconducting gravimeters in the Black Forest Observatory in South-Western Germany and at Djougou, Northern Benin in Western Africa we set an upper limit on the maximum mass of any dark object orbiting the Earth as a function of orbital radius. For semi-major axis less than two earth radii we exclude all black holes or CDOs with masses larger than 6.7x10^{13} kg. Lower mass primordial black holes may be strongly constrained by Hawking radiation. We conclude that near Earth black holes are extremely unlikely.

The fundamental nature and origin of dark energy are one of the premier mysteries of theoretical physics. In General Relativity Theory, the cosmological constant $\Lambda$ is the simplest explanation for dark energy. On the other hand, the cosmological constant $\Lambda$ suffers from a delicate issue so-called fine-tuning problem. This motivates one to modify the spacetime geometry of Einstein's GR. The $f(Q)$ gravity is a recently proposed modified theory of gravity in which the non-metricity scalar $Q$ drives the gravitational interaction. In this article, we consider a linear $f(Q)$ model, specifically $f(Q)=\alpha Q + \beta$, where $\alpha$ and $\beta$ are free parameters. Then we estimate the best fit values of model parameters that would be in agreement with the recent observational data sets. We use 57 points of the updated $H(z)$ data sets, 6 points of the BAO data sets, and 1048 points from the Pantheon supernovae samples. We apply the Bayesian analysis and likelihood function along with the Markov Chain Monte Carlo (MCMC) method. Further, we analyse the physical behaviour of cosmological parameters such as density, deceleration, and the EoS parameters corresponding to the constraint values of the model parameters. The evolution of deceleration parameter predicts a transition from decelerated to accelerated phases of the universe. Further, the evolution of equation of state parameter depicts quintessence type behaviour of the dark energy fluid part. We found that our $f(Q)$ cosmological model can effectively describe the late time cosmic acceleration without invoking any dark energy component in the matter part.

Raghuveer Garani, Michele Redi, Andrea Tesi

27 pages, 8 figures

We study the imprints of secluded dark sectors with a mass gap and self-interactions on the matter power spectrum. When Dark Matter (DM) is sufficiently light, in the ballpark of a few KeV, and self-interacting we find qualitative difference with respect to $\Lambda$CDM and also to free streaming DM. In order to emphasize the role of interactions for the evolution of the primordial perturbations we discuss various regimes: ranging from the ideal case of a tightly coupled perfect fluid to the free case of Warm Dark Matter, including the realistic case of small but non-vanishing self-interactions. We compute the matter power spectrum in all these regimes with the aid of Boltzmann solvers. Light dark sectors with self-interactions are efficiently constrained by Lyman-$\alpha$ data and we find that the presence of self-interactions relaxes the bound on the DM mass. As a concrete realization we study models with dark QCD-like sectors, where DM is made of light dark-pions.

Guey-Lin Lin, Thi Thuy Linh Nguyen, Martin Spinrath, Thi Dieu Hien Van, Tse-Chun Wang

34 pages, 15 figures

In this paper we discuss the prospects to take a picture of an extended neutrino source, i.e., resolving its angular neutrino luminosity distribution. This is challenging since neutrino directions cannot be directly measured but only estimated from the directions of charged particles they interact with in the detector material. This leads to an intrinsic blurring effect. We first discuss the problem in general terms and then apply our insights to solar neutrinos scattering elastically with electrons. Despite the aforementioned blurring we show how with high statistics and precision the original neutrino distributions could be reconstructed.

The holographic Ricci dark energy can be treated as a running vacuum due to its analogy in the energy density, which is a combination of $H$ and $\dot H$, the model can predict either eternal acceleration or eternal deceleration. In the earlier works, we have shown that the presence of additive constant in the energy density or by considering possible interaction between dark sectors through a phenomenological term, the model can predict a transition from a prior decelerated to a late accelerated epoch. This paper analyses the cosmic evolution of holographic Ricci dark energy as a running vacuum with a nonlinear interaction between dark sectors in a flat FLRW universe. We consider three possible nonlinear interaction forms which give analytically feasible solutions. We have constrained the model using the Type1a Supernova(Pantheon)+CMB(Planck 2018)+BAO(SDSS) data and evaluated the best-estimated values of all the model parameters. We have analyzed the evolution of the Hubble parameter and deceleration parameter of all three cases. We perform state finder analysis of the model, which implies the quintessence nature of the model and found that it is distinguishably different from the standard $\Lambda$CDM model. The dynamical system analysis of all three cases confirms the evolution of the universe from an unstable prior matter-dominated era to a stable end de Sitter phase.

We formulate the theory of first-order dissipative magnetohydrodynamics in an arbitrary hydrodynamic frame under the assumption of parity-invariance and discrete charge symmetry. We study the mode spectrum of Alfv\'en and magnetosonic waves as well as the spectrum of gapped excitations and derive constraints on the transport coefficients such that generic equilibrium states with constant magnetic fields are stable and causal under linearised perturbations. We solve these constraints for a specific equation of state and show that there exists a large family of hydrodynamic frames that renders the linear fluctuations stable and causal. This theory does not require introducing new dynamical degrees of freedom and therefore is a promising and simpler alternative to M\"{u}ller-Israel-Stewart-type theories. Together with a detailed analysis of transport, entropy production and Kubo formulae, the theory presented here is well suited for studying dissipative effects in various contexts ranging from heavy-ion collisions to astrophysics.

Shu-Yu Ho, Pyungwon Ko, Chih-Ting Lu

30 pages, 10 figures, 3 tables, and 2 appendices

In this paper, we construct for the first time a two-component strongly interacting massive particles (SIMP) dark matter (DM) model, where a complex scalar and a vector-like fermion play the role of the SIMP DM candidates. These two particles are stable due to an accidental $\mathbb{Z}^{}_4$ symmetry after the breaking of a $\text{U}(1)^{}_\textsf{D}$ gauge symmetry. By introducing one extra complex scalar as a mediator between the SIMP particles, this model can have $3 \to 2$ processes that determine the DM relic density. On the other hand, the SIMP DM particles can maintain kinetic equilibrium with the thermal bath until the DM freeze-out temperature via the $\text{U}(1)^{}_\textsf{D}$ gauge couplings. Most importantly, we find an unavoidable two-loop induced $2 \to 2$ process tightly connecting to the $3 \to 2$ process that would redistribute the SIMP DM number densities after the chemical freeze-out of DM. Moreover, this redistribution would significantly modify the predictions of the self-interacting cross section of DM compared with other SIMP models. It is crucial to include the two-loop induced $2 \to 2$ annihilations to obtain the correct DM phenomenology.

Joas Zapata, Rodrigo Negreiros, Thiago Sales, Prashanth Jaikumar

7 pages, 8 figures

In this article, we explore the cooling of isolated quark stars. These objects are structured of a homogeneous quark matter core and crusted by matter. To do this, we adopt two kinds of crust: (i) a crust made of purely nuclear matter following the Baym-Pethick-Sutherland (BPS) equation of state (EoS) and (ii) a crust made of nuggets of strange quark matter (strangelets). Both models have the same quark matter core described by the MIT bag model EoS. Our main purpose is to quantify the effects of a strangelet crust on the cooling and relaxation times of these strange stars. We also perform a thorough study of the thermal relaxation of quark stars, in which we have found that objects with a strangelet crust have a significantly different thermal relaxation time. Our study also includes the possible effects of color superconductivity in the quark core.

Tilt-to-length coupling is a technical term for the cross-coupling of angular or lateral jitter into an interferometric phase signal. It is an important noise source in precision interferometers and originates either from changes in the optical path lengths or from wavefront and clipping effects. Within this paper, we focus on geometric TTL coupling and categorize it into a number of different mechanisms for which we give analytic expressions. We then show that this geometric description is not always sufficient to predict the TTL coupling noise within an interferometer. We, therefore, discuss how understanding the geometric effects allows TTL noise reduction already by smart design choices. Additionally, they can be used to counteract the total measured TTL noise in a system. The presented content applies to a large variety of precision interferometers, including space gravitational wave detectors like LISA.

The holographic space-time approach to inflation provides a well defined and self-contained framework to study the early universe. Based on a quasi-local quantum gravity theory generalizing string theory beyond AdS backgrounds, it addresses fundamental questions like the arrow of time and the entropy of the initial cosmological state. It was recently argued that it also provides a naturel explanation for dark matter in the form of primordial black holes. The orignal idea however suffers from troubles that can be cured by considering instead Planck relics. This possibility is investigated and paths for observational confirmation are pointed out.