Papers for Wednesday, Apr 06 2022

The Einstein Telescope (ET) is going to bring a revolution for the future of multi-messenger astrophysics. In order to detect the counterparts of binary neutron star (BNS) mergers at high redshift, the high-energy observations will play a crucial role. Here, we explore the perspectives of ET, as single observatory and in a network of gravitational-wave (GW) detectors, operating in synergy with future $\gamma$-ray and X-ray satellites. We predict the high-energy emission of BNS mergers and its detectability in a theoretical framework which is able to reproduce the properties of the current sample of observed short GRBs (SGRB). We estimate the joint GW and high-energy detection rate for both the prompt and afterglow emissions, testing several combinations of instruments and observational strategies. We find that the vast majority of SGRBs detected in $\gamma$-rays will have a detectable GW counterpart; the joint detection efficiency approaches $100\%$ considering a network of third generation GW observatories. The probability of identifying the electromagnetic counterpart of BNS mergers is significantly enhanced if the sky localisation provided by GW instruments is observed by wide field X-ray monitors. We emphasize that the role of the future X-ray observatories will be very crucial for the detection of the fainter emission outside the jet core, which will allow us to probe the yet unexplored population of low-luminosity SGRBs in the nearby Universe, as well as to unveil the nature of the jet structure and the connections with the progenitor properties.

In this paper we use an RXTE library of spectral models from 10 black-hole and 9 pulsar X-ray binaries, as well as model spectra available in the literature from 13 extra-galactic Ultra-luminous X-ray sources (ULXs). We compute average bolometric corrections (BC=$\mathrm{L_{band}/L_{bol}}$) for our sample as a function of different accretion rates. We notice the same behaviour between black-hole and pulsars BCs only when ULX pulsars are included. These measurements provide a picture of the energetics of the accretion flow for an X-ray binary based solely on its observed luminosity in a given band. Moreover it can be a powerful tool in X-ray binary population synthesis models. Furthermore we calculate the X-ray (2-10 keV) to optical (V-band) flux ratios at different Eddington ratios for the black-hole X-ray binaries in our sample. This provides a metric of the maximum contribution of the disk to the optical emission of a binary system and better constraints on its nature (donor type etc). We find that the optical to X-ray flux ratio shows very little variation as a function of accretion rate, but testing for different disk geometries scenarios we find that the optical contribution of the disk increases as the $p$ value decreases ($T(r)\sim r^{-p}$). Moreover observational data are in agreement with a thicker disk scenario ($p<0.65$), which could also possibly explain the lack of observed high-inclination systems.

Accounting for all the relativistic effects, we have developed the fully nonlinear gauge-invariant formalism for describing the cosmological observables and presented the second-order perturbative expressions associated with light propagation and observations without choosing a gauge condition. For the first time, we have performed a complete verification of the validity of our second-order expressions by comparing their gauge-transformation properties from two independent methods: one directly obtained from their expressions in terms of metric perturbations and the other expected from their nonlinear relations. The expressions for the cosmological observables such as galaxy clustering and the luminosity distance are invariant under diffeomorphism and gauge-invariant at the observed position. We compare our results to the previous work and discuss the differences in the perturbative expressions. Our second-order gauge-invariant formalism constitutes a major step forward in the era of precision cosmology and its applications in the future will play a crucial role for going beyond the power spectrum and probing the early universe.

It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the $\alpha$-Fe plane. We employ Sloan Digital Sky Survey (SDSS-IV) data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE 2) to test this hypothesis. Adopting a random sampling technique, we contrast the abundances of 10 elements in Terzan 5 stars with those of their bulge field counterparts with comparable atmospheric parameters, finding that they differ at statistically significant levels. Abundances between the two groups differ by more than 1$\sigma$ in Ca, Mn, C, O, and Al, and more than 2$\sigma$ in Si and Mg. Terzan 5 stars have lower [$\alpha$/Fe] and higher [Mn/Fe] than their bulge counterparts. Given those differences, we conclude that Terzan 5 is not the remnant of a $major$ building block of the bulge. We also estimate the stellar mass of the Terzan 5 progenitor based on predictions by the Evolution and Assembly of GaLaxies and their Environments (EAGLE) suite of cosmological numerical simulations, concluding that it may have been as low as $\sim3\times10^8$ M$_\odot$ so that it was likely unable to significantly influence the mean chemistry of the bulge/inner disk, which is significantly more massive ($\sim10^{10}$ M$_\odot$). We briefly discuss existing scenarios for the nature of Terzan 5 and propose an observational test that may help elucidate its origin.

A wide array of astrophysical systems can only be accurately modelled when the behaviour of their baryonic gas components is well understood. The residual distribution (RD) family of partial differential equation (PDE) solvers can be used to solve the fluid equations, which govern this gas. We present a new implementation of the RD method to do just this. The solver efficiently and accurately calculates the evolution of the fluid, to second order accuracy in both time and space, across an unstructured Delaunay triangulation, built from an arbitrary distribution of vertices, in either 2D and 3D. We implement a novel new variable time stepping routine, which applies a drifting mechanism to greatly improve the computational efficiency of the method. We conduct extensive testing of the new implementation, demonstrating its innate ability to resolve complex fluid structures, even at very low resolution. Our implementation can resolve complex fluid structures with as few as 3-5 resolution elements, demonstrated by Kelvin-Helmholtz and Sedov blast tests. It includes three residual calculation modes, the LDA, N and blended schemes, each designed for different scenarios. These range from smooth flows (LDA), to extreme shocks (N), and scenarios where either may be encountered (blended). We compare our RD solver results to state-of-the-art solvers used in other astrophysical codes, demonstrating the competitiveness of the new approach, particularly at low resolution. This is of particular interest in large scale astrophysical simulations, where important structures, such as star forming gas clouds, are often resolved by small numbers of fluid elements.

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales.

Heating from active galactic nuclei (AGN) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magneto-hydrodynamic (MHD) simulations to model light AGN jet feedback with different accretion modes (Bondi-Hoyle-Lyttleton and cold accretion) in an idealised Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet density and resolution, the cluster self-regulates with central entropies and cooling times consistent with observed cool-core clusters in this non-cosmological setting. We find that increased jet efficiencies lead to more intermittent jet powers and enhanced star formation rates. Our fiducial low-density jets can easily be deflected by orbiting cold gaseous filaments, which redistributes angular momentum and leads to more extended cold gas distributions and isotropic bubble distributions. In comparison to our fiducial low momentum-density jets, high momentum-density jet heats less efficiently and enables the formation of a persistent cold-gas disc perpendicular to the jet that is centrally confined. Cavity luminosities measured from our simulations generally reflect the cooling luminosities of the intracluster medium (ICM) and correspond to averaged jet powers that are relatively insensitive to short periods of low-luminosity jet injection. Cold gas structures in our MHD simulations with low momentum-density jets generally show a variety of morphologies ranging from discy to very extended filamentary structures. In particular, magnetic fields are crucial to inhibit the formation of unrealistically massive cold gas discs by redistributing angular momentum between the hot and cold phases and by fostering the formation of elongated cold filaments that are supported by magnetic pressure.

In this work, we study the effect of differential rotation, finite temperature and strangeness on the quasi stationary sequences of hyper massive neutron stars (HMNS). We generate constant rest mass sequences of differentially rotating and uniformly rotating stars. The nucleonic matter relevant to the star interior is described within the framework of the relativistic mean field model with the DD2 parameter set. We also consider the strange $\Lambda$ hyperons using the BHB$\Lambda\phi$ equation of state (EoS). Additionally, we probe the behaviour of neutron stars (NS) with these compositions at different temperatures. We report that the addition of hyperons to the EoS produces a significant boost to the spin-up phenomenon. Moreover, increasing the temperature can make the spin-up more robust. We also study the impact of strangeness and thermal effects on the T/W instability. Finally, we analyse equilibrium sequences of a NS following a stable transition from differential rotation to uniform rotation. The decrease in frequency relative to angular momentum loss during this transition is significantly smaller for EoS containing hyperons, compared to nucleonic EoS.

We present predictions, derived from the EAGLE $\Lambda$CDM cosmological hydrodynamical simulations, for the abundance and of galaxies expected to be detected at high redshift by the {\it James Webb Space Telescope} (\JWST). We consider the galaxy as a whole and focus on the sub-population of progenitors of Milky Way (MW) analogues, defined to be galaxies with accretion histories similar to the MW's, that is, galaxies that underwent a merger the Gaia-Enceladus-Sausage (GES) event and that contain an analogue of the Large Magellanic Cloud (LMC) satellite today. We the luminosity function of all EAGLE galaxies in \JWST/NIRCam passbands, in the redshift range $z=2-8$, taking into account dust obscuration and different exposure times. For an exposure time of $T=10^5$s, average MW progenitors are observable as far back as $z\sim6$ in most bands, and this changes to $z\sim5$ and $z\sim4$ for the GES and LMC progenitors, respectively. The progenitors of GES and LMC analogues are, on average, $\sim 2$ and $\sim 1$ mag fainter than the MW progenitors at most redshifts. They lie, on average, within $\sim 60$ and $30$ arcsec, respectively, of their future MW host at all times, and thus will appear within the field-of-view of \JWST/NIRCam. We conclude that galaxies resembling the main progenitor of the MW and its major accreted components should be observable with \JWST\ beyond redshift $2$, providing a new and unique window in studying the formation history of our own galaxy.

The cosmic microwave background (CMB), carrying the inhomogeneous information of the very early universe, is of great significance for understanding the origin and evolution of our universe. However, observational CMB maps contain serious foreground contaminations from several sources, such as galactic synchrotron and thermal dust emissions. Here, we build a deep convolutional neural network (CNN) to recover the tiny CMB signal from various huge foreground contaminations. Focusing on the CMB temperature fluctuations, we find that the CNN model can successfully recover the CMB temperature maps with high accuracy, and that the deviation of the recovered power spectrum $C_\ell$ is smaller than the cosmic variance at $\ell>10$. We then apply this method to the current Planck observation, and find that the recovered CMB is quite consistent with that disclosed by the Planck collaboration, which indicates that the CNN method can provide a promising approach to the component separation of CMB observations. Furthermore, we test the CNN method with simulated CMB polarization maps based on the CMB-S4 experiment. The result shows that both the EE and BB power spectra can be recovered with high accuracy. Therefore, this method will be helpful for the detection of primordial gravitational waves in current and future CMB experiments. The CNN is designed to analyze two-dimensional images, thus this method is not only able to process full-sky maps, but also partial-sky maps. Therefore, it can also be used for other similar experiments, such as radio surveys like the Square Kilometer Array.

This work is part of a multi-wavelength program to study the effects of X-ray/UV/optical stellar radiation on the chemistry of the circumbinary disk around the young high-eccentricity binary DQ Tau. ALMA observations for near/around December 5, 2021 periastron were postponed due to bad weather, but supporting Swift-XRT-UVOT TOO observations were successful. These Swift observations along with previous X-ray-optical-mm data show that DQ Tau keeps exhibiting powerful flares near periastron, offering a unique laboratory for studies of flare effects on the gas-phase ion chemistry in protoplanetary disks.

Celestial bodies approximated with rigid triaxial ellipsoids in a two-body system can rotate chaotically due to the time-varying gravitational torque from the central mass. At small orbital eccentricity values, rotation is short-term orderly and predictable within the commensurate spin-orbit resonances, while at eccentricity approaching unity, chaos completely takes over. Here, we present the full 3D rotational equations of motion around all three principle axes for triaxial minor planets and two independent methods of numerical solution based on Euler rotations and quaternion algebra. The domains of chaotic rotation are numerically investigated over the entire range of eccentricity with a combination of trial integrations of Euler's equations of motion and the GALI($k$) method. We quantify the dependence of the order--chaos boundaries on shape by changing a prolateness parameter, and find that the main 1:1 spin-orbit resonance disappears for specific moderately prolate shapes already at eccentricities as low as 0.3. The island of short-term stability around the main 1:1 resonance shrinks with increasing eccentricity at a fixed low degree of prolateness and completely vanishes at approximately 0.8. This island is also encroached by chaos on longer time scales indicating longer Lyapunov exponents. Trajectories in the close vicinity of the 3:2 spin-orbit resonance become chaotic at smaller eccentricities, but separated enclaves of orderly rotation emerge at eccentricities as high as 0.8. Initial perturbations of rotational velocity in latitude away from the exact equilibrium result in a spectrum of free libration, nutation, and polar wander, which is not well matched by the linearized analysis omitting the inertial terms.

This study forges new tools to discriminate between dynamical and static dark energy. It provides accurate evolutionary templates of dynamical cosmological parameters and fundamental constants as analytic functions of the scale factor. They are designed to replace the commonly used parameterizations in likelihood calculations with evolutionary templates based on the physics of specific dynamical cosmologies and dark energy potentials. Thus they are termed Cosmology and Potential Specific, CPS, templates. A suite of CPS templates are calculated for a flat quintessence cosmology with a dark energy potential of the same mathematical form as the Higgs potential. This Higgs inspired, HI, polynomial potential produces a rich set of evolutions unlike most monomial potentials. The study produces CPS templates that are analytic functions of the scale factor. It uses a recently developed beta function formalism that provides a differential function for the scalar in terms of the scale factor. This establishes a methodology for easily producing CPS templates for other dark energy potentials and cosmologies to determine the likelihoods of dynamical cosmologies relative to Lambda CDM. The study also examines the evolution of fundamental constants such as the proton to electron mass ratio and the fine structure constant involving an intersection between particle physics and cosmology. Appendix A displays an abridged suite of CPS templates for flat quintessence and the HI dark energy potential.

Application Specific Integrated Circuits (ASICs) are used in space-borne instruments for signal processing and detector readout. The electrical interface of these ASICs to frontend printed circuit boards (PCBs) is commonly accomplished with wire bonds. Through Silicon Via (TSV) technology has been proposed as an alternative interconnect technique that will reduce assembly complexity of ASIC packaging by replacing wire bonding with flip-chip bonding. TSV technology is advantageous in large detector arrays where TSVs enable close detector tiling on all sides. Wafer-level probe card testing of TSV ASICs is frustrated by solder balls introduced onto the ASIC surface for flip-chip bonding that hamper alignment. Therefore, we developed the ASIC Test Stand (ATS) to enable rapid screening and characterization of individual ASIC die. We successfully demonstrated ATS operation on ASICs originally developed for CdZnTe detectors on the Nuclear Spectroscopic and Telescope Array (\textit{NuSTAR}) mission that were later modified with TSVs in a via-last process. We tested both back-side blind-TSVs and front-side through-TSVs, with results from internal test pulser measurements that demonstrate performance equal to or exceeding the probe card wafer-level testing data. The ATS can easily be expanded or duplicated in order to parallelize ASIC screening for large area imaging detectors of future space programs.

We perform maximum likelihood estimates (MLEs) for single and double flavor ultralight dark matter (ULDM) models using the Spitzer Photometry and Accurate Rotation Curves (SPARC) database. These estimates are compared to MLEs for several commonly used cold dark matter (CDM) models. By comparing various CDM models we find, in agreement with previous studies, that the Burkert and Einasto models tend to perform better than other commonly used CDM models. We focus on comparisons between the Einasto and ULDM models and analyze cases for which the ULDM particle masses are: free to vary; and fixed. For each of these analyses, we perform fits assuming the soliton and halo profiles are: summed together; and matched at a given radius. When we let the particle masses vary, we find a negligible preference for any particular range of particle masses, within $10^{-25}\,\text{eV}\leq m\leq10^{-19}\,\text{eV}$, when assuming the summed models. For the matched models, however, we find that almost all galaxies prefer particles masses in the range $10^{-23}\,\text{eV}\lesssim m\lesssim10^{-20}\,\text{eV}$. For both double flavor models we find that most galaxies prefer approximately equal particle masses. We find that the summed models give much larger variances with respect to the soliton-halo (SH) relation than the matched models. When the particle masses are fixed, the matched models give median and mean soliton and halo values that fall within the SH relation bounds, for most masses scanned. When the particle masses are fixed in the fitting procedure, we find the best fit results for the particle mass $m=10^{-20.5}\,\text{eV}$ (for the single flavor models) and $m_1=10^{-20.5}\,\text{eV}$, $m_2=10^{-20.2}\,\text{eV}$ for the double flavor, matched model. We discuss how our study will be furthered using a reinforcement learning algorithm.

Planck's Cosmic Microwave Background temperature and polarization observations are the premier dataset for constraining cosmological models. Cosmic variance limited temperature at large and intermediate scales today dominates the constraints; polarization provides additional constraining power and further scrutiny of the models. To complete this picture from Planck, ground-based experiments, such as the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) continue to add temperature and polarization measurements at small scales, allowing for the extraction of competitive cosmological constraints from the $TE$ and $EE$ power spectra. Matching at the same time all these stringent probes is a key challenge and validation step for any cosmological model. In particular, $\Lambda$CDM requires a tight consistency between the temperature and polarization measurements. In this paper, we present a number of methods to identify and quantify possible inconsistencies between temperature and polarization, we apply them to the latest Planck, ACT and SPT data and find no evidence for a deviation from $\Lambda$CDM. Application of these methods will have increased importance for future, more constraining CMB data.

We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Peg III, we find an elliptical half-light radius of $a_h{{=}}1.88^{+0.42}_{-0.33}$ arcminutes ($118^{+31}_{-30}$ pc) and $M_V{=}{-4.17}^{+0.19}_{-0.22}$; for Psc II, we measure $a_h{=}1.31^{+0.10}_{-0.09}$ arcminutes ($69\pm8$ pc) and $M_V{=}{-4.28}^{+0.19}_{-0.16}$. We do not find any morphological features that indicate a significant interaction between the two has occurred, despite their close separation of only $\sim$40 kpc. Using proper motions (PMs) from Gaia early Data Release 3, we investigate the possibility of any past association by integrating orbits for the two UFDs in a MW-only and a combined MW and Large Magellanic Cloud (LMC) potential. We find that including the gravitational influence of the LMC is crucial, even for these outer-halo satellites, and that a possible orbital history exists where Peg III and Psc II experienced a close ($\sim$10-20 kpc) passage about each other just over $\sim$1 Gyr ago, followed by a collective passage around the LMC ($\sim$30-60 kpc) just under $\sim$1 Gyr ago. Considering the large uncertainties on the PMs and the restrictive priors imposed to derive them, improved PM measurements for Peg III and Psc II will be necessary to clarify their relationship. This would add to the rare findings of confirmed pairs of satellites within the Local Group.

The radius valley, a dip in the radius distribution of exoplanets at ~1.9 Earth radii separates compact rocky Super-Earths and Sub-Neptunes with lower density. Various hypotheses have been put forward to explain the radius valley. Characterizing the radius valley morphology and its correlation to stellar properties will provide crucial observation constraints on its origin mechanism and deepen the understanding of planet formation and evolution. In this paper, the third part of the Planets Across the Space and Time (PAST) series, using the LAMOST-Gaia-Kepler catalog, we perform a systematical investigation into how the radius valley morphology varies in the Galactic context, i.e., thin/thick galactic disks, stellar age and metallicity abundance ([Fe/H] and [alpha/Fe]). We find that (1) The valley becomes more prominent with the increase of both age and [Fe/H]. (2) The number ratio of super-Earths to sub-Neptunes monotonically increases with age but decreases with [Fe/H] and [alpha/Fe]. (3) The average radius of planets above the valley (2.1-6 Earth radii) decreases with age but increases with [Fe/H]. (4) In contrast, the average radius of planets below the valley (R < 1.7 Earth radii) is broadly independent on age and metallicity. Our results demonstrate that the valley morphology as well as the whole planetary radius distribution evolves on a long timescale of giga-years, and metallicities (not only Fe but also other metal elements, e.g., Mg, Si, Ca, Ti) play important roles in planet formation and in the long term planetary evolution.

We report the mean metallicity and absolute magnitude of RR Lyrae stars in a sample of 37 globular clusters, calculated via the Fourier decomposition of their light curves and ad hoc semiempirical calibrations, in an unprecedented homogeneous approach. This enabled a new discussion of the metallicity dependence of the horizontal branch (HB) luminosity, as a fundamental distance indicator. The calibration for the RRab and RRc stars should be treated separately. For the RRab the dispersion is larger and non-linear. For the RRc stars the correlation is less steep, very tight and linear. The relevance of the HB structural parameter L, is highlighted and offer a non-linear calibration of the form MV ([Fe/H],L). Excellent agreement is found between values of [Fe/H] and MV from the light curve decomposition with spectroscopic values and distances obtained via Gaia-DR3 and HST. The variable stars census in 35 clusters includes 326 stars found by our program.

The nature of the long period radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627 for short) is discussed. We try to understand both its radio emission and pulsation in the neutron star scenario. We think that: (1) From the radio emission point of view, GLEAM-X J1627 can be a radio-loud magnetar. (2) From the rotational evolution point of view, GLEAM-X J1627 is unlikely to be an isolated magnetar. (3) The 1091s period is unlikely to be the precession period. (4) GLEAM-X J1627 may be a radio-loud magnetar spun-down by a fallback disk. The pulsar death line is modified due to the presence of a fallback disk. This may explain why GLEAM-X J1627 is still radio active with such a long pulsation period. General constraint on the neutron star magnetic field and initial disk mass are given analytically. Possible ways to discriminate between different modelings are also discussed.

For transient sources with timescales of 1-100 seconds, standardized imaging for all observations at each time step become impossible as large modern interferometers produce significantly large data volumes in this observation time frame. Here we propose a method based on machine learning and using interferometric closure products as input features to detect transient source candidates directly from the spatial frequency domain without imaging. We train a simple neural network classifier on a synthetic dataset of Noise/Transient/RFI events, which we construct to tackle the lack of labelled data. We also use the hyper-parameter dropout rate of the model to allow the model to approximate Bayesian inference, and select the optimal dropout rate to match the posterior prediction to the actual underlying probability distribution of the detected events. The overall F1-score of the classifier on the simulated dataset is greater than 85\%, with the signal-to-noise at 7$\sigma$. The performance of the trained neural network with Monte Carlo dropout is evaluated on semi-real data, which includes a simulated transient source and real noise. This classifier accurately identifies the presence of transient signals in the detectable signal-to-noise levels (above 4$\sigma$) with the optimal variance. Our findings suggest that a feasible radio transient classifier can be built up with only simulated data for applying to the prediction of real observation, even in the absence of annotated real samples for the purpose of training.

The time series of light reflected from exoplanets by future direct imaging can provide spatial information with respect to the planetary surface. We apply sparse modeling to the retrieval method that disentangles the spatial and spectral information from multi-band reflected light curves termed as spin-orbit unmixing. We use the $\ell_1$-norm and the Total Squared Variation norm as regularization terms for the surface distribution. Applying our technique to a toy model of cloudless Earth, we show that our method can infer sparse and continuous surface distributions and also unmixed spectra without prior knowledge of the planet surface. We also apply the technique to the real Earth data as observed by DSCOVR/EPIC. We determined the representative components that can be interpreted as cloud and ocean. Additionally, we found two components that resembled the distribution of land. One of the components captures the Sahara Desert, and the other roughly corresponds to vegetation although their spectra are still contaminated by clouds. Sparse modeling significantly improves the geographic retrieval, in particular, of cloud and leads to higher resolutions for other components when compared with spin-orbit unmixing using Tikhonov regularization.

We report the complete statistical planetary sample from the prime fields ($\Gamma \geq 2~{\rm hr}^{-1}$) of the 2019 Korea Microlensing Telescope Network (KMTNet) microlensing survey. We develop the optimized KMTNet AnomalyFinder algorithm and apply it to the 2019 KMTNet prime fields. We find a total of 14 homogeneously selected planets and report the analysis of three planetary events, KMT-2019-BLG-(1042,1552,2974). The planet-host mass ratios, $q$, for the three planetary events are $6.34 \times 10^{-4}, 4.89 \times 10^{-3}$ and $6.18 \times 10^{-4}$, respectively. A Bayesian analysis indicates the three planets are all cold giant planets beyond the snow line of their host stars. The 14 planets are basically uniform in $\log q$ over the range $-5.0 < \log q < -1.5$. This result suggests that the planets below $q_{\rm break} = 1.7 \times 10^{-4}$ proposed by the MOA-II survey may be more common than previously believed. This work is an early component of a large project to determine the KMTNet mass-ratio function, and the whole sample of 2016--2019 KMTNet events should contain about 120 planets.

The search for rotating radio transients (RRAT) at declination from -9o to +42o was carried out in the semi-annual monitoring data obtained on the Large Phased Array (LPA) radio telescope at the frequency of 111 MHz. A neural network was used to search for candidates. 4 new RRATs were detected, having dispersion measures (DM) 5-16 pc/cm3. A comparison with an earlier RRAT search conducted using the same data shows that the neural network reduced the amount of interference by 80 times, down to 1.3% of the initial amount of interferences. The loss of real pulsar pulses does not exceed 6% of their total number.

We have investigated the optical spectral behavior of a large sample of Fermi blazars (40 FSRQs and 13 BL Lacs), and found two new universal optical spectral behaviors. In the low state the optical spectrum gradually becomes softer (steeper) or harder (flatter) but more and more slowly when the brightness increases, and then tends to stable in the high state, which are briefly named the redder-stable-when-brighter (RSWB) and bluer-stable-when-brighter (BSWB) behaviors, respectively. 34 FSRQs and 7 BL Lacs exhibit clear RSWB behavior, and 2 FSRQs and 5 BL Lacs show distinct BSWB behavior, which mean that FSRQs favor more RSWB than BSWB behavior, while BL Lacs have no clear preference among both behaviors. We have put forward a unified nonlinear formula to quantitatively characterize the optical spectral behaviors of FSRQs and BL Lacs, which can fit both kinds of behaviors very well. We argue that the RSWB and BSWB behaviors originate from the same mechanism, and they are the universal optical spectral behaviors for blazars. The frequently observed redder-when-brighter (RWB) and bluer-when-brighter (BWB) trends can be considered to be the approximations of the behaviors of RSWB and BSWB, respectively. The rarely observed stable-when-brighter (SWB) trend can also be viewed as an approximation or a special case of the RSWB or BSWB behavior. We have developed a model with two constant-spectral-index components which can not only explain well both two kinds of optical spectral behaviors, but also successfully interpret the differential behaviors between FSRQs and BL Lacs.

Moving structures have been detected in coronal bright points and in a solar flare in active regions, which were bi-directional, symmetrical, simultaneous, and quasi-periodic (Ning & Guo 2014; Ning 2016; Li et al. 2016a). They could be regarded as observational evidence of plasma outflows via magnetic reconnection. In this article, we explored pairs of moving structures in fifteen ultraviolet bright points (UBPs), which were observed in the quiet Sun or inside a small active region on 19 November 2020, and measured by the High Resolution (HRI) Telescopes of the Extreme Ultraviolet Imager (EUI) on board the Solar Orbiter (SolO) in two passbands, HRIEUV 174 {\AA} and HRILy{\alpha} 1216 {\AA}. Moving structures observed in ten UBPs as starting from their bright cores and propagating toward two ends, are interpreted as diverging motions of bi-directional moving structures. These moving structures are also characterized by simultaneity and symmetry and in the case of seven UBPs they exhibit quasi-periodicity. They could be generated by outflows after magnetic reconnections. Moving structures seen in another five UBPs as originating from double ends and moving closer, and merging together, are manifested as converging motions, which might be caused by inflows through the magnetic reconnection, or might be interpreted as upflows driven by the chromospheric evaporation.

The interstellar medium (ISM) is a turbulent, highly structured multi-phase medium. State-of-the-art cosmological simulations of the formation of galactic discs usually lack the resolution to accurately resolve those multi-phase structures. However, small-scale density structures play an important role in the life cycle of the ISM, and determine the fraction of cold, dense gas, the amount of star formation and the amount of radiation and momentum leakage from cloud-embedded sources. Here, we derive a $statistical\, model$ to calculate the unresolved small-scale ISM density structure from coarse-grained, volume-averaged quantities such as the $gas\, clumping\, factor$, $\mathcal{C}$, and mean density $\left<\rho\right>_V$. Assuming that the large-scale ISM density is statistically isotropic, we derive a relation between the three-dimensional clumping factor, $\mathcal{C}_\rho$, and the clumping factor of the $4\pi$ column density distribution on the cloud surface, $\mathcal{C}_\Sigma$, and find $\mathcal{C}_\Sigma=\mathcal{C}_\rho^{2/3}$. Applying our model to calculate the covering fraction, i.e., the $4\pi$ sky distribution of optically thick sight-lines around sources inside interstellar gas clouds, we demonstrate that small-scale density structures lead to significant differences at fixed physical ISM density. Our model predicts that gas clumping increases the covering fraction by up to 30 per cent at low ISM densities compared to a uniform medium. On the other hand, at larger ISM densities, gas clumping suppresses the covering fraction and leads to increased scatter such that covering fractions can span a range from 20 to 100 per cent at fixed ISM density. All data and example code is publicly available at GitHub.

We present new CO(2-1) observations of a representative sample of 24 local (z$<$0.02) luminous infrared galaxies (LIRGs) obtained at high spatial resolution ($<$100 pc) from ALMA. We derive the effective radii of the CO(2-1) and the 1.3 mm continuum emissions using the curve-of-growth method. LIRGs show an extremely compact cold molecular gas distribution (median R$_{CO}$ $\sim$0.7 kpc), which is a factor 2 smaller than the ionized gas, and 3.5 times smaller than the stellar size. The molecular size of LIRGs is similar to that of early-type galaxies (R$_{CO}\sim$1 kpc) and about a factor of 6 more compact than local Spirals of similar stellar mass. Only the CO emission in low-z ULIRGs is more compact than these local LIRGs by a factor of 2. Compared to high-z (1$<$z$<$6) systems, the stellar sizes and masses of local LIRGs are similar to those of high-z MS star-forming galaxies (SFG) and about a factor of 2-3 lower than sub-mm galaxies (SMG). The molecular sizes of high-z MS SFGs and SMGs are larger than those derived for LIRGs by a factor of $\sim$3 and $\sim$8, respectively. These results indicate that while low-z LIRGs and high-z MS-SFGs have similar stellar masses and sizes, the regions of current star formation (traced by the ionized gas) and of potential star-formation (traced by the molecular gas) are substantially smaller in LIRGs, and constrained to the central kpc region. High-z galaxies represent a wider population but their star-forming regions are more extended, even covering the overall size of the host galaxy. High-z galaxies have larger fractions of gas than low-z LIRGs, and therefore the formation of stars could be induced by interactions and mergers in extended disks or filaments with large enough molecular gas surface density involving physical mechanisms similar to those identified in the central kpc of LIRGs.

We present SHEEP, a new machine learning approach to the classic problem of astronomical source classification, which combines the outputs from the XGBoost, LightGBM, and CatBoost learning algorithms to create stronger classifiers. A novel step in our pipeline is that prior to performing the classification, SHEEP first estimates photometric redshifts, which are then placed into the dataset as an additional feature for classification model training; this results in significant improvement in the subsequent classification performance. SHEEP contains two distinct classification methodologies: (i) Multi-class; (ii) one vs all with correction by a meta-learner. We demonstrate the performance of SHEEP for the classification of stars, galaxies and quasars using a dataset composed of SDSS and WISE photometry of 3.5 million astronomical sources. The resulting F1-scores are as follows: (i) 0.992 for galaxies; (ii) 0.967 for quasars; (iii) and 0.985 for stars. In terms of the F1-scores for the three classes, SHEEP is found to outperform the recent RandomForest-based classification approach of Clarke et al. (2020) using an essentially identical dataset. Our methodology also facilitates model and dataset explainability via feature importances; it also allows the selection of sources whose uncertain classifications may make them interesting sources for follow-up observations.

To understand the role of the galaxy group environment on galaxy evolution, we present a study of radio luminosity functions (RLFs) of group galaxies based on the Karl G. Jansky Very Large Array-COSMOS 3 GHz Large Project. The radio-selected sample of 7826 COSMOS galaxies with robust optical/near-infrared counterparts, excellent photometric coverage, and the COSMOS X-ray galaxy groups (M_200c > 10^13.3 M_0) enables us to construct the RLF of group galaxies (GGs) and their contribution to the total RLF since z ~ 2.3. Using the Markov chain Monte Carlo algorithm, we fit a redshift-dependent pure luminosity evolution model and a linear and power-law model to the luminosity functions. We compare it with past RLF studies from VLA-COSMOS on individual populations of radio-selected star-forming galaxies (SFGs) and galaxies hosting active galactic nuclei (AGN). These populations are classified based on the presence or absence of a radio excess concerning the star-formation rates derived from the infrared emission. We find that the fraction of radio group galaxies evolves by a factor of ~ 3 from z ~ 2 to the present day. The increase in the galaxy group contribution is due to the radio activity in groups being nearly constant at z < 1, while it is declining in the field. We show that massive galaxies inside galaxy groups remain radio active below redshift 1, contrary to the ones in the field. This evolution in the GG RLF is driven mainly by satellite galaxies in groups. Group galaxies associated with SFGs dominate the GG RLF at z_med = 0.3, while at z_med = 0.8, the peak in the RLF, coinciding with a known overdensity in COSMOS, is mainly driven by AGN. The study provides an observational probe for the accuracy of the numerical predictions of the radio emission in galaxies in a group environment.

Context. The central molecular zone at the Galactic center is currently being studied intensively to understand how star formation proceeds under the extreme conditions of a galactic nucleus. Knowing the position of molecular clouds along the line of sight toward the Galactic center has had important implications in our understanding of the physics of the gas and star formation in the central molecular zone.It was recently claimed that the dense molecular cloud G0.253 + 0.016 (the Brick) has a distance of $\sim$7.20 kpc from the Sun. That would place it outside of the central molecular zone, and therefore of the nuclear stellar disk, but still inside the Bulge. Aims. Theoretical considerations as well as observational studies show that stars that belong to the nuclear stellar disk have different kinematics from those that belong to the inner Bulge. Therefore, we aim to constrain the distance to the Brick by studying the proper motions of the stars in the area. Results. The stellar population seen toward the nuclear stellar disk shows the following three kinematic components: 1) Bulge stars with an isotropic velocity dispersion of $\sim$3.5 micro-arc second per year; 2) eastward moving stars on the near side of the nuclear stellar disk; and 3) westward moving stars on the far side of the nuclear stellar disk. We clearly see all three components toward the comparison field. However, toward the Brick, which blocks the light from stars behind it, we can only see kinematic components 1) and 2). Conclusions. While the Brick blocks the light from stars on the far side of the nuclear stellar disk, the detection of a significant component of eastward streaming stars implies that the Brick must be located inside the nuclear stellar disk and, therefore, that it forms part of the central molecular zone.

We present the global solutions of low angular momentum, inviscid, advective accretion flow around Kerr-Taub-NUT (KTN) black hole in presence and absence of shock waves. These solutions are obtained by solving the governing equations that describe the relativistic accretion flow in KTN spacetime which is characterized by the Kerr parameter ($a_{\rm k}$) and NUT parameter ($n$). During accretion, rotating flow experiences centrifugal barrier that eventually triggers the discontinuous shock transition provided the relativistic shock conditions are satisfied. In reality, the viability of shocked accretion solution appears more generic over the shock free solution as the former possesses high entropy content at the inner edge of the disc. Due to shock compression, the post-shock flow (equivalently post-shock corona, hereafter PSC) becomes hot and dense, and therefore, can produce high energy radiations after reprocessing the soft photons from the pre-shock flow via inverse Comptonization. In general, PSC is characterized by the shock properties, namely shock location ($r_s$), compression ratio ($R$) and shock strength ($S$), and we examine their dependencies on the energy (${\cal E}$) and angular momentum ($\lambda$) of the flow as well as black hole parameters. We identify the effective domain of the parameter space in $\lambda-{\cal E}$ plane for shock and observe that shock continues to form for wide range of flow parameters. We also find that $a_{\rm k}$ and $n$ act oppositely in determining the shock properties and shock parameter space. Finally, we calculate the disc luminosity ($L$) considering free-free emissions and observe that accretion flows containing shocks are more luminous compared to the shock free solutions.

High-redshift radio sources provide plentiful opportunities for studying the formation and evolution of early galaxies and supermassive black holes. However, the number of known radio-loud active galactic nuclei (AGN) above redshift 4 is rather limited. At high redshifts, it appears that blazars, with relativistically beamed jets pointing towards the observer, are in majority compared to radio-loud sources with jets misaligned with respect to the line of sight. To find more of these misaligned AGN, milliarcsec-scale imaging studies carried out with very long baseline interferometry (VLBI) are needed, as they allow us to distinguish between compact core--jet radio sources and those with more extended emission. Previous high-resolution VLBI studies revealed that some of the radio sources among blazar candidates in fact show unbeamed radio emission on milliarcsecond scales. The most accurate optical coordinates determined with the Gaia astrometric space mission are also useful in the classification process. Here, we report on dual-frequency imaging observations of 13 high-redshift (4 < z < 4.5) quasars at 1.7 and 5 GHz with the European VLBI Network. This sample increases the number of z>4 radio sources for which VLBI observations are available by about a quarter. Using structural and physical properties, such as radio morphology, spectral index, variability, brightness temperature, as well as optical coordinates, we identified six blazars and six misaligned radio AGNs, with the remaining one tentatively identified as blazar.

To unravel the star formation process, we present a multi-scale and multi-wavelength study of the filamentary infrared dark cloud (IRDC) G333.73+0.37, which hosts previously known two HII regions located at its center. Each HII region is associated with a mid-infrared source, and is excited by a massive OB star. Two filamentary structures and a hub-filament system (HFS) associated with one HII region are investigated in absorption using the Spitzer 8.0 $\mu$m image. The $^{13}$CO(J = 2-1) and C$^{18}$O(J = 2-1) line data reveal two velocity components (around $-$35.5 and $-$33.5 km s$^{-1}$) toward the IRDC, favouring the presence of two filamentary clouds at different velocities. Nonthermal (or turbulent) motions are depicted in the IRDC using the C$^{18}$O line data. The spatial distribution of young stellar objects (YSOs) identified using the VVV near-infrared data traces star formation activities in the IRDC. Low-mass cores are identified toward both the HII regions using the ALMA 1.38 mm continuum map. The VLT/NACO adaptive-optics L$^{\prime}$-band images show the presence of at least three point-like sources and the absence of small-scale features in the inner 4000 AU around YSOs NIR31 and MIR 16 located toward the HII regions. The HII regions and groups of YSO are observed toward the central part of the IRDC, where the two filamentary clouds intersect. A scenario of cloud-cloud collision or converging flows in the IRDC seems to be applicable, which may explain star formation activities including HFS and massive stars.

The X-ray emission in BL Lac objects is believed to be dominated by synchrotron emission from their relativistic jets. However, when the jet emission is not strong, one could expect signatures of X-ray emission from inverse Compton scattering of accretion disc photons by the corona. Moreover, the observed X-ray variability can also originate in the disc, and gets propagated and amplified by the jet. Here, we present results on the BL Lac object Mrk 421 using the NuSTAR data acquired during 2017 when the source was in a moderate X-ray brightness state. For comparison with high jet activity, we also considered one epoch data in April 2013 during a very high X-ray brightness state. Our aim is to explore the possibility of the signature of accretion disc emission in the overall X-ray emission from Mrk 421. The spectral fitting of the data using the two component advective flow model shows gives (a) the size of the dynamic corona at the base of the jet from ~28 to 10 r$_s$, (b) the disc mass accretion rate from 0.021 to 0.051 $\dot M_{\rm Edd}$, (c) the halo mass accretion rate from 0.22 to 0.35 $\dot M_{\rm Edd}$, and (d) the viscosity parameter of the Keplerian accretion disc from 0.18 $-$ 0.25. In the assumed model, the total flux, disc and jet flux correlate with the radio flux observed during these epochs. We conclude that the spectra of all the epochs of Mrk 421 in 2017 are well described by the accretion disc based two component advective flow model. The estimated disc and jet flux relations with radio flux show that accretion disc can contribute to the observed X-ray emission, when X-ray data (that covers a small portion of the broad band spectral energy distribution of Mrk 421) is considered in isolation. However, the present disc based models are disfavoured with respect to the relativistic jet models when considering the X-ray data in conjunction with data at other wavelengths.

Changing look active galactic nuclei (CLAGNs) show complex nature in their X-ray spectral shape and line of sight column density variation. The physical mechanisms responsible for these variations are unclear. Here, we study the spectral properties of a CLAGN, NGC\,1365 using combined {\it XMM-Newton} and {\it NuSTAR} observations to understand the CL behavior. The model fitted mass accretion rate varied between $0.003\pm 0.001$ and $0.009\pm0.002$ $\dot M_{\rm Edd}$ and the dynamic corona changed from $28\pm 3$ to $10\pm1$ $r_g$. We found that the variable absorption column density correlates with the mass accretion rate and the geometry of the corona. The derived wind velocity was sufficiently low compared to the escape velocity to drive the wind away from the disc for the epochs when column densities were high. This suggests that the high and variable absorption can be due to failed winds from the disc. Our estimated ratio of mass outflow to inflow rate from the inner region of the disc lies between $0.019\pm0.006$ and $0.12\pm0.04$. From spectral fitting of the combined data, we found the mass of the central black hole to be constant $4.38\pm0.34 - 4.51\pm0.29 \times10^{6} M_\odot$, consistent with earlier findings. The confidence contours of $N_H$ with other model parameters show that the model fitted parameters are robust and non-degenerate. Our study construed that the changing accretion rate, which is a fundamental physical quantity and the geometry of the corona driving the CL phenomena in NGC\,1365. The physical picture considered in this work connects both variable continuum and variable absorbing medium scenarios.

Several studies have been presented in the last few years applying some kind of automatic processing of data to estimate the fundamental parameters of open clusters. These parameters are later on employed in larger scale analyses, for example the structure of the Galaxy's spiral arms. The distance is one of the more straightforward parameters to estimate, yet enormous differences can still be found among published data. This is particularly true for open clusters located more than a few kpc away. We cross-matched several published catalogues and selected the twenty-five most distant open clusters ($>$9000 pc). We then performed a detailed analysis of their fundamental parameters, with emphasis on their distances, to determine the agreement between catalogues and our estimates.} Photometric and astrometric data from the Gaia EDR3 survey was employed. The data was processed with our own membership analysis code (pyUPMASK), and our package for automatic fundamental cluster's parameters estimation ($\texttt{ASteCA}$). We find differences in the estimated distances of up to several kpc between our results and those catalogued, even for the catalogues that show the best matches with $\texttt{ASteCA}$ values. Large differences are also found for the age estimates. As a by-product of the analysis we find that vd Bergh-Hagen 176 could be the open cluster with the largest heliocentric distance catalogued to date. Caution is thus strongly recommended when using catalogued parameters of open clusters to infer large-scale properties of the Galaxy, particularly for those located more than a few kpc away.

The baryon acoustic oscillation (BAO) peak, seen in the cosmic matter distribution at redshifts up to ~3.5, reflects the continued expansion of the sonic horizon first identified in temperature anisotropies of the cosmic microwave background. The BAO peak position can now be measured to better than ~1% accuracy using galaxies, and ~1.4-1.6% precision with Ly-alpha forrests and the clustering of quasars. In conjunction with the Alcock-Paczy\'nski (AP) effect, which arises from the changing ratio of angular to spatial/redshift size of (presumed) spherically-symmetric source distributions with distance, the BAO measurement is viewed as one of the most powerful tools to use in assessing the geometry of the Universe. In this paper, we employ five BAO peak measurements from the final release of the Sloan Digital Sky Survey IV, at average redshifts <z>=0.38, 0.51, 0.70, 1.48 and 2.33, to carry out a direct head-to-head comparison of the standard model, Lambda-CDM, and one of its principal competitors, known as the R_h=ct universe. For completeness, we complement the AP diagnostic with a volume-averaged distance probe that assumes a constant comoving distance scale r_d. Both probes are free of uncertain parameters, such as the Hubble constant, and are therefore ideally suited for this kind of model selection. We find that R_h=ct is favored by these measurements over the standard model based solely on the AP effect, with a likelihood ~75% versus ~25%, while Planck-Lambda-CDM is favored over R_h=ct based solely on the volume-averaged distance probe, with a likelihood ~80% versus ~20%. A joint analysis using both probes produces an inconclusive outcome, yielding comparable likelihoods to both models. We are therefore not able to confirm with this work that the BAO data, on their own, support an accelerating Universe.

Using high-resolution hydrodynamical simulations of galaxy clusters, we study the interaction between the brightest cluster galaxy, its supermassive black hole (BH) and the intracluster medium (ICM). We create initial conditions for which the ICM is in hydrostatic equilibrium within the gravitational potential from the galaxy and an NFW dark matter halo. Two free parameters associated with the thermodynamic profiles determine the cluster gas fraction and the central temperature, where the latter can be used to create cool-core or non-cool-core systems. Our simulations include radiative cooling, star formation, BH accretion, and stellar and active galactic nucleus (AGN) feedback. Even though the energy of AGN feedback is injected thermally and isotropically, it leads to anisotropic outflows and buoyantly rising bubbles. We find that the BH accretion rate (BHAR) is highly variable and only correlates strongly with the star formation rate (SFR) and the ICM when it is averaged over more than $1~\rm Myr$. We generally find good agreement with the theoretical precipitation framework. In $10^{13}~\rm M_\odot$ haloes, AGN feedback quenches the central galaxy and converts cool-core systems into non-cool-core systems. In contrast, higher-mass, cool-core clusters evolve cyclically. Episodes of high BHAR raise the entropy of the ICM out to the radius where the ratio of the cooling time and the local dynamical time $t_{\rm cool}/t_{\rm dyn} > 10$, thus suppressing condensation and, after a delay, the BHAR. The corresponding reduction in AGN feedback allows the ICM to cool and become unstable to precipitation, thus initiating a new episode of high SFR and BHAR.

We present weak lensing mass estimates for a sample of 458 galaxy clusters from the redMaPPer Sloan Digital Sky Survey DR8 catalogue using Hyper Suprime-Cam weak lensing data. We develop a method to quickly estimate cluster masses from weak lensing shear and use this method to estimate weak lensing masses for each of the galaxy clusters in our sample. Subsequently, we constrain the mass-richness relation as well as the intrinsic scatter between the cluster richness and the measured weak lensing masses. When calculating the mass-richness relation for all clusters with a richness $\lambda>20$, we find a tension in the slope of the mass-richness relation with the Dark Energy Survey Year 1 stacked weak lensing analysis. For a reduced sample of clusters with a richness $\lambda>40$, our results are consistent with the Dark Energy Survey Year 1 analysis.

NASA's Transiting Exoplanet Survey Satellite (TESS) is performing a homogeneous survey of the sky from space in search of transiting exoplanets. The collected data are also being used for detecting passing Solar system objects, including 17 new outer Solar system body candidates located at geocentric distances in the range 80-200 au, that need follow-up observations with ground-based telescope resources for confirmation. Here, we present results of a proof-of-concept mini-survey aimed at recovering two of these candidates that was carried out with the 4.2 m William Herschel Telescope and a QHY600L CMOS camera mounted at its prime focus. For each candidate attempted, we surveyed a patch of sky of over one square degree around its expected coordinates in Sloan r'. The same location was revisited in five consecutive or nearly consecutive nights, reaching S/N > 4 at r' < 23 mag. We focused on the areas of sky around the circumpolar TESS candidates located at equatorial coordinates (07:00:15, +86:55:19), 202.8 au from Earth, and (06:39:47, +83:43:54) at 162.1 au, but we could not recover either of them at r' < 23 mag. Based on the detailed analysis of the acquired images, we confirm that either both candidates are much fainter than predicted or that they are false positives.

We present deep optical and near-infrared photometry of UID 30901, a superluminous supernova (SLSN) discovered during the UltraVISTA survey. The observations were obtained with VIRCAM ($YJHK_{s}$) mounted on the VISTA telescope, DECam ($griz$) on the Blanco telescope, and SUBARU Hyper Suprime-Cam (HSC; $grizy$). These multi-band observations comprise +700 days making UID 30901 one of the best photometrically followed SLSNe to date. The host galaxy of UID 30901 is detected in a deep HST F814W image with an AB magnitude of $27.3 \pm 0.2$. While no spectra exist for the SN or its host galaxy, we perform our analysis assuming $z = 0.37$, based on the photometric redshift of a possible host galaxy found at a projected distance of 7 kpc. Fitting a blackbody to the observations, the radius, temperature, and bolometric light curve are computed. We find a maximum bolometric luminosity of $5.4 \pm 0.34 \times 10^{43}$ erg s$^{-1}$. A flattening in the light curve beyond 600 days is observed and several possible causes are discussed. We find the observations to clearly favour a SLSN type I, and plausible power sources such as the radioactive decay of $^{56}$Ni and the spin-down of a magnetar are compared to the data. We find that the magnetar model yields a good fit to the observations with the following parameters: a magnetic field $B = 1.4 \pm 0.3 \times 10^{14} \ G$, spin period of $P = 6.0 \pm 0.1 \ ms$ and ejecta mass $M_{ej} = 11.9^{+4.8}_{-6.4} M_{\odot}$.

We study the evolution and final dispersal of protoplanetary discs that evolve under the action of internal and external photoevaporation, and different degrees of viscous transport. We identify five distinct dispersal pathways, which are i) very long lived discs ($>20\,$Myr), ii) inside-out dispersal where internal photoevaporation dominates and opens inner holes, iii) outside-in dispersal where external photoevaporation dominates through disc truncation and two intermediate regimes characterised by lingering material in the inner disc with the outer disc dispersed predominantly by either internal or external photoevaporation. We determine how the lifetime, relative impact of internal and external winds and clearing pathway varies over a wide, plausible, parameter space of stellar/disc/radiation properties. There are a number of implications, for example in high UV environments because the outer disc lifetime is shorter than the time-scale for clearing the inner disc we do not expect transition discs to be common, which appears to be reflected in the location of transition disc populations towards the Orion Nebular Cluster. Irrespective of environment, we find that ongoing star formation is required to reproduce observed disc fractions as a function of stellar cluster age. This work demonstrates the importance of including both internal and external winds for understanding protoplanetary disc evolution.

Reverse shock (RS) emission can be used to probe the properties of the relativistic ejecta, especially the degree of magnetization $\sigma$, in gamma-ray burst (GRB) afterglows. However, there has been confusion in the literature regarding the physical condition for the RS formation, and the role of magnetic fields in the RS dynamics in the Poynting-flux-dominated regime is not fully understood. Exploiting the shock jump conditions, we characterize the properties of a magnetized RS. We compare the RS dynamics and forming conditions from different theories and numerical simulations, and reconcile the discrepancies among them. The strict RS forming condition is found to be $\sigma < \sigma_\mathrm{cr}=(8/3)\gamma_4^2(n_1/n_4)$, where $n_4$ and $n_1$ are the rest-frame number densities of the ejecta and the ambient medium, respectively, $\gamma_4$ is the bulk Lorentz factor, and $\sigma_\mathrm{cr}$ is the critical magnetization. Contrary to previous claims, we prove that this condition agrees with other theoretical and simulated results, which can be further applied to the setup and consistency check of future numerical simulations. Using this condition, we propose a characteristic radius for RS formation, and categorize the magnetized shell into three regimes: 'thick shell' (relativistic RS), 'thin shell' (trans-relativistic RS), and 'no RS' regimes. The critical magnetization $\sigma_\mathrm{cr}$ is generally below unity for thin shells, but can potentially reaches $\sim 100-1000$ in the 'thick shell' regime. Our results could be applied to the dynamical evolution of Poynting-flux-dominated ejecta, with potential applications to self-consistent lightcurve modelling of magnetized relativistic outflows.

We propose a phenomenological generalisation of the standard model with only one extra degree of freedom that parametrises the evolution of a scalar field responsible for the cosmic acceleration. The model also foresees an additional parameter in the form of a coupling between dark energy and dark matter. This model captures a large diversity of dark energy evolutions at low redshift and could usefully complement common CPL parametrisations widely used. In this context, we have been constraining the parametrisation with data from Planck and KiDS, bringing different results between the early and late universe observations.

While large numbers of supermassive black holes have been detected at z>6, their origin is still essentially unclear. Numerical simulations have shown that the conditions for the classical direct collapse scenario are very restrictive and fragmentation is very difficult to be avoided. We thus consider here a more general case of a dense massive protostar cluster at low metallicity (<~ 10^{-3} Z_solar) embedded in gas. We estimate the mass of the central massive object, formed via collisions and gas accretion, considering the extreme cases of a logarithmically flat and a Salpeter-type initial mass function. Objects with masses of at least 10^4 solar could be formed for inefficient radiative feedback, whereas ~10^3 solar mass objects could be formed when the accretion time is limited via feedback. These masses will vary depending on the environment and could be considerably larger, particularly due to the continuous infall of gas into the cloud. As a result, one may form intermediate mass black holes of ~ 10^4 solar masses or more. Upcoming observations with the James Webb Space Telescope (JWST) and other observatories may help to detect such massive black holes and their environment, thereby shedding additional light on such a formation channel.

We present ALMA data on the 3mm continuum emission, CO isotopologues (12CO, 13CO, C18O), and high-density molecular tracers (HCN, HCO+, HNC, HNCO, CS, CN, and CH3OH) in NGC4526. These data enable a detailed study of the physical properties of the molecular gas in a longtime resident of the Virgo Cluster; comparisons to more commonly-studied spiral galaxies offer intriguing hints into the processing of molecular gas in the cluster environment. Many molecular line ratios in NGC4526, along with our inferred abundances and CO/H2 conversion factors, are similar to those found in nearby spirals. One striking exception is the very low observed 12CO/13CO(1-0) line ratio, $3.4\pm0.3$, which is unusually low for spirals though not for Virgo Cluster early-type galaxies. We carry out radiative transfer modeling of the CO isotopologues with some archival (2-1) data, and we use Bayesian analysis with Markov chain Monte Carlo techniques to infer the physical properties of the CO-emitting gas. We find surprisingly low [12CO/13CO] abundance ratios of $7.8^{+2.7}_{-1.5}$ and $6.5^{+3.0}_{-1.3}$ at radii of 0.4 kpc and 1 kpc. The emission from the high-density tracers HCN, HCO+, HNC, CS and CN is also relatively bright, and CN is unusually optically thick in the inner parts of NGC4526. These features hint that processing in the cluster environment may have removed much of the galaxy's relatively diffuse, optically thinner molecular gas along with its atomic gas. Angular momentum transfer to the surrounding intracluster medium may also have caused contraction of the disk, magnifying radial gradients such as we find in [13CO/C18O]. More detailed chemical evolution modeling would be interesting in order to explore whether the unusual [12CO/13CO] abundance ratio is entirely an environmental effect or whether it also reflects the relatively old stellar population in this early-type galaxy.

We present a cosmic shear analysis with an improved redshift calibration for the fourth data release of the Kilo-Degree Survey (KiDS-1000) using self-organising maps (SOMs). Compared to the previous analysis of the KiDS-1000 data, we expand the redshift calibration sample to more than twice its size, now consisting of data of 17 spectroscopic redshift campaigns, and significantly extending the fraction of KiDS galaxies we are able to calibrate with our SOM redshift methodology. We then enhance the calibration sample with precision photometric redshifts from COSMOS2015 and the Physics of the Accelerated Universe Survey (PAUS), allowing us to fill gaps in the spectroscopic coverage of the KiDS data. Finally we perform a Complete Orthogonal Sets of E/B-Integrals (COSEBIs) cosmic shear analysis of the newly calibrated KiDS sample. We find $S_8 = 0.748_{-0.025}^{+0.021}$, which is in good agreement with previous KiDS studies and increases the tension with measurements of the cosmic microwave background to 3.4{\sigma}. We repeat the redshift calibration with different subsets of the full calibration sample and obtain, in all cases, agreement within at most 0.5{\sigma} in $S_8$ compared to our fiducial analysis. Including additional photometric redshifts allows us to calibrate an additional 6 % of the source galaxy sample. Even though further systematic testing with simulated data is necessary to quantify the impact of redshift outliers, precision photometric redshifts can be beneficial at high redshifts and to mitigate selection effects commonly found in spectroscopically selected calibration samples.

We present fBLS -- a novel fast-folding technique to search for transiting planets, based on the fast-folding algorithm (FFA), which is extensively used in pulsar astronomy. For a given lightcurve with $N$ data points, fBLS simultaneously produces all the binned phase-folded lightcurves for an array of $N_p$ trial periods. For each folded lightcurve produced by fBLS, the algorithm generates the standard BLS periodogram and statistics. The number of performed arithmetic operations is $\mathcal{O}\big(N_p\cdot\log N_p \big)$, while regular BLS requires $\mathcal{O}\big(N_p\cdot N\big)$ operations. fBLS can be used to detect small rocky transiting planets, with periods shorter than one day, a period range for which the computation is extensive. We demonstrate the capabilities of the new algorithm by performing a preliminary fBLS search for planets with ultra-short periods in the Kepler main-sequence lightcurves. In addition, we developed a simplistic signal validation scheme for vetting the planet candidates. This two-stage preliminary search identified all known ultra-short planet candidates and found three new ones.

We report a new black hole scalarization mechanism and disclose novel dynamical critical phenomena in the process of the nonlinear accretion of the scalar field into black holes. The accretion process can transform a seed black hole into a final scalarized or bald black hole, depending on the initial parameter of the scalar field $p$. There is a critical parameter $p_{\ast}$ and near it all intermediate solutions are attracted to a critical solution and stay there for a time scaling as $T\propto-\gamma\ln|p-p_{\ast}|$. At late times, the solutions evolve into scalarized black holes if $p>p_{\ast}$, or bald black holes if $p<p_{\ast}$. The final masses of the resulting scalarized/bald black holes satisfy power-laws $M_{p}-M_{\pm}\propto|p-p_{\ast}|^{\gamma_{\pm}}$ where $M_{\pm}$ are the masses of the scalarized/bald black holes when $p\to p_\ast$ from above/below, and $\gamma_{\pm}$ the corresponding exponents.

We consider the classically scale invariant Higgs-dilaton model of dynamical symmetry breaking extended with an extra scalar field that plays the role of dark matter. The Higgs boson is light near a critical boundary between different symmetry breaking phases, where quantum corrections beyond the usual Gildener-Weinberg approximation become relevant. This implies a tighter connection between dark matter and Higgs phenomenology. The model has only three free parameters, yet it allows for the observed relic abundance of dark matter while respecting all constraints. The direct detection cross section mediated by the Higgs boson is determined by the dark matter mass alone and is testable at future experiments.

We propose the infrared gauge symmetry of our sector includes an unbroken discrete gauged subgroup of baryon minus lepton number of order $2 \times 3 \text{ colors} \times 3 \text{ generations}$, the inclusion of which does not modify local physics. We UV complete this at a scale $\Lambda$ as the familiar $U(1)_{B-N_cL}$ Abelian Higgs theory, and the early universe phase transition forms cosmic strings which are charged under an emergent higher-form gauge symmetry. These topological defects catalyze interactions which turn $3$ baryons into $3$ leptons at strong scale rates in an analogue of the Callan-Rubakov effect. The cosmological lithium problem -- that the observed primordial abundance is lower than theoretical expectations by a factor of a few -- is perhaps the most statistically significant anomaly of $\text{SM}+\Lambda\text{CDM}$, and has resisted decades of attempts by cosmologists, nuclear physicists, and astronomers alike to root out systematics. We write down a model in which $B-N_cL$ strings superconduct bosonic global baryon plus lepton number currents and catalyze solely $3p^+ \rightarrow 3e^+$. We suggest that such cosmic strings have disintegrated $\mathcal{O}(1)$ of the lithium nuclei formed during Big Bang Nucleosynthesis and estimate the rate, with our benchmark model finding $\Lambda \sim 10^8 \text{ GeV}$ gives the right number density of strings.

We define compactness of a gravitational lens as the scaled closest distance of approach (i.e., $r_0/M$) of the null geodesic giving rise to an image. We model forty supermassive dark objects as Schwarzschild lenses and compute compactness of lenses (determined by the formation of the first order relativistic image). We then obtain a novel formula for the compactness of a lens as a function of mass to the distance ratio ($M/D_d$) and the ratio of lens-source to the observer-source distances ($D_{ds}/D_s$). This formula yields an interesting result: Just an observation of a relativistic image would give an incredibly accurate upper bound to the compactness of the lens without having any knowledge of mass of the lens, angular source position, and observer-source and lens-source distances. Similarly, we show that the observation of the second order relativistic image would give a lower value of upper bound to the compactness. These results, though obtained for supermassive dark objects at galactic centers, are valid for any object compact enough to give rise to relativistic images.

We present the facilities of the Aragats Space Environmental Center in Armenia used during multi-year observations of the thunderstorm ground enhancements (TGEs) and corresponding environmental parameters. We analyze the characteristics of the scintillation detectors, operated on Aragats, and describe the coordinated detection of TGEs by the network of scintillation detectors, field meters, and environmental parameters. By using a fast synchronized data acquisition system we reveal correlations of the multivariate data on time scales from second to nanosecond which allow us to gain insight into the TGE and lightning origin and their interrelations.

We present a calculation of the radiative capture cross section $p(n,\gamma )d$ in the low-energy range, where the $M1$ reaction channel dominates. Employing the LENPIC nucleon-nucleon interaction up to the fifth order (N4LO) that is regularized by the semi-local coordinate space regulators, we obtain the initial and final state wave functions, and evaluate the phase shifts of the scattering state and deuteron properties. We derive the transition operator from the chiral effective field theory up to the next-to-next-to leading order (N2LO), where we also regularize the transition operator using regulators consistent with those of the interactions. We compute the capture cross sections and the results show a converging pattern with the chiral-order expansion of the nucleon-nucleon interaction, where the regulator dependence of the results is weak when higher-order nucleon-nucleon interactions are employed. We quantify the uncertainties of the cross-section results due to the chiral-order truncation. The chirally complete and consistent cross-section results are performed up to N2LO and they compare well with the experiments and other theoretical predictions.

We model the supermassive dark object $M87^*$ as a Schwarzschild lens and study the variations in tangential, radial, and total (the product of tangential and radial) magnifications of images (primary, secondary, and relativistic) against the changes in angular source position and the ratio of lens-source to the observer-source distance. Further, we study the behavior of partial derivatives (with respect to the angular source position) of total magnifications of images against the angular source position. Finally, we model supermassive dark objects at centers of forty galaxies as Schwarzschild lenses and study the variations in tangential, radial, and total magnifications of images against the change in the ratio of mass of the lens to its distance. These studies yield many non-intuitive results which are likely to be significant for ngEHT observations. We {\em hypothesize} that there exists a distortion parameter such that their signed sum of all images of singular gravitational lensing of a source identically vanishes. We test this with images of Schwarzschild lensing in weak and strong gravitational fields and find that this aesthetically appealing hypothesis succeeds with flying colors.

Universal fault-tolerant quantum computers, which promise to revolutionize computing, are currently limited by excessive noise in their constituent superconducting qubits. Determining the dominant sources of this excess noise will lead to a clearer understanding of how to mitigate it in future superconducting systems. Superconducting Nanowire Single-Photon Detectors (SNSPDs) are devices that do not appear to suffer from such effects and have extremely low dark-count backgrounds. We propose to use SNSPDs as low-background laboratories to study noise accumulation processes in superconducting systems with the purpose of explaining and mitigating noise in related quantum information systems. Through these studies we also aim to increase the sensitive wavelengths of SNSPDs above the current limits of 10 microns, which would open new regimes for dark matter detection, biology, space sciences, and quantum sensing.

In the minimal supersymmetric standard model (MSSM) extended by a singlet superfield, when the coupling between the singlet sector and the MSSM sector is tiny, the singlet sector can be a quasi dark sector with supersymmetry (SUSY). We investigate the cosmological phenomena in this scenario and obtain the following observations: (i) In the parameter space solving the small cosmological scale anomalies via self-interacting singlino dark matter (SIDM), a first-order phase transition (FOPT) can readily happen but requires rather light dark matter below MeV; (ii) The corresponding parameter space indicated by FOPT and SIDM can be partially covered by detecting the phase-transition gravitational waves (GWs) at the near-future projects, such as LISA, TianQin and Taiji. Therefore, the recently developed GW astronomy could be a novel probe to such a SUSY scenario.

The Korea-IBS-Daegu-SKKU energy density functional (KIDS-EDF) models, derived from the universal Skyrme functional, have been successfully and widely applied in describing the properties of finite nuclei and infinite nuclear matter. In the present work, we extend the applications of the KIDS-EDF models to investigate the implications of the nucleon effective mass and nuclear symmetry energy obtained from the KIDS-EDF models on the properties of neutron star (NS) and neutrino interaction with the NS constituents matter in the linear response approximation (LRA). We then analyze the total differential cross-section of neutrino, neutrino mean free path (NMFP), and the NS mass-radius (M-R) relations. We find that the NS M-R relations predictions for all KIDS-EDF models are in excellent agreement with the recent observations as well as the NICER result. Remarkable prediction results on the NMFPs are given by the KIDS0-m*77 and KIDS0-m*99 models with $M_n^* /M \lesssim 1$ which are quite higher in comparison with those obtained for the KIDS0, KIDS-A, and KIDS-B models with $M_n^*/M \gtrsim 1$. For the KIDS0, KIDS-A, and KIDS-B models, we obtain the $\lambda \lesssim R_{\textrm{NS}}$, indicating that these models support the slow NS cooling and neutrino trapping in NS. On the contrary, both KIDS0-m*77 and KIDS0-m*99 models support faster NS cooling and a small possibility of neutrino trapping within NS, predicting $\lambda \gtrsim R_{\textrm{NS}}$. More interestingly the NMFP decreases as the density and neutrino energy increase, which is consistent with those obtained in the Brussels-Montreal Skyrme (BSk17 and BSk18) models at saturation density.

Trajectories and radiation of the accelerating electrons are studied in the pulsar magnetosphere approximated as the electromagnetic field of the Deutsch's solutions. Because the electrons are accelerated rapidly to ultra-relativistic velocity near the neutron star surface, the electron velocity vector (and then its trajectory) is derived from the balance between Lorentz force and radiation reaction force, which makes the pitch angle between electron trajectories and magnetic field lines nonzero in most part of the magnetosphere. In such a case, the spectral energy distributions (SEDs) of synchro-curvature radiation for the accelerating electrons with a mono-energetic form are calculated. Our results indicate that: (i) the pitch angle is the function of electron position ($r, \theta, \phi$) in the open field line regions, and increases with increasing $r$ and $\theta$ as well as increasing the inclination angle; (ii) the radius of curvature becomes large along the particle trajectory, and (iii) the SED appears a double peak structure depending on the emission position, where the synchrotron radiation plays an important role in X-ray band and curvature radiation mainly works in GeV band, which is only determined by parameters $\alpha$ and $\zeta$

In this paper, we consider the hypothesis that fractions of dark matter could be constituted by primordial black holes (PBHs). To test this possibility, we work out the observational properties of a static black hole embedded in the dark matter envelope made of a PBH source. The corresponding modifications of geometry due to such a physical system are investigated, with particular focus on the accretion disk luminosity in spiral galaxies. The impact of the PBH presence is analyzed through modification of the disk luminosity and kinematic quantities. Thus, we discuss possible constraints on the PBH abundance in view of the most recent theoretical bounds. The results of our study indicate that suitable PBH masses are $M_\text{PBH}\in[10^6,10^{12}]M_\odot$ for PBH fractions $f_\text{PBH}\in[10^{-3},1]$. In particular, a comparison with the predictions of the exponential sphere density profile for dark matter suggests that the best-matching configuration is achieved for $f_\text{PBH}=1$ and $M_\text{PBH}=10^6 M_\odot$. Consequences with respect to the current knowledge on primordial black hole physics are discussed.

We use neutron star mass and radius measurements to constrain the spontaneous scalarization phenomenon in scalar-tensor theories using Bayesian analysis. Neutron star structures in this scenario can be significantly different from the case of general relativity, which can be used to constrain the theory parameters. We utilize this idea to obtain lower bounds on the coupling parameter $\beta$ for the case of massless scalars. These constraints are currently weaker than the ones coming from binary observations, and they have relatively low precision due to the approximations in our method. Nevertheless, our results clearly demonstrate the power of the mass-radius data in testing gravity, and can be further improved with future observations. The picture is different for massive scalars, for which the same data is not able to effectively constrain the theory parameters in an unexpected manner. We identify the main reason for this to be a large high-likelihood region in the parameter space where deviations from general relativity are relatively small. We hope this initial study to be an invitation to use neutron star structure measurements more effectively to test alternative theories in general.

We heal the cosmological constant problem by means of a \emph{cancellation mechanism} that adopts a phase transition during which quantum fluctuations are eliminated. To this purpose, we propose that a generalized scalar (dark) matter field with a non-vanishing pressure term can remove the vacuum energy contribution, if its corresponding thermodynamics is written in terms of a \emph{quasi-quintessence} representation. In such a picture, pressure differs from quintessence as it shows a zero kinetic contribution. Using this field, we investigate a metastable transition phase, in which the universe naturally passes through an inflationary phase. To reach this target, we single out a double exponential potential, describing the metastable inflationary dynamics by considering suitable boundaries and thermodynamic conditions. We analyze stability investigating saddle, stable and unstable points and we thus predict a chaotic inflation that mimics the Starobinsky exponential potential. Consequently, the role of the proposed dark matter field is investigated throughout the overall universe evolution. To do so, we provide a physical explanation on unifying the dark sector with inflation by healing the cosmological constant problem.

The production of light nuclei and hypernuclei together with heavy baryons, both hyperons and $\Delta$-baryons, in low density matter as found in stellar environments such as supernova or binary mergers is studied within relativistic mean-field models. Five light nuclei were considered together with three light hypernuclei. The presence of both hyperons and $\Delta$-baryons shift the dissolution of clusters to larger densities and increase the abundance of clusters. This effect is larger the smaller the charge fraction and the higher the temperature. The couplings of the $\Delta$-baryons were chosen imposing that the nucleon effective mass remains finite inside neutron stars.

Spherical harmonic modes of gravitational waveforms for inspiraling compact binaries in eccentric orbits from post-Newtonian (PN) theory accurate to third post-Newtonian order, and those extracted from Numerical Relativity (NR) simulations for binary black holes (BBHs) are compared. We combine results from the two approaches (PN and NR) to construct time-domain hybrid waveforms that describe the complete evolution of BBH mergers through inspiral(I), merger(M) and the ringdown(R) stages. These hybrids are then used in constructing a fully analytical dominant mode (l=2, |m|=2) eccentric IMR model. Overlaps with quasi-circular IMR waveform models including the effect of higher modes, maximized over a time- and phase-shift, hint at the importance (mismatches > 1%) of including eccentricity in gravitational waveforms when analysing BBHs lighter than ~ 80 M_sun, irrespective of the binary's eccentricity (as it enters the LIGO bands), or mass-ratio. Combined impact of eccentricity and higher modes seems to become more apparent through smaller overlaps with increasing inclination angles and mass ratios. Finally, we show that the state-of-the-art quasi-circular models including the effect of higher modes will not be adequate in extracting source properties for signals with initial eccentricities e0 >= 0.1.