Papers for Wednesday, Jun 23 2021

Searches for counterparts to multimessenger events with optical imagers use difference imaging to detect new transient sources. However, even with existing artifact detection algorithms, this process simultaneously returns several classes of false positives: false detections from poor quality image subtractions, false detections from low signal-to-noise images, and detections of pre-existing variable sources. Currently, human visual inspection to remove remaining false positives is a central part of multimessenger follow-up observations, but when next generation gravitational wave and neutrino detectors come online and increase the rate of multimessenger events, the visual inspection process will be prohibitively expensive. We approach this problem with two convolutional neural networks operating on the difference imaging outputs. The first network focuses on removing false detections and demonstrates an accuracy above 95\%. The second network focuses on sorting all real detections by the probability of being a transient source within a host galaxy and distinguishes between various classes of images that previously required additional human inspection. We find the number of images requiring human inspection will decrease by a factor of 1.5 using our approach alone and a factor of 3.6 using our approach in combination with existing algorithms, facilitating rapid multimessenger counterpart identification by the astronomical community.

We set out to rederive the 8 Hubble parameter values obtained from estimated relative galaxy ages in Simon et al, PRD 71 (2005). We find that to obtain the level of precision claimed in $H(z)$, unrealistically small galaxy age uncertainties have to be assumed. Also, some parameter values will be correlated. In our analysis we find that the uncertainties in the Hubble parameter values are significantly larger when 8 independent $H(z)$ are obtained using Monte Carlo sampling. Smaller uncertainties can be obtained using Gaussian processes, but at the cost of strongly correlated results. We do not obtain any useful constraints on the Hubble parameter from the galaxy data employed.

Image Processing algorithms for vision-based navigation require reliable image simulation capacities. In this paper we explain why traditional rendering engines may present limitations that are potentially critical for space applications. We introduce Airbus SurRender software v7 and provide details on features that make it a very powerful space image simulator. We show how SurRender is at the heart of the development processes of our computer vision solutions and we provide a series of illustrations of rendered images for various use cases ranging from Moon and Solar System exploration, to in orbit rendezvous and planetary robotics.

We investigate the contribution of sub-Chandrasekhar mass Type Ia supernovae to the chemical enrichment of the Gaia Sausage galaxy, the progenitor of a significant merger event in the early life of the Milky Way. Using a combination of data from Nissen & Schuster (2010), the 3rd GALAH data release (with 1D NLTE abundance corrections) and APOGEE data release 16, we fit analytic chemical evolution models to a 9-dimensional chemical abundance space (Fe, Mg, Si, Ca, Cr, Mn, Ni, Cu, Zn) in particular focusing on the iron-peak elements, Mn and Ni. We find that low [Mn/Fe] $\sim-0.15\,\mathrm{dex}$ and low [Ni/Fe] $\sim-0.3\,\mathrm{dex}$ Type Ia yields are required to explain the observed trends beyond the [$\alpha$/Fe] knee of the Gaia Sausage (approximately at [Fe/H] $=-1.4\,\mathrm{dex}$). Comparison to theoretical yield calculations indicate a significant contribution from sub-Chandrasekhar mass Type Ia supernovae in this system (from $\sim60$% to $100$% depending on the theoretical model with an additional $\pm10$% systematic from NLTE corrections). We compare to results from other Local Group environments including dwarf spheroidal galaxies, the Magellanic Clouds and the Milky Way's bulge, finding the Type Ia [Mn/Fe] yield must be metallicity-dependent. Our results suggest that sub-Chandrasekhar mass channels are a significant, perhaps even dominant, contribution to Type Ia supernovae in metal-poor systems, whilst more metal-rich systems could be explained by metallicity-dependent sub-Chandrasekhar mass yields, possibly with additional progenitor mass variation related to star formation history, or an increased contribution from Chandrasekhar mass channels at higher metallicity.

We use ZFIRE and ZFOURGE observations with the Spectral Energy Distribution (SED) fitting tool Prospector to reconstruct the star formation histories (SFHs) of proto-cluster and field galaxies at $z\sim 2 $ and compare our results to the TNG100 run of the IllustrisTNG cosmological simulation suite. In the observations, we find that massive proto-cluster galaxies ($\log[{\rm M}_{\ast}/{\rm M}_{\odot}]>$10.5) form $45 \pm 8 \%$ of their total stellar mass in the first $2$ Gyr of the Universe compared to $31 \pm 2 \%$ formed in the field galaxies. In both observations and simulations, massive proto-cluster galaxies have a flat/declining SFH with decreasing redshift compared to rising SFH in their field counterparts. Using IllustrisTNG, we find that massive galaxies ($\log[{\rm M}_{\ast}/{\rm M}_{\odot}] \geq 10.5$) in both environments are on average $\approx190$ Myr older than low mass galaxies ($\log[{\rm M}_{\ast}/{\rm M}_{\odot}]= 9-9.5$). However, the difference in mean stellar ages of cluster and field galaxies is minimal when considering the full range in stellar mass ($\log[{\rm M}_{\ast}/{\rm M}_{\odot}] \geq 9$). We explore the role of mergers in driving the SFH in IllustrisTNG and find that massive cluster galaxies consistently experience mergers with low gas fraction compared to other galaxies after 1 Gyr from the Big Bang. We hypothesize that the low gas fraction in the progenitors of massive cluster galaxies is responsible for the reduced star formation.

We present resolved thermodynamic profiles out to 500 kpc, about $r_{500}$, of the $z=1.75$ galaxy cluster IDCS J1426.5+3508 with 40 kpc resolution. Thanks to the combination of Sunyaev-Zel'dovich and X-ray datasets, IDCS J1426.5+3508 becomes the most distant cluster with resolved thermodynamic profiles. These are derived assuming a non-parametric pressure profile and a very flexible model for the electron density profile. The shape of the pressure profile is flatter than the universal pressure profile. The IDCS J1426.5+3508 temperature profile is increasing radially out to 500 kpc. To identify the possible future evolution of IDCS J1426.5+3508 , we compared it with its local descendants that numerical simulations show to be $0.65\pm0.12$ dex more massive. We found no evolution at 30 kpc, indicating a fine tuning between cooling and heating at small radii. At $30<r<300$ kpc, our observations show that entropy and heat must be deposited with little net gas transfer, while at 500 kpc the gas need to be replaced by a large amount of cold, lower entropy gas, consistent with theoretical expectation of a filamentary gas stream, which brings low entropy gas to 500 kpc and energy at even smaller radii. At $r \gtrsim 400$ kpc the polytropic index takes a low value, which indicates the presence of a large amount of non-thermal pressure. Our work also introduces a new definition of the evolutionary rate, which uses unscaled radii, unscaled thermodynamic quantities, and different masses at different redshifts to compare ancestors and descendants. It has the advantage of separating cluster evolution, dependence on mass, pseudo-evolution and returns a number with unique interpretation, unlike other definitions used in literature.

We present 17 transit light curves of seven known warm-Jupiters observed with the CHaracterising ExOPlanet Satellite (CHEOPS). The light curves have been collected as part of the CHEOPS Guaranteed Time Observation (GTO) program that searches for transit-timing variation (TTV) of warm-Jupiters induced by a possible external perturber to shed light on the evolution path of such planetary systems. We describe the CHEOPS observation process, from the planning to the data analysis. In this work we focused on the timing performance of CHEOPS, the impact of the sampling of the transit phases, and the improvement we can obtain combining multiple transits together. We reached the highest precision on the transit time of about 13-16 s for the brightest target (WASP-38, G = 9.2) in our sample. From the combined analysis of multiple transits of fainter targets with G >= 11 we obtained a timing precision of about 2 min. Additional observations with CHEOPS, covering a longer temporal baseline, will further improve the precision on the transit times and will allow us to detect possible TTV signals induced by an external perturber.

We continue our program of investigation of the proper motions of spectrally separated structures in the Herbig-Haro outflows with the aid of Fabry-Perot scanning interferometry. This work mainly focuses on the physical nature of various structures in the jets. The aim of the present study is to measure the proper motions of the previously discovered kinematically separated structures in the working surface of the HH 83 collimated outflow. We used observations from two epochs separated by 15 years, which were performed on the 6m telescope with Fabry-Perot scanning interferometer. We obtained images corresponding to different radial velocities for the two separate epochs, and used them to measure proper motions. In the course of our data analysis, we discovered a counter bow-shock of HH 83 flow with positive radial velocity, which makes this flow a relatively symmetric bipolar system. The second epoch observations confirm that the working surface of the flow is split into two structures with an exceptionally large (250 km\ s$^{-1}$) difference in radial velocity. The proper motions of these structures are almost equal, which suggests that they are physically connected. The asymmetry of the bow shock and the turning of proper motion vectors suggests a collision between the outflow and a dense cloud. The profile of the H$\alpha$ line for the directly invisible infrared source HH 83 IRS, obtained by integration of the data within the reflection nebula, suggests it to be of P Cyg type with a broad absorption component characteristic of the FU Ori like objects. If this object underwent an FU Ori type outburst, which created the HH 83 working surfaces, its eruption took place about 1500 years ago according to the kinematical age of the outflow.

We present polarimetric observations of the Keyhole Nebula in the Carina Nebula Complex carried out using the Stratospheric Observatory for Infrared Astronomy. The Keyhole Nebula located to the west of $\eta$ Carinae is believed to be disturbed by the stellar winds from the star. We observed the Keyhole Nebula at 89 $\mu$m wavelength with the HAWC+ instrument. The observations cover the entire Keyhole Nebula spanning 8$'$ by 5$'$ with central position RA = 10:44:43 and Dec = -59:38:04. The typical uncertainty of polarization measurement is less than 0.5\% in the region with intensity above 5,500 MJy sr$^{-1}$. The polarization has a mean of 2.4\% with a standard deviation of 1.6\% in the region above this intensity, similar to values in other high--mass star--forming regions. The magnetic field orientation in the bar--shaped structure is similar to the large--scale magnetic field orientation. On the other hand, the magnetic field direction in the loop is not aligned with the large--scale magnetic fields but has tight alignment with the loop itself. Analysis of the magnetic field angles and the gas turbulence suggests that the field strength is $\sim$70 $\mu$G in the loop. A simple comparison of the magnetic field tension to the ram pressure of $\eta$ Carinae's stellar wind suggests that the magnetic fields in the Keyhole Nebula are not strong enough to maintain the current structure against the impact of the stellar wind, and that the role of the magnetic field in resisting stellar feedback in the Keyhole Nebula is limited.

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

We apply a friends-of-friends (FoF) algorithm to identify galaxy clusters and we use the catalog to explore the evolutionary synergy between BCGs and their host clusters. We base the cluster catalog on the dense HectoMAP redshift survey (2000 redshifts deg$^{-2}$). The HectoMAP FoF catalog includes 346 clusters with 10 or more spectroscopic members. We list these clusters and their members (5992 galaxies with a spectroscopic redshift). We also include central velocity dispersions ($\sigma_{*, BCG}$) for all of the FoF cluster BCGs, a distinctive feature of the HectoMAP FoF catalog. HectoMAP clusters with higher galaxy number density (80 systems) are all genuine clusters with a strong concentration and a prominent BCG in Subaru/Hyper Suprime-Cam images. The phase-space diagrams show the expected elongation along the line-of-sight. Lower-density systems include some false positives. We establish a connection between BCGs and their host clusters by demonstrating that $\sigma_{*,BCG}/\sigma_{cl}$ decreases as a function of cluster velocity dispersion ($\sigma_{cl}$), in contrast, numerical simulations predict a constant $\sigma_{*, BCG}/\sigma_{cl}$. Sets of clusters at two different redshifts show that BCG evolution in massive systems is slow over the redshift range $z < 0.4$. The data strongly suggest that minor mergers may play an important role in BCG evolution in these clusters ($\sigma_{cl} \gtrsim 300$ km s$^{-1}$). For systems of lower mass ($\sigma_{cl} < 300$ km s$^{-1}$), the data indicate that major mergers may play a significant role. The coordinated evolution of BCGs and their host clusters provides an interesting test of simulations in high density regions of the universe.

We compare Baryonic Acoustic Oscillation (BAO) and Redshift Space Distortion (RSD) measurements from recent galaxy surveys with their Fisher matrix based predictions. Measurements of the position of the BAO signal lead to constraints on the comoving angular diameter distance $D_{M}$ and the Hubble distance $D_{H}$ that agree well with their Fisher matrix based expectations. However, RSD-based measurements of the growth rate $f \sigma_{8}$ do not agree with the predictions made before the surveys were undertaken, even when repeating those predictions using the actual survey parameters. We show that this is due to a combination of effects including degeneracies with the geometric parameters $D_{M}$ and $D_{H}$, and optimistic assumptions about the scale to which the linear signal can be extracted. We show that measurements using current data and large-scale modelling techniques extract an equivalent amount of signal to that in the linear regime for $k\simlt 0.08 \,h\,{\rm Mpc}^{-1}$, remarkably independent of the sample properties and redshifts covered.

We obtained high-resolution spectra of the ultra-cool M-dwarf TRAPPIST-1 during the transit of its planet `b' using two high dispersion near-infrared spectrographs, IRD instrument on the Subaru 8.2m telescope and HPF instrument on the 10m Hobby-Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet `b' using the APO/ARCTIC, which assisted us to capture the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets `b', `e' and `f' from IRD, we attempt to constrain the atmospheric escape of the planet using the He I triplet 10830 {\AA} absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet related absorption, we place an upper limit on the equivalent widths of <7.754 m{\AA} for planet `b', <10.458 m{\AA} for planet `e', and <4.143 m{\AA} for planet `f' at 95% confidence from the IRD data, and <3.467 m{\AA} for planet `b' at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This non-detection of extended atmospheres having low mean-molecular weight in all three planets aids in further constraining their atmospheric composition by steering the focus towards the search of high molecular weight species in their atmospheres.

We present a new approach for fast calculation of gravitational lensing properties, including the lens potential, deflection angles, convergence, and shear, of elliptical Navarro-Frenk-White (NFW) and Hernquist density profiles, by approximating them by superpositions of elliptical density profiles for which simple analytic expressions of gravitational lensing properties are available. This model achieves high fractional accuracy better than $10^{-4}$ in the range of the radius normalized by the scale radius of $10^{-4}-10^3$. These new approximations are $\sim 300$ times faster in solving the lens equation for a point source compared with the traditional approach resorting to expensive numerical integrations, and are implemented in {\tt glafic} software.

Molecular cations are present in various astronomical environments, most notably in cometary atmospheres and tails where sunlight produces exceptionally bright near-UV to visible transitions. Such cations typically have longer-wavelength and brighter electronic emission than their corresponding neutrals. A robust understanding of their near-UV to visible properties would allow these cations to be used as tools for probing the local plasma environments or as tracers of neutral gas in cometary environments. However, full spectral models are not possible for characterization of small, oxygen containing molecular cations given the body of molecular data currently available. The five simplest such species (H2O+, CO+2 , CO+, OH+, and O+2 ) are well characterized in some spectral regions but are lacking robust reference data in others. Such knowledge gaps hinder fully quantitative models of cometary spectra, specifically, hindering accurate estimates of physical-chemical processes originating with the most common molecules in comets. Herein the existing spectral data are collected for these molecules and highlight the places where future work is needed, specifically where the lack of such data would greatly enhance the understanding of cometary evolution.

AGBs and YSOs often share the same domains in IR color-magnitude or color-color diagrams leading to potential mis-classification. We extracted a list of AGB interlopers from the published YSO catalogues using the periodogram analysis on NEOWISE time series data. YSO IR variability is typically stochastic and linked to episodic mass accretion. Furthermore, most variable YSOs are at an early evolutionary stage, with significant surrounding envelope and/or disk material. In contrast, AGBs are often identified by a well defined sinusoidal variability with periods of a few hundreds days. From our periodogram analysis of all known low mass YSOs in the Gould Belt, we find 85 AGB candidates, out of which 62 were previously classified as late-stage Class III YSOs. Most of these new AGB candidates have similar IR colors to O-rich AGBs. We observed 73 of these AGB candidates in the H2O, CH3OH and SiO maser lines to further reveal their nature. The SiO maser emission was detected in 10 sources, confirming them as AGBs since low mass YSOs, especially Class III YSOs, do not show such maser emission. The H2O and CH3OH maser lines were detected in none of our targets.

Millimetre continuum observations of debris discs can provide insights into the physical and dynamical properties of the unseen planetesimals that these discs host. The material properties and collisional models of planetesimals leave their signature on the grain size distribution, which can be traced through the millimetre spectral index. We present 8.8 mm observations of the debris discs HD 48370, CPD 72 2713, HD 131488, and HD 32297 using the Australian Telescope Compact Array (ATCA) as part of the PLanetesimals Around TYpicalPre-main seqUence Stars (PLATYPUS) survey. We detect all four targets with a characteristic beam size of 5 arcseconds and derive a grain size distribution parameter that is consistent with collisional cascade models and theoretical predictions for parent planetesimal bodies where binding is dominated by self-gravity. We combine our sample with 19 other millimetre-wavelength detected debris discs from the literature and calculate a weighted mean grain size power law index which is close to analytical predictions for a classical steady state collisional cascade model. We suggest the possibility of two distributions of q in our debris disc sample; a broad distribution (where q is approximately 3.2 to 3.7) for "typical" debris discs (gas-poor/non-detection), and a narrow distribution (where q is less than 3.2) for bright gas-rich discs. Or alternatively, we suggest that there exists an observational bias between the grain size distribution parameter and absolute flux which may be attributed to the detection rates of faint debris discs at cm wavelengths.

We present analyses of host galaxy properties of type 1 and type 2 X-ray selected AGNs in the XMM-XXL field, which have available optical spectroscopic classification. We model their optical to far-infrared spectral energy distributions (SEDs) using the X-CIGALE code. X-CIGALE allows the fitting of X-ray flux and accounts for the viewing angle of dusty torus and the attenuation from polar dust. By selecting matched type 1 and 2 subsamples in the X-ray luminosity and redshift parameter space, we find that both types live in galaxies with similar star formation. However, type 2 AGN tend to reside in more massive systems ($10.87^{+0.06}_{-0.12}\,\rm M_\odot$) compared to their type 1 counterparts ($10.57^{+0.20}_{-0.12}\,\rm M_\odot$). In the second part of our analysis, we compare the spectroscopic classification with that from the SED fitting. X-CIGALE successfully identifies all spectroscopic type 2 sources either by estimating an inclination angle that corresponds to edge on viewing of the source or by measuring increased polar dust in these systems. $\sim 85\%$ of spectroscopic type 1 AGN are also identified as such, based on the SED fitting analysis. There is a small number of sources ($\sim 15\%$ of the sample), that present broad lines in their spectra, but show strong indications of obscuration, based on SED analysis. These, could be systems that are viewed face on and have an extended dust component along the polar direction. The performance of X-CIGALE in classifying AGN is similar at low and high redshifts, under the condition that there is sufficient photometric coverage. Finally, usage of optical/mid-IR colour criteria to identify optical red AGN (${\it{u}}-\rm W3$), suggests that these criteria are better suited for IR selected AGN and their efficiency drops for the low to moderate luminosity sources included in X-ray samples.

Aims.Molecular masers, including methanol and hydroxyl masers, and in particular the ones in excited rotational states (ex-OHmasers), are one of the most informative tools for studying star-forming regions. So, the discovery, of new maser sources in theseregions is of great importance. Many studies and surveys of ex-OH maser sources have been carried out in the southern celestialhemisphere, but only a few have been done in the northern hemisphere. The specific aim of this work is to close this gap.Methods.The star-forming regions in the northern hemisphere with known active methanol masers were observed to search for newex-OH maser sources with the 32 m and 16 m radio telescopes of the Ventspils International Radio Astronomy Centre (VIRAC).Results.Three OH maser lines in the excited state at the 6035 MHz in three northern hemisphere star-forming regions are detected.The maser 189.030+0.783 was previously known, but we suggest this maser is a possible variable. We confirm recent detections ofthe ex-OH masers 85.41+0.00 and 90.92+1.49 by other authors. The magnetic field strength in the masering regions is estimated byusing right circular polarization (RCP) and left circular polarization (LCP) pair splitting. The high-velocity resolution provides uswith an estimation of a comparatively small magnetic field strength for the 189.030+0.783 and 90.92+1.49 star-forming regions

Recent high-resolution interferometric images of submillimetre galaxies (SMGs) reveal fascinatingly complex morphologies. This raises a number of questions: how does the relative orientation of a galaxy affect its observed submillimetre emission, and does this result in an `orientation bias' in the selection and analysis of such galaxies in flux-limited cosmological surveys? We investigate these questions using the Simba cosmological simulation paired with the dust radiative transfer code Powderday. We select eight simulated SMGs ($S_{850}\gtrsim2$ mJy) at $z = 2$, and measure the variance of their `observed' emission over 50 random orientations. Each galaxy exhibits significant scatter in its emission close to the peak of the thermal dust emission, with variation in flux density of up to $\sim$50 mJy at the peak. This results in an appreciable dispersion in the inferred dust temperatures and infrared luminosities ($16^{\mathrm{th}}-84^{\mathrm{th}}$ percentile ranges of 5 K and 0.1 dex, respectively) and therefore a fundamental uncertainty in derived parameters such as dust mass and star formation rate ($\sim$30% for the latter using simple calibrations). Using a Monte Carlo simulation we also assess the impact of orientation on flux-limited surveys, finding a bias in the selection of SMGs towards those with face-on orientations, as well as those at lower redshifts. We predict that the orientation bias will affect flux-limited single-dish surveys, most significantly at THz frequencies, and this bias should be taken into account when placing the results of targeted follow-up studies in a statistical context.

We conduct three-dimensional hydrodynamical simulations of eccentric common envelope jets supernova (CEJSN) impostors, i.e., a neutron star (NS) that crosses through the envelope of a red supergiant star on a highly eccentric orbit and launches jets as it accretes mass from the envelope. Because of numerical limitations we apply a simple prescription where we inject the assumed jets' power into two opposite conical regions inside the envelope. We find the outflow morphology to be very complicated, clumpy, and non-spherical, having a large-scale symmetry only about the equatorial plane. The outflow morphology can substantially differ between simulations that differ by their jets' power. We estimate by simple means the light curve to be very bumpy, to have a rise time of one to a few months, and to slowly decay in about a year to several years. These eccentric CEJSN impostors will be classified as `gap' objects, i.e., having a luminosity between those of classical novae and typical supernovae (termed also ILOTs for intermediate luminosity optical transients). We strengthen a previous conclusion that CEJSN impostors might account for some peculiar ILOTs, in particular those that might repeat over timescales of months to years.

We exploit the unique characteristics of a sample of open clusters (OCs) and field stars for which high-precision 7Li abundances and stellar parameters are homogeneously derived by the Gaia-ESO Survey (GES). We derive possibly undepleted 7Li abundances for 26 OCs and star forming regions with ages from young to old spanning a large range of Galactocentric distances, which allows us to reconstruct the local late Galactic evolution of lithium as well as its current abundance gradient along the disc. Field stars are added to look further back in time and to constrain 7Li evolution in other Galactic components. The data are then compared to theoretical tracks from chemical evolution models that implement different 7Li forges. We find that the upper envelope of the 7Li abundances measured in field stars of nearly solar metallicities traces very well the level of lithium enrichment attained by the ISM as inferred from observations of cluster stars. We confirm previous findings that the abundance of 7Li in the solar neighbourhood does not decrease at supersolar metallicity. The comparison of the data with the chemical evolution model predictions favours a scenario in which the majority of the 7Li abundance in meteorites comes from novae. Current data also seem to suggest that the nova rate flattens out at later times. This requirement might have implications for the masses of the white dwarf nova progenitors and deserves further investigation. Neutrino-induced reactions taking place in core-collapse supernovae also produce some fresh lithium. This likely makes a negligible contribution to the meteoritic abundance, but could be responsible for a mild increase of the 7Li abundance in the ISM of low-metallicity systems that would counterbalance the astration processes.

The tilt angle of sunspot groups is crucial in the BL type dynamo. Some studies have shown that the tilt coefficient is anti-correlated with the cycle strength. If the anti-correlation exists, it will be shown to act as an effective nonlinearity of the BL-type dynamo to modulate the solar cycle. However, some studies have shown that the anti-correlation has no statistical significance. We aim to investigate the causes behind the controversial results of tilt angle studies and to establish whether the tilt coefficient is indeed anti-correlated with the cycle strength. We first analyzed the tilt angles from DPD. Based on the methods applied in previous studies, we took two criteria to select the data, along with the linear and square-root functions to describe Joy's law, and three methods to derive the tilt coefficients for cycles 21-24. This allowed us to evaluate different methods based on comparisons of the differences among the tilt coefficients and the tilt coefficient uncertainties. Then we utilized Monte Carlo experiments to verify the results. Finally, we extended these methods to analyze the separate hemispheric DPD data and the tilt angle data from Kodaikanal and Mount Wilson. The tilt angles exhibit an extremely wide scatter due to both the intrinsic mechanism for its generation and measurement errors, for instance, the unipolar regions included in data sets. Different methods to deal with the uncertainties are mainly responsible for the controversial character of the previous results. The linear fit to the tilt-latitude relation of sunspot groups with $\Delta s>2.5$ of a cycle carried out without binning the data can minimize the effect of the tilt scatter on the uncertainty of the tilt coefficient. Based on this method the tilt angle coefficient is anti-correlated with the cycle strength with strong statistical significance.

We present the results of ALMA band 3 observations of a nearby type Ic supernova (SN) 2020oi. Under the standard assumptions on the SN-circumstellar medium (CSM) interaction and the synchrotron emission, the data indicate that the CSM structure deviates from a smooth distribution expected from the steady-state mass loss in the very vicinity of the SN (~10^{15} cm), which is then connected to the outer smooth distribution (~10^{16} cm). This structure is further confirmed through the light curve modeling of the whole radio data set as combined with data at lower frequency previously reported. Being an explosion of a bare carbon-oxygen (C+O) star having a fast wind, we can trace the mass-loss history of the progenitor of SN 2020oi in the final year. The inferred non-smooth CSM distribution corresponds to fluctuations on the sub-year time scale in the mass-loss history toward the SN explosion. Our finding suggests that the pre-SN activity is likely driven by the accelerated change in the nuclear burning stage in the last moments just before the massive star's demise. The structure of the CSM derived in this study is beyond the applicability of the other methods at optical wavelengths, highlighting an importance and uniqueness of quick follow-up observations of SNe by ALMA and other radio facilities.

In the context of tomographic cosmic shear surveys, there exists a nulling transformation of weak lensing observations (also called BNT transform) that allows us to simplify the correlation structure of tomographic cosmic shear observations, as well as to build observables that depend only on a localised range of redshifts and thus independent from the low-redshift/small-scale modes. This procedure renders possible accurate, and from-first-principles, predictions of the convergence and aperture mass one-point distributions (PDF). We here explore other consequences of this transformation on the (reduced) numerical complexity of the estimation of the joint PDF between nulled bins and demonstrate how to use these results to make theoretical predictions.

The magnetic nature of the formation of solar active regions lies at the heart of understanding solar activity and, in particular, solar eruptions. A widespread model, used in many theoretical studies, simulations and the interpretation of observations is that the basic structure of an active region is created by the emergence of a large tube of pre-twisted magnetic field. Despite plausible reasons and the availability of various proxies suggesting the veracity of this model, there has not yet been any direct observational evidence of the emergence of large twisted magnetic flux tubes. Thus, the fundamental question, "are active regions formed by large twisted flux tubes?" has remained open. In this work, we answer this question in the affirmative and provide direct evidence to support this. We do this by investigating a robust topological quantity, called magnetic winding, in solar observations. This quantity, combined with other signatures that are currently available, provides the first direct evidence that large twisted flux tubes do emerge to create active regions.

Context. Extragalactic Planetary Nebulae (PNe) are useful distance indicators and are often used to trace the dark-matter content in external galaxies. At the same time, PNe can also be used as probes of their host galaxy stellar populations and to help understanding the later stages of stellar evolution. Previous works have indicated that specific number of PNe per stellar luminosity can vary across different galaxies and as a function of stellar-population properties, for instance increasing with decreasing stellar metallicity. Aims. In this study we further explore the importance of stellar metallicity in driving the properties of the PNe population in early-type galaxies, using three edge-on galaxies in the Fornax cluster offering a clear view into their predominantly metal-rich and metal-poor regions near the equatorial plane or both below and above it, respectively . Methods. Using VLT-MUSE integral-field observations and dedicated PNe detection procedures, we construct the PNe luminosity function and compute the luminosity-specific number of PNe alpha in both in- and off-plane regions of our edge-on systems. Results. Comparing these alpha values with metallicity measurements also based on the same MUSE data, we find no evidence for an increase in the specific abundance of PNe when transitioning between metal-rich and metal-poor regions. Conclusions. Our analysis highlights the importance of ensuring spatial consistency to avoid misleading results when investigating the link between PNe and their parent stellar populations and suggest that in passively-evolving systems variations in the specific number of PNe may pertain to rather extreme metallicity regimes found either in the innermost or outermost regions of galaxies.

Astrophysical black holes are thought to be the Kerr black holes predicted by general relativity, but macroscopic deviations from the Kerr solution can be expected from a number of scenarios involving new physics. In Paper I, we studied the reflection features in NuSTAR and XMM-Newton spectra of the supermassive black hole at the center of the galaxy MCG-06-30-15 and we constrained a set of deformation parameters proposed by Konoplya, Rezzolla & Zhidenko (Phys. Rev. D93, 064015, 2016). In the present work, we analyze the X-ray data of a stellar-mass black hole within the same theoretical framework in order to probe a different curvature regime. We consider a NuSTAR observation of the X-ray binary EXO 1846-031 during its outburst in 2019. As in the case of Paper I, all our fits are consistent with the Kerr black hole hypothesis, but some deformation parameters cannot be constrained well.

Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM using these RT schemes to hot Jupiter simulation results from the literature. We perform a HD 209458b-like simulation with the three schemes and compare their results. These simulations are then post-processed to compare their observable differences. The semi-grey scheme results show qualitative agreement with previous studies in line with variations seen between GCM models. The real gas model reproduces well the temperature and dynamical structures from other studies. After post-processing our non-grey picket fence scheme compares very favourably with the real gas model, producing similar transmission spectra, emission spectra and phase curve behaviours. Exo-FMS is able to reliably reproduce the essential features of contemporary GCM models in the hot gas giant regime. Our results suggest the picket fence approach offers a simple way to improve upon RT realism beyond semi-grey schemes.

We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: MOND and Verlinde's emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6\sigma$ between the RARs of early- and late-type galaxies (split by S\'{e}rsic index and $u-r$ colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour. The difference might be explained if only the early-type galaxies have significant ($M_{gas} \approx M_*$) circumgalactic gaseous haloes. The observed behaviour is also expected in $\Lambda$CDM models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a $\Lambda$CDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys will be able to further distinguish between MG and $\Lambda$CDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced.

We propose and apply a new test of Einstein's Equivalence Principle (EEP) based on the gravitational redshift induced by the central super massive black hole of quasars in the surrounding accretion disk. Specifically, we compare the observed gravitational redshift of the Fe III$\lambda\lambda$2039-2113 emission line blend in quasars with the predicted values in a wide, uncharted, cosmic territory ($0 \lesssim z_{cosm}\lesssim3$). For the first time we measure, with statistical uncertainties comparable or better than those of other classical methods outside the Solar System, the ratio between the observed gravitational redshifts and the theoretical predictions in 10 independent cosmological redshift bins in the $1 \lesssim z_{cosm}\lesssim3$ range. The average of the measured over predicted gravitational redshifts ratio in this cosmological redshift interval is $\langle z^m_g/z_g^p\rangle=1.05\pm 0.06$ with scatter $0.13\pm 0.05$ showing no cosmological evolution of EEP within these limits. This method can benefit from larger samples of measurements with better S/N ratios, paving the way for high precision tests (below 1\%) of EEP on cosmological scales.

The possibility that primordial black holes (PBHs) form a part of dark matter has been considered over a wide mass range from the Planck mass ($10^{-5}~\rm g$) to the level of the supermassive black hole in the center of the galaxy. Primordial origin might be one of the most important formation channel of massive black holes. We propose the lensing effect of very long baseline interferometer observations of compact radio sources with extremely high angular resolution as a promising probe for the presence of intergalactic PBHs in the mass range $\sim10^2$-$10^9~M_{\odot}$. For a sample of well-measured 543 compact radio sources, no millilensing multiple images are found with angular separations between $0.2$ milliarcsecond and $50$ milliarcseconds. From this null search result, we derive that the fraction of dark matter made up of PBHs in the mass range $\sim10^4$-$10^8~M_{\odot}$ is $\lesssim0.56\%$ at $68\%$ confidence level.

The covariance matrix $\boldsymbol{\Sigma}$ of non-linear clustering statistics that are measured in current and upcoming surveys is of fundamental interest for comparing cosmological theory and data and a crucial ingredient for the likelihood approximations underlying widely used parameter inference and forecasting methods. The extreme number of simulations needed to estimate $\boldsymbol{\Sigma}$ to sufficient accuracy poses a severe challenge. Approximating $\boldsymbol{\Sigma}$ using inexpensive but biased surrogates introduces model error with respect to full simulations, especially in the non-linear regime of structure growth. To address this problem we develop a matrix generalization of Convergence Acceleration by Regression and Pooling (CARPool) to combine a small number of simulations with fast surrogates and obtain low-noise estimates of $\boldsymbol{\Sigma}$ that are unbiased by construction. Our numerical examples use CARPool to combine GADGET-III $N$-body simulations with fast surrogates computed using COmoving Lagrangian Acceleration (COLA). Even at the challenging redshift $z=0.5$, we find variance reductions of at least $\mathcal{O}(10^1)$ and up to $\mathcal{O}(10^4)$ for the elements of the matter power spectrum covariance matrix on scales $8.9\times 10^{-3}<k_\mathrm{max} <1.0$ $h {\rm Mpc^{-1}}$. We demonstrate comparable performance for the covariance of the matter bispectrum, the matter correlation function and probability density function of the matter density field. We compare eigenvalues, likelihoods, and Fisher matrices computed using the CARPool covariance estimate with the standard sample covariance estimators and generally find considerable improvement except in cases where $\Sigma$ is severely ill-conditioned.

Nine point sources appeared within half an hour on a region within $\sim$ 10 arcmin of a red-sensitive photographic plate taken in April 1950 as part of the historic Palomar Sky Survey. All nine sources are absent on both previous and later photographic images, and absent in modern surveys with CCD detectors which go several magnitudes deeper. We present deep CCD images with the 10.4-meter Gran Telescopio Canarias (GTC), reaching brightness $r \sim 26$ mag, that reveal possible optical counterparts, although these counterparts could equally well be just chance projections. The incidence of transients in the investigated photographic plate is far higher than expected from known detection rates of optical counterparts to e.g.\ flaring dwarf stars, Fast Radio Bursts (FRBs), Gamma Ray Bursts (GRBs) or microlensing events. One possible explanation is that the plates have been subjected to an unknown type of contamination producing mainly point sources with of varying intensities along with some mechanism of concentration within a radius of $\sim$ 10 arcmin on the plate. If contamination as an explanation can be fully excluded, another possibility is fast (t $<0.5$ s) solar reflections from objects near geosynchronous orbits. An alternative route to confirm the latter scenario is by looking for images from the First Palomar Sky Survey where multiple transients follow a line.

Amino acids are the essential keys in chemistry that contribute to the study of the formation of life. The complex organic molecule glycine (NH$_{2}$CH$_{2}$COOH) is the simplest amino acid that has been investigated in the interstellar medium for a long period to search for a potential connection between the Universe and the origin of life. Several attempts have failed to search for glycine in the last forty years, which made the researcher look for some glycine precursor in the interstellar medium as an alternative approach. We report the successful detection of the rotational emission lines of interstellar glycine with confirmer I and II in the hot molecular core G10.47+0.03 between the frequency range of $\nu$ = 158.6$-$160.4 GHz with Atacama Large Millimeter/Submillimeter Array (ALMA) observation. In hot molecular core G10.47+0.03, the fractional abundance of glycine is found between the range of (4.01$-$4.61)$\times$10$^{-10}$ which refers to the "medium warm-up" case. The detection of glycine in the interstellar medium is very complicated but many theoretical and laboratory studies indicated the possibilities of the presence of glycine and its precursors in hot molecular cores. We also detected the emission lines of complex organic molecules CHOCHOHCH$_{2}$OH, $^{13}$CH$_{2}$OHCHO, CHD(OH)CHO, CH$_{2}$OH$^{13}$CHO, cis-CH$_{2}$OHCHO, G$^{\prime}$Gg$^{\prime}$-CH$_{2}$(OH)CH(OH)CH$_{2}$OH, and CH$_{2}$DOH in the hot molecular core G10.47+0.03.

Parametric resonance in a single-field inflationary model with a periodic structure on the potential gives rise to curvature perturbations with large amplitudes on small scales, which could result in observable primordial black holes (PBHs) and concomitant gravitational waves (GWs) induced by curvature perturbations in the radiation-dominated era. In such a model, besides GWs associated with the PBH formation, there exists a stochastic GW background sourced by inflaton perturbations resonantly amplified during inflation. We compute the energy spectra of the induced GWs produced both during inflation and in the radiation-dominated era, and find that the peak of the energy spectrum of the former is much higher than that of the latter, but is located at a lower frequency. Moreover, the energy spectrum of the induced GWs produced during inflation exhibits a unique oscillating character in the ultraviolet region. Both the stochastic GW backgrounds are expected to be detected by future space-based laser interferometers.

Mapping the evolution of galaxy colours, from blue star-forming to red passive systems, is fundamental to understand the processes involved in galaxy evolution. To this end, we reconstruct the colour evolution of low-redshift galaxies, combining stellar templates with star formation and metallicity histories of galaxies from the Galaxy And Mass Assembly survey and Shark semi-analytic model. We use these colour histories to robustly characterise the evolution of red and blue galaxy populations over cosmic time. Using a Gaussian Mixture Model to characterise the colour distribution at any given epoch and stellar mass, we find both observations and simulations strongly favour a model with only two populations (blue and red), with no evidence for a third "green" population. We map the evolution of mean, weight, and scatter of the blue and red populations as a function of both stellar mass and lookback time. Using our simulated galaxy catalogue as a testbed, we find that we can accurately recover galaxies colour histories up to a lookback time of $\sim6$ Gyr. We find that both populations show little change in the mean colour for low-mass galaxies, while the colours at the massive end become significantly redder with time. The stellar mass above which the galaxy population is predominantly red decreases by $0.3$ dex in the last $5$ Gyrs. We find a good agreement between observations and simulations, with the largest tension being that massive galaxies from Shark are too blue (a known issue with many galaxy evolution models).

{The Galactic centre (GC) is a unique astrophysical laboratory to study the stellar population of galactic nuclei because it is the only galactic nucleus whose stars can be resolved down to milliparsec scales. However, the extreme and spatially highly variable interstellar extinction towards the GC poses a serious obstacle to photometric stellar classification.} {Our goal is to identify hot, massive stars in the nuclear stellar disc (NSD) region through combining near-infrared (NIR) and mid-infrared (MIR) photometry, and thus to demonstrate the feasibility of this technique, which may gain great importance with the arrival of the James Webb Space Telescope (JWST).} {We combined the GALACTICNUCLEUS NIR survey with the IRAC/Spitzer MIR survey of the GC. We applied the so-called Rayleigh-Jeans colour excess (RJCE) de-reddening method to our combined NIR-MIR data to identify potential hot stars in colour-magnitude diagrams (CMDs).} {Despite the very low angular resolution of IRAC we find 12 clear candidates for young massive stars among the $1\,065$ sources that meet our selection criteria. Seven out of these 12 stars are previously known hot, massive stars belonging to the Arches and Quintuplet clusters, as well as sources detected by the Hubble Space Telescope/NICMOS Paschen-$\alpha$ survey. Five of our massive star candidates have not been previously reported in the literature.} {We show that the RJCE method is a valuable tool to identify hot stars in the GC using photometry alone. Upcoming instruments with high angular resolution MIR imaging capabilities such as the JWST could surely make more substantial use of this de-reddening method and help establish a far more complete census of hot, young stars in the GC area than what is possible at the moment.}

We propose that fourteen co-moving groups of stars uncovered by Kounkel & Covey (2019) may be related to known nearby moving groups and bridge those and nearby open clusters with similar ages and space velocities. This indicates that known nearby moving groups may be spatially much more extended than previously though, and some of them might be parts of tidal tails around the cores of known open clusters, reminiscent of those recently found around the Hyades and a handful of other nearby clusters. For example, we find that both the nearby Carina and Columba associations may be linked to Theia 208 from Kounkel & Covey (2019) and together form parts of a large tidal tail around the Platais 8 open cluster. The AB Doradus moving group and Theia 301 may form a trailing tidal tail behind the Pleiades open cluster, with hints of a possible leading tidal tail in Theia 369. We similarly find that IC 2391 and its tidal tails identified by Meingast et al. (2021) may be extended by the nearby Argus association and are possibly further extended by Theia 115. The nearby Octans and Octans-Near associations, as well as Theia 94 and 95, may form a large tidal tail leading the poorly studied Platais 5 open cluster candidate. While a preliminary analysis of Gaia color-magnitude sequences hint that these structures are plausibly related, more observational evidence is still required to corroborate their consistent ages and space velocities. These observations may change our current understanding of nearby moving groups and the different pathways through which they can form. While some moving groups may have formed loosely in extended star-formation events with rich spatial structure, others may in fact correspond to the tidal tails of nearby open clusters.

A filament, a dense cool plasma supported by the magnetic fields in the solar corona, often becomes unstable and erupts. It is empirically known that the filament often demonstrates some activations such as a turbulent motion prior to eruption. In our previous study (Seki et al. 2017), we analysed the Doppler velocity of an H{\alpha} filament and found that the standard deviation of the line-of-sight-velocity (LOSV) distribution in a filament, which indicates the increasing amplitude of the small-scale motions, increased prior to the onset of the eruption. Here, we present a further analysis on this filament eruption, which initiated approximately at 03:40UT on 2016 November 5 in the vicinity of NOAA AR 12605. It includes a coronal line observation and the extrapolation of the surrounding magnetic fields. We found that both the spatially averaged micro-turbulence inside the filament and the nearby coronal line emission increased 6 and 10 hours prior to eruption, respectively. In this event, we did not find any significant changes in the global potential-field configuration preceding the eruption for the past 2 days, which indicates that there is a case in which it is difficult to predict the eruption only by tracking the extrapolated global magnetic fields. In terms of space weather prediction, our result on the turbulent motions in a filament could be used as the useful precursor of a filament eruption.

We report our analyses of the multi-epoch (2015-2017) ALMA archival data of the Class II binary system XZ Tau at Bands 3, 4 and 6. The millimeter dust continuum images show compact, unresolved (r <~ 15 au) circumstellar disks (CSDs) around the individual binary stars; XZ Tau A and B, with a projected separation of ~ 39 au. The 12CO (2-1) emission associated with those CSDs traces the Keplerian rotations, whose rotational axes are misaligned with each other (P.A. ~ -5 deg for XZ Tau A and ~ 130 deg for XZ Tau B). The similar systemic velocities of the two CSDs (VLSR ~ 6.0 km s-1) suggest that the orbital plane of the binary stars is close to the plane of the sky. From the multi-epoch ALMA data, we have also identified the relative orbital motion of the binary. Along with the previous NIR data, we found that the elliptical orbit (e = 0.742+0.025-0.034, a = 0''.172+0''.002-0''.003, and {\omega} = -54.2+2.0-4.7 deg) is preferable to the circular orbit. Our results suggest that the two CSDs and the orbital plane of the XZ Tau system are all misaligned with each other, and possible mechanisms to produce such a configuration are discussed. Our analyses of the multi-epoch ALMA archival data demonstrate the feasibility of time-domain science with ALMA.

The traditional clustering analyses of galaxy redshift surveys compress the clustering data into a set of late-time physical variables in a model-independent way. This approach has recently been extended by an additional shape variable encoding early-time physics information. We apply this new technique, ShapeFit, to SDSS-III BOSS data and show that it matches the constraining power of alternative, model-dependent approaches, which directly constrain the model's parameters adopting a cosmological model ab-initio. ShapeFit is $\sim30$ times faster, model-independent, naturally splits early- and late-time variables, and enables a better control of observational systematics.

Finite density effects can destabilize the metastable vacua in relaxion models. Focusing on stars as nucleation seeds, we derive the conditions that lead to the formation and runaway of a relaxion bubble of a lower energy minimum than in vacuum. The resulting late-time phase transition in the universe allows us to set new constraints on the parameter space of relaxion models. We also find that similar instabilities can be triggered by the large electromagnetic fields around rotating neutron stars.

In this paper we study the effect of electroweak sphaleron transition or electroweak phase transition (EWPT) in balancing baryon excess to the excess stable quarks of $4^{th}$ generation. Considering the conservation of all quantum number and charges, sphaleron transition between between baryons, leptons and $4^{th}$ family of leptons and quarks is possible. We have tried to established a possible definite relationship between the value and sign of the $4^{th}$ family excess relative to baryon asymmetry under the framework of second order EWPT. In passing by we show the small, yet negligible dilution in the pre-existing dark matter density due the sphaleron transition.

In this article, we reevaluate supernovae (SN) constraints on the diffusion time of neutrinos for a family of extensions of the Standard Model that incorporate new light scalar and vector mediators. We compute the neutrino mean free path, taking into account medium effects in the neutrino-nucleon scattering cross-section, and a radial dependence of the density, energy, and temperature inside the proto-neutron star to determine the coupling strengths compatible with SN1987A constraints on the time duration signal of diffusing neutrinos. We show that medium effects can induce an order of magnitude enhancement in the neutrino mean free path with respect to the vacuum calculation. The increase is more significant when new physics terms dominate over the Standard Model contribution (that is, for small mediator mass and large couplings). Finally, we interpret these results as bounds on the parameter space of a vector $U(1)_{B-L}$ model and scalar lepton number conserving and lepton number violating scenarios, improving on previous results in the literature where medium effects were ignored. We show that SN constraints on the neutrino diffusion time lie within regions of the parameter space that are already ruled out by other experimental constraints. We also comment on potential limits due to changes in the SN equation of state or right-handed neutrino free-streaming, but argue that detailed numerical simulations are needed to improve the reliability of these limits.

One-loop renormalised quantum effective action for gravity contains quadratic in curvature terms. We have found an approximate analytic black hole solution in quadratic gravity by keeping only the radial spherically symmetric fluctuations and dimensionally reducing the 4-dimensional (4D) theory down to the 2-dimensional (2D) dilaton gravity with a potential. The solution reduces to the Schwarzschild black hole in the limit of Einstein's gravity, but otherwise admits non-negative Arnowitt-Deser-Misner (ADM) and positive quasi-local Misner-Sharp masses that can differ significantly. We then study the thermodynamics of such quantum corrected black holes and compute their lifetime under Hawking evaporation. We note that for some range of parameters, black holes increase in mass while emitting Hawking radiation. This pathological behaviour is related to the negative energy states that are present in quadratic gravity. We also find that the micro-lensing of non-Schwarzschild black holes could significantly deviate from the micro-lensing of their Schwarzschild counterparts. These findings have important ramifications for the phenomenology of primordial black holes (PBHs) as dark matter. In particular, the quoted constraints on PBH dark matter from micro-lensing data can be completely evaded, thus making PBHs in the mass range $\sim 10^{-12} - 10~M_{\odot}$ viable dark matter candidates.