Abstract visibility
Abstract text size

Papers for Thursday, Mar 03 2022

Papers with local authors

Leander Thiele, Miles Cranmer, William Coulton, Shirley Ho, David N. Spergel

11 pages, 5 figures; condensed version accepted at the Fourth Workshop on Machine Learning and the Physical Sciences (NeurIPS 2021) as "Equivariant and Modular DeepSets with Applications in Cluster Cosmology"

1 vote
Paper 8 — arXiv:2203.00026
1 vote
Paper 8 — arXiv:2203.00026

Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continuous electron pressure field in galaxy clusters from gravity-only simulations. Modeling clusters is challenging for neural networks as most of the gas pressure is concentrated in a handful of voxels and even the largest hydrodynamical simulations contain only a few hundred clusters that can be used for training. Instead of conventional convolutional neural net (CNN) architectures, we choose to employ a rotationally equivariant DeepSets architecture to operate directly on the set of dark matter particles. We argue that set-based architectures provide distinct advantages over CNNs. For example, we can enforce exact rotational and permutation equivariance, incorporate existing knowledge on the tSZ field, and work with sparse fields as are standard in cosmology. We compose our architecture with separate, physically meaningful modules, making it amenable to interpretation. For example, we can separately study the influence of local and cluster-scale environment, determine that cluster triaxiality has negligible impact, and train a module that corrects for mis-centering. Our model improves by 70 % on analytic profiles fit to the same simulation data. We argue that the electron pressure field, viewed as a function of a gravity-only simulation, has inherent stochasticity, and model this property through a conditional-VAE extension to the network. This modification yields further improvement by 7 %, it is limited by our small training set however. (abridged)

Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

26 pages, 12 figures. Submitted to MNRAS. Comments welcome

1 vote
Paper 12 — arXiv:2203.00040
1 vote
Paper 12 — arXiv:2203.00040

Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback.

Scott G. Carlsten, Jenny E. Greene, Rachael L. Beaton, Shany Danieli, Johnny P. Greco

Submitted to ApJ. Comments welcome. 54 pages with the results starting on page 17 and conclusions on page 27

0 votes
Paper 2 — arXiv:2203.00014
0 votes
Paper 2 — arXiv:2203.00014

We present the final results of the Exploration of Local VolumE Satellites (ELVES) Survey, a survey of the dwarf satellites of a nearly volume-limited sample of Milky Way (MW)-like hosts in the Local Volume. Hosts are selected simply via a cut in luminosity ($M_{K_s}<-22.1$ mag) and distance ($D<12$ Mpc). We have cataloged the satellites of 25 of the 31 such hosts, with another five taken from the literature. All hosts are surveyed out to at least 150 projected kpc ($\sim R_\mathrm{vir}/2$) with the majority surveyed to 300 kpc ($\sim R_\mathrm{vir}$). Satellites are detected using a consistent semi-automated algorithm specialized for low surface brightness dwarfs. As shown through extensive tests with injected galaxies, the catalogs are complete to $M_V\sim-9$ mag and $\mu_{0,V}\sim26.5$ mag arcsec$^{-2}$. Candidates are confirmed to be real satellites through distance measurements including redshift, tip of the red giant branch, and surface brightness fluctuations. Across all 30 surveyed hosts, there are 338 confirmed satellites with a further 105 candidates awaiting distance measurement. For the vast majority of these, we provide consistent multi-band S\'{e}rsic photometry. We show that satellite abundance correlates with host mass, with the MW being quite typical amongst comparable systems, and that satellite quenched fraction rises steeply with decreasing satellite mass, mirroring the quenched fraction for the MW and M31. The ELVES survey represents a massive increase in the statistics of surveyed systems with known completeness, and the provided catalogs are a unique dataset to explore various aspects of small-scale structure and dwarf galaxy evolution.

Yingtian Chen, Oleg Y. Gnedin

20 pages, 17 figures, submitted to MNRAS

0 votes
Paper 42 — arXiv:2203.00599
0 votes
Paper 42 — arXiv:2203.00599

Globular clusters (GCs) are old massive star clusters that serve as "fossils" of galaxy formation. The advent of Gaia observatory has enabled kinematics studies of the Galactic GCs with great detail and revolutionized our understanding of the connections between GC properties and galaxy assembly. However, lack of kinematic measurements of extragalactic GCs limits the sample size of GC systems that we can fully study observationally. In this work, we present a model for GC formation and evolution, which includes positional and kinematic information of individual GCs by assigning them to particles in the Illustris TNG50-1 simulation based on the age and location. The model predictions are consistent with many observed properties of GC systems, including the distributions of position, systemic velocity, velocity dispersion, anisotropy parameter, orbital actions, and metallicities. We also analyze the properties of GCs from different origins. In the outer galaxy, ex-situ clusters are more dominant than the clusters formed in-situ. This leads to the GC metallicities decreasing outwards due to the increasing abundance of accreted, metal-poor clusters. We also find the ex-situ GCs to have greater velocity dispersions and orbital actions, in agreement with their accretion origin.

All other papers

Electromagnetic follow-up of gravitational wave detections is very resource intensive, taking up hours of limited observation time on dozens of telescopes. Creating more efficient schedules for follow-up will lead to a commensurate increase in counterpart location efficiency without using more telescope time. Widely used in operations research and telescope scheduling, mixed integer linear programming (MILP) is a strong candidate to produce these higher-efficiency schedules, as it can make use of powerful commercial solvers that quickly and efficiently find globally optimal solutions to provided problems. We detail a new target of opportunity scheduling algorithm designed with Zwicky Transient Facility in mind that uses mixed integer linear programming. We compare its performance to \texttt{gwemopt}, the tuned heuristic scheduler used by the Zwicky Transient Facility and other facilities during the third LIGO-Virgo gravitational wave observing run. This new algorithm uses variable-length observing blocks to enforce cadence requirements and ensure field observability, along with having a secondary optimization step to minimize slew time. Our algorithm shows improvement over \texttt{gwemopt} in successfully scheduling observations for a simulated binary neutron star merger data set consistent with LIGO-Virgo's third observing run. We demonstrate the ability to achieve better efficiencies with this straightforward algorithm, highlighting the potential of mixed integer target of opportunity schedulers for future multimessenger follow-up surveys.

Seoyoung Lyla Jung, Douglas Rennehan, Vida Saeedzadeh, Arif Babul, Michael Tremmel, Thomas R. Quinn, S. Ilani Loubser, E. O'Sullivan, Sukyoung K. Yi

25 pages, 17 figures, submitted to MNRAS and under revision

Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central galaxies in group-scale halos at $z=0$ to observed BGGs. The comparison encompasses the stellar mass-halo mass relation, various kinematic properties and scaling relations, morphologies, and the star formation rates. Generally, we find that Romulus reproduces the full spectrum of diversity in the properties of the BGGs very well, albeit with a tendency toward lower than the observed fraction of quenched BGGs. We find both early-type S0 and elliptical galaxies as well as late-type disk galaxies; we find Romulus galaxies that are fast-rotators as well as slow-rotators; and we observe galaxies transforming from late-type to early-type following strong dynamical interactions with satellites. We also carry out case studies of selected Romulus galaxies to explore the link between their properties, and the recent evolution of the stellar system as well as the surrounding intragroup/circumgalactic medium. In general, mergers/strong interactions quench star-forming activity and disrupt the stellar disk structure. Sometimes, however, such interactions can also trigger star-formation and galaxy rejuvenation. Black hole feedback can also lead to a decline of the star formation rate but by itself, it does not typically lead to complete quenching of the star formation activity in the BGGs.

Kyuseok Oh, Michael J. Koss, Yoshihiro Ueda, Daniel Stern, Claudio Ricci, Benny Trakhtenbrot, Meredith C. Powell, Jakob S. Den Brok, Isabella Lamperti, Richard Mushotzky, Federica Ricci, Rudolf E. Bär, Alejandra F. Rojas, Kohei Ichikawa, Rogerio Riffel, Ezequiel Treister, Fiona Harrison, C. Megan Urry, Franz E. Bauer, Kevin Schawinski

30 pages, 17 figures, 11 tables, accepted for publication in ApJS; part of BASS DR2 special issue

We present the second catalog and data release of optical spectral line measurements and AGN demographics of the BAT AGN Spectroscopic Survey, which focuses on the of Swift-BAT hard X-ray detected AGNs. We use spectra from dedicated campaigns and publicly available archives to investigate spectral properties of most of the AGNs listed in the 70-month Swift-BAT all-sky catalog; specifically, 743 of the 746 unbeamed and unlensed AGNs (99.6%). We find a good correspondence between the optical emission line widths and the hydrogen column density distributions using the X-ray spectra, with a clear dichotomy of AGN types for NH = 10^22 cm-2. Based on optical emission-line diagnostics, we show that 48%-75% of BAT AGNs are classified as Seyfert, depending on the choice of emission lines used in the diagnostics. The fraction of objects with upper limits on line emission varies from 6% to 20%. Roughly 4% of the BAT AGNs have lines too weak to be placed on the most commonly used diagnostic diagram, [O III]{\lambda}5007/H\b{eta} versus [N II]{\lambda}6584/H{\alpha}, despite the high signal-to-noise ratio (S/N) of their spectra. This value increases to 35% in the [O III]{\lambda}5007/[O II]{\lambda}3727 diagram, owing to difficulties in line detection. Compared to optically-selected narrow-line AGNs in the Sloan Digital Sky Survey, the BAT narrow-line AGNs have a higher rate of reddening/extinction, with H{\alpha}/H\b{eta} > 5 (~ 36%), indicating that hard X-ray selection more effectively detects obscured AGNs from the underlying AGN population. Finally, we present a subpopulation of AGNs that feature complex broad-lines (34%, 250/743) or double-peaked narrow emission lines (2%, 17/743).

S. J. Molyneux, R. Smit, D. Schaerer, R. J. Bouwens, L. Bradley, J. A. Hodge, S. N. Longmore, S. Schouws, P. van der Werf, A. Zitrin, S. Phillips

9 pages, 15 figures. Accepted for publication in MNRAS

We present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman break galaxy in the Epoch of Reionisation, A1703-zD1, through the detection of [C II] at a redshift of z = 6.8269 +/- 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic L$_{UV}$ ~ L$^*$$_{z=7}$ luminosity and a very blue Spitzer/IRAC [3.6]-[4.5] colour, implying high equivalent width line emission of [O III]+H$\beta$. [C II] is reliably detected at 6.1$\sigma$ co-spatial with the rest-frame UV counterpart, showing similar spatial extent. Correcting for the lensing magnification, the [C II] luminosity in A1703-zD1 is broadly consistent with the local L$_{[CII]}$ - SFR relation. We find a clear velocity gradient of 103 +/- 22 km/s across the source which possibly indicates rotation or an ongoing merger. We furthermore present spectral scans with no detected [C II] above 4.6$\sigma$ in two unlensed Lyman break galaxies in the EGS-CANDELS field at z ~ 6.6 - 6.9. This is the first time that NOEMA has been successfully used to observe [C II] in a 'normal' star-forming galaxy at z > 6, and our results demonstrate its capability to complement ALMA in confirming galaxies in the Epoch of Reionisation.

Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically-dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully-kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and lightcurves which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, $F_\nu\propto {\rm const}$; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by $\sim 0.8\, {\rm rad}$; (3) the peak flux exceeds the one at inferior conjunction by a factor of ten. We demonstrate that the double-peaked signature in the lightcurve is due to Doppler boosting in the post-shock flow.

A. Sollima, V. D'Orazi, R. Gratton, R. Carini, E. Carretta, A. Bragaglia, S. Lucatello

11 pages, 10 figures, acepted for publication by A&A

NGC1850 is the nearest Young Massive Cluster of the Local Group with a mass similar to those of Galactic globular clusters. Recent studies have revealed an extended morphology of its MSTO, which can be interpreted as a spread in either age or internal rotation. An accurate spectroscopic determination of its chemical properties is still missing. We analyse spectra obtained with MUSE in adaptive optics mode of 1167 stars in both components of this cluster (NGC1850A and NGC1850B). Thanks to this dataset, we measure an average metallicity of <[M/H]>=-0.31 +/- 0.01, a mean Ba abundance of <[Ba/Fe]>=+0.40 +/- 0.02 and a systemic radial velocity of <v_{LOS}>=251.1 +/- 0.3 km/s. The dispersion of the radial velocities suggests a dynamical mass of log(M/Ms)=4.84 +/- 0.1, while no significant systemic rotation is detected. We detect a significant bimodality in OI line strength among the TO stars of NGC1850A with ~66% of stars with [O/Fe]~-0.16 and the rest with no detectable line. The majority of O-weak stars populate preferentially the red side of the MSTO and show H lines in emission, suggesting that they are Be stars rotating close to their critical velocity. Among normal MSTO stars, red stars have on average broader line profiles than blue ones, suggesting a correlation between colour and rotational velocity. The mean metallicity of this cluster lies at the metal-rich side of the metallicity distribution of the LMC following its age-metallicity relation. The Ba and O abundances agree with those measured in the bar of this galaxy. The observed spread in OI line width among its MS stars can be interpreted as an effect of rotational mixing occurring in the envelopes of O-weak stars. The correlation between line broadening and colour suggests that the observed colour spread among turn-off stars can be due to a wide range in rotational velocity covered by these stars.

Patrick D. Sheehan, John J. Tobin, Leslie L. Looney, S. Thomas Megeath

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

We present Markov Chain Monte Carlo radiative transfer modeling of a joint ALMA 345 GHz and spectral energy distribution dataset for a sample of 97 protostellar disks from the VLA and ALMA Nascent Disk and Multiplicity Survey of Orion Protostars. From this modeling, we derive disk and envelope properties for each protostar, allowing us to examine the bulk properties of a population of young protostars. We find that disks are small, with a median dust radius of $ 29.4^{+ 4.1}_{- 2.7}$ au and a median dust mass of $ 5.8^{+ 4.6}_{- 2.7}$ M$_{\oplus}$. We find no statistically significant difference between most properties of Class 0, I, and Flat Spectrum sources with the exception of envelope dust mass and inclination. The distinction between inclination is an indication that the Class 0/I/Flat Spectrum system may be difficult to tie uniquely to the evolutionary state of protostars. When comparing with Class II disk dust masses in Taurus from similar radiative transfer modeling, we further find that the trend of disk dust mass decreasing from Class 0 to Class II disks is no longer present, though it remains unclear whether such a comparison is fair due to differences in star forming region and modeling techniques. Moreover, the disks we model are broadly gravitationally stable. Finally, we compare disk masses and radii with simulations of disk formation and find that magnetohydrodynamical effects may be important for reproducing the observed properties of disks.

Lucas Hellström, Abbas Askar, Alessandro A. Trani, Mirek Giersz, Ross P. Church, Johan Samsing

14 pages, 6 figures, 12 tables, submitted to MNRAS

In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission (using post-Newtonian (PN) corrections) and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH-star mergers by up to 75 per cent but it does not affect the number of BH-BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH-BH mergers during binary-single encounters. However, BH-star mergers are largely unaffected by the inclusion of these terms.

Akshay Suresh, James M. Cordes, Shami Chatterjee, Vishal Gajjar, Karen I. Perez, Andrew P. V. Siemion, Matt Lebofsky, David H. E. MacMahon, Cherry Ng

20 pages, 10 figures, submitted to ApJ

The Galactic Center (GC), with its high density of massive stars, is a promising target for radio transient searches. In particular, the discovery and timing of a pulsar orbiting the central supermassive black hole (SMBH) of our Galaxy will enable stringent strong-field tests of gravity and accurate measurements of SMBH properties. We performed multi-epoch 4-8 GHz observations of the inner $\approx$ 15 pc of our Galaxy using the Robert C. Byrd Green Bank Telescope in 2019 August-September. Our investigations constitute the most sensitive 4-8 GHz GC pulsar survey conducted to date, reaching down to a 6.1 GHz pseudo-luminosity threshold of $\approx$ 1 mJy kpc$^2$ for a pulse duty cycle of 2.5$\%$. We searched our data in the Fourier domain for periodic signals incorporating a constant or linearly changing line-of-sight pulsar acceleration. We report the successful detection of the GC magnetar PSR J1745$-$2900 in our data. Our pulsar searches yielded a non-detection of novel periodic astrophysical emissions above a 6$\sigma$ detection threshold in harmonic-summed power spectra. We reconcile our non-detection of GC pulsars with inadequate sensitivity to a likely GC pulsar population dominated by millisecond pulsars. Alternatively, close encounters with compact objects in the dense GC environment may scatter pulsars away from the GC. The dense central interstellar medium may also favorably produce magnetars over pulsars.

Sidney Lower, Desika Narayanan, Joel Leja, Benjamin D. Johnson, Charlie Conroy, Romeel Dave

Submitted to ApJ. 15 pages, 10 figures

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.

Spencer James Gibson, Aditya Narendra, Maria Giovanna Dainotti, Malgorzata Bogdan, Agniezska Pollo, Artem Poliszczuk, Enrico Rinaldi, Ioannis Liodakis

21 Pages, 11 Figures, 3 Tables

Redshift measurement of active galactic nuclei (AGNs) remains a time-consuming and challenging task, as it requires follow up spectroscopic observations and detailed analysis. Hence, there exists an urgent requirement for alternative redshift estimation techniques. The use of machine learning (ML) for this purpose has been growing over the last few years, primarily due to the availability of large-scale galactic surveys. However, due to observational errors, a significant fraction of these data sets often have missing entries, rendering that fraction unusable for ML regression applications. In this study, we demonstrate the performance of an imputation technique called Multivariate Imputation by Chained Equations (MICE), which rectifies the issue of missing data entries by imputing them using the available information in the catalog. We use the Fermi-LAT Fourth Data Release Catalog (4LAC) and impute 24% of the catalog. Subsequently, we follow the methodology described in Dainotti et al. (2021) and create an ML model for estimating the redshift of 4LAC AGNs. We present results which highlight positive impact of MICE imputation technique on the machine learning models performance and obtained redshift estimation accuracy.

A. G. de Wijn, R. Casini, A. Carlile, A. R. Lecinski, S. Sewell, P. Zmarzly, A. D. Eigenbrot, C. Beck, F. Wöger, M. Knölker

The Daniel K. Inouye Solar Telescope (DKIST) Visible Spectro-Polarimeter (ViSP) is a traditional slit-scanning spectrograph, with the ability to observe solar regions up to a $120\times78~\mathrm{arcsec}^2$ area. The design implements dual-beam polarimetry, a polychromatic polarization modulator, a high-dispersion echelle grating, and three spectral channels that can be automatically positioned. A defining feature of the instrument is its capability to tune anywhere within the 380-900~nm range of the solar spectrum, allowing for a virtually infinite number of combinations of three wavelengths to be observed simultaneously. This enables the ViSP user to pursue well-established spectro-polarimetric studies of the magnetic structure and plasma dynamics of the solar atmosphere, as well as completely novel investigations of the solar spectrum. Within the suite of first-generation instruments at the DKIST, ViSP is the only wavelength-versatile spectro-polarimeter available to the scientific community. It was specifically designed to be a discovery instrument, for the exploration of new spectroscopic and polarimetric diagnostics, and to test improved models of polarized line formation, through high spatial-, spectral-, and temporal-resolution observations of the Sun's polarized spectrum. In this instrument article, we describe the science requirements and design drivers of ViSP, and we present preliminary science data collected during the commissioning of the instrument.

C.L. Carilli (NRAO), F. Walter (MPIA), R. Decarli (INAF), M. Aravena (UDP), Dominik A. Riechers (Cologne), J. Gonzalez-Lopez (UDP), Yali Shao (MPIfR), L. Boogaard (Leiden), R. Bouwens (Leiden), M. Neeleman (MPIA)

5 pages. 4 figures. to appear in Radio Science

We review the evolution of the cosmic average molecular gas density to large look-back times, using observations of rotational transitions of CO. Molecular gas is the fuel for star formation in galaxies. Deep searches for CO emission from distant galaxies have delineated the density of molecular gas back to $z \sim 5$, or within 1~Gyr of the Big Bang. The results show a rise and fall in the gas density that parallels, and likely drives, the rise and fall of the cosmic star formation rate density. We present the potential for the next generation Very Large Array to image the distribution and dynamics of the molecular gas in early galaxies, and to make a precise measurement of the dense gas history of the Universe.

Long-Fei Chen, Di Li, Donghui Quan, Xia Zhang, Qiang Chang, Xiaohu Li, Lin Xiao

25 pages, 11 figures, accepted for publication in ApJ

We investigated the chemical evolutions of gas phase and grain surface species across the Taurus molecular cloud-1 (TMC-1) filament from translucent phase to dense phase. By comparing observations with modeling results from an up-to-date chemical network, we examined the conversion processes for the carbon-, oxygen-, nitrogen- and sulfur-bearing species, i.e.from their initial atomic form to their main molecular reservoir form both in the gas phase and on the grain surface. The conversion processes were found to depend on the species and A$_V$. The effect of initial carbon to oxygen elemental abundances ratio (C/O) by varying O on the chemistry was explored, and an initial carbon elemental abundance of 2.5 $\times$ 10$^{-4}$ and a C/O ratio of 0.5 could best reproduce the abundances of most observed molecules at TMC-1 CP, where more than 90 molecules have been identified. Based on the TMC-1 condition, we predicted a varied grain ice composition during the evolutions of molecular clouds, with H$_2$O ice as the dominant ice composition at A$_V$ $>$ 4 mag, CO$_2$ ice as the dominant ice composition at A$_V$ $<$ 4 mag, while CO ice severely decreased at A$_V$ around 4--5 mag.

Abdolhosein Hashemizadeh, Simon P. Driver, Luke J. M. Davies, Aaron S. G. Robotham, Sabine Bellstedt, Rogier A. Windhorst, Matt Jarvis, Benne W. Holwerda, Malgorzata Siudek, Caroline Foster, Steven Phillipps, Jessica E. Thorne, Christian Wolf

The majority of the stellar mass in the Universe today resides in galaxies with two primary structural components (bulge and disk). In this work, we use the largest contiguous HST imaging region (COSMOS), the Deep Extragalactic VIsible Legacy Survey (DEVILS), and the galaxy structural fitting code {\sc ProFit}, to deconstruct $\sim 35,000$ galaxies into their sub-structures. We use this sample to determine the stellar mass density (SMD) sub-divided by structural components and its evolution since $z = 1$. We find that the majority of stellar mass at all epochs lies in disk-like structures. The SMD in the disk population increases in $z = 1-0.35$ and stabilizes/decreases to $z = 0$ (contributing $\sim 60\%$ to the total SMD at $z = 1$ and declining to $\sim 30\%$ at $z = 0$). This decline is countered by a rapid rise of the SMD in pseudo-bulge population and a consistent growth of spheroidal structures (classical bulges and ellipticals) with significant stellar mass growth in $z = 1-0.35$ and a somewhat flattened trend to $z = 0$. While the physical mechanisms for this are not obvious, the results are consistent with a transition from a Universe dominated by disk growth, to a Universe in which pseudo-bulges are emerging and spheroids are growing. Further study is required before this can be apportioned to internal secular processes, gas accretion, minor and major mergers. However, it is clear that since $z = 0.35$ the processes mentioned above now dominate over quiescent disk growth as the cosmic star-formation history declines.

M.B. Naik, D.K. Ojha, S.K. Ghosh, P. Manoj, J.P. Ninan, S. Ghosh, S.L.A. D'Costa, S.S. Poojary, S.B. Bhagat, P.R. Sandimani, H. Shah, R.B. Jadhav, S.M. Gharat, G.S. Meshram, B.G. Bagade

19 pages, 28 figures, 2 tables, Accepted for publication in The Journal of Astronomical Instrumentation (JAI)

The Infrared Astronomy Group (Department of Astronomy and Astrophysics) at Tata Institute of Fundamental Research (TIFR) is presently developing controllers for the Teledyne HxRG Focal Plane Arrays (FPAs) to be used on board the Infrared Spectroscopic Imaging Survey (IRSIS) satellite payload. In this manuscript we discuss the results of our tests with different FPA controllers like the Astronomical Research Cameras (ARC) controller, Teledyne's SIDECAR ASIC as well as our new in-house designed Array controller. As part of the development phase of the IRSIS instrument, which is an optical fibre based Integral Field Unit (IFU) Near-Infrared (NIR) Spectrometer, a laboratory model with limited NIR bandwidth was built which consisted of various subsystems like a Ritchey-Chretien (RC) 30 cm telescope, optical fibre IFU, spectrometer optics, and the Teledyne H2RG detector module. We discuss the various developments during the building and testing of the IRSIS laboratory model and the technical aspects of the prototype in-house H2RG controller.

Katherine L. Hill, Colin Littlefield, Peter Garnavich, Simone Scaringi, Paula Szkody, Paul A. Mason, Mark R. Kennedy, Aarran W. Shaw, Ava E. Covington

12 pages, 4 figures, submitted to the Astronomical Journal and reviewed favorably

We present the Transiting Exoplanet Surveying Satellite (TESS) light curve of the intermediate polar YY Draconis (YY Dra, also known as DO Dra). The power spectrum indicates that while there is stream-fed accretion for most of the observational period, there is a day-long, flat-bottomed low state at the beginning of 2020 during which the only periodic signal is ellipsoidal variation and there is no appreciable flickering. We interpret this low state to be a complete cessation of accretion, a phenomenon that has been observed only once before in an intermediate polar. Simultaneous ground-based observations of this faint state establish that when accretion is negligible, YY Dra fades to $g=17.37\pm0.12$, which we infer to be the magnitude of the combined photospheric contributions of the white dwarf and its red dwarf companion. Using survey photometry, we identify additional low states in 2018-2019 during which YY Dra repeatedly fades to -- but never below -- this threshold. This implies relatively frequent cessations in accretion. Spectroscopic observations during future episodes of negligible accretion can be used to directly measure the field strength of the white dwarf by Zeeman splitting. Separately, we search newly available catalogs of variable stars in an attempt to resolve the long-standing dispute over the proper identifier of this system.

pmbonney@uark.edu, (2) Department of Physics, University of Arkansas)

17 pages, 6 figures, 3 tables. Submitted to The Planetary Science Journal (2/28/22)

A main determinant of the habitability of exoplanets is the presence of stable liquid surface water. In an era of abundant possible targets, the potential to find a habitable world remains a driving force in prioritization. We present here a data-forward method to investigate the likelihood of a stable hydrosphere on the timescales of the formation of life, 1 Gyr, and beyond. As our primary application, we use this method to examine the potential hydrospheres of TESS Objects of Interest 700 d, 256 b, and 203 b. We first present our selection criteria, which is based on an implementation of the Earth Similarity Index, as well as the results of an initial investigation into the desiccation of the targets, which reveals that TOI 203 b is almost certainly desiccated based on TESS observations. We then describe the characterization of the remaining targets and their host stars from 2MASS, Gaia, and TESS data and the derivation of sampled probability distributions for their parameters. Following this, we describe our process of simulating the desiccation of the targets hydrospheres using the Virtual Planet Simulator, VPLanet, with inputs directly linked to the previously derived probability distributions. We find that 50.86 percent of the likely cases for TOI 700 d are desiccated and no modeled cases for TOI 256 b are without water. In addition, we calculate the remaining water inventory for the targets, the percentage of cases that are continuing to lose water, and the rate at which these cases are losing water.

Hong Wu, Zhongzu Wu, Yu. Sotnikova, Yongjun Chen, Bo Zhang, T. Mufakharov, Zhiqiang Shen, Xi Chen, A. Mikhailov, M. Mingaliev, Xianming L. Han, Prabhakar Misra

11 pages,9 figures, Accepted by A&A

Based on the two epochs EVN archive data from OH line observations of IIZw 096, we confirm that the high-resolution OH emission in this source mainly comes from two spots (OH1 and OH2) of comp D1 of this merging system. We found no significant variations in the OH line emission. The OH 1665 MHz line emission is detected at about 6 $\sigma$ level in the OH1 region by combining two epoch EVN observations. We found that the comp D1 shows the brightest CO, HCO+ line emission, as well as multi-band radio continuum emission. The environment around D1 shows no clear velocity structure associated with circular motions, making it different from most other OHMs in the literature, which might have been caused by an effect during the merger stage. Meanwhile, we found that the CO emission shows three velocity structures around D1, including the central broad FWHM region, the double peak region where the CO line profile shows two separated peaks, and the region of the high-velocity clouds where the CO line peaks at a high velocity ($\sim$ 11000 \kms). \HI in absorption also show high-velocity clouds around the D1 region, which might be due to inflows caused by the merging of two or more galaxy components. Based on the high-resolution K-band VLA and L-band VLBA observations of the radio continuum emission, we derived the brightness temperature in the range $10^{5}$ K to $10^{6}$ K, which is consistent with other starburst dominant OHM sources in the literature. The multi-band VLA observations show that the radio continuum emission of comp D might also have contributions from free-free emission, besides synchrotron emission. As a concenquence, these results support a starburst origin for the OHMs, without the presence of an AGN.

Haomin Sun, Hui Deng, Feng Wang, Ying Mei, Tingting Xu, Oleg Smirnov, Linhua Deng, Shoulin Wei

10 Pages, 9 figures, Accepted by MNRAS

The rapid development of new generation radio interferometers such as the Square Kilometer Array (SKA) has opened up unprecedented opportunities for astronomical research. However, anthropogenic Radio Frequency Interference (RFI) from communication technologies and other human activities severely affects the fidelity of observational data. It also significantly reduces the sensitivity of the telescopes. We proposed a robust Convolutional Neural Network (CNN) model to identify RFI based on machine learning methods. We overlaid RFI on the simulation data of SKA1-LOW to construct three visibility function datasets. One dataset was used for modeling, and the other two were used for validating the model's usability. The experimental results show that the Area Under the Curve (AUC) reaches 0.93, with satisfactory accuracy and precision. We then further investigated the effectiveness of the model by identifying the RFI in the actual observational data from LOFAR and MeerKAT. The results show that the model performs well. The overall effectiveness is comparable to AOFlagger software and provides an improvement over existing methods in some instances.

Ö. Baştürk, E. M. Esmer, S. Yalçınkaya, Ş. Torun, L. Mancini, F. Helweh, E. Karamanlı, J. Southworth, S. Aliş, A. Wünsche, F. Tezcan, Y. Aladağ, N. Aksaker, E. Tunç, F. Davoudi, S. Fişek, M. Bretton, D. F. Evans, C. Yeşilyaprak, M. Yılmaz, C. T. Tezcan, K. Yelkenci

16 pages, 5 figures, accepted for publication by the MNRAS

We study the transit timings of 10 exoplanets in order to investigate potential Transit Timing Variations (TTVs) in them. We model their available ground-based light curves, some presented here and others taken from the literature, and homogeneously measure the mid-transit times. We statistically compare our results with published values and find that the measurement errors agree. However, in terms of recovering the possible frequencies, homogeneous sets can be found to be more useful, of which no statistically relevant example has been found for the planets in our study. We corrected the ephemeris information of all ten planets we studied and provide these most precise light elements as references for future transit observations with space-borne and ground-based instruments. We found no evidence for secular or periodic changes in the orbital periods of the planets in our sample, including the ultra-short period WASP-103 b, whose orbit is expected to decay on an observable timescale. Therefore, we derive the lower limits for the reduced tidal quality factors (Q$^{\prime}_{\star}$) for the host stars based on best fitting quadratic functions to their timing data. We also present a global model of all available data for WASP-74 b, which has a Gaia parallax-based distance value ~25% larger than the published value.

The magnetic interaction between a classical T Tauri star and its surrounding accretion disk is thought to influence its rotational evolution. We use 2.5D magnetohydrodynamic, axisymmetric simulations of star-disk interaction, computed via the PLUTO code, to calculate the net torque acting on these stars. We divide the net torque into three contributions: accretion (spin-up), stellar winds (spin-down), and magnetospheric ejections (MEs) (spin-up or down). In Paper I, we explored interaction regimes in which the stellar magnetosphere truncates the inner disk at a location spinning faster than the star, resulting in a strong net spin-up contribution from accretion and MEs ("steady accretion" regime). In this paper, we investigate interaction regimes in which the truncation radius gets closer to and even exceeds corotation, where it is possible for the disk material to gain angular momentum and be periodically ejected by the centrifugal barrier ("propeller" regime). This reduces the accretion torque, can change the sign of the ME torque, and can result in a net stellar spin-down configuration. These results suggest it is possible to have a net spin-down stellar torque even for truncation radii within the corotation radius ($R_\text{t} \gtrsim 0.7 R_\text{co}$). We fit semi-analytic functions for the truncation radius, and the torque associated with star-disk interaction (i.e., the sum of accretion and ME torques) and stellar wind, allowing for the prediction of the net stellar torque for a parameter regime covering both net spin-up and spin-down configurations, as well as the possibility of investigating rotational evolution via 1D stellar evolution codes.

Allard Jan van Marle, Artem Bohdan, Paul J. Morris, Martin Pohl, Alexandre Marcowith

Accepted for publication in ApJ

The current paradigm of cosmic ray (CR) origin states that the most part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. Then, with this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate the upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely triggered at high Mach number shocks.

S. Zieba, M. Zilinskas, L. Kreidberg, T. G. Nguyen, Y. Miguel, N. B. Cowan, R. Pierrehumbert, L. Carone, L. Dang, M. Hammond, T. Louden, R. Lupu, L. Malavolta, K. B. Stevenson

31 pages, 27 figures, 12 tables, accepted for publication in A&A

K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of the planet. We measure an infrared eclipse depth of 143 +/- 39 ppm and a peak to trough amplitude variation of 121 +/- 43 ppm. The best fit model to the Spitzer data shows no significant thermal hotspot offset, in contrast to the previously observed offset for the well-studied USP planet 55 Cnc e. We also jointly analyze the new Spitzer observations with the photometry collected by Kepler during two separate K2 campaigns. We model the planetary emission with a range of toy models that include a reflective and a thermal contribution. With a two-temperature model, we measure a dayside temperature of 2049 +/- 361 K and a night-side temperature that is consistent with zero (Tp,n < 1712 K at 2 sigma). Models with a steep dayside temperature gradient provide a better fit to the data than a uniform dayside temperature (DeltaBIC = 22.2). We also find evidence for a non-zero geometric albedo Ag = 0.28 +/- 0.07. We also compare the data to a physically motivated, pseudo-2D rock vapor model and a 1D turbulent boundary layer model. Both models fit the data well. Notably, we find that the optical eclipse depth can be explained by thermal emission from a hot inversion layer, rather than reflected light. A thermal inversion may also be responsible for the deep optical eclipse observed for another USP, Kepler-10 b. Finally, we significantly improve the ephemerides for K2-141 b and c, which will facilitate further follow-up observations of this interesting system with state-of-the-art observatories like JWST.

S.V.Jeffers (MPS, Goettingen), J.R.Barnes, P.Scheofer, A.Quirrenbach, M.Zechmeister, P.J.Amado, J.A.Caballero, M.Fernandez, E.Rodriguez, I.Ribas, A.Reiners, C.Cardona Guillen, C.Cifuentes, S.Czesla, A.P.Hatzes, M.Kurster, D.Montes, J.C.Morales, S.Pedraz, S.Sadegi

accepted by A&A

Current exoplanet surveys using the radial velocity (RV) technique are targeting M dwarfs because any habitable zone terrestrial-mass planets will induce a high RV and orbit on shorter periods than for more massive stars. One of the main caveats is that M dwarfs show a wide range of activity levels from inactive to very active, which can induce an asymmetry in the line profiles and, consequently, a spurious RV measurement. We aim to benchmark the impact of stellar activity on high-precision RV measurements using regular-cadence CARMENES visible and near-infrared observations of the active M3.5 dwarf EV Lac. We used the newly developed technique of low-resolution Doppler imaging to determine the centre-of-light, or spot-induced RV component, for eight observational epochs. We confirm a statistically significant and strong correlation between the independently measured centre-of-light and the chromatic index, which is a measure of the amplitude variation with wavelength of the RVs. We also find circular closed-loop relations of several activity indices with RV for a subset of data that covers only several rotation periods. We also investigate the implications of large phase gaps in the periodograms of activity indicators. Finally, by removing the spot-induced RV component we improve the planet-mass sensitivity by a factor of at least three. We conclude that for active M stars, a regular-cadence observing strategy is the most efficient way to identify and eliminate sources of correlated noise.

Marino Yamamoto, Jun Kataoka, Yoshiaki Sofue

11 pages, 5 figs, accepted in MNRAS

We investigated the detailed plasma condition of the North Polar Spur (NPS)/Loop I using archival $Suzaku$ data. In previous research collisional ionization equilibrium (CIE) have been assumed for X-ray plasma state, but we also assume non-equilibrium ionization (NEI) to check the plasma condition in more detail. We found that most of the plasma in the NPS/Loop I favors the state of NEI, and has the density-weighted ionization timescale of $n_e t\sim10^{11-12}$ s cm$^{-3}$ and the electron number density $n_e\sim$ a few $\times$ 10$^{-3}$ cm$^{-3}$. The plasma shock age, $t$, or the time elapsed after the shock front passed through the plasma, is estimated to be on the order of a few $\rm{Myr}$ for the NPS/Loop I, which puts a strict lower limit to the age of the whole NPS/Loop I structure. We found that NEI results in significantly higher temperature and lower emission measure than those currently derived under CIE assumption. The electron temperature under NEI is estimated to be as high as 0.5 keV toward the brightest X-ray NPS ridge at $\Delta\theta=-20^\circ$, which decreases to 0.3 keV at $-10^\circ$, and again increases to $\sim 0.6$ keV towards the outer edge of Loop I at $\Delta\theta\sim0^\circ$, about twice the currently estimated temperatures. Here, $\Delta \theta$ is the angular distance from the outer edge of Loop I. We discuss the implication of introducing NEI for the research in plasma states in astrophysical phenomena.

Felix Dannert, Maurice Ottiger, Sascha P. Quanz, Romain Laugier, Emile Fontanet, Adrian Gheorghe, Olivier Absil, Colin Dandumont, Denis Defrère, Carlos Gascón, Adrian M. Glauser, Jens Kammerer, Tim Lichtenberg, Hendrik Linz, Jerôme Loicq, the LIFE collaboration

Accepted for publication in A&A; 15 pages (main text incl. 13 figures and 1 table) + appendix; comments are welcome

The Large Interferometer For Exoplanets (LIFE) initiative is developing the science and a technology roadmap for an ambitious space mission featuring a space-based mid-infrared (MIR) nulling interferometer in order to detect the thermal emission of hundreds of exoplanets and characterize their atmospheres. In order to quantify the science potential of such a mission, in particular in the context of technical trade-offs, an instrument simulator is required. In addition, signal extraction algorithms are needed to verify that exoplanet properties (e.g., angular separation, spectral flux) contained in simulated exoplanet datasets can be accurately retrieved. We present LIFEsim, a software tool developed for simulating observations of exoplanetary systems with an MIR space-based nulling interferometer. It includes astrophysical noise sources (i.e., stellar leakage and thermal emission from local zodiacal and exo-zodiacal dust) and offers the flexibility to include instrumental noise terms in the future. LIFEsim provides an accessible way for predicting the expected SNR of future observations as a function of various key instrument and target parameters. The SNRs of the extracted spectra are photon-noise dominated, as expected from our current simulations. From single epoch observations in our mock survey of small ($R < 1.5 R_\mathrm{Earth}$) planets orbiting within the habitable zones of their stars, we find that typical uncertainties in the estimated effective temperature of the exoplanets are $\lesssim$10%, for the exoplanet radius $\lesssim$20%, and for the separation from the host star $\lesssim$2%. SNR values obtained in the signal extraction process deviate less than 10% from purely photon-counting statistics based SNRs. (abridged)

Alexandra Kozyreva, Hans-Thomas Janka, Daniel Kresse, Stefan Taubenberger

10 pages, 6 figures, submitted to MNRAS

SN 2020cxd is a representative of the family of low-energy, underluminous Type IIP supernovae (SNe), whose observations and analysis were recently reported by Yang et al. (2021). Here we re-evaluate the observational data for the diagnostic SN properties by employing the hydrodynamic explosion model of a 9 MSun red supergiant progenitor with an iron core and a pre-collapse mass of 8.75 Msun. The explosion of the star was obtained by the neutrino-driven mechanism in a fully self-consistent simulation in three dimensions (3D). Multi-band light curves and photospheric velocities for the plateau phase are computed with the one-dimensional radiation-hydrodynamics code STELLA, applied to the spherically averaged 3D explosion model as well as spherisized radial profiles in different directions of the 3D model. We find that the overall evolution of the bolometric light curve, duration of the plateau phase, and basic properties of the multi-band emission can be well reproduced by our SN model with its explosion energy of only 0.7x10^50 erg and an ejecta mass of 7.4 Msun. These values are considerably lower than the previously reported numbers, but they are compatible with those needed to explain the fundamental observational properties of the prototype low-luminosity SN 2005cs. Because of the good compatibility of our photospheric velocities with line velocities determined for SN 2005cs, we conclude that the line velocities of SN 2020cxd are probably overestimated by up to a factor of about 3. The evolution of the line velocities of SN 2005cs compared to photospheric velocities in different explosion directions might point to intrinsic asymmetries in the SN ejecta.

Joanna K. Barstow, Quentin Changeat, Katy L.Chubb, Patricio E. Cubillos, Billy Edwards, Ryan J. MacDonald, Michiel Min, Ingo P. Waldmann

28 pages, 14 figures. Accepted in Experimental Astronomy (2022)

The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a critical preparatory step for the mission. Here, we present the results of a retrieval challenge. We use five different exoplanet retrieval codes to analyse the same synthetic datasets, and test a) the ability of each to recover the correct input solution and b) the consistency of the results. We find that generally there is very good agreement between the five codes, and in the majority of cases the correct solutions are recovered. This demonstrates the reproducibility of retrievals for transit spectra of exoplanets, even when codes are not previously benchmarked against each other.

D. Maurin, E. Ferronato Bueno, Y. Génolini, L. Derome, M. Vecchi

20 pages, 16 figures, 2 tables (submitted to A&A)

Li, Be, and B high-precision data from AMS-02 provide the best constraints on Galactic cosmic-ray transport parameters. We re-evaluate the impact of Fe fragmentation on the Li, Be, and B modelling. We discuss the consequences on the transport parameter determination and reassess whether a primary source of Li is needed to match AMS-02 data. We renormalise several cross-section parametrisations to existing data for the most important reactions producing Li, Be, and B. We use the \usine{} code with these new cross-section sets to re-analyse Li/C, Be/C, and B/C AMS-02 data. We build three equally plausible cross-section sets. Compared to the use of the initial cross-section sets, they lead to an average enhanced production of Li ($\sim20-50\%$) and Be ($\sim5-15\%$), while leaving the B flux mostly unchanged. In particular, Fe fragmentation is found to contribute to up to 10\% of the Li and Be fluxes. Used in the combined analysis of AMS-02 Li/C, Be/C, and B/C data, the fit is significantly improved, with an enhanced diffusion coefficient ($\sim 20\%)$. The three updated cross-section sets are found to either slightly undershoot or overshoot the Li/C and B/C ratios: this strongly disfavours evidences for a primary source of Li in cosmic rays. We stress that isotopic cosmic-ray ratios of Li (and to a lesser extent Be), soon to be released by AMS-02, are also impacted by the use of these updated sets. Almost no nuclear data exist for the production of Li and B isotopes from Ne, Mg, Si, and Fe, whereas these reactions are estimated to account for $\sim 20\%$ of the total production. Some new nuclear measurements are clearly desired to better exploit the high-precision AMS-02 cosmic-ray data.

We study the evolution of random hypermagnetic fields (HMFs) in the symmetric phase of the early universe before the electroweak phase transition. The behavior of HMFs is driven by the analog of the chiral magnetic effect accounting for the asymmetries of leptons and Higgs bosons. These asymmetries are also dynamical variables of the model and evolve together with HMFs. Moreover, we account for the contribution of the hyper-MHD turbulence in the effective diffusion coefficient and the $\alpha$-dynamo parameter. The realistic spectrum of seed HMFs consists of two branches: Batchelor and Kolmogorov ones. The impact of HMFs on the production of relic gravitational waves (GWs) and the baryon asymmetry of the universe (BAU), as well as flavor oscillations of supernova neutrinos in the GWs generated are considered. We establish the constraint on the strength of the seed HMF comparing the spectral density of produced GWs with the observations of the LIGO-Virgo-KAGRA collaborations. The stronger upper bound on the seed HMF is obtained from the condition of not exceeding the observed value of BAU.

Izaskun Jimenez-Serra (1), Jesus Martin-Pintado (1), Aran Insausti (2, 3), Elena R. Alonso (2, 3), Emilio J. Cocinero (2, 3), Tyler L. Bourke (4) ((1) Centro de Astrobiologia (CSIC/INTA), Torrejon de Ardoz, Spain, (2) Departamento de Quimica Fisica, Facultad de Ciencia y Tecnologia, Universidad del Pais Vasco (UPV-EHU), Bilbao, Spain, (3) Biofisika Institute (CSIC, UPV/EHU), Leioa, Spain, (4) SKA Observatory, Jodrell Bank, Macclesfield, United Kingdom)

12 pages, 1 figure, 1 table, accepted for Frontiers in Astronomy and Space Sciences, section Astrochemistry

One of the theories for the origin of life proposes that a significant fraction of prebiotic material could have arrived to Earth from outer space between 4.1 and 3.8 billion years ago. This suggests that those prebiotic compounds could have originated in interstellar space, to be later on incorporated to small Solar-system bodies and planetesimals. The recent discovery of prebiotic molecules such as hydroxylamine and ethanolamine in the interstellar medium, strongly supports this hypothesis. However, some species such as sugars, key for the synthesis of ribonucleotides and for metabolic processes, remain to be discovered in space. The unmatched sensitivity of the Square Kilometer Array (SKA) at centimeter wavelengths will be able to detect even more complex and heavier prebiotic molecules than existing instrumentation. In this contribution, we illustrate the potential of the SKA to detect simple sugars with three and four carbon atoms, using a moderate investment of observing time.

M.C. Bezuidenhout, E. Barr, M. Caleb, L.N. Driessen, F. Jankowski, M. Kramer, M. Malenta, V. Morello, K. Rajwade, S. Sanidas, B.W. Stappers, M. Surnis

16 pages, 14 figures

MeerTRAP is a real-time untargeted search project using the MeerKAT telescope to find single pulses from fast radio transients and pulsars. It is performed commensally with the MeerKAT large survey projects (LSPs), using data from up to 64 of MeerKAT's 13.96~m dishes to form hundreds of coherent beams on sky, each of which is processed in real time to search for millisecond-duration pulses. We present the first twelve Galactic sources discovered by MeerTRAP, with DMs in the range of 33--381~pc~cm$^{-3}$. One source may be Galactic or extragalactic depending on the Galactic electron density model assumed. Follow-up observations performed with the MeerKAT, Lovell, and Parkes radio telescopes have detected repeat pulses from seven of the twelve sources. Pulse periods have been determined for four sources. Another four sources could be localised to the arcsecond-level using a novel implementation of the tied-array beam localisation method.

We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis with the goal of understanding how each of the updates affects star cluster populations. Two key parameters are the momentum boost of supernova feedback $f_{\mathrm{boost}}$ and star formation efficiency per freefall time $\epsilon_{\mathrm{ff}}$. We find that $f_{\mathrm{boost}}$ has a strong influence on the galactic star formation rate, with higher values leading to less star formation. The efficiency $\epsilon_{\mathrm{ff}}$ does not have a significant impact on the global star formation rate, but dramatically changes cluster properties, with increasing $\epsilon_{\mathrm{ff}}$ leading to a higher maximum cluster mass, shorter age spread of stars within clusters, and higher integrated star formation efficiencies. We also explore the redshift evolution of the observable cluster mass function, finding that most massive clusters have formed at high redshift $z>4$. Extrapolation of cluster disruption to $z=0$ produces good agreement with both the Galactic globular cluster mass function and age-metallicity relation. Our results emphasize the importance of using small-scale properties of galaxies to calibrate subgrid models of star cluster formation and feedback.

Thomas Pähtz, Orencio Duŕan, Francesco Comola

Sediment transport by atmospheric winds shapes the surface and affects the climates of planetary bodies. Reliably predicting the occurrence and rate of sediment transport in the Solar System has been notoriously difficult because fluid density, grain size and soil cohesiveness vary across many orders of magnitude. Here, we use recent advances in analytical and numerical sediment transport modeling to derive general scaling relations for planetary transport. In particular, we show that the equations of motion of rebounding grains predict that the minimum threshold fluid shear velocity needed to sustain transport (transport cessation threshold) scales with the particle-fluid-density ratio ($s$) as $s^{1/3}$, in contrast to the $s^{1/2}$-scaling exhibited by the threshold for transport initiation. The grain size corresponding to this minimum is in the range $80{-}290~\mu\mathrm{m}$ for Solar System bodies. Our results, summarized in phase diagrams for the cessation threshold, mean transport rate and dust emission potential, explain the observed eastward propagation of Titan's dunes, in spite of a predominantly westward wind circulation, indicate active dust cycles on Earth, Mars and Titan, and suggest marginal but active atmospheric transport on Venus, Triton and Pluto.

Henry J. Leath, Michael A. Beasley, Alexandre Vazdekis, Nuria Salvador-Rusiñol, Anastasia Gvozdenko

15 pages, 13 figures, accepted for publication in MNRAS

The leading method for the determination of relevant stellar population parameters of unresolved extragalactic Globular Clusters is through the study of their integrated spectroscopy, where Balmer line-strength indices are considered to be age sensitive. Previously, a splitting in the highly optimised spectral line-strength index H$\beta_o$ was observed in a sample of Galactic globular clusters at all metallicities resulting in an apparent "upper branch" and "lower branch" of globular clusters in the H$\beta_o$ - [MgFe] diagram. This was suggested to be caused by the presence of hot Blue straggler stars (BSSs), resulting in an underestimation of 'spectroscopic' ages in the upper branch. Over a decade on, we look to re-evaluate these findings. We make use of new, large Galactic Globular Cluster integrated spectroscopy datasets. To produce a large, homogeneously combined sample we have considered a number of factors including the radial dependence of Balmer and metal lines. Using this new sample, in disagreement with previous work, we find the splitting in H$\beta_o$ only occurs at intermediate to high metallicities ([M/H]$>-1$), and is not the result of an increased fraction of BSSs, but rather is due to an increased Helium abundance. We explore the possible impact of varying Helium on simple stellar population models to provide a theoretical basis for our hypothesis and then use the relationship between upper branch candidacy and enhanced Helium to predict the Helium content of three M31 clusters. We discuss what this can tell us about their mass and fraction of first generation stars.

Giovanni Gandolfi, Andrea Lapi, Stefano Liberati

39 pages, 22 figures. Accepted by ApJ

We look for empirical evidence of a non-minimal coupling (NMC) between dark matter (DM) and gravity in the dynamics of local spiral galaxies. In particular, we consider a theoretically motivated NMC that may arise dynamically from the collective behavior of the coarse-grained DM field (e.g., via Bose-Einstein condensation) with averaging/coherence length $L$. In the Newtonian limit, this NMC amounts to modify the Poisson equation by a term $L^{2} \nabla^{2} \rho$ proportional to the Laplacian of the DM density itself. We show that such a term, when acting as a perturbation over the standard Navarro-Frenk-White (NFW) profile of cold DM particles, can substantially alter the dynamical properties of galaxies, in terms of their total radial acceleration within the disk and rotation velocity. Specifically, we find that this NMC model can properly fit the stacked rotation curves of local spiral galaxies with different velocities at the optical radius, including dwarfs and low-surface brightness systems, at a level of precision comparable to, and in some instances even better than, the phenomenological Burkert profile. Finally we show that, extrapolating down to smaller masses the scaling of $L$ vs. halo mass found from the above rotation curve analysis, the NMC model can adequately reproduce the radial acceleration relation (or RAR) in shape and normalization down to the dwarf spheroidal galaxy range, a task which constitutes a serious challenge for alternative DM models even inclusive of baryonic effects.

A. Andrianov, S. Chernov, I. Girin, S. Likhachev, A. Lyakhovets, Yu. Shchekinov

Accepted to publish in Physical Review D

Photon rings near the edge of a black hole shadow is supposed to be a unique tool to validate general relativity and provide reliable measurements of principal black hole parameters: spin and mass. Such measurements are possible though only for nearby supermassive black holes (SMBH) with Space-Earth Very Long Baseline Interferometry (S-VLBI) in the submillimeter wavelength range. For subrings to be distinguished S-EVLBI observations with long baselines at the Lagrangian Sun-Earth L2 libration point are needed. However, the average fluxes of nearby SMBH: Sagittarius A$^\ast$ (Sgr A$^\ast$) and M87$^\ast$ -- $F_\nu\sim 1$ Jy, are still insufficient to detect the signal from the photon rings with even such long baselines. We argue that only manifestations of flares in the submillimeter waveband in their accretion disks can reveal observable signals from the photon rings with the S-EVLBI at L2. Such observations will become possible within the planned join program of the {\it Event Horizon Telescope} (EHT) and {\it Millimetron Space Observatory} (MSO), and within the planned {\it next generation EHT} (ngEHT) project. Two different observational tests for photons rings are discussed. The first one involves observations of a time series of responds from subsequent subrings as can be seen in a 1D visibility function within the join EHT-MSO configuration, the second one -- measurements of an increase of the angle between subsequent subrings in the 2D VLBI image which can be obtained within the ngEHT project.

M. Serylak, V. Venkatraman Krishnan, P. C. C. Freire, T. M. Tauris, M. Kramer, M. Geyer, A. Parthasarathy, M. Bailes, M. C. i Bernadich, S. Buchner, M. Burgay, F. Camilo, A. Karastergiou, M. E. Lower, A. Possenti, D. J. Reardon, R. M. Shannon, R. Spiewak, I. H. Stairs, W. van Straten

20 pages, 9 figures, 5 tables, accepted for publication by the A&A

PSR J0955$-$6150 is a member of a class of eccentric MSP+He WD systems (eMSPs), whose binary evolution is poorly understood and believed to be different to that of traditional MSP+He WD systems. Measuring the masses of the stars in this system is important for testing hypotheses for the formation of eMSPs. We have carried out observations of this pulsar with the Parkes and MeerKAT radio telescopes. Our observations reveal a strong frequency evolution of this pulsar's intensity, with a spectral index ($\alpha$) of $-3.13(2)$. The sensitivity of MeerKAT has resulted in a $>10$-fold improvement in the timing precision compared to older Parkes observations. Combined with the 8-year timing baseline, it has allowed precise measurements of a proper motion and three orbital "post-Keplerian" parameters: the rate of advance of periastron, $\dot{\omega} = 0.00152(1) \, {\rm deg} \, yr^{-1}$ and the orthometric Shapiro delay parameters, $h_3 = 0.89(7) \, \mu$s and $\varsigma = 0.88(2)$. Assuming general relativity, we obtain $M_{p} = 1.71(2) \, M_{\odot}$ for the mass of the pulsar and $M_{c} = 0.254(2) \, M_{\odot}$ for the mass of the companion; the orbital inclination is 83.2(4) degrees. We find that the spin axis has a misalignment relative to the orbital angular momentum of $> 4.8$ degrees at 99% CI. While the value of $M_{\rm p}$ falls within the wide range observed in eMSPs, $M_{\rm c}$ is significantly smaller than expected, allowing several formation hypotheses being ruled out. $M_{\rm c}$ is also significantly different from the expected value for an ideal low mass X-ray binary evolution scenario. The putative misalignment between the spin axis of the pulsar and the orbital angular momentum suggests that the unknown process that created the orbital eccentricity of the binary was also capable of changing its orbital orientation, an important evidence for understanding the origin of eMSPs.

A.V. Smith Castelli, A.R. Lopes, A. Cortesi, P. Dimauro, R. Cid Fernandes, G. Lucatelli, C. Mendes de Oliveira, F. Almeida-Fernandes, J.T.S.C. Batista, D. Brambila, S. Dhiwar, P. Lopes, K. Saha

3 pages, 5 figures. Accepted to be published in Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia

The photometric and morphological analysis of galaxies in clusters provides invaluable information regarding the evolutionary stage of the cluster itself. In addition, it helps to understand how the environment affects the properties of the galaxies and, as a consequence, their evolutionary path. In this contribution we present the first steps on the photometric and morphological analysis of galaxies in the Fornax cluster using S-PLUS data. We expect that the S-PLUS novel filter set and wide field coverage allow us to obtain new information about Fornax and its galaxy population.

Tidal effects in gravitational-wave (GW) observations from binary neutron star mergers have the potential to probe ultra-dense matter and shed light on the unknown nuclear equation of state of neutron stars. Tidal effects in inspiralling neutron star binaries become relevant at GW frequencies of a few hundred Hz and require detectors with exquisite high-frequency sensitivity. Third generation GW detectors such as the Einstein Telescope or Cosmic Explorer will be particularly sensitive in this high-frequency regime, allowing us to probe neutron star tides beyond the adiabatic approximation. Here we assess whether dynamical tides can be measured from a neutron star inspiral. We find that the measurability of dynamical tides depends strongly on the neutron star mass and equation of state. For a semi-realistic population of 10,000 inspiralling binary neutron stars, we conservatively estimate that on average $\mathcal{O}(50)$ binaries will have measurable dynamical tides. As dynamical tides are characterised not only by the star's tidal deformability but also by its fundamental ($f$-) mode frequency, they present a possibility of probing higher-order tidal effects and test consistency with quasi-universal relations. For a GW170817-like signal in a third generation detector network, we find that the stars' $f$-mode frequencies can be measured to within a few hundred Hz.

F. Haberl, C. Maitra, G. Vasilopoulos, P. Maggi, A. Udalski, I. M. Monageng, D. A. H. Buckley

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

The Magellanic Clouds host a large population of high-mass X-ray binary (HMXB) systems, but although the Large Magellanic Cloud (LMC) is an order of magnitude more massive than the Small Magellanic Cloud, significantly fewer HMXBs are known. We conducted a search for new HMXBs in XMM-Newton observations, which we performed to investigate supernova remnant candidates in the supergiant shells LMC5 and LMC7. The three observed fields are located in regions, which were little explored in X-rays before. We analysed the XMM-Newton data to look for sources with hard X-ray spectrum and counterparts with optical colours and brightness typical for HMXBs. We report the discovery of three new Be/X-ray binaries, two of them showing pulsations in their X-ray flux. With a luminosity of 6.5e34 erg/s, XMMU J045315.1-693242 in LMC7 was relatively X-ray faint. The long-term OGLE I-band light curve of the V = 15.5 mag counterpart suggests a 49.6 day or 24.8 day orbital period for the binary system. XMMU J045736.9-692727, also located in LMC7 was brighter with a luminosity of 5.6e35 erg/s and hard spectrum with a power-law photon index of 0.63. The X-ray flux revealed clear pulsations with a period of 317.7 s. We obtained optical high resolution spectra from the V = 14.2 mag counterpart using the SALT-HRS spectrograph. Halpha and Hbeta were observed in emission with complex line profiles and equivalent widths of -8.0 A and -1.3 A, respectively. The I-band light curve obtained from OGLE shows a series of four strong outbursts followed by a sudden drop in brightness by more than 1 mag within 73-165 days and a recovery to the level before the outbursts. RX J0524.2-6620, previously classified as X-ray binary candidate, is located at the eastern part of LMC5. We report the discovery of 360.7 s pulsations. (abridged)

Rayhane Sadeghi, Ehsan Tavabi

8 pages, 10 figs., accepted in MNRAS

In this investigation oscillation periods in Mg II k line intensity, brightness temperature, and Doppler velocity obtained above Magnetic Bright Points (MBPs) are investigated. For that purpose, data from the Interface Region Imaging Spectrometer (IRIS) observing the higher chromosphere and transition region (TR) were analysed together with imaging and magnetogram data obtained by Solar Dynamics Observatory (SDO). The MBPs were identified in combining Si IV 1403 A Slit Jaw Images (SJIs) with the magnetogram information from the Heliospheric and Magnetic Imager (HMI). A time-slice analysis followed by a wavelet inspection were carried out on the Mg II k (2796 A and 10,000 K) resonance lines for the detection of the oscillation period. Finally, a power spectrum analysis was performed to characterise the oscillations with the result that network points feature a typical intensity, temperature, and velocity oscillation period of about 300 seconds; The internetwork points have a mean intensity oscillation period of about 180 seconds, mean temperature oscillation period of about 202, and mean velocity oscillation period of about 202 seconds. In addition, one BP was analysed in detail, which demonstrates intensity oscillation periods with a value of 500 seconds, obviously not related to the common 3- or 5-minute oscillations found typically elsewhere in chromospheric/photospheric structures.

We derive the spectrum and analyse the detectability prospects of secondary gravity waves (GWs) associated to primordial black hole (PBH) production in a class of string inflationary models called Fibre Inflation. The inflationary potential features a near inflection point that induces a period of ultra slow-roll responsible for an enhancement of the scalar perturbations which can lead to PBHs with different masses and contributions to dark matter (DM) in agreement with current observational bounds, including CMB constraints on the scalar spectral index and the tensor-to-scalar ratio. This enhancement of the curvature perturbations sources secondary GWs which can be detected by either LISA or DECIGO and BBO, depending on the GW frequency but regardless of the amount of PBH DM since secondary GWs remain detectable even if the PBH contribution to DM is exponentially suppressed. The possibility to see a secondary GW signal is instead due to the presence of an ultra slow-roll epoch between CMB horizon exit and the end of inflation.

Federica Pompa, Francesco Capozzi, Olga Mena, Michel Sorel

5 pages, 3 figures

Time of flight delay in the supernova neutrino signal offers a unique tool to set model-independent constraints on the absolute neutrino mass. The presence of a sharp time structure during a first emission phase, the so-called neutronization burst in the electron neutrino flavor time distribution, makes this channel a very powerful one. Large liquid argon underground detectors will provide precision measurements of the time dependence of the electron neutrino fluxes. We derive here a new $\nu$ mass sensitivity attainable at the future DUNE far detector from a future supernova collapse in our galactic neighborhood, finding a sub-eV reach under favorable scenarios. These values are competitive with those expected for laboratory direct neutrino mass searches.

We analyze the effects of chiral symmetry restoration in hadronic matter, including the lowest-lying baryonic resonance $\Delta$ based on the parity doublet model. We study the role of $\Delta$ and its chiral partner on the equation of state (EoS) of dense matter under neutron star (NS) conditions of $\beta$-equilibrium and charge neutrality. We find that the softening of the EoS driven by the early onset of $\Delta$ matter due to partial restoration of chiral symmetry allows accommodating the modern multi-messenger astrophysical constraints on the mass, radius, and tidal deformability. The softening above the saturation density is accompanied by subsequent stiffening at high densities. We also find that the matter composition in the NS cores may be different upon variations of the repulsive interactions of $\Delta$ baryons in hadronic matter.

Z. Feizi Mangoudehi

32 pager, 15 figures accepted in Astrophysics and Space Science

The purpose of this paper is to study the Tsallis agegraphic dark energy with an interaction term between dark energy and dark matter in the DGP brane-world scenario. For this, we assume some initial conditions to obtain the dark energy density, deceleration, dark energy EoS, and total EoS parameters. Then, we analyze the statefinder parameters, $\omega'{}_{DE}-\omega_{DE}$ plots, and classical stability features of the model. The results state that the deceleration parameter provides the phase transition from decelerated to accelerated phase. The $\omega_{DE}$ graphs show the phantom behavior, while the $\omega_{tot}$ exhibits the quintessence and phantom during the evolution of the Universe. Following the graphs, the Statefinder analysis shows the quintessence behavior of the model for the past and present. However, it tends to the $\Lambda CDM$ in the following era. The $\omega'{}_{DE}-\omega_{DE}$ plot indicates the thawing or freezing area depending on the type of era and different values of $b^{2}$, $\delta$, and $m$. By the square of the sound speed, we see the model is stable in the past, stable or unstable at the current time, and unstable in the future for selected values of $b^{2}$, $\delta$, and $m$. To test the model, we use the recent Hubble data. We also employ Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) to compare the model with the $\Lambda CDM$ as the reference model. In addition, we test the model using the $H-z$ plot, and we see a turning point in the future time. The results from the best fit values for the $\omega_{tot}$ plot emphasize that the Universe is in the quintessence region in the current time. It will enter the phantom phase, and then it will approach the $\Lambda$ state in the future. But, the $\omega_{DE}$ always stays on the phantom region. The model is unstable in the present and progressive era.

Tessa Baker, Gianluca Calcagni, Anson Chen, Matteo Fasiello, Lucas Lombriser, Katarina Martinovic, Mauro Pieroni, Mairi Sakellariadou, Gianmassimo Tasinato, Daniele Bertacca, Ippocratis D. Saltas

38+20 pages, 22 figures

The propagation speed of gravitational waves, $c_T$, has been tightly constrained by the binary neutron star merger GW170817 and its electromagnetic counterpart, under the assumption of a frequency-independent $c_T$. Drawing upon arguments from Effective Field Theory and quantum gravity, we discuss the possibility that modifications of General Relativity allow for transient deviations of $c_T$ from the speed of light at frequencies well below the band of current ground-based detectors. We motivate two representative Ans\"atze for $c_T(f)$, and study their impact upon the gravitational waveforms of massive black hole binary mergers detectable by the LISA mission. We forecast the constraints on $c_T(f)$ obtainable from individual systems and a population of sources, from both inspiral and a full inspiral-merger-ringdown waveform. We show that LISA will enable us to place stringent independent bounds on departures from General Relativity in unexplored low-frequency regimes, even in the absence of an electromagnetic counterpart.

The characteristics of the cosmic microwave background provide circumstantial evidence that the hot radiation-dominated epoch in the early universe was preceded by a period of inflationary expansion. Here, we show how a measurement of the stochastic gravitational wave background can reveal the cosmic history and the physical conditions during inflation, subsequent pre- and reheating, and the beginning of the hot big bang era. This is exemplified with a particularly well-motivated and predictive minimal extension of the Standard Model which is known to provide a complete model for particle physics -- up to the Planck scale, and for cosmology -- back to inflation.

We investigate in scale-invariant $B-L$ scenario where the Standard Model (SM) is supplemented with a dark scalar $\phi$ associated with gauge \& Yukawa interactions, described by the couplings $g_{BL}$ and $y$ respectively, leading to radiative plateau inflation at scale $\phi=M$ in the ultraviolet (UV), while dynamically generating the Electroweak and Seesaw scales via Coleman-Weinberg in the infrared (IR). This is particularly achieved implementing threshold correction at an energy scale $\mu_T$ arising due to the presence of vector-like fermions. We show that implementing the inflationary observables makes the couplings solely dependent on the plateau scale or inflection-point scale $M$, leaving us with only two independent parameters $M$ and $\mu_T$. Within the theoretically consistent parameter space regions defined by $m_{Z_{BL}} > 850~\rm GeV$, from the assumption of independent evolution of the dark sector couplings from the SM couplings and $M < 5.67$ times the Planck scale ($M_P$) required for the realisation of inflationary plateau-like behaviour of the potential around $\phi=M$, we identify the parameter space that is excluded by current LHC results involving searches for the heavy $Z_{BL}$ boson. For typical benchmark points in the viable parameter regions, we estimate the reheating temperature to be $\sim\mathcal{O}(TeV)$ thus consistent with the standard Big Bang Nucleosynthesis (BBN) constraints. For typical benchmark points ($M=5.67,~1,~0.1~M_P$) we predict the scales of inflation to be $\mathcal{H}_{inf}=2.79\times10^{12}$ GeV, $1.53\times10^{10}$ GeV and $1.53\times10^7$ GeV respectively.