Papers for Wednesday, Jul 28 2021

The three-body problem is a fundamental long-standing open problem, with applications in all branches of physics, including astrophysics, nuclear physics and particle physics. In general, conserved quantities allow to reduce the formulation of a mechanical problem to fewer degrees of freedom, a process known as dynamical reduction. However, extant reductions are either non-general, or hide the problem's symmetry or include unexplained definitions. This paper presents a dynamical reduction that avoids these issues, and hence is general and natural. Any three-body configuration defines a triangle, and its orientation in space. Accordingly, we decompose the dynamical variables into the geometry (shape + size) and orientation of the triangle. The geometry variables are shown to describe the motion of an abstract point in a curved 3d space, subject to a potential-derived force and a magnetic-like force with a monopole charge. The orientation variables are shown to obey a dynamics analogous to the Euler equations for a rotating rigid body, only here the moments of inertia depend on the geometry variables, rather than being constant. The reduction rests on a novel symmetric solution to the center of mass constraint inspired by Lagrange's solution to the cubic. The formulation of the orientation variables is novel and rests on a little known generalization of the Euler-Lagrange equations to non-coordinate velocities. Applications to special exact solutions and to the statistical solution are described or discussed. Moreover, a generalization to the four-body problem is presented.

Dark matter interactions with electrons or protons during the early Universe leave imprints on the cosmic microwave background and the matter power spectrum, and can be probed through cosmological and astrophysical observations. We explore these interactions using a diverse suite of data: cosmic microwave background anisotropies, baryon acoustic oscillations, the Lyman-$\alpha$ forest, and the abundance of Milky-Way subhalos. We derive constraints using model-independent parameterizations of the dark matter--electron and dark matter--proton interaction cross sections and map these constraints onto concrete dark matter models. Our constraints are complementary to other probes of dark matter interactions with ordinary matter, such as direct detection, big bang nucleosynthesis, various astrophysical systems, and accelerator-based experiments.

We present new observational constraints on the elastic scattering of dark matter with electrons for dark matter masses between 10 keV and 1 TeV. We consider scenarios in which the momentum-transfer cross section has a power-law dependence on the relative particle velocity, with a power-law index $n \in \{-4,-2,0,2,4,6\}$. We search for evidence of dark matter scattering through its suppression of structure formation. Measurements of the cosmic microwave background temperature, polarization, and lensing anisotropy from \textit{Planck} 2018 data and of the Milky Way satellite abundance measurements from the Dark Energy Survey and Pan-STARRS1 show no evidence of interactions. We use these data sets to obtain upper limits on the scattering cross section, comparing them with exclusion bounds from electronic recoil data in direct detection experiments. Our results provide the strongest bounds available for dark matter--electron scattering derived from the distribution of matter in the Universe, extending down to sub-MeV dark matter masses, where current direct detection experiments lose sensitivity.

Young open clusters (t<200 Myr) have been observed to exhibit several peculiarities in their chemical compositions, from a slightly sub-solar iron content, super-solar abundances of some atomic species (e.g. ionised chromium), and atypical enhancements of [Ba/Fe], with values up to +0.7 dex. Regarding the behaviour of the other $s$-process elements like yttrium, zirconium, lanthanum, and cerium, there is general disagreement in the literature. In this work we expand upon our previous analysis of a sample of five young open clusters (IC2391, IC2602, IC4665, NGC2516, and NGC2547) and one star-forming region (NGC2264), with the aim of determining abundances of different neutron-capture elements, mainly CuI, SrI, SrII, YII, ZrII, BaII, LaII, and CeII. We analysed high-resolution, high signal-to-noise spectra of 23 solar-type stars observed within the \textit{Gaia}-ESO survey. We find that our clusters have solar [Cu/Fe] within the uncertainties, while we confirm the super-solar [Ba/Fe] values (from +0.22 to +0.64 dex). Our analysis also points to mildly enhanced [Y/Fe] values (from 0 and +0.3 dex). For the other $s$-process elements we find that [X/Fe] ratios are solar at all ages. It is not possible to reconcile the anomalous behaviour of Ba and Y at young ages with standard stellar yields and Galactic chemical evolution model predictions. Thus, we explore different possible scenarios related to the behaviour of spectral lines, from the sensitivity to the presence of magnetic fields to the first ionisation potential effect. We also investigate the possibility that they may arise from alterations of the structure of the stellar photosphere due to higher levels of activity in such young stars. We are still unable to explain these enhancements, but we suggest that other elements (i.e. La) might be more reliable tracer of the $s$-process at young ages and encourage further observations.

X-ray binaries are long-standing source candidates of Galactic cosmic rays and neutrinos. The compact object in a binary system can be the site for cosmic-ray acceleration, while high-energy neutrinos can be produced by the interactions of cosmic rays in the jet of the compact object, the stellar wind, or the atmosphere of the companion star. We report a time-dependent study of high-energy neutrinos from X-ray binaries with IceCube using 7.5 years of muon neutrino data and X-ray observations. In the absence of significant correlation, we report upper limits on the neutrino fluxes from these sources and provide a comparison with theoretical predictions.

We consider the Blandford-Znajek (BZ) mechanism for extracting black hole spin energy to drive astrophysical jets. In analyses of the BZ mechanism it is always assumed that the electric charge of the black hole remains zero. But, as noted by Wald and others, if the medium surrounding the black hole is an ionised plasma with mobile charges, then a spinning hole quickly acquires an electric charge. The effect of this charge is to nullify the electric field structures which drive the BZ mechanism -- the electric and magnetic fields then obey ${\bf E\cdot B= 0}$ everywhere. Since jets are now observed in a wide variety of classes of accreting objects, most of which do not contain a central black hole, it seems likely that the jet driving mechanism in all astrophysical objects uses energy directly from the accretion disc, rather than black hole spin.

The dominant uncertainty in the current measurement of the Hubble constant ($H_0$) with strong gravitational lensing time delays is attributed to uncertainties in the mass profiles of the main deflector galaxies. Strongly lensed supernovae (glSNe) can provide, in addition to measurable time delays, lensing magnification constraints when knowledge about the unlensed apparent brightness of the explosion is imposed. We present a hierarchical Bayesian framework to combine a dataset of SNe that are not strongly lensed and a dataset of strongly lensed SNe with measured time delays. We jointly constrain (i) $H_0$ using the time delays as an absolute distance indicator, (ii) the lens model profiles using the magnification ratio of lensed and unlensed fluxes on the population level and (iii) the unlensed apparent magnitude distribution of the SNe population and the redshift-luminosity relation of the relative expansion history of the Universe. We apply our joint inference framework on a future expected data set of glSNe, and forecast that a sample of 144 glSNe of Type~Ia with well measured time series and imaging data will measure $H_0$ to 1.5\%. We discuss strategies to mitigate systematics associated with using absolute flux measurements of glSNe to constrain the mass density profiles. Using the magnification of SNe images is a promising and complementary alternative to using stellar kinematics. Future surveys, such as the Rubin and \textit{Roman} observatories, will be able to discover the necessary number of glSNe, and with additional follow-up observations this methodology will provide precise constraints on mass profiles and $H_0$.

The interpretation of data from indirect detection experiments searching for dark matter annihilations requires computationally expensive simulations of cosmic-ray propagation. In this work we present a new method based on Recurrent Neural Networks that significantly accelerates simulations of secondary and dark matter Galactic cosmic ray antiprotons while achieving excellent accuracy. This approach allows for an efficient profiling or marginalisation over the nuisance parameters of a cosmic ray propagation model in order to perform parameter scans for a wide range of dark matter models. We identify importance sampling as particularly suitable for ensuring that the network is only evaluated in well-trained parameter regions. We present resulting constraints using the most recent AMS-02 antiproton data on several models of Weakly Interacting Massive Particles. The fully trained networks are released as DarkRayNet together with this work and achieve a speed-up of the runtime by at least two orders of magnitude compared to conventional approaches.

Numerous spherical ``shells" have been observed in young star-forming environments that host low- and intermediate-mass stars. These observations suggest that these shells may be produced by isotropic stellar wind feedback from young main-sequence stars. However, the driving mechanism for these shells remains uncertain because the momentum injected by winds is too low to explain their sizes and dynamics due to their low mass-loss rates. However, these studies neglect how the wind kinetic energy is transferred to the ISM and instead assume it is instantly lost via radiation, suggesting that these shells are momentum-driven. Intermediate-mass stars have fast ($v_w \gtrsim 1000$ km/s) stellar winds and therefore the energy injected by winds should produce energy-driven adiabatic wind bubbles that are larger than momentum-driven wind bubbles. Here, we explore if energy-driven wind feedback can produce the observed shells by performing a series of 3D magneto-hydrodynamic simulations of wind feedback from intermediate-mass and high-mass stars that are placed in a magnetized, turbulent molecular cloud. We find that, for the high-mass stars modeled, energy-driven wind feedback produces $\sim$pc scale wind bubbles in molecular clouds that agree with the observed shell sizes but winds from intermediate-mass stars can not produce similar shells because of their lower mass-loss rates and velocities. Therefore, such shells must be driven by other feedback processes inherent to low- and intermediate-mass star formation.

We present a spatio-temporal AI framework that concurrently exploits both the spatial and time-variable features of gravitationally lensed supernovae in optical images to ultimately aid in the discovery of such exotic transients in wide-field surveys. Our spatio-temporal engine is designed using recurrent convolutional layers, while drawing from recent advances in variational inference to quantify approximate Bayesian uncertainties via a confidence score. Using simulated Young Supernova Experiment (YSE) images as a showcase, we find that the use of time-series images yields a substantial gain of nearly 20 per cent in classification accuracy over single-epoch observations, with a preliminary application to mock observations from the Legacy Survey of Space and Time (LSST) yielding around 99 per cent accuracy. Our innovative deep learning machinery adds an extra dimension in the search for gravitationally lensed supernovae from current and future astrophysical transient surveys.

We report the analysis of the deep (270 ks) X-ray Chandra data of one of the most radio-loud, Seyfert 2 galaxies in the nearby Universe (z=0.01135), IC 5063. The alignment of the radio structure with the galactic disk and ionized bi-cone, enables us to study the effects of both radio jet and nuclear irradiation on the interstellar medium (ISM). The nuclear and bi-cone spectra suggest a low photoionization phase mixed with a more ionized or thermal gas component, while the cross-cone spectrum is dominated by shocked and collisionally ionized gas emission. The clumpy morphology of the soft (<3 keV) X-ray emission along the jet trails, and the large (~2.4 kpc) filamentary structure perpendicular to the radio jets at softer energies (<1.5 keV), suggest a large contribution of the jet-ISM interaction to the circumnuclear gas emission. The hard X-ray continuum (>3 keV) and the Fe K-alpha 6.4 keV emission are both extended to kpc size along the bi-cone direction, suggesting an interaction of nuclear photons with dense clouds in the galaxy disk, as observed in other Compton Thick (CT) active nuclei. The north-west cone spectrum also exhibits an Fe XXV emission line, which appears spatially extended and spatially correlated with the most intense radio hot-spot, suggesting jet-ISM interaction.

Spectrographs nominally contain a degree of quasi-static optical aberrations resulting from the quality of manufactured component surfaces, imperfect alignment, design residuals, thermal effects, and other other associated phenomena involved in the design and construction process. Aberrations that change over time can mimic the line centroid motion of a Doppler shift, introducing radial velocity (RV) uncertainty that increases time-series variability. Even when instrument drifts are tracked using a precise wavelength calibration source, barycentric motion of the Earth leads to a wavelength shift of stellar light causing a translation of the spectrum across the focal plane array by many pixels. The wavelength shift allows absorption lines to experience different optical propagation paths and aberrations over observing epochs. We use physical optics propagation simulations to study the impact of aberrations on precise Doppler measurements made by diffraction-limited, high-resolution spectrographs. We quantify the uncertainties that cross-correlation techniques introduce in the presence of aberrations and barycentric RV shifts. We find that aberrations which shift the PSF photo-center in the dispersion direction, in particular primary horizontal coma and trefoil, are the most concerning. To maintain aberration-induced RV errors less than 10 cm/s, phase errors for these particular aberrations must be held well below 0.05 waves at the instrument operating wavelength. Our simulations further show that wavelength calibration only partially compensates for instrumental drifts, owing to a behavioral difference between how cross-correlation techniques handle aberrations between starlight versus calibration light. Identifying subtle physical effects that influence RV errors will help ensure that diffraction-limited planet-finding spectrographs are able to reach their full scientific potential.

In this work we perform an observational data analysis on Einsteinian cubic gravity and $f(P)$ gravity with the objective of constraining the parameter space of the theories. We use the 30 point $z-H(z)$ cosmic chronometer data as the observational tool for our analysis along with the BAO and the CMB peak parameters. The $\chi^2$ statistic is used for the fitting analysis and it is minimized to obtain the best fit values for the free model parameters. We have used the Markov chain Monte Carlo algorithm to obtain bounds for the free parameters. To achieve this we used the publicly available \textit{CosmoMC} code to put parameter bounds and subsequently generate contour plots for them with different confidence intervals. Besides finding the Hubble parameter $H$ in terms of the redshift $z$ theoretically from our gravity models, we have exercised correlation coefficients and two \textit{machine learning} models, namely the linear regression (LR) and artificial neural network (ANN), for the estimation of $H(z)$. For this purpose, we have developed a \textit{Python} package for finding the parameter space, performing the subsequent statistical analysis and prediction analysis using machine learning. We compared both of our theoretical and estimated values of $H(z)$ with the observations. It is seen that our theoretical and estimated models from machine learning performed significantly well when compared with the observations.

Non-spherical rapid acceleration of mass (or energy) to a relativistic velocity is a natural source of gravitational radiation. Such conditions arise in both long and short gamma-ray bursts whose central engine ejects relativistic jets. The resulting gravitational wave signal is of a memory type, rising to a finite level (of order 4 G E/r) over a duration that corresponds to the longer of either the injection time and the acceleration time of the jet. We explore the properties of such signals and their potential detectability. Unfortunately, the expected signals are below the frequency band of Advanced LIGO-Virgo-Kagra, and above LISA. However, they fall within the range of the planned BBO and DECIGO. While current sensitivity is marginal for the detection of jet gravitational wave signals from GRBs, hidden relativistic jets that exist within some core collapse SNe could be detected. Such a detection would reveal the acceleration mechanism and the activity of the central engine, which cannot be explored directly in any other way

Ammonia (NH3) in a terrestrial planet atmosphere is generally a good biosignature gas, primarily because terrestrial planets have no significant known abiotic NH3 source. The conditions required for NH3 to accumulate in the atmosphere are, however, stringent. NH3's high water solubility and high bio-useability likely prevent NH3 from accumulating in the atmosphere to detectable levels unless life is a net source of NH3 and produces enough NH3 to saturate the surface sinks. Only then can NH3 accumulate in the atmosphere with a reasonable surface production flux. For the highly favorable planetary scenario of terrestrial planets with H2-dominated atmospheres orbiting M dwarf stars (M5V), we find a minimum of about 5 ppm column-averaged mixing ratio is needed for NH3 to be detectable with JWST, considering a 10 ppm JWST systematic noise floor. When the surface is saturated with NH3 (i.e., there are no NH3-removal reactions on the surface), the required biological surface flux to reach 5 ppm is on the order of 10^10 molecules cm-2 s-1, comparable to the terrestrial biological production of CH4. However, when the surface is unsaturated with NH3, due to additional sinks present on the surface, life would have to produce NH3 at surface flux levels on the order of 10^15 molecules cm-2 s-1 (approx. 4.5x10^6 Tg year-1). This value is roughly 20,000 times greater than the biological production of NH3 on Earth and about 10,000 times greater than Earth's CH4 biological production. Volatile amines have similar solubilities and reactivities to NH3 and hence share NH3's weaknesses and strengths as a biosignature. Finally, to establish NH3 as a biosignature gas, we must rule out mini-Neptunes with deep atmospheres, where temperatures and pressures are high enough for NH3's atmospheric production.

The InSight mission has operated on the surface of Mars for nearly two Earth years, returning detections of the first Marsquakes. The lander also deployed a meteorological instrument package and cameras to monitor local surface activity. These instruments have detected boundary layer phenomena, including small-scale vortices. These vortices register as short-lived, negative pressure excursions and closely resemble those that could generate dust devils. Although our analysis shows InSight encountered more than 900 vortices and collected more than 1000 images of the martian surface, no active dust devils were imaged. In spite of the lack of dust devil detections, we can leverage the vortex detections and InSight's daily wind speed measurements to learn about the boundary layer processes that create dust devils. We discuss our analysis of InSight's meteorological data to assess the statistics of vortex and dust devil activity. We also infer encounter distances for the vortices and, therefrom, the maximum vortex wind speeds. Surveying the available imagery, we place upper limits on what fraction of vortices carry dust (i.e., how many are bonafide dust devils) and estimate threshold wind speeds for dust lifting. Comparing our results to detections of dust devil tracks seen in space-based observations of the InSight landing site, we can also infer thresholds and frequency of track formation by vortices. Comparing vortex encounters and parameters with advective wind speeds, we find evidence that high wind speeds at InSight may have suppressed the formation of dust devils, explaining the lack of imaged dust devils.

The AMS-02 experiment has reported precise measurements of energy spectra of several cosmic-ray species in the range of ~(0.5-2000) GeV/n. An intriguing finding is the differences in the spectral shape between the different species. Protons exhibit the steepest spectrum of all the species, and helium, carbon, oxygen and iron spectra are found to be harder than that of neon, magnesium and silicon. These observations are difficult to explain as diffusive shock acceleration, the currently most plausible theory for cosmic particle acceleration at high energies, expects independence of the spectral index from mass and charge of the accelerated particle. Moreover, propagation in the Galaxy has been shown to not being able to compensate for this discrepancy. In this work, we present an explanation based on two-component model for the origin of cosmic rays in the Galaxy -- the first component originating from regular supernova remnants in the interstellar medium and the second component from Wolf-Rayet supernovae. Using recent results on cosmic-ray injection enhancement at supernova shocks in the uniform interstellar medium and in the wind environment of Wolf-Rayet stars, we show that the combination of the two components may explain most of the behavior observed by the AMS-02 experiment.

Varying fundamental constants (VFC) [e.g., the fine-structure constant, $\alpha_{\rm EM}$] can arise in numerous extended cosmologies. Through their effect on the decoupling of baryons and photons during last scattering and reionisation, these models can be directly constrained using measurements of the cosmic microwave background (CMB) temperature and polarization anisotropies. Previous investigations focused mainly on time-independent changes to the values of fundamental constants. Here we generalize to time-dependent variations. Instead of directly studying various VFC parameterizations, we perform a model-independent principal component analysis (PCA), directly using an eigenmode decomposition of the varying constant during recombination. After developing the formalism, we use Planck 2018 data to obtain new VFC limits, showing that three independent VFC modes can be constrained at present. No indications for significant departures from the standard model are found with Planck data. Cosmic variance limited modes are also compared and simple forecasts for The Simons Observatory are carried out, showing that in the future improvements of the current constraints by a factor of $\simeq 3$ can be anticipated. Our modes focus solely on VFC at redshifts $z\geq 300$. This implies that they do not capture some of the degrees of freedom relating to the reionisation era. This aspect provides important new insights into the possible origin of the Hubble tension, hinting that indeed a combined modification of recombination and reionisation physics could be at work. An extended PCA, covering both recombination and reionisation simultaneously, could shed more light on this question, as we emphasize here.

The persistently bright ultra-compact neutron star low-mass X-ray binary 4U 1820$-$30 displays a $\sim$170 d accretion cycle, evolving between phases of high and low X-ray modes, where the 3 -- 10 keV X-ray flux changes by a factor of up to $\approx 8$. The source is generally in a soft X-ray spectral state, but may transition to a harder state in the low X-ray mode. Here, we present new and archival radio observations of 4U 1820$-$30 during its high and low X-ray modes. For radio observations taken within a low mode, we observed a flat radio spectrum consistent with 4U 1820$-$30 launching a compact radio jet. However, during the high X-ray modes the compact jet was quenched and the radio spectrum was steep, consistent with optically-thin synchrotron emission. The jet emission appeared to transition at an X-ray luminosity of $L_{\rm X (3-10 keV)} \sim 3.5 \times 10^{37} (D/\rm{7.6 kpc})^{2}$ erg s$^{-1}$. We also find that the low-state radio spectrum appeared consistent regardless of X-ray hardness, implying a connection between jet quenching and mass accretion rate in 4U 1820$-$30, possibly related to the properties of the inner accretion disk or boundary layer.

The nature and geometry of the accretion flow in the low/hard state of black hole binaries is currently controversial. While most properties are generally explained in the truncated disc/hot inner flow model, the detection of a broad residual around the iron line argues for strong relativistic effects from an untruncated disc. Since spectral fitting alone is somewhat degenerate, we combine it with the additional information in the fast X-ray variability and perform a full spectral-timing analysis for NICER and NuSTAR data on a bright low/hard state of MAXI J1820+070. For the first time, we model the variability with propagating mass accretion rate fluctuations by combining two separate current insights: that the hot flow is spectrally inhomogeneous, and that there is a discontinuous jump in viscous time-scale between the hot flow and variable disc. Our model naturally gives the 'double hump' shape of the power spectra, and the increasing high frequency variability with energy in the second hump. Including reflection quantitatively reproduces the switch in the lag-frequency spectra, from hard lagging soft at low frequencies (propagation through the variable flow) to the soft lagging hard at the high frequencies (reverberation from the hard X-ray continuum illuminating the disc). The light travel time derived from the model corresponds to a distance of $\sim$ 45 gravitational radii, supporting the truncated disc model geometry for the low/hard state. The propagation lags allow us to measure the viscous time-scale in the hot flow, and the results favour SANE rather than MAD models for this source.

Planets are born from disks of gas and dust, and observations of protoplanetary disks are used to constrain the initial conditions of planet formation. However, dust mass measurements of Class II disks with ALMA have called into question whether they contain enough solids to build the exoplanets that have been detected to date. In this paper, we calculate the mass and spatial scale of solid material around Sun-like stars probed by transit and radial velocity exoplanet surveys, and compare those to the observed dust masses and sizes of Class II disks in the same stellar mass regime. We show that the apparent mass discrepancy disappears when accounting for observational selection and detection biases. We find a discrepancy only when the planet formation efficiency is below 100%, or if there is a population of undetected exoplanets that significantly contributes to the mass in solids. We identify a positive correlation between the masses of planetary systems and their respective orbital periods, which is consistent with the trend between the masses and the outer radii of Class II dust disks. This implies that, despite a factor 100 difference in spatial scale, the properties of protoplanetary disks seem to be imprinted on the exoplanet population.

We present the results of surveying [CI] $^3P_1-^3P_0$, $^{12}$CO $J=4-3$, and 630 $\mu$m dust continuum emission for 36 nearby ultra/luminous infrared galaxies (U/LIRGs) using the Band 8 receiver mounted on the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array. We describe the survey, observations, data reduction, and results; the main results are as follows. (i) We confirmed that [CI] $^3P_1-^3P_0$ has a linear relationship with both the $^{12}$CO $J=4-3$and 630 $\mu$m continuum. (ii) In NGC 6052 and NGC 7679, $^{12}$CO $J=4-3$ was detected but [CI] $^3P_1-^3P_0$ was not detected with a [CI] $^3P_1-^3P_0$/ $^{12}$CO $J=4-3$ ratio of $\lesssim0.08$. Two possible scenarios of weak [CI] $^3P_1-^3P_0$ emission are C$^0$-poor/CO-rich environments or an environment with an extremely large [CI] $^3P_1-^3P_0$ missing flux. (iii) There is no clear evidence showing that galaxy mergers, AGNs, and dust temperatures control the ratios of [CI] $^3P_1-^3P_0$/ $^{12}$CO $J=4-3$ and $L'_{\rm [CI](1-0)}/L_{\rm 630\mu m}$. (iv) We compare our nearby U/LIRGs with high-z galaxies, such as galaxies on the star formation main sequence (MS) at z$\sim1$ and submillimeter galaxies (SMGs) at $z=2-4$. We found that the mean value for the [CII] $^3P_1$--$^3P_0$/ $^{12}$CO $J=4-3$ ratio of U/LIRGs is similar to that of SMGs but smaller than that of galaxies on the MS.

Fitting half-integer generalized Laguerre functions to the evolved, real-space dark matter and halo correlation functions provides a simple way to reconstruct their initial shapes. We show that this methodology also works well in a wide variety of realistic, assembly biased, velocity biased and redshift-space distorted mock galaxy catalogs. We use the linear point feature in the monopole of the redshift-space distorted correlation function to quantify the accuracy of our approach. We find that the linear point estimated from the mock galaxy catalogs is insensitive to the details of the biasing scheme at the sub-percent level. However, the linear point scale in the nonlinear, biased, and redshift-space distorted field is systematically offset from its scale in the unbiased linear density fluctuation field by more than 1%. In the Laguerre reconstructed correlation function, this is reduced to sub-percent values, so it provides comparable accuracy and precision to methods that reconstruct the full density field before estimating the distance scale. The linear point in the reconstructed density fields provided by these other methods is likewise precise, accurate, and insensitive to galaxy bias. All reconstructions depend on some input parameters, and marginalizing over uncertainties in the input parameters required for reconstruction can degrade both accuracy and precision. The linear point simplifies the marginalization process, enabling more realistic estimates of the precision of the distance scale estimate for negligible additional computational cost. We show this explicitly for Laguerre reconstruction.

X-ray flares in gamma-ray bursts (GRBs) are believed to be generated by the late activities of central engine, and thus provide an useful tool to diagnose the properties of central objects. In this paper, we work on a GRB X-ray flare sample whose bulk Lorentz factors are constrained by two different methods and the jet opening angles are determined by the jet breaks in afterglow lightcurves. Considering a hyperaccreting stellar-mass black hole (BH) as the central engine of GRBs and the Blandford \& Znajek process (BZ) as the jet production mechanism, we constrain the parameters of central engine by using the X-ray flare data. We find that the BZ mechanism is so powerful making it possible to interpret both GRB prompt emissions and bright X-ray flares. The wind parameter ($p$) and accreted mass ($M_d$) fall into reasonable ranges. Our result is also applied to GRB 170817A. The late X-ray flare in GRB 170817A, if it is true, might not be a BH origin.

The Keck Planet Imager and Characterizer (KPIC) is a purpose-built instrument to demonstrate new technological and instrumental concepts initially developed for the exoplanet direct imaging field. Located downstream of the current Keck II adaptive optic system, KPIC contains a fiber injection unit (FIU) capable of combining the high-contrast imaging capability of the adaptive optics system with the high dispersion spectroscopy capability of the current Keck high resolution infrared spectrograph (NIRSPEC). Deployed at Keck in September 2018, this instrument has already been used to acquire high resolution spectra ($R > 30,000$) of multiple targets of interest. In the near term, it will be used to spectrally characterize known directly imaged exoplanets and low-mass brown dwarf companions visible in the northern hemisphere with a spectral resolution high enough to enable spin and planetary radial velocity measurements as well as Doppler imaging of atmospheric weather phenomena. Here we present the design of the FIU, the unique calibration procedures needed to operate a single-mode fiber instrument and the system performance.

Current state-of-the-art computational modeling makes it possible to build realistic models of stellar convection zones and atmospheres that take into account chemical composition, radiative effects, ionization, and turbulence. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of the stellar interior dynamics. Mixing-length models provide an initial approximation of stellar structure that can be used to initialize 3D radiative hydrodynamics simulations which include realistic modeling of turbulence, radiation, and other phenomena. In this paper, we present 3D radiative hydrodynamic simulations of an F-type main-sequence star with 1.47 solar mass. The computational domain includes the upper layers of the radiation zone, the entire convection zone, and the photosphere. The effects of stellar rotation is modeled in the f-plane approximation. These simulations provide new insight into the properties of the convective overshoot region, the dynamics of the near-surface, highly turbulent layer, and the structure and dynamics of granulation. They reveal anti-solar-type differential rotation and latitudinal dependence of the tachocline location.

Carbonaceous chondrite meteorites are so far the only available samples representing carbon-rich asteroids and in order to allow future comparison with samples returned by missions such as Hayabusa 2 and OSIRIS-Rex, is important to understand their physical properties. Future characterization of asteroid primitive classes, some of them targeted by sample-return missions, requires a better understanding of their mineralogy, the consequences of the exposure to space weathering, and how both affect the reflectance behavior of these objects. In this paper, the reflectance spectra of two chemically-related carbonaceous chondrites groups, precisely the Vigrano (CVs) and Karoonda (CKs), are measured and compared. The available sample suite includes polished sections exhibiting different petrologic types: from 3 (very low degree of thermal metamorphism) to 5 (high degree of thermal metamorphism). We found that the reflective properties and the comparison with the Cg asteroid reflectance class point toward a common chondritic reservoir from which the CV-CK asteroids collisionally evolved. In that scenario the CV and CK chondrites could be originated from 221 Eos asteroid family, but because of its collisional disruption, both chondrite groups evolved separately, experiencing different stages of thermal metamorphism, annealing and space weathering.

In this work, we study the statistical properties of soft gamma-/hard X-ray bursts from the soft gamma-ray repeater (SGR) 1806--20 and SGR J1935+2154 and of radio bursts from the repeating fast radio burst (FRB) 121102. For SGRs, we show that the probability density functions for the differences of fluences, fluxes, and durations at different times have fat tails with a $q$-Gaussian form. The $q$ values in the $q$-Gaussian distributions are approximately steady and independent of the temporal interval scale adopted, implying a scale-invariant structure of SGRs. These features indicate that SGR bursts may be governed by a self-organizing criticality (SOC) process. The recent discovery of a Galactic FRB associated with a hard X-ray burst from the Galactic magnetar SGR J1935+2154 has established the magnetar origin of at least some FRBs. Very recently, 1652 independent bursts from the repeating FRB 121102 have been detected by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Here we also investigate the scale-invariant structure of FRB 121102 based on the latest observation of FAST, and show that FRB 121102 and SGRs share similar statistical properties. Scale invariance in both FRB 121102 and SGRs can be well explained within the same physical framework of fractal-diffusive SOC systems.

Ultra-high-energy (UHE) cosmic rays (CRs) of energies $\sim (10^{18}-10^{20})~{\rm eV}$, accelerated in violent astrophysical environments, interact with cosmic background radiation fields via photo-hadronic processes, leading to strong attenuation. Typically, the Universe would become `opaque' to UHE CRs after several tens of Mpc, setting the boundary of the Greisen-Zatsepin-Kuz'min (GZK) horizon. In this work, we investigate the contribution of sources beyond the conventional GZK horizon to the UHE CR flux observed on Earth, when photo-spallation of the heavy nuclear CRs is taken into account. We demonstrate this contribution is substantial, despite the strong attenuation of UHE CRs. A significant consequence is the emergence of an isotropic background component in the observed flux of UHE CRs, coexisting with the anisotropic foreground component that are associated with nearby sources. Multi-particle CR horizons, which evolve over redshift, are determined by the CR nuclear composition. Thus, they are dependent on the source populations and source evolutionary histories.

IGR J00370+6122 is a high-mass X-ray binary, of which the primary is a B1 Ib star, whereas the companion is suggested to be a neutron star by the detection of 346-s pulsation in one-off 4-ks observation. To better understand the nature of the compact companion, the present work performs timing and spectral studies of the X-ray data of this object, taken with XMM-Newton, Swift, Suzaku, RXTE, and INTEGRAL. In the XMM-Newton data, a sign of coherent 674 s pulsation was detected, for which the previous 346-s period may be the 2nd harmonic. The spectra exhibited the "harder when brighter" trend in the 1$-$10 keV range, and a flat continuum without clear cutoff in the 10$-$80 keV range. These properties are both similar to those observed from several low-luminosity accreting pulsars, including X Persei in particular. Thus, the compact object in IGR J00370+6122 is considered to be a magnetized neutron star with a rather low luminosity. The orbital period was refined to $15.6649 \pm 0.0014$ d. Along the orbit, the luminosity changes by 3 orders of magnitude, involving a sudden drop from $\sim 4 \times 10^{33}$ to $\sim 1\times10^{32}$ erg s$^{-1}$ at an orbital phase of 0.3 (and probably vice verse at 0.95). Although these phenomena cannot be explained by a simple Hoyle-Lyttleton accretion from the primary's stellar winds, they can be explained when incorporating the propeller effect with a strong dipole magnetic field of $\sim 5 \times10^{13}$ G. Therefore, the neutron star in IGR J00370+6122 may have a stronger magnetic field compared to ordinary X-ray pulsars.

We present statistical analyses of a large homogenous data sample of Fermi-detected blazars thoroughly studied in order to reassess the relationship between flat-spectrum radio quasars (FSRQs) and subclasses of BL Lacertae objects (BL Lacs) blazar populations. We discovered from the average values of gamma-ray and X-ray spectral indices that the sequence of distribution is indicative of the blazar orientation scheme. Analyses of FSRQs and BL Lacs data show difference in the shape of gamma-ray and X-ray indices: significant anti-correlation (r~ greater than - 0.79) exists between gamma-ray and X-ray spectral indices. The spectral energy distributions of the blazar subclasses show that FSRQs and BL Lacs have similar spectral properties which can be unified through an evolutionary sequence. Nevertheless, there is a significant difference between the shapes of X-ray and gamma-ray spectra of blazars suggesting that different mechanisms are responsible for spectral variations in the two energy bands. All these results suggest that there is a form of a unified scheme for all blazars.

The CoMET is an R$\&$D project aiming to design a very-high-energy (VHE) gamma-ray observatory sensitive to energies above $\sim$ 200 GeV. The science goals include continuous observation of soft-spectrum VHE gamma-ray sources such as Active Galactic Nuclei (AGNs) and transients like Gamma-Ray Bursts (GRBs). With these objectives, CoMET is designed to have a low energy threshold with a wide field-of-view of about 2 sr, at a high altitude, and combines ALTO particle detectors with CLiC air-Cherenkov detectors. In this contribution, we focus on the ALTO particle detector array performance only. Water Cherenkov detectors are used for the detection of secondary particles in atmospheric air showers while scintillators serve as muon counters. A detailed study is presented through air-shower, detector and trigger simulations, followed by the reconstruction of the event parameters and the extraction of the signal (gamma-rays) from the background (cosmic-rays). We present the sensitivity of the ALTO detectors to a list of astrophysical sources using two SEMLA analysis configurations.

The formation of planetesimals was a key step in the assemblage of planetary bodies, yet many aspects of their formation remain poorly constrained. Notably, the mechanism by which chondrules -- sub-millimetric spheroids that dominate primitive meteorites -- were incorporated into planetesimals remains poorly understood. Here we classify and analyze particle-size distributions in various CO carbonaceous chondrites found in the Atacama Desert. Our results show that the average circle-equivalent diameters of chondrules define a positive trend with the petrographic grade, which reflects the progressive role of thermal metamorphism within the CO parent body. We show that this relationship could not have been established by thermal metamorphism alone but rather by aerodynamic sorting during accretion. By modeling the self-gravitational contraction of clumps of chondrules, we show that (i) the accretion of the CO parent body(ies) would have generated a gradual change of chondrule size with depth in the parent body, with larger chondrules being more centrally concentrated than smaller ones, and (ii) any subsequent growth by pebble accretion would have been insignificant. These findings give substantial support to the view that planetesimals formed via gravitational collapse.

Aims. We present a solution to determine the actual or physical relative positions between CCD chips, rather than determine their relative positions directly from the distortion-free coordinate. We verify the applicability of this solution from two types of astrometry. One is to make use of the astrometric catalog -- Gaia DR2 (Gaia Collaboration et al. 2018) to compute stars' celestial positions. We refer to the practice as photographic astrometry in this paper. And the other practice is to compute stars' relative positions by only using their pixel coordinates. We refer to it as differential astrometry. Meanwhile, we also aim to test the consistency between two types of astrometry from the results. Methods. By taking advantage of the GD solutions derived from observations, we relate the physical positions of the adjacent pixel edges of two CCD chips to each other and estimate the actual relative positions between chips. We implement the technique for two epochs of observations, which are taken from the CCD mosaic chips of the Bok 2.3-m telescope at Kitt Peak. Results. There is a good agreement between the two types of astrometry for the relative positions between chips. And the results provide us with more confidence in the differential astrometry of the planned Chinese Space Station Telescope (CSST). We achieved a 0.04 pixel uncertainty for the relative positions between chips on average, at least the factor of two improvement over the method which is used to monitor the progress of the interchip offset of CCD mosaic chips in the Hubble Space Telescope (HST). For the first-epoch observations, although suffered from possible variations in the non-linear parts of distortion, we still achieved a better uncertainty, 0.02 pixel on average.

Gamma-ray bursts (GRBs), as a possible probe to extend the Hubble diagram to high redshifts, have attracted much attention recently. In this paper, we select two samples of GRBs that have a plateau phase in X-ray afterglow. One is short GRBs with plateau phases dominated by magnetic dipole (MD) radiations. The other is long GRBs with gravitational-wave (GW) dominated plateau phases. These GRBs can be well standardized using the correlation between the plateau luminosity $L_0$ and the end time of plateau $t_b$. The so-called circularity problem is mitigated by using the observational Hubble parameter data and Gaussian process method. The calibrated \ltb ~correlations are also used to constrain $\Lambda$CDM and $w(z)$ = $w_{0}$ models. Combining the MD-LGRBs sample from Wang et al. (2021) and the MD-SGRBs sample, we find $\Omega_{m} = 0.33_{-0.09}^{+0.06}$ and $\Omega_{\Lambda}$ = $1.06_{-0.34}^{+0.15}$ excluding systematic uncertainties in the nonflat $\Lambda$CDM model. Adding type Ia supernovae from Pantheon sample, the best-fitting results are $w_{0}$ = $-1.11_{-0.15}^{+0.11}$ and $\Omega_{m}$ = $0.34_{-0.04}^{+0.05}$ in the $w=w_0$ model. These results are in agreement with the $\Lambda$CDM model. Our result supports that selection of GRBs from the same physical mechanism is crucial for cosmological purposes.

We present our analysis of the narrow emission lines produced in the plasma regions within the bright active galactic nucleus of NGC 4151, from an ORBYTS research-with-schools public engagement project. Our goal was to test whether the properties of these plasma regions changed between XMM- Newton observations spanning 15 years from 2000 to 2015, by measuring the outflow velocities and distances. From this study, we found that NGC 4151 has at least two to three plasma regions. There is no evidence of the outflowing wind properties changing as the velocities and distances are consistent throughout the observations.

We present the HR-pyPopStar model, which provides a complete set (in ages) of high resolution (HR) Spectral Energy Distributions of Single Stellar Populations. The model uses the most recent high wavelength-resolution theoretical atmosphere libraries for main sequence, post-AGB/planetary nebulae and Wolf-Rayet stars. The Spectral Energy Distributions are given for more than a hundred ages ranging from 0.1 Myr to 13.8 Gyr, at four different values of the metallicity (Z = 0.004, 0.008, 0.019 and 0.05), considering four different IMFs. The wavelength range goes from 91 to 24 000 {\AA} in linear steps {\delta}{\lambda} = 0.1 {\AA}, giving a theoretical resolving power R_{th,5000} ~ 50 000 at 5000 {\AA}. This is the main novelty of these spectra, unique for their age and wavelength ranges. The models include the ionising stellar populations that are relevant both at young (massive hot stars) as well as old (planetary nebulae) ages. We have tested the results with some examples of HR spectra recently observed with MEGARA at GTC. We highlight the importance of wavelength-resolution in reproducing and interpreting the observational data from the last and forthcoming generations of astronomical instruments operating at 8-10m class telescopes, with higher spectral resolution than their predecessors.

Determining the heavy-element accretion rate of growing giant planets is crucial for understanding their formation and bulk composition. The solid (heavy-element) accretion rate should be carefully modeled during the various stages of giant planet formation and therefore, the planetary capture radius must be determined. In some simulations that model the heavy-element accretion rate, such as in N-body simulations, the presence of the gaseous envelope is either neglected, or treated in an over-simplified manner. In this paper, we present an approximation for the capture radius that does not require the numerical solution of the stellar structure equations. Our approximation for the capture radius works extremely well for various planetesimal sizes and compositions. We show that the commonly assumed constant density assumption for inferring the capture radius leads to a large error in the calculated capture radius and we therefore suggest that our approximation should be implemented in future simulations.

In a collaboration between astroparticle physicists, animation artists from the award-winning Science Communication Lab, and musician Carsten Nicolai (a.k.a. Alva Noto), two cosmic particle accelerators have been brought to life: the massive binary star Eta Carinae, and the exploding star, which resulted in the gamma-ray burst GRB190829A. For Eta Carinae, the computer-generated images are close to reality because the measured orbital, stellar and wind parameters were used for this purpose. Particle acceleration in the jet of GRB190829A has also been animated at a level of detail not seen before. The internationally acclaimed multimedia artist Carsten Nicolai, who uses the pseudonym Alva Noto for his musical works, exclusively composed the sound for the animations. The multimedia projects aim at making the discoveries more accessible to the general public, and to mediate scientific results and their reference to reality from an artistic point of view.

The measurement of the expansion history of the Universe from the redshift unknown gravitational wave (GW) sources (dark GW sources) detectable from the network of LIGO-Virgo-KAGRA (LVK) detectors depends on the synergy with the galaxy surveys having accurate redshift measurements over a broad redshift range, large sky coverage, and detectability of fainter galaxies. In this work, we explore the possible synergy of the LVK with the spectroscopic galaxy surveys such as DESI and SPHEREx to measure the cosmological parameters which are related to the cosmic expansion history and the GW bias parameters. We show that by using the three-dimensional spatial cross-correlation between the dark GW sources and the spectroscopic galaxy samples, we can measure the value of Hubble constant with about $2\%$ and $1.5\%$ precision from LVK+DESI and LVK+SPHEREx respectively from the five years of observation with $50\%$ duty-cycle for the GW merger rates driven by the star formation history. Similarly, the dark energy equation of state can be measured with about $10\%$ and $8\%$ precision from LVK+DESI and LVK+SPHEREx respectively. We find that due to the larger sky coverage of SPHEREx than DESI, the performance in constraining the cosmological parameters is better from the former than the latter. By combining Euclid along with DESI, and SPHEREx a marginal gain in the measurability of the cosmological parameters is possible from the sources at high redshift ($z\geq 0.9$).

Up to 2020, the Chandra ACIS gain has been calibrated using the External Calibration Source (ECS). The ECS consists of an Fe-55 radioactive source and is placed in the ACIS housing such that all chips are fully illuminated. Since the radioactive source decays over time with a half-life of 2.7 years, count rates are becoming too low for gain calibration. Instead, astrophysical calibration sources will be needed, which do not fill and illuminate the entire field of view. Here, we determine the dominant spatial components of the gain maps through principal component analysis (PCA). We find that, given the noise levels observed today, all ACIS gain maps can be sufficiently described by just a few (often only one) spatial components. We conclude that illuminating a small area is sufficient for gain calibration. We apply this to observations of the astrophysical source Cassiopeia A. The resulting calibration is found to be accurate to 0.6% in at least 68% of the chip area, following the same definition for the calibration accuracy that has been used since launch.

Aims. We provide new insights into the nature of HESS J1857+026, a very-high-energy {\gamma}-ray source whose complex morphology in the TeV band was attributed to the superposition of two distinct sources. Methods. We performed radio continuum observations to look for the pulsar wind nebula and the supernova remnant associated with the pulsar PSR J1856+0245, which might be powering part of the {\gamma}-ray emission. We observed HESS J1857+026 with the Karl G. Jansky Very Large Array (VLA) at 1.5 GHz in the C configuration. In addition, using the same array configuration, we observed a region of $0.4^\circ \times 0.4^\circ$ towards PSR J1856+0245 at 6.0 GHz. We obtained complementary data for the neutral hydrogen and molecular gas emission from public surveys in order to investigate the properties of the interstellar medium in the direction of HESS J1857+026. Results. The new observations at 1.5GHz do not show evidence of emission above the noise level of $0.7\,\rm{mJy\,beam^{-1}}$ that could be associated with either HESS J1857+026 or PSR J1856+0245. Also, in the new image at 6.0GHz we do not detect radio emission from a pulsar wind nebula powered by PSR J1856+0245. The neutral gas analysis shows the existence of a superbubble in the direction of the {\gamma}-ray source. We suggest that this structure is located at $\sim5.5\,\rm{kpc}$, compatible with the distance to the pulsar PSR J1856+0245. Conclusions. We conclude that TeV emission from HESS J1857+026 originates in a superbubble, arguing in favour of a single {\gamma}-ray source rather than the superposition of two distinct sources. The pulsar PSR J1856+0245 could also be contributing as a source of {\gamma}-rays within the bubble.

We analyze the slow periodicities identified in burst sequences from FRB 121102 and FRB 180916 with periods of about 16 and 160 d, respectively, while also addressing the absence of any fast periodicity that might be associated with the spin of an underlying compact object. Both phenomena can be accounted for by a young, highly magnetized, precessing neutron star that emits beamed radiation with significant imposed phase jitter. Sporadic narrow-beam emission into an overall wide solid angle can account for the necessary phase jitter, but the slow periodicities with 25 to 55\% duty cycles constrain beam traversals to be significantly smaller. Instead, phase jitter may result from variable emission altitudes that yield large retardation and aberration delays. A detailed arrival-time analysis for triaxial precession includes wobble of the radio beam and the likely larger, cyclical torque resulting from the changes in the spin-magnetic moment angle. These effects will confound identification of the fast periodicity in sparse data sets longer than about a quarter of a precession cycle unless fitted for and removed as with orbital fitting. Stochastic spin noise, likely to be much larger than in radio pulsars, may hinder detection of any fast-periodicity in data spans longer than a few days. These decoherence effects will dissipate as FRB sources age, so they may evolve into objects with properties similar to Galactic magnetars.

Differential flows among different ion species are often observed in the solar wind, and such ion differential flows can provide the free energy to drive Alfv\'en/ion-cyclotron and fast-magnetosonic/whistler instabilities. Previous works mainly focused on the ion beam instability under the parameters representative of the solar wind nearby 1 au. In this paper we further study the proton beam instability using the radial models of the magnetic field and plasma parameters in the inner heliosphere. We explore a comprehensive distribution of the proton beam instability as functions of the heliocentric distance and the beam speed. We also perform a detailed analysis of the energy transfer between unstable waves and particles and quantify how much the free energy of the proton beam flows into unstable waves and other kinds of particle species (i.e., proton core, alpha particle and electron). This work clarifies that both parallel and perpendicular electric field are responsible for the excitation of oblique Alfv\'en/ion-cyclotron and oblique fast-magnetosonic/whistler instabilities. Moreover, this work proposes an effective growth length to estimate whether the instability is efficiently excited or not. It shows that the oblique Alfv\'en/ion-cyclotron instability, oblique fast-magnetosonic/whistler instability and oblique Alfv\'en/ion-beam instability can be efficiently driven by proton beams drifting at the speed $\sim 600-1300$ km/s in the solar atmosphere. In particular, oblique Alfv\'en/ion-cyclotron waves driven in the solar atmosphere can be significantly damped therein, leading to the solar corona heating. These results are helpful for understanding the proton beam dynamics in the inner heliosphere and can be verified through in situ satellite measurements.

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5<b<+5 deg. and 30<l<215 deg. From the beginning, the incorporation of narrowband H-alpha imaging has been a unique and distinctive feature of this effort. Alongside a focused discussion of the nature and application of the H-alpha data, we present the IGAPS world-accessible database of images for all 5 survey filters, i, r, g, U-RGO and narrowband H-alpha, observed on a pixel scale of 0.33 arcsec and at an effective (median) angular resolution of 1.1 to 1.3 arcsec. The background, noise, and sensitivity characteristics of the narrowband H-alpha filter images are outlined. Typical noise levels in this band correspond to a surface brightness at full one-arcsec resolution of around 2e-16 erg/cm2/s/arcsec2. Illustrative applications of the H-alpha data to planetary nebulae and Herbig-Haro objects are outlined and, as part of a discussion of mosaicking technique, we present a very large background-subtracted narrowband mosaic of the supernova remnant, Simeis 147. Finally we lay out a method that exploits the database via an automated selection of bright ionized diffuse interstellar emission targets for the coming generation of wide-field massive-multiplex spectrographs. Two examples of the diffuse H-alpha maps output from this selection process are presented and compared with previously published data.

Competition between decay and re-interaction of charged pions and kaons depends on the temperature/density profile of the upper atmosphere. The amplitude and phase of the variations depend on the minimum muon energy required to reach the detector and on muon multiplicity in the detector. Here we compare different methods for characterizing the muon production profile and the corresponding effective temperature. A muon production profile based on a parameterization of simulations of muons as a function of primary energy is compared with approximate analytic solutions of the cascade equation integrated over primary energy. In both cases, we compare two definitions of effective temperature. We emphasize applications to compact underground detectors like MINOS and OPERA, while indicating how they relate to extended detectors like IceCube.

To aid the understanding of the non-linear relationship between galaxy properties and predicted spectral energy distributions (SED), we present a new interactive graphical user interface (GUI) tool pipes_vis based on Bagpipes \citep{arXiv:1712.04452,arXiv:1903.11082}. It allows for real-time manipulation of a model galaxy's star formation history, dust and other relevant properties through sliders and text boxes, with each change's effect on the predicted SED reflected instantaneously. We hope the tool will assist in building intuition about what affects the SED of galaxies, potentially helping to speed up fitting stages such as prior construction, and aid in undergraduate and graduate teaching. pipes_vis is available online (pipes_vis is maintained and documented online at https://github.com/HinLeung622/pipes_vis, or version 0.4.1 is archived in Zenodo and also available for installation through pip install pipes_vis).

According to radiative models, radio galaxies and quasars are predicted to produce gamma rays from the earliest stages of their evolution. Exploring their high-energy emission is crucial for providing information on the most energetic processes, the origin and the structure of the newly born radio jets. Taking advantage of more than 11 years of \textit{Fermi}-LAT data, we investigate the gamma-ray emission of 162 young radio sources (103 galaxies and 59 quasars), the largest sample of young radio sources used so far for such a gamma-ray study. We separately analyze each source and perform the first stacking analysis of this class of sources to investigate the gamma-ray emission of the undetected sources. We detect significant gamma-ray emission from 11 young radio sources, four galaxies and seven quasars, including the discovery of significant gamma-ray emission from the compact radio galaxy PKS 1007+142 (z=0.213). The cumulative signal of below-threshold young radio sources is not significantly detected. However, it is about one order of magnitude below than those derived from the individual sources, providing stringent upper limits on the gamma-ray emission from young radio galaxies ($F_{\gamma}< 4.6 \times 10^{-11}$ ph cm$^{-2}$ s$^{-1}$) and quasars ($F_{\gamma}< 10.1 \times 10^{-11}$ ph cm$^{-2}$ s$^{-1}$), and enabling a comparison with the models proposed. With this analysis of more than a decade of \textit{Fermi}-LAT observations, we can conclude that while individual young radio sources can be bright gamma-ray emitters, the collective gamma-ray emission of this class of sources is not bright enough to be detected by \textit{Fermi}-LAT.

Magnetic flares create hot relativistic shocks outside the light cylinder radius of a magnetised star. Radio emission produced in such a shock or at a radius smaller than the shock undergoes free-free absorption while passing through the shocked medium. In this work, we demonstrate that this free-free absorption can lead to a negative drift in the frequency-time spectra. Whether it is related to the downward drift pattern observed in fast radio bursts (FRBs) is unclear. However, if the FRB down drifting is due to this mechanism then it will be pronounced in those shocks that have isotropic kinetic energies $\gtrsim10^{44}$ erg. In this model, for an internal shock with a Lorentz factor $\sim100$, the normalised drift rate $|{\rm DR_{\rm obs}}|/\nu_{\rm mean}$ is $\sim10^{-2}$ per ms, where $\nu_{\rm mean}$ is the central frequency of the radio pulses. The corresponding radius of the shocked shell is, therefore, in the range of $10^{10}$ cm and $10^{11}$ cm. This implies that, for an outflow consisting of hydrogen ion, the upper limit on the mass of the relativistic shocks is a few $\times~10^{-10}~M_\odot$, which is considerably low compared to that ejected from SGR 1806-20 during the 2004 outburst.

The late-time modifications of the standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmological model can be parameterized by three time-dependent functions describing the expansion history of the Universe and gravitational effects on light and matter in the Large Scale Structure. In this Letter, we present the first joint reconstruction of these three functions performed in a non-parametric way from a combination of recent cosmological observations. The reconstruction is performed with a theory-informed prior, built on the general Horndeski class of scalar-tensor theories. We find that current data can constrain 15 combined modes of these three functions with respect to the prior. Our methodology enables us to identify the phenomenological features that alternative theories would need to have in order to ease some of the tensions between datasets within $\Lambda$CDM, and deduce important constraints on broad classes of modified gravity models.

We present a novel scenario for the growth of dust grains in galaxies at high-redshift ($z\sim 6$). In our model, the mechanical feedback from massive star clusters evolving within high-density pre-enriched media allows to pile-up a large amount of matter into massive supershells. If the gas metallicity ($\geq$ Z$_{\odot}$), number density ($\geq 10^6$ cm$^{-3}$) and dust-to-gas mass ratio ($\sim 1/150 \times Z$) within the supershell are sufficiently large, such supershells may become optically thick to the starlight emerging from their host star clusters and even to radiation from the Cosmic Microwave Background (CMB). Based on semi-analytic models, we argue that this mechanism, occurring in the case of massive ($\geq 10^7$ M$_{\odot}$) molecular clouds hosting $\geq 10^6$ M$_{\odot}$ star clusters, allows a large mass of gas and dust to acquire a temperature below that of the CMB, whereupon dust grain growth may occur with ease. In galaxies with total stellar mass $M_{*}$, grain growth within supershells may increase the dust mass by $\sim 10^6$ M$_{\odot}$ $(M_{*}/10^{8}$ M$_{\odot}$).

In type IIB Fibre Inflation models the inflaton is a Kaehler modulus which is kinetically coupled to the corresponding axion. In this setup the curvature of the field space induces tachyonic isocurvature perturbations normal to the background inflationary trajectory. However we argue that the associated instability is unphysical since it is due to the use of ill-defined entropy variables. In fact, upon using the correct relative entropy perturbation, we show that in Fibre Inflation axionic isocurvature perturbations decay during inflation and the dynamics is essentially single-field.

Axion-Like Particles (ALPs) coupled with electrons would be produced in a Supernova (SN) via electron-proton bremsstrahlung and electron-positron fusion. We evaluate the ALP emissivity from these processes by taking into account the ALP mass and thermal effects on electrons in the strongly degenerate and relativistic SN plasma. Using a state-of-the-art SN simulation, we evaluate the SN 1987A cooling bound on ALPs for masses in the range $1-200$ MeV, which excludes currently unprobed regions down to $g_{ae}\sim 2.5\times 10^{-10}$ at $m_a\sim 120$ MeV.

We study static neutron stars in the context of a class of non-minimally coupled inflationary potentials, the universal attractors. Universal attractors are known to generate a viable inflationary era, and they fall into the same category of inflationary phenomenology as the $R^2$ model and other well-known cosmological attractors. We present the essential features of universal attractors in both the Einstein and Jordan frame, and we extract the Tolman-Oppenheimer-Volkoff equations in the Einstein frame using the usual notation of theoretical astrophysics. We use a python 3 based double shooting numerical code for our numerical analysis and we construct the $M-R$ graphs for the universal attractor potential, using piecewise polytropic equation of state the small density part of which is the WFF1 or the APR or the SLy equation of state. As we show, all the studied cases predict larger maximum masses for the neutron stars, and all the results are compatible with the GW170817 constraints imposed on the radii of the neutron stars.

We review analytical solutions of the Einstein equations which are expressed in terms of elementary functions and describe Friedmann-Lema\^itre-Robertson-Walker universes sourced by multiple (real or effective) perfect fluids with constant equations of state. Effective fluids include spatial curvature, the cosmological constant, and scalar fields. We provide a description with unified notation, explicit and parametric forms of the solutions, and relations between different expressions present in the literature. Interesting solutions from a modern point of view include interacting fluids and scalar fields. Old solutions, integrability conditions, and solution methods keep being rediscovered, which motivates a review with modern eyes.

We demonstrate strong negative electrothermal feedback accelerating and linearizing the response of a thermal kinetic inductance detector (TKID). TKIDs are a proposed highly multiplexable replacement to transition-edge sensors and measure power through the temperature-dependent resonant frequency of a superconducting microresonator bolometer. At high readout probe power and probe frequency detuned from the TKID resonant frequency, we observe electrothermal feedback loop gain up to $\mathcal L$ $\approx$ 16 through measuring the reduction of settling time. We also show that the detector response has no detectable non-linearity over a 38% range of incident power and that the noise-equivalent power is below the design photon noise.

Since gravitational and electromagnetic waves from a compact binary coalescence carry independent information about the source, the joint observation is important for understanding the physical mechanisms of the emissions. Rapid detection and source localization of a gravitational wave signal are crucial for the joint observation to be successful. For a signal with a high signal-to-noise ratio, it is even possible to detect it before the merger, which is called early warning. In this letter, we estimate the performances of the early warning for neutron-star black-hole binaries, considering the precession effect of a binary orbit, with the near-future detectors such as A+, AdV+, KAGRA+, and Voyager. We find that a gravitational wave source can be localized in $100 \,\mathrm{deg^2}$ on the sky before $\sim 10$--$40 \,\mathrm{s}$ of time to merger once per year.

The effect of the electroweak sphaleron transition in balance between baryon excess and and the excess of stable quarks of 4th generation is studied in this paper. Considering the non-violation of $SU(2)$ symmetry and the conservation of electroweak and new charges and quantum numbers of the new family, it makes possible sphaleron transitions between baryons, leptons and 4th family of leptons and quarks. In this paper, we have tried to established a possible definite relationship between the value and sign of the 4th family excess relative to baryon asymmetry. If $U$-type quarks are the lightest quarks of the 4th family and sphaleron transitions provide excessive $\bar U$ antiquarks, asymmetric dark matter in the form of dark atom bound state of ($\bar{U} \bar{U} \bar{U}$) with primordial He nuclei is balanced with baryon asymmetry.

We calculate the production of ultra-light axion-like particles (ALPs) in a nearby supernova progenitor. Once produced, ALPs escape from the star and a part of them is converted into photons during propagation in the Galactic magnetic field. It is found that the MeV photon flux that reaches Earth may be detectable by gamma ray telescopes for ALPs lighter than ~1 neV when Betelgeuse undergoes oxygen and silicon burning. (Non-)detection of gamma rays from a supernova progenitor with next-generation gamma ray telescopes just after pre-supernova neutrino alerts would lead to an independent constraint on ALP parameters as stringent as a SN 1987A limit.

We study the mass density distribution of Newtonian self-gravitating systems. Modeling the system as a fluid in hydrostatical equilibrium, we obtain from first principle the field equation and its solution of correlation function $\xi(r)$ of the mass density fluctuation itself. We apply thid to studies of the large-scale structure of the Universe within a small redshift range. The equation tells that $\xi(r)$ depends on the point mass $m$ and the Jeans wavelength scale $\lambda_{0}$, which are different for galaxies and clusters. It explains several longstanding, prominent features of the observed clustering : that the profile of $\xi_{cc}(r)$ of clusters is similar to $\xi_{gg}(r)$ of galaxies but with a higher amplitude and a longer correlation length, and that the correlation length increases with the mean separation between clusters as a universal scaling $r_0\simeq 0.4d$. Our solution $\xi(r)$ also yields the observed power-law correlation function of galaxies $\xi_{gg}(r)\simeq (r_0/r)^{1.7}$ valid only in a range $1<r<10 h^{-1}$Mpc. At larger scales the solution $\xi(r)$ breaks below the power law and goes to zero around $\sim 50h^{-1}$Mpc, just as the observational data have demonstrated. With a set of fixed model parameters, the solutions $\xi_{gg}(r)$ for galaxies, the corresponding power spectrum, and $\xi_{cc}(r)$ for clusters, simultaneously, agree with the observational data from the major surveys of galaxies, and of clusters.

We present an application of anomaly detection techniques based on deep recurrent autoencoders to the problem of detecting gravitational wave signals in laser interferometers. Trained on noise data, this class of algorithms could detect signals using an unsupervised strategy, i.e., without targeting a specific kind of source. We develop a custom architecture to analyze the data from two interferometers. We compare the obtained performance to that obtained with other autoencoder architectures and with a convolutional classifier. The unsupervised nature of the proposed strategy comes with a cost in terms of accuracy, when compared to more traditional supervised techniques. On the other hand, there is a qualitative gain in generalizing the experimental sensitivity beyond the ensemble of pre-computed signal templates. The recurrent autoencoder outperforms other autoencoders based on different architectures. The class of recurrent autoencoders presented in this paper could complement the search strategy employed for gravitational wave detection and extend the reach of the ongoing detection campaigns.

The KM3NeT Collaboration is constructing two deep-sea Cherenkov detectors in the Mediterranean Sea. The ARCA detector aims at TeV-PeV neutrino astronomy, while the ORCA detector is optimised for atmospheric neutrino oscillation studies at energies of a few GeV. In this contribution, an analysis of the data collected with the first deployed detection units of the ARCA detector is presented. A high-purity sample of atmospheric neutrinos is selected demonstrating the capability of the ARCA detector.

In this paper we develop the formalism for the stochastic approach to inflation at all order in slow-roll parameters. This is done by including the momentum and Hamiltonian constraints into the stochastic equations. We then specialise to the widely used Starobinski approximation where interactions between IR and UV modes are neglected. We show that, whenever this approximation holds, no significant deviations are observed when comparing the two-point correlation functions (power spectrum) calculated with stochastic methods, to the ones calculated with the QFT approach to linear theory. As a byproduct, we argue that: a) the approaches based on the Starobinski approximation, generically, do not capture any loop effects of the quantum scalar-gravity system; b) correlations functions can only be calculated in the linear theory regimes, thus, no non-perturbative statistics can be extracted within this approximation, as commonly claimed.

We study an inevitable cosmological consequence in PeV scale SUSY-breaking scenarios. We focus on the SUSY-breaking scale corresponding to the gravitino mass $m_{3/2}=100{\rm eV}-1{\rm keV}$. We argue that the presence of an early matter-dominated era and the resulting entropy production are requisite for the Universe with this gravitino mass. We infer the model-independent minimum amount of the entropy production $\Delta$ by requiring that the number of dwarf satellite galaxies $N_{\rm sat}$ in the Milky Way exceed the currently observed value, i.e. $N_{\rm sat}\gtrsim63$. This entropy production is inevitably imprinted on the primordial gravitational waves (pGWs) produced during the inflationary era. We study how the information on the value of $\Delta$ and the time of entropy production are encoded in the pGW spectrum $\Omega_{\rm GW}$. If the future GW surveys observe a suppression feature in the pGW spectrum for the frequency range $\mathcal{O}(10^{-10}){\rm Hz}\lesssim f_{\rm GW}\lesssim\mathcal{O}(10^{-5}){\rm Hz}$, it works as a smoking gun for PeV SUSY-breaking scenarios. Even if they do not, our study can be used to rule out all such scenarios.

In relativistic magnetized plasmas, asymmetry in the number densities of left- and right-handed fermions, i.e., a non-zero chiral chemical potential mu_5, leads to an electric current along the magnetic field. This causes a chiral dynamo instability for a uniform mu_5, but our simulations reveal dynamos even for fluctuating mu_5 with zero mean. This generates small-scale magnetic helicity and turbulence. A large-scale mu_5 emerges due to chirality conservation. These effects amplify a mean magnetic field via the magnetic alpha effect and produce a universal scale-invariant mu_5 spectrum.