Papers for Tuesday, Jan 26 2021

Dynamically formed black hole (BH) binaries (BBHs) are important sources of gravitational waves (GWs). Globular clusters (GCs) are one of the major environments to produce such BBHs, but the total mass of the known GCs is small compared to that in the Galaxy, thus the fraction of BBHs formed in GCs is also small. However, this assumes that GCs contain a canonical initial mass function like that in the field stars. This might not be true because several studies suggest that extreme dense and metal-poor environment can result in top-heavy IMFs, where GCs may come from. Although GCs with top-heavy IMFs were easily disrupted or have become dark clusters, the contribution to the GW sources can be significant. Using a high-performance and accurate $N$-body code, \textsc{petar}, we investigate the effect of varying IMFs by carrying out four star-by-star simulations of dense GCs with the initial mass of $5\times10^5 M_\odot$ and the half-mass radius of $2$~pc. We find that the BBH merger rate does not monotonically correlate with the slope of IMFs. Due to a rapid expansion, top-heavy IMFs lead to less efficient formation of merging BBHs. The formation rate continuously decreases as the cluster expands due to the dynamical heating caused by BHs. However, in star clusters with a top-heavier IMF, the total number of BHs is larger, and therefore the final contribution to merging BBHs can still be more than clusters with the standard IMF, if the initial cluster mass and density is higher than the model we used.

The mean and the scatter of the HI content of a dark-matter halo as a function of the halo mass are useful statistics that can be used to test models of structure and galaxy formation. We investigate the possibility of constraining this HI-halo mass relation (HIHMR) from intensity maps of the redshifted 21 cm line. In particular, we use the geometry and topology of the brightness-temperature isocontours in a single frequency channel as quantified by the Minkowski functionals. First, we generate mock maps from a large N-body simulation considering the impact of thermal noise and foreground removal. We then use the Fisher information formalism to forecast constraints on a parametric model for the HIHMR. We consider a 20,000 deg$^2$ survey (originally proposed for dark-energy science) conducted with the Square Kilometre Array Phase 1 (SKA-1) MID observatory operating in single-dish mode. For a channel bandwidth of 2 MHz, we show that an integration time of a few$\,\times\,10^4$ s per pointing is sufficient to image the smoothed HI distribution at redshift $z \simeq 1$ and to measure the HIHMR in a nearly optimal way from the Minkowski functionals. Tighter constraints on some of the parameters can be obtained by using also an independent measurement of the mean HI density. Combining the results from different frequency channels provides exquisite constraints on the evolution of the HIHMR, especially in the central frequency range of the data cube.

We present analytic estimates of the fractional uncertainties on the mass, radius, surface gravity, and density of a transiting planet, using only empirical or semi-empirical measurements. We first express these parameters in terms of transit photometry and radial velocity (RV) observables, as well as the stellar radius $R_{\star}$, if required. In agreement with previous results, we find that, assuming a circular orbit, the surface gravity of the planet ($g_p$) depends only on empirical transit and RV parameters; namely, the planet period $P$, the transit depth $\delta$, the RV semi-amplitude $K_{\star}$, the transit duration $T$, and the ingress/egress duration $\tau$. However, the planet mass and density depend on all these quantities, plus $R_{\star}$. Thus, an inference about the planet mass, radius, and density must rely upon an external constraint such as the stellar radius. For bright stars, stellar radii can now be measured nearly empirically by using measurements of the stellar bolometric flux, the effective temperature, and the distance to the star via its parallax, with the extinction $A_V$ being the only free parameter. For any given system, there is a hierarchy of achievable precisions on the planetary parameters, such that the planetary surface gravity is more accurately measured than the density, which in turn is more accurately measured than the mass. We find that surface gravity provides a strong constraint on the core mass fraction of terrestrial planets. This is useful, given that the surface gravity may be one of the best measured properties of a terrestrial planet.

If a single line of sight (LOS) intercepts multiple dust clouds of different spectral energy distributions and magnetic field orientations, the frequency scaling of each of the Stokes $Q$ and $U$ parameters of thermal dust emission may be different ("LOS frequency decorrelation"). We present first evidence for LOS frequency decorrelation in $Planck$ data. We use independent, neutral-hydrogen--measurements of the number of clouds per LOS and the magnetic field orientation in each cloud to select two sets of sightlines: (i) a target sample (pixels likely to exhibit LOS frequency decorrelation); (ii) a control sample (pixels lacking complex LOS structure). We test the null hypothesis that LOS frequency decorrelation is not detectable in $Planck$ 353 and 217~GHz polarization data at high Galactic latitudes. The data reject this hypothesis at high significance. The detection is robust against choice of CMB map and map-making pipeline. The observed change in polarization angle due to LOS frequency decorrelation is detectable above the $Planck$ noise level. The probability that the detected effect is due to noise alone ranges from $5\times 10^{-2}$ to $4\times 10^{-7}$, depending on the CMB subtraction algorithm and treatment of residual systematics; correcting for residual systematics increases the significance of the effect. The LOS decorrelation effect is stronger for sightlines with more misaligned magnetic fields, as expected. We estimate that an intrinsic variation of $\sim15\%$ in the ratio of 353 to 217~GHz polarized emission between clouds is sufficient to reproduce the measured effect. Our finding underlines the importance of ongoing studies to map the 3D structure of the magnetized dusty ISM that could help component separation methods to account for frequency decorrelation effects in CMB polarization studies.

Magnetic fields are ubiquitous in galaxy clusters, yet their radial profile, power spectrum, and connection to host cluster properties are poorly known. Merging galaxy clusters hosting diffuse polarized emission in the form of radio relics offer a unique possibility to study the magnetic fields in these complex systems. In this paper, we investigate the intra-cluster magnetic field in Abell 2345. This cluster hosts two radio relics that we detected in polarization with 1-2 GHz JVLA observations. X-ray XMM-Newton images show a very disturbed morphology. We derived the Rotation Measure (RM) of five polarized sources within $\sim$ 1 Mpc from the cluster center applying the RM synthesis. Both, the average RM and the RM dispersion radial profiles probe the presence of intra-cluster magnetic fields. Using the thermal electron density profile derived from X-ray analysis and simulating a 3D magnetic field with fluctuations following a power spectrum derived from magneto-hydrodynamical cosmological simulations, we build mock RM images of the cluster. We constrained the magnetic field profile in the eastern radio relic sector by comparing simulated and observed RM images. We find that, within the framework of our model, the data require a magnetic field scaling with thermal electron density as $B(r)\propto n_e(r)$. The best model has a central magnetic field (within a 200 kpc radius) of $2.8\pm0.1$ $\mu$G. The average magnetic field at the position of the eastern relic is $\sim$0.3 $\mu$G, a factor 2.7 lower than the equipartition estimate.

Diffuse filaments connect galaxy clusters to form the cosmic web. Detecting these filaments could yield information on the magnetic field strength, cosmic ray population and temperature of intercluster gas, yet, the faint and large-scale nature of these bridges makes direct detections very challenging. Using multiple independent all-sky radio and X-ray maps we stack pairs of luminous red galaxies as tracers for cluster pairs. For the first time, we detect an average surface brightness between the clusters from synchrotron (radio) and thermal (X-ray) emission with $\gtrsim 5\sigma$ significance, on physical scales larger than observed to date ($\geq 3\,$Mpc). We obtain a synchrotron spectral index of $\alpha \simeq -1.0$ and estimates of the average magnetic field strength of $ 30 \leq B \leq 60 \,$nG, derived from both equipartition and Inverse Compton arguments, implying a 5 to 15$\,$per cent degree of field regularity when compared with Faraday rotation measure estimates. While the X-ray detection is inline with predictions, the average radio signal comes out higher than predicted by cosmological simulations and dark matter annihilation and decay models. This discovery demonstrates that there are connective structures between mass concentrations that are significantly magnetised, and the presence of sufficient cosmic rays to produce detectable synchrotron radiation.

We present a theoretical calculation of the influence of ultraviolet radiative pumping on the excitation of the rotational levels of the ground vibrational state for HD molecules under conditions of the cold diffuse interstellar medium (ISM). Two main excitation mechanisms have been taken into account in our analysis: (i) collisions with atoms and molecules and (ii) radiative pumping by the interstellar ultraviolet (UV) radiation field. The calculation of the radiative pumping rate coefficients $\Gamma_{\rm ij}$ corresponding to Drane's model of the field of interstellar UV radiation, taking into account the self-shielding of HD molecules, is performed. We found that the population of the first HD rotational level ($J = 1$) is determined mainly by radiative pumping rather than by collisions if the thermal gas pressure $p_{\rm th}\le10^4\left(\frac{I_{\rm{UV}}}{1}\right)\,\mbox{K\,cm}^{-3}$ and the column density of HD is lower than $\log N({\rm{HD}})<15$. Under this constraint the populations of rotational levels of HD turns out to be as well a more sensitive indicator of the UV radiation intensity than the fine-structure levels of atomic carbon. We suggest that taking into account radiative pumping of HD rotational levels may be important for the problem of the cooling of primordial gas at high redshift: ultraviolet radiation from first stars can increase the rate of HD cooling of the primordial gas in the early Universe.

We measured the Rossiter-McLaughlin effect of WASP-107b during a single transit with Keck/HIRES. We found the sky-projected inclination of WASP-107b's orbit, relative to its host star's rotation axis, to be $|\lambda| = {118}^{+38}_{-19}$ degrees. This confirms the misaligned/polar orbit that was previously suggested from spot-crossing events and adds WASP-107b to the growing population of hot Neptunes in polar orbits around cool stars. WASP-107b is also the fourth such planet to have a known distant planetary companion. We examined several dynamical pathways by which this companion could have induced such an obliquity in WASP-107b. We find that nodal precession and disk dispersal-driven tilting can both explain the current orbital geometry while Kozai-Lidov cycles are suppressed by general relativity. While each hypothesis requires a mutual inclination between the two planets, nodal precession requires a much larger angle which for WASP-107 is on the threshold of detectability with future Gaia astrometric data. As nodal precession has no stellar type dependence, but disk dispersal-driven tilting does, distinguishing between these two models is best done on the population level. Finding and characterizing more extrasolar systems like WASP-107 will additionally help distinguish whether the distribution of hot-Neptune obliquities is a dichotomy of aligned and polar orbits or if we are uniformly sampling obliquities during nodal precession cycles.

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

We use a sample of 14 massive, dynamically relaxed galaxy clusters to constrain the Hubble Constant, $H_0$, by combining X-ray and Sunyaev-Zel'dovich (SZ) effect signals measured with Chandra, Planck and Bolocam. This is the first such analysis to marginalize over an empirical, data-driven prior on the overall accuracy of X-ray temperature measurements, while our restriction to the most relaxed, massive clusters also minimizes astrophysical systematics. For a cosmological-constant model with $\Omega_m = 0.3$ and $\Omega_{\Lambda} = 0.7$, we find $H_0 = 67.3^{+21.3}_{-13.3}$ km/s/Mpc, limited by the temperature calibration uncertainty (compared to the statistically limited constraint of $H_0 = 72.3^{+7.6}_{-7.6}$ km/s/Mpc). The intrinsic scatter in the X-ray/SZ pressure ratio is found to be $13 \pm 4$ per cent ($10 \pm 3$ per cent when two clusters with significant galactic dust emission are removed from the sample), consistent with being primarily due to triaxiality and projection. We discuss the prospects for reducing the dominant systematic limitation to this analysis, with improved X-ray calibration and/or precise measurements of the relativistic SZ effect providing a plausible route to per cent level constraints on $H_0$.

Characterising the atmospheres of exoplanets is key to understanding their nature and provides hints about their formation and evolution. High-resolution measurements of the helium triplet, He(2$^{3}$S), absorption of highly irradiated planets have been recently reported, which provide a new mean to study their atmospheric escape. In this work, we study the escape of the upper atmospheres of HD 189733 b and GJ 3470 b by analysing high-resolution He(2$^{3}$S) absorption measurements and using a 1D hydrodynamic model coupled with a non-LTE model for the He(2$^{3}$S) state. We also use the H density derived from Ly$\alpha$ observations to further constrain their temperatures, T, mass-loss rates,$\dot M$, and H/He ratios. We have significantly improved our knowledge of the upper atmospheres of these planets. While HD 189733 b has a rather compressed atmosphere and small gas radial velocities, GJ 3470 b, with a gravitational potential ten times smaller, exhibits a very extended atmosphere and large radial outflow velocities. Hence, although GJ 3470 b is much less irradiated in the XUV, and its upper atmosphere is much cooler, it evaporates at a comparable rate. In particular, we find that the upper atmosphere of HD 189733 b is compact and hot, with a maximum T of 12400$^{+400}_{-300}$ K, with very low mean molecular mass (H/He=(99.2/0.8)$\pm0.1$), almost fully ionised above 1.1 R$_p$, and with $\dot M$=(1.1$\pm0.1$)$\times$10$^{11}$ g/s. In contrast, the upper atmosphere of GJ 3470 b is highly extended and relatively cold, with a maximum T of 5100$\pm900$ K, also with very low mean molecular mass (H/He=(98.5/1.5)$^{+1.0}_{-1.5}$), not strongly ionised and with $\dot M$=(1.9$\pm1.1$)$\times$10$^{11}$ g/s. Furthermore, our results suggest that the upper atmospheres of giant planets undergoing hydrodynamic escape tend to have very low mean molecular mass (H/He$\gtrsim$97/3).

We perform a series of high-resolution 2D hydrodynamical simulations of equal-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) accretion disks to study whether these binaries can be driven to merger by the surrounding gas. We find that the gravitational softening adopted for the BBH has a profound impact on this result. When the softening is less than ten percent of the binary separation, we show that, in agreement with recent simulations of isolated equal-mass binaries, prograde BBHs expand in time rather than contract. Eventually, however, the binary separation becomes large enough that the tidal force of the central AGN disrupts them. Only when the softening is relatively large do we find that prograde BBHs harden. We determine through detailed analysis of the binary torque, that this dichotomy is due to a loss of spiral structure in the circum-single disks orbiting each BH when the softening is a significant fraction of the binary separation. Properly resolving these spirals -- both with high resolution and small softening -- results in a significant source of binary angular momentum. Only for retrograde BBHs do we find consistent hardening, regardless of softening, as these BBHs lack the important spiral structure in their circum-single disks. This suggests that the gas-driven inspiral of retrograde binaries can produce a population of compact BBHs in the gravitational-wave-emitting regime in AGN disks, which may contribute a large fraction to the observed BBH mergers.

Using parsec scale resolution hydrodynamical adaptive mesh refinement simulations we have studied the mass transport process throughout a galactic merger. The aim of such study is to connect both the peaks of mass accretion rate onto the BHs and star formation bursts with both gravitational and hydrodynamic torques acting on the galactic gaseous component. Our merger initial conditions were chosen to mimic a realistic system. The simulations include gas cooling, star formation, supernovae feedback, and AGN feedback. Gravitational and hydrodynamic torques near pericenter passes trigger gas funneling to the nuclei which is associated with bursts of star formation and black hole growth. Such episodes are intimately related with both kinds of torques acting on the galactic gas. Pericenters trigger both star formation and mass accretion rates of $\sim$ few $(1-10)$ $M_\odot$/yr. Such episodes last $\sim$ $(50-75)$ Myrs. Close passes also can produce black hole accretion that approaches and reaches the Eddington rate, lasting $\sim$ few Myrs. Our simulation shows that both gravitational and hydrodynamic torques are enhanced at pericenter passes with gravitational torques tending to have higher values than the hydrodynamic torques throughout the merger. We also find that in the closest encounters, hydrodynamic and gravitational torques can be comparable in their effect on the gas, the two helping in the redistribution of both angular momentum and mass in the galactic disc. Such phenomena allow inward mass transport onto the BH influence radius, fueling the compact object and lighting up the galactic nuclei.

We present the detailed optical evolution of a type Ib SN 2015dj in NGC 7371, using data spanning up to $\sim$ 170 days after discovery. SN 2015dj shares similarity in light curve shape with SN 2007gr and peaks at M$_{V}$ = $-17.37\pm$0.02 mag. Analytical modelling of the quasi bolometric light curve yields 0.06$\pm$0.01 M$_{\odot}$ of $^{56}$Ni, ejecta mass $M_{\rm ej} = 1.4^{+1.3}_{-0.5}$ \msol\, and kinetic energy $E_{\rm k} = 0.7^{+0.6}_{-0.3} \times 10^{51}$ erg. The spectral features show a fast evolution and resemble those of spherically symmetric ejecta. The analysis of nebular phase spectral lines indicate a progenitor mass between 13-20 M$_{\odot}$ suggesting a binary scenario.

We present comprehensive characterization of the Galactic open cluster M 36. Some two hundred member candidates, with an estimated contamination rate of $\sim$8%, have been identified on the basis of proper motion and parallax measured by the $Gaia$ DR2. The cluster has a proper motion grouping around ($\mu_{\alpha} \cos\delta = -$0.15 $\pm$ 0.01 mas yr$^{-1}$, and $\mu_{\delta} = -$3.35 $\pm$ 0.02 mas yr$^{-1}$), distinctly separated from the field population. Most member candidates have parallax values 0.7$-$0.9 mas, with a median value of 0.82 $\pm$ 0.07 mas (distance $\sim$1.20 $\pm$ 0.13 kpc). The angular diameter of 27$'$ $\pm$ $0\farcm4$ determined from the radial density profile then corresponds to a linear extent of 9.42 $\pm$ 0.14 pc. With an estimated age of $\sim$15 Myr, M 36 is free of nebulosity. To the south-west of the cluster, we discover a highly obscured ($A_{V}$ up to $\sim$23 mag), compact ($\sim$ $1\farcm9 \times 1\farcm2$) dense cloud, within which three young stellar objects in their infancy (ages $\lesssim$ 0.2 Myr) are identified. The molecular gas, 3.6 pc in extent, contains a total mass of (2$-$3)$\times$10$^{2}$ M$_{\odot}$, and has a uniform velocity continuity across the cloud, with a velocity range of $-$20 to $-$22 km s$^{-1}$, consistent with the radial velocities of known star members. In addition, the cloud has a derived kinematic distance marginally in agreement with that of the star cluster. If physical association between M 36 and the young stellar population can be unambiguously established, this manifests a convincing example of prolonged star formation activity spanning up to tens of Myrs in molecular clouds.

The galactic diffuse $\gamma$-ray emission, as seen by Fermi Large Area Telescope (LAT), shows a sharp peak in the region around 4 kpc from the Galactic center, which can be interpreted either as due to an enhanced density of cosmic-ray accelerators or to a modification of the particle diffusion in that region. Observations of $\gamma$-rays originating in molecular clouds are a unique tool to infer the cosmic-ray density point by point, in distant regions of the Galaxy. We report here the analysis of 11 yr Fermi-LAT data, obtained in the direction of nine molecular clouds located in the 1.5--4.5 kpc region. The cosmic-ray density measured at the locations of these clouds is compatible with the locally measured one. We demonstrate that the cosmic-ray density gradient inferred from the diffuse gamma-ray emission is the result of the presence of cosmic-ray accelerators rather than a global change of the sea of Galactic cosmic rays due to their propagation.

We have derived the stellar atmospheric parameters, the effective temperature T$_{eff}$, the microturbulent velocity $\zeta$, the surface gravity log g, and the metallicity [Fe/H] for HE 0017+0055, HE 2144-1832, HE 2339-0837, HD 145777, and CD-27 14351 from local thermodynamic equilibrium analyses using model atmospheres. Elemental abundances of C, N, $\alpha$-elements, iron-peak elements, and several neutron-capture elements are estimated using the equivalent width measurement technique as well as spectrum synthesis calculations in some cases. In the context of the double enhancement observed in four of the programme stars, we have critically examined whether the literature i-process model yields ([X/Fe]) of heavy elements can explain the observed abundance distribution. The estimated metallicity [Fe/H] of the programme stars ranges from -1.63 to -2.74. All five stars show enhanced abundance for Ba, and four of them exhibit enhanced abundance for Eu. Based on our analysis, HE 0017+0055, HE 2144-1832, and HE 2339-0837 are found to be CEMP-r/s stars, whereas HD 145777 and CD-27 14351 show characteristic properties of CEMP-s stars. From a detailed analysis of different classifiers of CEMP stars, using a large sample of similar stars from the literature, we have identified the one which best describes the CEMP-s and CEMP-r/s stars. We have also examined if [hs/ls] alone can be used as a classifier, and if there are any limiting values for [hs/ls] ratio that can be used to distinguish CEMP-s and CEMP-r/s stars. In spite of peaking at different values of [hs/ls], CEMP-s and CEMP-r/s stars show an overlap in the range 0.0 < [hs/ls] < 1.5 and hence this ratio alone can not be used to distinguish CEMP-s and CEMP-r/s stars.

One of the most important early results from the Parker Solar Probe (PSP) is the ubiquitous presence of magnetic switchbacks, whose origin is under debate. Using a three-dimensional direct numerical simulation of the equations of compressible magnetohydrodynamics from the corona to 40 solar radii, we investigate whether magnetic switchbacks emerge from granulation-driven Alfv\'en waves and turbulence in the solar wind. The simulated solar wind is an Alfv\'enic slow-solar-wind stream with a radial profile consistent with various observations, including observations from PSP. As a natural consequence of Alfv\'en-wave turbulence, the simulation reproduced magnetic switchbacks with many of the same properties as observed switchbacks, including Alfv\'enic v-b correlation, spherical polarization (low magnetic compressibility), and a volume filling fraction that increases with radial distance. The analysis of propagation speed and scale length shows that the magnetic switchbacks are large-amplitude (nonlinear) Alfv\'en waves with discontinuities in the magnetic field direction. We directly compare our simulation with observations using a virtual flyby of PSP in our simulation domain. We conclude that at least some of the switchbacks observed by PSP are a natural consequence of the growth in amplitude of spherically polarized Alfv\'en waves as they propagate away from the Sun.

The porosity of an asteroid is important when studying the evolution of our solar system through small bodies and for planning mitigation strategies to avoid disasters due to asteroid impacts. Our knowledge of asteroid porosity largely relies on meteorites sampled on Earth. However, chondrites sampled on Earth are suggested to be sorted by strength. In this study, we obtained an estimate of the most porous structure of primordial "granular" chondrite parent bodies based on measurements of the compaction behavior of chondrite component analogs. We measured compaction curves of dust and dust-beads mixture samples. The dust sample consisted of various spherical and irregular particles with diameters in the order of 10^0-10^1 $\mu$m. The mixture sample consisted of dust and beads with different dust volume fractions (~0.2-1). We used 1.5 and 4.8 $\mu$m particles as dust as a first step although the typical size of materials in matrix may be much smaller. We approximated the compaction curve of each sample with a power-law form and calculated the porosity structure of the primordial chondrite parent bodies using the experimental results. Our results show that the primordial parent bodies are likely to have higher porosity than the chondrites. Moreover, the relatively higher volume fraction of the matrix may be one of the reasons why most meteorites with high porosity are carbonaceous chondrites.

IceCube has observed a flux of cosmic neutrinos, with a "bump" in the energy range $10 \lesssim E/{\rm TeV} \lesssim 100$ that creates a $3\sigma$ tension with gamma-ray data from the Fermi satellite. This has been interpreted as evidence for a population of hidden cosmic-ray accelerators. We propose an alternative explanation of this conundrum on the basis of cold dark matter which decays into sterile neutrinos that after oscillations produce the bump in the cosmic neutrino spectrum.

The scaling relations between the black hole (BH) mass and soft lag properties for both active galactic nuclei (AGN) and BH X-ray binaries (BHXRBs) suggest the same underlying physical mechanism at work in accreting BH systems spanning a broad range of mass. However, the low-mass end of AGN has never been explored in detail. In this work, we extend the existing scaling relations to lower-mass AGN, which serve as anchors between the normal-mass AGN and BHXRBs. For this purpose, we construct a sample of low-mass AGN ($M_{\rm BH}<3\times 10^{6} M_{\rm \odot}$) from the XMM-Newton archive and measure frequency-resolved time delays between the soft (0.3-1 keV) and hard (1-4 keV) X-ray emission. We report that the soft band lags behind the hard band emission at high frequencies $\sim[1.3-2.6]\times 10^{-3}$ Hz, which is interpreted as a sign of reverberation from the inner accretion disc in response to the direct coronal emission. At low frequencies ($\sim[3-8]\times 10^{-4}$ Hz), the hard band lags behind the soft band variations, which we explain in the context of the inward propagation of luminosity fluctuations through the corona. We find that the X-ray source for the sample extends at an average radius of around $8r_{\rm g}$ and a median height of around $15r_{\rm g}$ above the disc plane, assuming a lamppost geometry for the corona. Our results confirm that the scaling relations between the BH mass and soft lag amplitude/frequency derived for higher-mass AGN can safely extrapolate to lower-mass AGN, and the accretion process is indeed independent of the BH mass.

Ali CMB Polarization Telescope (AliCPT-1) is the first CMB degree-scale polarimeter to be deployed on the Tibetan plateau at 5,250m above sea level. AliCPT-1 is a 90/150 GHz 72 cm aperture, two-lens refracting telescope cooled down to 4 K. Alumina lenses, 800mm in diameter, image the CMB in a 33.4{\deg} field of view on a 636mm wide focal plane. The modularized focal plane consists of dichroic polarization-sensitive Transition-Edge Sensors (TESes). Each module includes 1,704 optically active TESes fabricated on a 150mm diameter silicon wafer. Each TES array is read out with a microwave multiplexing readout system capable of a multiplexing factor up to 2,048. Such a large multiplexing factor has allowed the practical deployment of tens of thousands of detectors, enabling the design of a receiver that can operate up to 19 TES arrays for a total of 32,376 TESes. AliCPT-1 leverages the technological advancements in the detector design from multiple generations of previously successful feedhorn-coupled polarimeters, and in the instrument design from BICEP-3, but applied on a larger scale. The cryostat receiver is currently under integration and testing. During the first deployment year, the focal plane will be populated with up to 4 TES arrays. Further TES arrays will be deployed in the following years, fully populating the focal plane with 19 arrays on the fourth deployment year. Here we present the AliCPT-1 receiver design, and how the design has been optimized to meet the experimental requirements.

The prototype of the Large-Sized Telescope (LST) of the Cherenkov Telescope Array (CTA) is presently going through its commissioning phase. A total of four LSTs, among others, will operate together at Observatorio del Roque de Los Muchachos, which will host the CTA North site. A computing center endowed with 1760 cores and several petabytes disk space is installed onsite. It is used to acquire, process, and analyze the data produced, at a rate of 3~TB/hour during operation. The LST On-site Analysis LSTOSA is a set of scripts written in Python which connects the different steps of lstchain, the analysis pipeline developed for the LST. It processes the data in a semiautomatic way producing high-level data and quality plots including detailed provenance logs. Data are analyzed before the next observation night to help in the commissioning procedure and debugging.

We present an independent examination of the parallax zero-point of the Third Gaia Early Data Release (hereafter EDR3), using the LAMOST primary red clump (PRC) stellar sample. A median parallax offset of around $26 \mu$as, slightly larger than that found by examination of distant quasars, is found for both the five- and six-parameter solutions in EDR3, based on samples of over 66,000 and 2000 PRC stars, respectively. Similar to the previous investigation of Lindegren et al., to which we compare our results, the parallax zero-point exhibits clear dependencies on the $G$ magnitudes, colors, and positions of the objects. Based on our analysis, the zero-point of the revised parallax can be reduced to a few $\mu$as, and some significant patterns, e.g., discontinuities with stellar magnitude, can be properly removed. However, relatively large offsets ($> 10 \mu$as) are still found for the revised parallaxes over different positions on the sky.

We present results of spectroscopic monitoring observations of the Ultra-Luminous Infra Red Galaxy F01004-2237. This galaxy was observed to undergo changes in its optical spectrum, detected by comparing a spectrum from 2015 with one from 2000. These changes were coincident with photometric brightening. The main changes detected in the optical spectrum are enhanced He II $\lambda$4686 emission and the appearance of He I $\lambda$3898,$\lambda$5876 emission lines. The favoured interpretation of these changes was that of a tidal disruption event (TDE) happening in 2010. However, subsequent work suggested that these changes are caused by another hitherto unknown reason related to variations in the accretion rate in the active galactic nucleus (AGN). Our optical spectroscopic monitoring observations show that the evolution of the He lines is in line with the evolution seen in TDEs and opposite of what observed from reverberation mapping studies of AGNs, renewing the discussion on the interpretation of the flare as a TDE.

We present \verb|piXedfit|, pixelized spectral energy distribution (SED) fitting, a Python package that provides tools for analyzing spatially resolved properties of galaxies using multiband imaging data alone or in combination with integral field spectroscopy (IFS) data. \verb|piXedfit| has six modules that can handle all tasks in the spatially resolved SED fitting. The SED fitting module uses the Bayesian inference technique with two kinds of posteriors sampling methods: Markov Chain Monte Carlo (MCMC) and random densely-sampling of parameter space (RDSPS). We test the performance of the SED fitting module using mock SEDs of simulated galaxies from IllustrisTNG. The SED fitting with both posteriors sampling methods can recover physical properties and star formation histories of the IllustrisTNG galaxies well. We further test the performance of \verb|piXedfit| modules by analyzing 20 galaxies observed by the CALIFA and MaNGA surveys. The data comprises of 12-band imaging data from GALEX, SDSS, 2MASS, and WISE, and the IFS data from CALIFA or MaNGA. \verb|piXedfit| can spatially match (in resolution and sampling) of the imaging and IFS data. By fitting only the photometric SEDs, \verb|piXedfit| can predict the spectral continuum, $\text{D}_{\rm n}4000$, $H_{\alpha}$, and $H_{\beta}$ well. The star formation rate (SFR) derived by \verb|piXedfit| is consistent with that derived from $H_{\alpha}$ emission. The RDSPS method gives equally good fitting results as the MCMC and it is much faster than the MCMC. \verb|piXedfit| is a versatile tool equipped with a parallel computing module for efficient analysis of large datasets, and will be made publicly available (https://github.com/aabdurrouf/piXedfit).

The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a $8.31$\,deg$^2$ sub-region of the GAMA~23~field which comprises of surveys covering the frequency range 0.1--9\,GHz. We present a sample of 104 radio galaxies compiled from observations conducted by the Murchison Wide-field Array (216\,MHz), the Australia Square Kilometer Array Pathfinder (887\,MHz), and the Australia Telescope Compact Array (5.5\,GHz). We adopt an `absent radio core' criterion to identify 10 radio galaxies showing no evidence for an active nucleus. We classify these as new candidate remnant radio galaxies. Seven of these objects still display compact emitting regions within the lobes at 5.5\,GHz; at this frequency the emission is short-lived, implying a recent jet switch-off. On the other hand, only three show evidence of aged lobe plasma by the presence of an ultra-steep spectrum ($\alpha<-1.2$) and a diffuse, low surface-brightness radio morphology. The predominant fraction of young remnants is consistent with a rapid fading during the remnant phase. Within our sample of radio galaxies, our observations constrain the remnant fraction to $4\%\lesssim f_{\mathrm{rem}} \lesssim 10\%$; the lower limit comes from the limiting case in which all remnant candidates with hotspots are simply active radio galaxies with faint, undetected radio cores. Finally, we model the synchrotron spectrum arising from a hotspot to show they can persist for 5--10\,Myr at 5.5\,GHz after the jets switch off -- radio emission arising from such hotspots can therefore be expected in an appreciable fraction of genuine remnants.

In the new era of Extremely Large Telescopes (ELTs) currently under construction, challenging requirements drive spectrograph designs towards techniques that efficiently use a facility's light collection power. Operating in the single-mode (SM) regime, close to the diffraction limit, reduces the footprint of the instrument compared to a conventional high-resolving power spectrograph. The custom built injection fiber system with 3D-printed micro-lenses on top of it for the replicable high-resolution exoplanet and asteroseismology spectrograph at Subaru in combination with extreme adaptive optics of SCExAO, proved its high efficiency in a lab environment, manifesting up to ~77% of the theoretical predicted performance.

The detonation of a helium shell on top of a carbon-oxygen white dwarf has been argued as a potential explosion mechanism for type Ia supernovae (SNe~Ia). The ash produced during helium shell burning can lead to light curves and spectra that are inconsistent with normal SNe~Ia, but may be viable for some objects showing a light curve bump within the days following explosion. We present a series of radiative transfer models designed to mimic predictions from double detonation explosion models. We consider a range of core and shell masses, and systematically explore multiple post-explosion compositions for the helium shell. We find that a variety of luminosities and timescales for early light curve bumps result from those models with shells containing $^{56}$Ni, $^{52}$Fe, or $^{48}$Cr. Comparing our models to SNe~Ia with light curve bumps, we find that these models can reproduce the shapes of almost all of the bumps observed, but only those objects with red colours around maximum light ($B-V \gtrsim 1$) are well matched throughout their evolution. Consistent with previous works, we also show that those models in which the shell does not contain iron-group elements provide good agreement with normal SNe~Ia of different luminosities from shortly after explosion up to maximum light. While our models do not amount to positive evidence in favour of the double detonation scenario, we show that provided the helium shell ash does not contain iron-group elements, it may be viable for a wide range of normal SNe~Ia.

Keck-telescope spectrophotometry of the companion of PSR J1810+1744 shows a flat, but asymmetric light-curve maximum and a deep, narrow minimum. The maximum indicates strong gravity darkening near the L_1 point, along with a heated pole and surface winds. The minimum indicates a low underlying temperature and substantial limb darkening. The gravity darkening is a consequence of extreme pulsar heating and the near-filling of the Roche lobe. Light-curve modeling gives a binary inclination i=65.7+/-0.4deg. With the Keck-measured radial-velocity amplitude K_c=462.3+/-2.2km/s, this gives an accurate neutron star mass M_NS=2.13+/-0.04M_o, with important implications for the dense-matter equation of state. A classic direct-heating model, ignoring the L_1 gravitational darkening, would predict an unphysical M_NS>3M_o. A few other ``spider" pulsar binaries have similar large heating and fill factor; thus, they should be checked for such effects.

Recently published ALMA observations suggest the presence of 20 ppb PH$_3$ in the upper clouds of Venus. This is an unexpected result, as PH$_3$ does not have a readily apparent source and should be rapidly photochemically destroyed according to our current understanding of Venus atmospheric chemistry. While the reported PH$_3$ spectral line at 266.94 GHz is nearly co-located with an SO$_2$ spectral line, the non-detection of stronger SO$_2$ lines in the wideband ALMA data is used to rule out SO$_2$ as the origin of the feature. We present a reassessment of wideband and narrowband datasets derived from these ALMA observations. The ALMA observations are re-reduced following both the initial and revised calibration procedures discussed by the authors of the original study. We also investigate the phenomenon of apparent spectral line dilution over varying spatial scales resulting from the ALMA antenna configuration. A 266.94 GHz spectral feature is apparent in the narrowband data using the initial calibration procedures, but this same feature can not be identified following calibration revisions. The feature is also not reproduced in the wideband data. While the SO$_2$ spectral line is not observed at 257.54 GHz in the ALMA wideband data, our dilution simulations suggest that SO$_2$ abundances greater than the previously suggested 10 ppb limit would also not be detected by ALMA. Additional millimeter, sub-millimeter, and infrared observations of Venus should be undertaken to further investigate the possibility of PH$_3$ in the Venus atmosphere.

IceCube has performed several all-sky searches for point-like neutrino sources using track-like events, including a recent time-integrated analysis using 10 years of IceCube data. This paper accompanies the public data release of these neutrino candidates detected by IceCube between April 6, 2008 and July 8, 2018. The selection includes through-going tracks, primarily due to muon neutrino candidates, that reach the detector from all directions, as well as neutrino track events that start within the instrumented volume. An updated selection and reconstruction for data taken after April 2012 slightly improves the sensitivity of the sample. While more than 80% of the sample overlaps between the old and new versions, differing events can lead to changes relative to the previous 7 year event selection. An a posteriori estimate of the significance of the 2014-2015 TXS flare is reported with an explanation of observed discrepancies with previous results. This public data release, which includes 10 years of data and binned detector response functions for muon neutrino signal events, shows improved sensitivity in generic time-integrated point source analyses and should be preferred over previous releases.

The observation of a 266.94 GHz feature in the Venus spectrum has been attributed to PH$_3$ in the Venus clouds, suggesting unexpected geological, chemical or even biological processes. Since both PH$_3$ and SO$_2$ are spectrally active near 266.94 GHz, the contribution to this line from SO$_2$ must be determined before it can be attributed, in whole or part, to PH$_3$. An undetected SO$_2$ reference line, interpreted as an unexpectedly low SO$_2$ abundance, suggested that the 266.94 GHz feature could be attributed primarily to PH$_3$. However, the low SO$_2$ and the inference that PH$_3$ was in the cloud deck posed an apparent contradiction. Here we use a radiative transfer model to analyze the PH$_3$ discovery, and explore the detectability of different vertical distributions of PH$_3$ and SO$_2$. We find that the 266.94 GHz line does not originate in the clouds, but above 80 km in the Venus mesosphere. This level of line formation is inconsistent with chemical modeling that assumes generation of PH$_3$ in the Venus clouds. Given the extremely short chemical lifetime of PH$_3$ in the Venus mesosphere, an implausibly high source flux would be needed to maintain the observed value of 20$\pm$10 ppb. We find that typical Venus SO$_2$ vertical distributions and abundances fit the JCMT 266.94 GHz feature, and the resulting SO$_2$ reference line at 267.54 GHz would have remained undetectable in the ALMA data due to line dilution. We conclude that nominal mesospheric SO$_2$ is a more plausible explanation for the JCMT and ALMA data than PH$_3$.

Multiple stellar populations (MPs) are a distinct characteristic of Globular Clusters (GCs). Their general properties have been widely studied among main sequence, red giant branch (RGB) and horizontal branch (HB) stars, but a common framework is still missing at later evolutionary stages. We studied the MP phenomenon along the AGB sequences in 58 GCs, observed with the Hubble Space Telescope in ultraviolet (UV) and optical filters. By using UV-optical color-magnitude diagrams, we selected the AGB members of each cluster and identified the AGB candidates of the metal-enhanced population in type II GCs. We studied the photometric properties of AGB stars and compared them to theoretical models derived from synthetic spectra analysis. We observe the following features: i) the spread of AGB stars in photometric indices sensitive to variations of light-elements and helium is typically larger than that expected from photometric errors; ii) the fraction of metal-enhanced stars in the AGB is lower than in the RGB in most of the type II GCs; iii) the fraction of 1G stars derived from the chromosome map of AGB stars in 15 GCs is larger than that of RGB stars; v) the AGB/HB frequency correlates with the average mass of the most helium-enriched population. These findings represent a clear evidence of the presence of MPs along the AGB of Galactic GCs and indicate that a significant fraction of helium-enriched stars, which have lower mass in the HB, does not evolve to the AGB phase, leaving the HB sequence towards higher effective temperatures, as predicted by the AGB-manqu\'e scenario.

The space gravitational wave (GW) detector Laser Interferometer Space Antenna (LISA) that was planed to launch in the early 2030s is to detect the low-frequency GW signals in the Galaxy. AM CVn stars were generally thought to be important low-frequency GW sources. Employing the MESA code, in this work we calculate the evolution of a great number of binary systems consisting of a white dwarf (WD) and a main sequence (MS) star, and diagnose whether their descendants-AM CVn stars can be visible by the LISA. The simulated results show that the progenitors of these LISA sources within a distance of 1 kpc are WD-MS binaries with a donor star of $1.0-1.4~M_\odot$ (for initial WD mass of $0.5~M_\odot$) or $1.0-2.0~M_\odot$ (for initial WD mass of $0.7~M_\odot$), and an initial orbital period slightly smaller than the bifurcation period. Our simulations also indicate ten verification AM CVn sources can be reproduced by the standard magnetic braking model, and are potential LISA sources. Based on the birthrate of AM CVns simulated by the population synthesis, the birthrate of AM CVn-LISA sources evolving from the evolved donor star channel within a distance of 1 kpc can be estimated to be $(0.6-1.4)\times10^{-6}~\rm yr^{-1}$, and the predicted number of AM CVn-LISA sources is about $340-810$. Therefore, the evolved donor star channel play an important role in forming AM CVn-LISA sources in the Galaxy.

Accurate modelling of redshift-space distortions (RSD) is challenging in the non-linear regime for two-point statistics e.g. the two-point correlation function (2PCF). We take a different perspective to split the galaxy density field according to the local density, and cross-correlate those densities with the entire galaxy field. Using mock galaxies, we demonstrate that combining a series of cross-correlation functions (CCFs) offers improvements over the 2PCF as follows: 1. The distribution of peculiar velocities in each split density is nearly Gaussian. This allows the Gaussian streaming model for RSD to perform accurately within the statistical errors of a ($1.5\,h^{-1}$Gpc)$^3$ volume for almost all scales and all split densities. 2. The PDF of the density field at small scales is non-Gaussian, but the CCFs of split densities capture the non-Gaussianity, leading to improved cosmological constraints over the 2PCF. We can obtain unbiased constraints on the growth parameter $f\sigma_{12}$ at the per-cent level, and Alcock-Paczynski (AP) parameters at the sub-per-cent level with the minimal scale of $15\,h^{-1}{\rm Mpc}$. This is a $\sim$30 per cent and $\sim$6 times improvement over the 2PCF, respectively. The diverse and steep slopes of the CCFs at small scales are likely to be responsible for the improved constraints of AP parameters. 3. Baryon acoustic oscillations (BAO) are contained in all CCFs of split densities. Including BAO scales helps to break the degeneracy between the line-of-sight and transverse AP parameters, allowing independent constraints on them. We discuss and compare models for RSD around spherical densities.

Both the absolute magnitude of type Ia supernovae (SNe Ia) and the luminosity distance of them are modified in the context of the minimally extended varying speed of light (meVSL) model compared to those of general relativity (GR). We have analyzed the likelihood of various dark energy models under meVSL by using the Pantheon SNe Ia data. Both $\omega$CDM and CPL parameterization dark energy models indicate a cosmic variation of the speed of light at the 1-$\sigma$ level. For $\Omega_{\text{m} 0} = 0.30, 0.31$, and 0.32 with $(\omega_0 \,, \omega_a) = (-1 \,, 0)$, 1-$\sigma$ range of $\dot{\tilde{c}}_0/\tilde{c}_0 \, (10^{-13} \, \text{yr}^{-1}) $ are (-8.76 \,, -0.89), (-11.8 \,, 3.93), and (-14.8 \,, -6.98), respectively. Meanwhile, 1-$\sigma$ range of $\dot{\tilde{c}}_0/\tilde{c}_0 (10^{-12} \, \text{yr}^{-1}) $ for the CPL dark energy models with $-1.05 \leq \omega_{0} \leq -0.95$ and $0.28 \leq \Omega_{\text{m} 0} \leq 0.32$, are (-6.31\,, -2.98). The value of $\tilde{c}$ at $z = 3$ can be larger than that of the present by $0.2 \sim 3$ \% for $\omega$CDM models and $5 \sim 13$ \% for CPL models. We also obtain $-25.6 \leq \dot{\tilde{G}}_0/\tilde{G}_0 \, (10^{-12} \, \text{yr}^{-1}) \leq -0.36$ for viable models except for CPL model for $\Omega_{\text{m} 0} = 0.28$. We obtain the increasing rate of the gravitational constant as $1.65 \leq \dot{\tilde{G}}_0/\tilde{G}_0 \, (10^{-12} \, \text{yr}^{-1}) \leq 3.79$ for that model.

Deep near-infrared $J$- and $K$-band photometry of three Local Group dwarf spheroidal galaxies: Fornax, Carina, and Sculptor, is made available for the community. Until now, these data have only been used by the Araucaria Project to determine distances using the tip of the red giant and RR Lyrae stars. Now, we present the entire data collection in a form of a database, consisting of accurate $J$- and $K$-band magnitudes, sky coordinates, ellipticity measurements, and timestamps of observations, complemented by stars' loci in their reference images. Depth of our photometry reaches about 22 mag at 5$\sigma$ level, and is comparable to NIR surveys, like the UKIRT Infrared Deep Sky Survey (UKIDSS) or the VISTA Hemisphere Survey (VHS). Small overlap with VHS and no overlap with UKIDSS makes our database a unique source of quality photometry.

We have examined the star formation history (SFH) of Andromeda VII (And VII), the brightest and most massive dwarf spheroidal (dSph) satellite of the Andromeda galaxy (M 31). Although M 31 is surrounded by several dSph companions with old stellar populations and low metallicity, it has a metal-rich stellar halo with an age of 6$-$8 Gyr. This indicates that any evolutionary association between the stellar halo of M 31 and its dSph system is frail. Therefore, the question is whether And VII (a high-metallicity dSph located $\sim$220 kpc from M 31), can be associated with M 31's young, metal-rich halo. Here, we perform the first reconstruction of the SFH of And VII employing long-period variable (LPV) stars. As the most-evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars, the birth mass of LPVs can be determined by connecting their near-infrared photometry to theoretical evolutionary tracks. We found 55 LPV candidates within two half-light radii, using multi-epoch imaging with the Isaac Newton Telescope in the $i$ and $V$ bands. Based on their birth mass function, the star-formation rate (SFR) of And VII was obtained as a function of cosmic time. The main epoch of star formation occurred $\simeq 6.2$ Gyr ago with a SFR of $0.006\pm0.002$ M$_\odot$ yr$^{-1}$. Over the past 6 Gyr, we find slow star formation, which continued until 500 Myr ago with a SFR $\sim0.0005\pm0.0002$ M$_\odot$ yr$^{-1}$. We determined And VII's stellar mass $M=(13.3\pm5.3)\times10^6$ M$_\odot$ within a half-light radius $r_{\frac{1}{2}}=3.8\pm0.3$ arcmin and metallicity $Z=0.0007$, and also derived its distance modulus of $\mu=24.38$ mag.

We present the Galactic Radio Explorer (GReX), an all-sky monitor to probe the brightest bursts in the radio sky. Building on the success of STARE2, we will search for fast radio bursts (FRBs) emitted from Galactic magnetars as well as bursts from nearby galaxies. GReX will search down to ~ten microseconds time resolution, allowing us to find new super giant radio pulses from Milky Way pulsars and study their broadband emission. The proposed instrument will employ ultra-wide band (0.7--2 GHz) feeds coupled to a high performance (receiver temperature <10 K) low noise amplifier (LNA) originally developed for the DSA-110 and DSA-2000 projects. In GReX Phase I (GReX-I), unit systems will be deployed at Owens Valley Radio Observatory (OVRO), NASA's Goldstone station, and at Telescope Array, Delta Utah. Phase II will expand the array, placing feeds in India, Australia, and elsewhere in order to build up to continuous coverage of nearly 4$\pi$ steradians and to increase our exposure to the Galactic plane. We model the local magnetar population to forecast for GReX, finding the improved sensitivity and increased exposure to the Galactic plane could lead to dozens of FRB-like bursts per year.

Lyman-$\alpha$ forest data are known to be a good probe of the small scale matter power. In this paper, we explore the redshift evolution of the observable effective optical depth $\tau_{\rm eff} (z)$ from the Lyman-$\alpha$ data as a discriminator between dark matter models that differ from the $\Lambda$CDM model on small scales. We consider the thermal warm dark matter (WDM) and the ultra-light axion (ULA) models for the following set of parameters: the mass of ULA, $m_a \simeq 10^{-24}\hbox{--}10^{-22} \, \rm eV$ and WDM mass, $m_{\rm wdm} = 0.1 \hbox{--} 4.6 \, \rm keV$. We simulate the line-of-sight HI density and velocity fields using semi-analytic methods. The simulated effective optical depth for the alternative dark matter models diverges from the $\Lambda$CDM model for $z \gtrsim 3$, which provides a meaningful probe of the matter power at small scales. Using likelihood analysis, we compare the simulated data with the high-resolution Lyman-$\alpha$ cloud data in the redshift range $2 < z < 4.2$. The one-dimensional posterior probabilities of the relevant masses peak at $m_a \simeq 5 \times 10^{-23} \, \rm eV$ and $m_{\rm wdm} \simeq 1.1 \, \rm keV$. The posterior probabilities remain flat for larger masses, which shows that the data is compatible with the $\Lambda$CDM model.

We report multi-wavelength study of a complex M-class solar eruptive flare that consists of three different sets of flare ribbons, viz. circular, parallel, and remote ribbons. Magnetic field modeling of source active region NOAA 12242 exhibits the presence of 3D null-point magnetic topology which encompasses an inner bipolar region. The event initiates with the faint signatures of the circular ribbon along with remote brightening right from the pre-flare phase which points toward the ongoing slow yet persistent null-point reconnection. We first detected flux cancellation and an associated brightening, which are likely signatures of tether-cutting reconnection that builds the flux rope near the polarity inversion line (PIL) of the inner bipolar region. In the next stage, with the onset of M8.7 flare, there is a substantial enhancement in the brightening of circular ribbon which essentially suggests an increase in the rate of ongoing null-point reconnection. Finally, the eruption of underlying flux rope triggers ``standard flare reconnection" beneath it producing an abrupt rise in the intensity of the parallel ribbons as well as enhancing the rate of null-point reconnection by external forcing. We show that within the the fan dome, the region with magnetic decay index n>1.5 borders the null-point QSL. Our analysis suggests that both the torus instability and the breakout model have played role toward the triggering mechanism for the eruptive flare. This event is a nice example of the dynamical evolution of a flux rope initially confined in a null-point topology, that subsequently activates and erupts with the progression of the circular -- cum -- parallel ribbon flare.

The sensitivity of the Advanced LIGO detectors to gravitational waves can be affected by environmental disturbances external to the detectors themselves. Since the transition from the former initial LIGO phase, many improvements have been made to the equipment and techniques used to investigate these environmental effects. These methods have aided in tracking down and mitigating noise sources throughout the first three observing runs of the advanced detector era, keeping the ambient contribution of environmental noise below the background noise levels of the detectors. In this paper we describe the methods used and how they have led to the mitigation of noise sources, the role that environmental monitoring has played in the validation of gravitational wave events, and plans for future observing runs.

We measure the scalar induced gravitational waves from the cosmic microwave background (CMB) observations and the gravitational wave observations. In the $\Lambda$CDM+$r$ model, we constrain the cosmological parameters within the evolution of the scalar induced gravitational waves by the additional scalar power spectrum. The two special cases called narrow power spectrum and wide power spectrum have influence on the cosmological parameters, especially the combinations of Planck18+BAO+BK15+LISA. We also compare these numerical results from four datasets within LIGO, LISA, IPTA and FAST projects, respectively. The constraints from FAST have a significant impact on tensor-to-scalar ratio.

Due to its proximity to Earth, Jupiter of the Solar System serves as a unique case study for gas-giant exoplanets. In the current study, we perform fits of ab initio, reflective, semi-infinite, homogeneous model atmospheres to 61 phase curves from 0.40 to 1.00 $\mu$m, obtained from the Cassini spacecraft, within a Bayesian framework. We reproduce the previous finding that atmospheric models using classic reflection laws (Lambertian, Rayleigh, single Henyey-Greenstein) provide poor fits to the data. Using the double Henyey-Greenstein reflection law, we extract posterior distributions of the single-scattering albedo and scattering asymmetry factors and tabulate their median values and uncertainties. We infer that the aerosols in the Jovian atmosphere are large, irregular, polydisperse particles that produce strong forward scattering together with a narrow backscattering lobe. The near-unity values of the single-scattering albedos imply that multiple scattering of radiation is an important effect. We speculate that the observed narrow backscattering lobe is caused by coherent backscattering of radiation, which is usually associated with Solar System bodies with solid surfaces and regolith. Our findings demonstrate that precise, multi-wavelength phase curves encode valuable information on the fundamental properties of cloud/haze particles. The method described in this Letter enables single-scattering albedos and scattering asymmetry factors to be retrieved from James Webb Space Telescope phase curves of exoplanets.

We analyze new XMM-Newton and archival Chandra observations of the middle-aged $\gamma$-ray radio-quiet pulsar J1957+5033. We detect, for the first time, X-ray pulsations with the pulsar spin period of the point-like source coinciding by position with the pulsar. This confirms the pulsar nature of the source. In the 0.15--0.5 keV band, there is a single pulse per period and the pulsed fraction is $\approx18\pm6$ per cent. In this band, the pulsar spectrum is dominated by a thermal emission component that likely comes from the entire surface of the neutron star, while at higher energies ($\gtrsim0.7$ keV) it is described by a power law with the photon index $\Gamma \approx 1.6$. We construct new hydrogen atmosphere models for neutron stars with dipole magnetic fields and non-uniform surface temperature distributions with relatively low effective temperatures. We use them in the spectral analysis and derive the pulsar average effective temperature of $\approx(2-3)\times10^5$ K. This makes J1957+5033 the coldest among all known thermally emitting neutron stars with ages below 1 Myr. Using the interstellar extinction--distance relation, we constrain the distance to the pulsar in the range of 0.1--1 kpc. We compare the obtained X-ray thermal luminosity with those for other neutron stars and various neutron star cooling models and set some constraints on latter. We observe a faint trail-like feature, elongated $\sim 8$ arcmin from J1957+5033. Its spectrum can be described by a power law with a photon index $\Gamma=1.9\pm0.5$ suggesting that it is likely a pulsar wind nebula powered by J1957+5033.

In this work, we study the potential of the Cherenkov Telescope Array (CTA) for the detection of Galactic dark matter (DM) subhalos. We focus on low-mass subhalos that do not host any baryonic content and therefore lack any multiwavelength counterpart. If the DM is made of weakly interacting massive particles (WIMPs), these dark subhalos may thus appear in the gamma-ray sky as unidentified sources. A detailed characterization of the instrumental response of CTA to dark subhalos is performed, for which we use the {\it ctools} analysis software and simulate CTA observations under different array configurations and pointing strategies, such as the scheduled extragalactic survey. This, together with information on the subhalo population as inferred from N-body cosmological simulations, allows us to predict the CTA detectability of dark subhalos, i.e., the expected number of subhalos in each of the considered observational scenarios. In the absence of detection, for each observation strategy we set competitive limits to the annihilation cross section as a function of the DM particle mass, that are at the level of $\langle\sigma v\rangle\sim4\times10^{-24}$ ($7\times10^{-25}$) $\mathrm{cm^3s^{-1}}$ for the $b\bar{b}$ ($\tau^+\tau^-$) annihilation channel in the best case scenario. Interestingly, we find the latter to be reached with no dedicated observations, as we obtain the best limits by just accumulating exposure time from all scheduled CTA programs and pointings over the first 10 years of operation. This way CTA will offer the most constraining limits from subhalo searches in the intermediate range between $\sim 1-3$ TeV, complementing previous results with \textit{Fermi}-LAT and HAWC at lower and higher energies, respectively.

We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics.

The mechanisms that bring galaxies to strongly reduce their star formation activity (star-formation quenching) is still poorly understood. To better study galaxy evolution, we propose a classification based on the maps of the ionised hydrogen distribution, traced by kpc-resolved, equivalent width of H$\alpha$ maps, and the nuclear activity of the galaxies using information from the BPT diagnostic diagrams. Using these tools, we group a sample of 238 galaxies from the CALIFA survey in six quenching stages (QS): objects dominated by recent star formation; systems that present a quiescent-nuclear-ring structure in their centre; galaxies that are centrally-quiescent; galaxies with no clear pattern in their ionisation gas distribution - mixed; systems that posses only a few star-forming regions - nearly-retired, or galaxies that are completely quiescent - fully-retired. Regarding their nuclear activity, we further divide the galaxies into two groups - active systems that host a weak or strong AGN in their centre, and non-active objects. Galaxies grouped into quenching stage classes occupy specific locations on the star-formation-rate versus stellar mass diagram. The "Blue cloud" is populated by the star-forming and the quiescent-nuclear-ring galaxies, the "Green valley" is populated by centrally-quiescent and mixed systems, "Red sequence" by the nearly- and fully-retired objects. Generally, galaxies that host a weak or strong AGN show properties, comparable to the non-active counterparts at the same quenching stages, except for the AGN-hosting star-forming systems. The degree of the star-formation quenching increases along the present emission-line pattern sequence from star-forming to fully-retired. The proposed emission-line classes reinforce the "inside-out" quenching scenario, which foresees that the suppression of the star-formation begins from the central regions of the galaxies.

All the studies of the interaction between tides and a convective flow assume that the large scale tides can be described as a mean shear flow which is damped by small scale fluctuating convective eddies. The convective Reynolds stress is calculated using mixing length theory, accounting for a sharp suppression of dissipation when the turnover timescale is larger than the tidal period. This yields tidal dissipation rates several orders of magnitude too small to account for the circularization periods of late-type binaries or the tidal dissipation factor of giant planets. Here, we argue that the above description is inconsistent, because fluctuations and mean flow should be identified based on the timescale, not on the spatial scale, on which they vary. Therefore, the standard picture should be reversed, with the fluctuations being the tidal oscillations and the mean shear flow provided by the largest convective eddies. We assume that energy is locally transferred from the tides to the convective flow. Using this assumption, we obtain values for the tidal $Q$ factor of Jupiter and Saturn and for the circularization periods of PMS binaries in good agreement with observations. The timescales obtained with the equilibrium tide approximation are however still 40 times too large to account for the circularization periods of late-type binaries. For these systems, shear in the tachocline or at the base of the convective zone may be the main cause of tidal dissipation.

We report microscopic, cathodoluminescence, chemical and O isotopic measurements of FeO-poor isolated olivine grains (IOG) in the carbonaceous chondrites Allende (CV3), Northwest Africa 5958 (C2-ung), Northwest Africa 11086 (CM2-an), Allan Hills 77307 (CO3.0). The general petrographic, chemical and isotopic similarity with bona fide type I chondrules confirms that the IOG derived from them. The concentric CL zoning, reflecting a decrease in refractory elements toward the margins, and frequent rimming by enstatite are taken as evidence of interaction of the IOG with the gas as stand-alone objects. This indicates that they were splashed out of chondrules when these were still partially molten. CaO-rich refractory forsterites, which are restricted to $\Delta^{17}O < -4\permil$ likely escaped equilibration at lower temperatures because of their large size and possibly quicker quenching. The IOG thus bear witness to frequent collisions in the chondrule-forming regions.

Pulsar timing is a technique that uses the highly stable spin periods of neutron stars to investigate a wide range of topics in physics and astrophysics. Pulsar timing arrays (PTAs) use sets of extremely well-timed pulsars as a Galaxy-scale detector with arms extending between Earth and each pulsar in the array. These challenging experiments look for correlated deviations in the pulsars' timing that are caused by low-frequency gravitational waves (GWs) traversing our Galaxy. PTAs are particularly sensitive to GWs at nanohertz frequencies, which makes them complementary to other space- and ground-based detectors. In this chapter, we will describe the methodology behind pulsar timing; provide an overview of the potential uses of PTAs; and summarise where current PTA-based detection efforts stand. Most predictions expect PTAs to successfully detect a cosmological background of GWs emitted by supermassive black-hole binaries and also potentially detect continuous-wave emission from binary supermassive black holes, within the next several years.

Compound chondrules, i.e. chondrules fused together, make a powerful probe of the density and compositional diversity in chondrule-forming environments, but their abundance among the dominating porphyritic textures may have been drastically underestimated. I report herein microscopic observations and LA-ICP-MS analyses of lobate chondrules in the CO3 chondrites Miller Range 07193 and 07342. Lobes in a given chondrule show correlated volatile and moderately volatile element abundances but refractory element concentrations are essentially independent. This indicates that they formed by the collision of preexisting droplets whose refractory elements behaved in closed system, while their more volatile elements were buffered by the same gaseous medium. The presence of lobes would otherwise be difficult to explain, as surface tension should have rapidly imposed a spherical shape at the temperature peak. In fact, since most chondrules across chondrite groups are nonspherical, a majority are probably compounds variously relaxed toward sphericity. The lack of correlation of refractory elements between conjoined compound chondrule components is inconsistent with derivation of chondrules from the disruption of homogenized melt bodies as in impact scenarios and evokes rather the melting of independent mm-size nebular aggregates. Yet a "nebular" setting for chondrule formation would need to involve not only increased solid concentration, e.g. by settling to the midplane, but also a boost in relative velocities between droplets during chondrule-forming events to account for observed compound chondrule frequencies .

Very steep reflection emissivity profiles in the inner part of accretion disks are commonly found in the analysis of X-ray observations of black hole binaries and AGNs, but there is some debate about their exact origin. While steep reflection emissivity profiles can be naturally produced by compact coronae close to black holes, the measured radial emissivity parameter can be further increased by the radial disk ionization profile when the theoretical model assumes a disk with constant ionization. In this paper, we implement the possibility of a radial disk ionization profile in the reflection model RELXILL_NK and we analyze a NuSTAR observation of the black hole binary EXO 1846-031, which was previously found to have a very high inner emissivity index. We find that the model with a radial disk ionization profile improves the fit, but the impact on the estimate of the black hole spin parameter and on the constraint of the deformation parameter is modest.

We analyzed the daily sunspot-group data reported by the Greenwich Photoheliographic Results (GPR) during the period 1874-1976 and Debrecen Photoheligraphic Data (DPD) during the period 1977-2017, and the revised Version-2 of ISSN during the period 1874-2017. We determined the amplitude of the Solar Cycles 12-24 and the 13-month smoothed monthly mean corrected areas of the sunspot groups in the Sun's whole-sphere (WSGA), northern hemisphere (NSGA), and southern hemisphere (SSGA) at the epochs of the maxima of the Solar Cycles 12-24. Using all these we obtained the relations similar to that found in our earlier analyzes -- i.e., the existence of a high correlation between the sum of the areas of sunspot groups in the southern-hemisphere near-equatorial band during a small interval (7-9 months) just after a maximum epoch of a solar cycle and the amplitude of next solar cycle -- separately for the Sun's whole-sphere and northern- and southern-hemispheres. By using these relations we predicted ~701 msh (millionth of solar hemisphere), ~429 msh, and ~366 msh for the values of WSGA, NSGA, and SSGA, respectively, at the maximum epoch of the Solar Cycle 25. We predicted 86 + or - 18 for the amplitude of Solar Cycle 25. The 13-month smoothed monthly mean sunspot-group area highly correlate with that of ISSN. Using this relation and the predicted values of WSGA, NSGA, and SSGA we have obtained 68 + or -11 for the amplitude of Solar Cycle 25, which is slightly lower than the aforementioned predicted value, and 39 + or - 4 and 31 + or - 6 for the values of northern- and southern-hemispheres' sunspot numbers at the maximum epoch of Solar Cycle 25. Overall, our results suggest that the amplitude of Solar Cycle 25 would be 25%-40% smaller, and the corresponding north-south asymmetry would be much smaller, than those of Solar Cycle 24.

In order to understand stellar evolution, it is crucial to efficiently determine stellar surface rotation periods. An efficient tool to automatically determine reliable rotation periods is needed when dealing with large samples of stellar photometric datasets. The objective of this work is to develop such a tool. Random forest learning abilities are exploited to automate the extraction of rotation periods in Kepler light curves. Rotation periods and complementary parameters are obtained from three different methods: a wavelet analysis, the autocorrelation function of the light curve, and the composite spectrum. We train three different classifiers: one to detect if rotational modulations are present in the light curve, one to flag close binary or classical pulsators candidates that can bias our rotation period determination, and finally one classifier to provide the final rotation period. We test our machine learning pipeline on 23,431 stars of the Kepler K and M dwarf reference rotation catalog of Santos et al. (2019) for which 60% of the stars have been visually inspected. For the sample of 21,707 stars where all the input parameters are provided to the algorithm, 94.2% of them are correctly classified (as rotating or not). Among the stars that have a rotation period in the reference catalog, the machine learning provides a period that agrees within 10% of the reference value for 95.3% of the stars. Moreover, the yield of correct rotation periods is raised to 99.5% after visually inspecting 25.2% of the stars. Over the two main analysis steps, rotation classification and period selection, the pipeline yields a global agreement with the reference values of 92.1% and 96.9% before and after visual inspection. Random forest classifiers are efficient tools to determine reliable rotation periods in large samples of stars. [abridged]

The abundance patterns observed in the Sun and in metal-poor stars show a clear odd-even effect. An important question is whether the odd-even effect in solar-metallicity stars is similar to the Sun, or if there are variations that can tell us about different chemical enrichment histories. In this work, we report for the first time observational evidence of a differential odd-even effect in the solar twin HIP 11915, relative to the solar odd-even abundance pattern. The spectra of this star were obtained with high resolving power (140 000) and signal-to-noise ratio ($\sim$420) using the ESPRESSO spectrograph and the VLT telescope. Thanks to the high spectral quality, we obtained extremely precise stellar parameters ($\sigma(T_{\rm eff})$ = 2 K, $\sigma(\rm{[Fe/H]})$ = 0.003 dex, and $\sigma(\log g)$ = 0.008 dex). We determine the chemical abundance of 20 elements ($Z\leq39$) with high precision ($\sim$0.01 dex), which shows a strong pattern of the odd-even effect even after performing Galactic Chemical Evolution corrections. The odd-even effect is reasonably well-reproduced by a core-collapse supernova of 13 $\rm{M_{\odot}}$ and metallicity Z = 0.001 diluted into a metal-poor gas of 1 $\rm{M_{\odot}}$. Our results indicate that HIP 11915 has an odd-even effect slightly different than the Sun, thus confirming a different supernova enrichment history.

Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million years old and could be experiencing rapid evolution of the neutral hydrogen fraction in the intergalactic medium. The protocluster is identified by an overdensity of 6 times the average galaxy density, and with 21 narrowband selected Lyman-$\alpha$ galaxies, among which 16 have been spectroscopically confirmed. At redshifts similar to or above this record, smaller protogroups with fewer members have been reported. LAGER-z7OD1 shows an elongated shape and consists of two subprotoclusters, which would have merged into one massive cluster with a present-day mass of $3.7 \times 10^{15}$ solar masses. The total volume of the ionized bubbles generated by its member galaxies is found to be comparable to the volume of the protocluster itself, indicating that we are witnessing the merging of the individual bubbles and that the intergalactic medium within the protocluster is almost fully ionized. LAGER-z7OD1 thus provides a unique natural laboratory to investigate the reionization process.

The Birmingham Solar Oscillations Network (BiSON) is a collection of ground-based automated telescopes observing oscillations of the Sun. The network has been operating since the early 1990s. We present development work on a prototype next generation observation platform, BiSON:NG, based almost entirely on inexpensive off-the-shelf components, and where the footprint is reduced to a size that can be inexpensively installed on the roof of an existing building. Continuous development is essential in ensuring that automated networks such as BiSON are well placed to observe the next solar cycle and beyond.

CONTEXT. Gaia EDR3 has produced parallaxes for 1.468x10^9 sources but there are calibration issues that require corrections to the published values and uncertainties. AIMS. We want to characterize the behavior of the uncertainties of the Gaia EDR3 parallaxes. We also aim to provide a procedure for the calculation of distances to stars and stellar clusters. METHODS. We reanalyze some of the data in the calibration papers for QSO and LMC parallaxes and combine those results with measurements for six bright GCs. We calculate the angular covariance of EDR3 parallaxes at small separations based on the LMC results and combine it with the results for larger angles using QSOs to obtain an analytical formula for the angular covariance over the whole sky. The results for the six GCs are used to validate the parallax bias correction as a function of magnitude, color, and ecliptic latitude and to determine the constant used to convert internal uncertainties to external ones. RESULTS. The angular covariance at zero separation is 106.2 muas^2, yielding a minimum uncertainty for EDR3 parallaxes of 10.3 muas for individual stars. That value can be only slightly reduced for GCs after considering the behavior of the angular covariance of the parallaxes for small separations. The Lindegren et al. parallax bias correction works quite well, except for the brighter magnitudes, suggesting improvements may be possible there. The value of k is 1.1-1.7 and depends on G. Stars with moderately large values of RUWE can still provide useful parallaxes albeit with larger values of k. We give accurate and precise Gaia EDR3 distances to the six GCs and for the specific case of 47 Tuc we are able to beat the angular covariance limit and derive a high-precision distance of 4.53+-0.06 kpc. Finally, a recipe for the derivation of distances to stars and stellar clusters using Gaia EDR3 parallaxes is given. [ABRIDGED]

Fluences of solar energetic particles (SEPs) are not easy to evaluate, especially for high-energy events (i.e. ground-level enhancements, GLEs). Earlier estimates of event-integrated SEP fluences for GLEs were based on partly outdated assumptions and data, and they required revisions. Here, we present the results of a full revision of the spectral fluences for most major SEP events (GLEs) for the period from 1956 -- 2017 using updated low-energy flux estimates along with greatly revisited high-energy flux data and applying the newly invented reconstruction method including an improved neutron-monitor yield function. Low- and high-energy parts of the SEP fluence were estimated using a revised space-borne/ionospheric data and ground-based neutron monitors, respectively. The measured data were fitted by the modified Band function spectral shape. The best-fit parameters and their uncertainties were assessed using a direct Monte Carlo method. As a result, a full reconstruction of the event-integrated spectral fluences was performed in the energy range above 30 MeV, parametrised, and tabulated for easy use along with estimates of the 68% confidence intervals. This forms a solid basis for more precise studies of the physics of solar eruptive events and the transport of energetic particles in the interplanetary medium, as well as the related applications.

The increase in the sensitivity of gravitational wave interferometers will bring additional detections of binary black hole and double neutron star mergers. It will also very likely add many merger events of black hole - neutron star binaries. Distinguishing mixed binaries from binary black holes mergers for high mass ratios could be challenging because in this situation the neutron star coalesces with the black hole without experiencing significant disruption. To investigate the transition of mixed binary mergers into those behaving more like binary black hole coalescences, we present results from merger simulations for different mass ratios. We show how the degree of deformation and disruption of the neutron star impacts the inspiral and merger dynamics, the properties of the final black hole, the accretion disk formed from the circularization of the tidal debris, the gravitational waves, and the strain spectrum and mismatches. The results also show the effectiveness of the initial data method that generalizes the Bowen-York initial data for black hole punctures to the case of binaries with neutron star companions.

Axion-like particles (ALPs) are predicted in some well-motivated theories beyond the Standard Model. The TeV gamma-rays from active galactic nuclei (AGN) suffers attenuation by the pair production interactions with the cosmic infrared background light (EBL/CMB) during its travel to the earth. The attenuation can be circumvented through photon-ALP conversions in the AGN and Galaxy magnetic-field, and a flux enhancement is expected to arise in the observed spectrum. In this work, we study the potential of the AGN gamma-ray spectrum for energy up to above 100 TeV to probe ALP-parameter space at around $\mu$eV, where the coupling $g_{a\gamma}$ is so far relatively weak constrained. We find the nearby and bright sources, Mrk 501, IC 310 and M 87, are suitable for our objective. Assuming an intrinsic spectrum exponential cutoff energy at $E_{\rm c}$=100 TeV, we extrapolate the observed spectra of these sources up to above 100 TeV by the models with/without ALPs. For $g_{a\gamma}\gtrsim 2\times$$10^{-11} \rm GeV^{-1}$ with $m_{a}\lesssim1\,\mu$eV, the flux at around 100 TeV predicted by the ALP model can be enhanced more than an order of magnitude than that from the standard absorption, and could be detected by LHAASO. Our result is subject to the uncertainty from the intrinsic spectrum above tens of TeV, which require further observations on these sources by the forthcoming CTA, LHAASO, SWGO and so on.

The Large-Scale Structure (LSS) of the universe has the potential to provide decisive answers to the remaining open questions in cosmology. Early attempts at modelling it analytically focused on using perturbation theory. However, small-scale effects introduced by gravitational collapse cannot be described perturbatively and this failure of perturbation theory is reflected even on the largest scales. The Effective Field Theory of Large Scale Structure (EFTofLSS) has emerged as a consistent method for describing LSS on large scales by introducing counterterms that account for the effects of small-scale dynamics. So far studies of the EFT have mostly focused on the two and three point functions with little attention devoted to the four point function or trispectrum. The trispectrum probes cubic interactions arising from non-linear clustering, biasing, and primordial non-Gaussianities, and constitutes a key element of the covariance matrix of the power spectrum. In this paper, we present explicit calibrations of the EFT counterterms for the one-loop trispectrum. Specifically, we find clear evidence for non-zero EFT corrections. We define two one-parameter ans\"{a}tze for the counterterm of the one-loop propagator and show that they provide a good correction to the residual at scales below k~0.07 h/Mpc. We then take the amplitudes of the linear and quadratic counterkernels calculated in our previous paper on the bispectrum and use them in the remaining counterterms, establishing consistency of the counterterms in the two, three and four point function. We also show that the commonly used EdS approximation for the growth of the density fields leads to errors that are of the same magnitude as loop corrections to the trispectrum on large scales.

The literature on precision differential abundances (PDAs) in stars is extensive. Surveys include sun-like stars in the solar neighborhood, binary systems, and Galactic clusters. Numerous references as well as a discussion of relevant mechanisms may be found in papers by Ramirez, et al. (2019) and Nagar, et al. (2020). A strong impetus for this work is the probability that the abundances have been influenced by exoplanetary systems and their evolution. We calculate the resulting differential abundances ([El/H]) assuming a given amount of material with the composition of the bulk earth (Wang, et al. 2018) was added to the stellar convection zone of a dwarf G-type star. The mass of the convection zone is uncertain and variable, depending on the spectral type. Here, we assume a mass of $5\cdot 10^{31}$ gm for the stellar convection zone (SCZ). This is 0.025 $M_\odot$ (Pinsonneault, et al. 2001, Chambers, 2010). For other SCZ masses, the parameters must be adjusted accordingly. In general, the sunlike star will not have exactly the solar composition. This contingency is roughly taken into account in our model. An appendix discusses issues of volatility and condensation temperature.

The Simons Observatory (SO) will be a cosmic microwave background (CMB) survey experiment with three small-aperture telescopes and one large-aperture telescope, which will observe from the Atacama Desert in Chile. In total, SO will field over 60,000 transition-edge sensor (TES) bolometers in six spectral bands centered between 27 and 280 GHz in order to achieve the sensitivity necessary to measure or constrain numerous cosmological quantities, as outlined in The Simons Observatory Collaboration et al. (2019). These telescopes require 33 highly transparent, large aperture, refracting optics. To this end, we developed mechanically robust, highly efficient, metamaterial anti-reflection (AR) coatings with octave bandwidth coverage for silicon optics up to 46 cm in diameter for the 22-55, 75-165, and 190-310 GHz bands. We detail the design, the manufacturing approach to fabricate the SO lenses, their performance, and possible extensions of metamaterial AR coatings to optical elements made of harder materials such as alumina.

Space-based telescopes offer unparalleled opportunities for characterising exoplanets, Solar System bodies and stellar objects. However, observatories in low Earth orbits (e.g. Hubble, CHEOPS, Twinkle and an ever increasing number of cubesats) cannot always be continuously pointed at a target due to Earth obscuration. For exoplanet observations consisting of transit, or eclipse, spectroscopy this causes gaps in the light curve, which reduces the information content and can diminish the science return of the observation. Terminus, a time-domain simulator, has been developed to model the occurrence of these gaps to predict the potential impact on future observations. The simulator is capable of radiometrically modelling exoplanet observations as well as producing light curves and spectra. Here, Terminus is baselined on the Twinkle mission but the model can be adapted for any space-based telescope and is especially applicable to those in a low-Earth orbit. Terminus also has the capability to model observations of other targets such as asteroids or brown dwarfs.

Spectra of stochastic gravitational waves (GW) generated in cosmological first-order phase transitions are computed within strongly correlated theories with a dual holographic description. The theories are mostly used as models of dark sectors. In particular, we consider the so-called Witten-Sakai-Sugimoto model, a $SU(N)$ gauge theory coupled to different matter fields in both the fundamental and the adjoint representations. The model has a well-known top-down holographic dual description which allows us to perform reliable calculations in the strongly coupled regime. We consider the GW spectra from bubble collisions and sound waves arising from two different kinds of first-order phase transitions: a confinement/deconfinement one and a chiral symmetry breaking/restoration one. Depending on the model parameters, we find that the GW spectra may fall within the sensibility region of ground-based and space-based interferometers, as well as of Pulsar Timing Arrays. In the latter case, the signal could be compatible with the recent potential observation by NANOGrav. When the two phase transitions happen at different critical temperatures, characteristic spectra with double frequency peaks show up. Moreover, in this case we explicitly show how to correct the redshift factors appearing in the formulae for the GW power spectra to account for the fact that adiabatic expansion from the first transition to the present times cannot be assumed anymore.

Existing searches for cosmic axions relics have relied heavily on the axion being non-relativistic and making up dark matter. However, light axions can be copiously produced in the early Universe and remain relativistic today, thereby constituting a Cosmic $\textit{axion}$ Background (C$a$B). As prototypical examples of axion sources, we consider thermal production, dark-matter decay, parametric resonance, and topological defect decay. Each of these has a characteristic frequency spectrum that can be searched for in axion direct detection experiments. We focus on the axion-photon coupling and study the sensitivity of current and future versions of ADMX, HAYSTAC, DMRadio, and ABRACADABRA to a C$a$B, finding that the data collected in search of dark matter can be repurposed to detect axion energy densities well below limits set by measurements of the energy budget of the Universe. In this way, direct detection of relativistic relics offers a powerful new opportunity to learn about the early Universe and, potentially, discover the axion.

The flux of cosmic rays (CRs) in the heliosphere is subjected to remarkable time variations caused by the 11-year cycle of solar activity. To help the study of this effect, we have developed a web application (Heliophysics Virtual Observatory) that collects real-time data on solar activity, interplanetary plasma, and charged radiation from several space missions or observatories. As we will show, our application can be used to visualize, manipulate, and download updated data on sunspots, heliospheric magnetic fields, solar wind, and neutron monitors counting rates. Data and calculations are automatically updated on daily basis. A nowcasting for the energy spectrum of CR protons near-Earth is also provided using calculations and real-time neutron monitor data as input.

We study effects of Lorentz-invariance violation on the rotation of neutron stars (NSs) in the minimal gravitational Standard-Model Extension framework, and calculate the quadrupole radiation generated by them. Aiming at testing Lorentz invariance with observations of continuous gravitational waves (GWs) from rotating NSs in the future, we compare the GW spectra of a rotating ellipsoidal NS under Lorentz-violating gravity with those of a Lorentz-invariant one. The former are found to possess frequency components higher than the second harmonic, which does not happen for the latter, indicating those higher frequency components to be potential signatures of Lorentz violation in continuous GW spectra of rotating NSs.

UltraNest is a general-purpose Bayesian inference package for parameter estimation and model comparison. It allows fitting arbitrary models specified as likelihood functions written in Python, C, C++, Fortran, Julia or R. With a focus on correctness and speed (in that order), UltraNest is especially useful for multi-modal or non-Gaussian parameter spaces, computational expensive models, in robust pipelines. Parallelisation to computing clusters and resuming incomplete runs is available.

Nested sampling (NS) computes parameter posterior distributions and makes Bayesian model comparison computationally feasible. Its strengths are the unsupervised navigation of complex, potentially multi-modal posteriors until a well-defined termination point. A systematic literature review of nested sampling algorithms and variants is presented. We focus on complete algorithms, including solutions to likelihood-restricted prior sampling. A new formulation of NS is presented, which casts the parameter space exploration as a search on a tree. Previously published ways of obtaining robust error estimates and dynamic variations of the number of live points are presented as special cases of this formulation.

A covariant, scalar-tensor gravity is constructed such that the static, spherically symmetric Rezzolla-Zhidenko metric is an exact solution to the theory. The equations describing gravitational perturbations of this spacetime, which represents a generic black hole possessing an arbitrary number of hairs, can then be derived. This allows for a self-consistent study of the associated quasi-normal modes. It is shown that mode spectra are tied to not only the non-Einstein parameters in the metric but also to those that appear at the level of the action, and that different branches of the exact theory can, in some cases, predict significantly different oscillation frequencies and damping times. For choices which make the theory appear more like general relativity in some precise sense, we find that a non-trivial Rezzolla-Zhidenko parameter space is permissible under current constraints on fundamental ringdown modes observed by Advanced LIGO.

The scaling of the turbulent spectra provides a key measurement that allows to discriminate between different theoretical predictions of turbulence. In the solar wind, this has driven a large number of studies dedicated to this issue using in-situ data from various orbiting spacecraft. While a semblance of consensus exists regarding the scaling in the MHD and dispersive ranges, the precise scaling in the transition range and the actual physical mechanisms that control it remain open questions. Using the high-resolution data in the inner heliosphere from Parker Solar Probe (PSP) mission, we find that the sub-ion scales (i.e., at the frequency f ~ [2, 9] Hz) follow a power-law spectrum f^a with a spectral index a varying between -3 and -5.7. Our results also show that there is a trend toward and anti-correlation between the spectral slopes and the power amplitudes at the MHD scales, in agreement with previous studies: the higher the power amplitude the steeper the spectrum at sub-ion scales. A similar trend toward an anti-correlation between steep spectra and increasing normalized cross helicity is found, in agreement with previous theoretical predictions about the imbalanced solar wind. We discuss the ubiquitous nature of the ion transition range in solar wind turbulence in the inner heliosphere.

Transport properties of dense quark matter are discussed in the strong magnetic field, B. B dependence as well as density dependence of the Hall conductivity is discussed in the inhomogeneous chiral phase. Anomalous Hall effect is intrinsic to the inhomogeneous chiral phase and resembles the one in Weyl semimetals in condensed matter physics. Some theoretical aspects inherent in anomalous Hall effect are revealed.

Muon detectors and neutron monitors were recently installed at Syowa Station, in the Antarctic, to observe different types of secondary particles resulting from cosmic ray interactions simultaneously from the same location. Continuing observations will give new insight into the response of muon detectors to atmospheric and geomagnetic effects. Operation began in February, 2018 and the system has been stable with a duty-cycle exceeding 94%. Muon data shows a clear seasonal variation, which is expected from the atmospheric temperature effect. We verified successful operation by showing that the muon and neutron data are consistent with those from other locations by comparing intensity variations during a space weather event. We have established a web page to make real time data available with interactive graphics (this http URL).

We consider the geodesic deviation equation, describing the relative accelerations of nearby particles, and the Raychaudhuri equation, giving the evolution of the kinematical quantities associated with deformations (expansion, shear and rotation) in the Weyl type $f(Q,T)$ gravity, in which the nonmetricity $Q$ is represented in the standard Weyl form, fully determined by the Weyl vector, while $T$ represents the trace of the matter energy-momentum tensor. The effects of the Weyl geometry and of the extra force induced by the nonmetricity-matter coupling are explicitly taken into account. The Newtonian limit of the theory is investigated, and the generalized Poisson equation, containing correction terms coming from the Weyl geometry, and from the geometry matter coupling, is derived. As a physical application of the geodesic deviation equation the modifications of the tidal forces, due to the nonmetricity-matter coupling, are obtained in the weak field approximation. The tidal motion of test particles is directly influenced by the gradients of the extra force, and of the Weyl vector. As a concrete astrophysical example we obtain the expression of the Roche limit (the orbital distance at which a satellite begins to be tidally torn apart by the body it orbits) in the Weyl type $f(Q,T)$ gravity.

Silicon detection is a mature technology for registering the passage of charged particles. At the same time it continues to evolve toward increasing radiation tolerance as well as precision and adaptability. For these reasons it is likely to remain a critical element of detection systems associated with extra-terrestrial exploration. Silicon sensor leakage current and depletion voltage depend upon the integrated fluence received by the sensor, and upon its thermal history during and after the irradiation process. For minimal assumptions on shielding and hence on particle energy spectrum, and using published data on Martian ground temperature, we predict the leakage current density and the depletion voltage, as a function of time, of silicon sensors deployed continuously on the Mars surface for a duration of up to 28 Earth-years, for several sensor geometries and a worst-case temperature scenario.

In this pedestrian approach I give my personal point of view on the various problems posed by dark matter in the universe. After a brief historical overview I discuss the various solutions stemming from high energy particle physics, and the current status of experimental research on candidate particles (WIMPS). In the absence of direct evidence, the theories can still be evaluated by comparing their implications for the formation of galaxies, clusters and superclusters of galaxies against astronomical observations. I conclude briefly with the attempts to circumvent the dark matter problem by modifying the laws of gravity.

We examine the impact of a faster expanding Universe on the phenomenology of scalar dark matter (DM) associated with $SU(2)_L$ multiplets. Earlier works with radiation dominated Universe have reported the presence of desert region for both inert $SU(2)_L$ doublet and triplet DM candidates where the DM is under abundant. We find that the existence of a faster expanding component before BBN can revive a major part of the desert parameter space consistent with relic density requirements and other direct and indirect search bounds. We also review the possible collider search prospects of the newly obtained parameter space and show that such region can be probed at the future colliders with improved sensitivity via a stable charged track.