Papers for Friday, May 14 2021

Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model FSPS, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model UNIVERSEMACHINE, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that UNIVERSEMACHINE predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with percent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.

One of the main predictions of excursion set theory is that the clustering of dark matter haloes only depends on halo mass. However, it has been long established that the clustering of haloes also depends on other properties, including formation time, concentration, and spin; this effect is commonly known as halo assembly bias. We use a suite of gravity-only simulations to study the dependence of halo assembly bias on cosmology; these simulations cover cosmological parameters spanning 10$\sigma$ around state-of-the-art best-fitting values, including standard extensions of the $\Lambda$CDM paradigm such as neutrino mass and dynamical dark energy. We find that the strength of halo assembly bias presents variations smaller than 0.05 dex across all cosmologies studied for concentration and spin selected haloes, letting us conclude that the dependence of halo assembly bias upon cosmology is negligible. We then study the dependence of galaxy assembly bias (i.e. the manifestation of halo assembly bias in galaxy clustering) on cosmology using subhalo abundance matching. We find that galaxy assembly bias also presents very small dependence upon cosmology ($\sim$ 2$\%$-4$\%$ of the total clustering); on the other hand, we find that the dependence of this signal on the galaxy formation parameters of our galaxy model is much stronger. Taken together, these results let us conclude that the dependence of halo and galaxy assembly bias on cosmology is practically negligible.

The concentration - virial mass relation is a well-defined trend that reflects the formation of structure in an expanding Universe. Numerical simulations reveal a marked correlation that depends on the collapse time of dark matter halos and their subsequent assembly history. However, observational constraints are mostly limited to the massive end via X-ray emission of the hot diffuse gas in clusters. An alternative approach, based on gravitational lensing over galaxy scales, revealed an intriguingly high concentration at Milky Way-sized halos. This letter focuses on the robustness of these results by adopting a bootstrapping approach that combines stellar and lensing mass profiles. We also apply the identical methodology to simulated halos from EAGLE to assess any systematic. We bypass several shortcomings of ensemble type lens reconstruction and conclude that the mismatch between observed and simulated concentration-to-virial-mass relations are robust, and need to be explained either invoking a lensing-related sample selection bias, or a careful investigation of the evolution of concentration with assembly history. For reference, at a halo mass of $10^{12} M_\odot$, the concentration of observed lenses is $c_{12}\sim 40\pm 5$, whereas simulations give $c_{12}\sim 15\pm1$.

Context: The energy-limited (EL) atmospheric escape approach is used to estimate mass-loss rates for a broad range of planets that host hydrogen-dominated atmospheres as well as for performing atmospheric evolution calculations. Aims: We aim to study the applicability range of the EL approximation. Methods: We revise the EL formalism and its assumptions. We also compare its results with those of hydrodynamic simulations, employing a grid covering planets with masses, radii, and equilibrium temperatures ranging between 1 $M_{\oplus}$ and 39 $M_{\oplus}$, 1 $R_{\oplus}$ and 10 $R_{\oplus}$, and 300 and 2000 K, respectively. Results: Within the grid boundaries, we find that the EL approximation gives a correct order of magnitude estimate for mass-loss rates for about 76% of the planets, but there can be departures from hydrodynamic simulations by up to three orders of magnitude in individual cases. Furthermore, we find that planets for which the mass-loss rates are correctly estimated by the EL approximation to within one order of magnitude have intermediate gravitational potentials as well as low-to-intermediate equilibrium temperatures and irradiation fluxes of extreme ultraviolet and X-ray radiation. However, for planets with low or high gravitational potentials, or high equilibrium temperatures and irradiation fluxes, the approximation fails in most cases. Conclusions: The EL approximation should not be used for planetary evolution calculations that require computing mass-loss rates for planets that cover a broad parameter space. In this case, it is very likely that the EL approximation would at times return mass-loss rates of up to several orders of magnitude above or below those predicted by hydrodynamic simulations. For planetary atmospheric evolution calculations, interpolation routines or approximations based on grids of hydrodynamic models should be used instead.

We present a new empirical model for the mass assembly of dark matter halos. We approximate the growth of individual halos as a simple power-law function of time, where the power-law index smoothly decreases as the halo transitions from the fast-accretion regime at early times, to the slow-accretion regime at late times. Using large samples of halo merger trees taken from high-resolution cosmological simulations, we demonstrate that our 3-parameter model can approximate halo growth with a typical accuracy of 0.1 dex for t > 1 Gyr for all halos of present-day mass greater than 10^11Msun, including subhalos and host halos in gravity-only simulations, as well as in the TNG hydrodynamical simulation. We additionally present a new model for the assembly of halo populations, which not only reproduces average mass growth across time, but also faithfully captures the diversity with which halos assemble their mass. Our python implementation is based on the autodiff library JAX, and so our model self-consistently captures the mean and variance of halo mass accretion rate across cosmic time. We show that the connection between halo assembly and the large-scale density field, known as halo assembly bias, is accurately captured by our model, and that residual errors in our approximations to halo assembly history exhibit a negligible residual correlation with the density field. Our publicly available source code can be used to generate Monte Carlo realizations of cosmologically representative halo histories; our differentiable implementation facilitates the incorporation of our model into existing analytical halo model frameworks.

General Relativity predicts that the passage of matter or radiation from an asymmetrically-emitting source should cause a permanent change in the local space-time metric. This phenomenon, called the \emph{gravitational memory effect}, has never been observed, however supernova neutrinos have long been considered a promising avenue for its detection in the future. With the advent of deci-Hertz gravitational wave interferometers, observing the supernova neutrino memory will be possible, with important implications for multimessenger astronomy and for tests of gravity. In this work, we develop a phenomenological (analytical) toy model for the supernova neutrino memory effect, which is overall consistent with the results of numerical simulations. This description is then generalized to several case studies of interest. We find that, for a galactic supernova, the dimensionless strain, $h(t)$, is of order $\sim 10^{-22} - 10^{-21}$, and develops over a typical time scale that varies between $\sim 0.1 - 10$ s, depending on the time-evolution of the anisotropy of the neutrino emission. The characteristic strain, $h_c(f)$, has a maximum at a frequency $f_{max} \sim {\mathcal O}(10^{-1}) - {\mathcal O}(1)$ Hz. The detailed features of the time- and frequency-structure of the memory strain will inform us of the matter dynamics near the collapsed core, and allow to distinguish between different stellar collapse scenarios. Next generation gravitational wave detectors like DECIGO and BBO will be sensitive to the neutrino memory effect for supernovae at typical galactic distances and beyond; with Ultimate DECIGO exceeding a detectability distance of 10 Mpc.

Observations of galactic white dwarfs with Gaia have allowed for unprecedented modeling of white dwarf cooling, resolving core crystallization and sedimentary heating from neutron rich nuclei. These cooling sequences are sensitive to the diffusion coefficients of nuclei in Coulomb plasmas which have order 10\% uncertainty and are often not valid across coupling regimes. Using large scale molecular dynamics simulations we calculate diffusion coefficients at high resolution in the regime relevant for white dwarf modeling. We present a physically motivated law for diffusion with a semi-empirical correction which is accurate at the percent level. Implemented along with linear mixing in stellar evolution codes, this law should reduce the error from diffusion coefficients by an order of magnitude.

Fast radio bursts (FRBs) are bright radio transient events with durations on the order of milliseconds. The majority of FRB sources discovered so far have a single peak, with the exception of a few showing multiple-peaked profiles, the origin of which is unknown. In this work, we show that the strong lensing effect of a point mass or a point mass $+$ external shear on a single-peak FRB can produce double peaks (i.e. lensed images). In particular, the leading peak will always be more magnified and hence brighter than the trailing peak for a point-mass lens model, while the point-mass $+$ external shear lens model can produce a less magnified leading peak. We find that, for a point-mass lens model, the combination of lens mass $M$ and redshift $z_l$ in the form of $M(1+z_l)$ can be directly computed from two observables -- the delayed time $\Delta t$ and the flux ratio of the leading peak to the trailing peak $R$. For a point-mass $+$ external shear lens model, upper and lower limits in $M(1+z_l)$ can also be obtained from $\Delta t$ and $R$ for a given external shear strength. In particular, tighter lens mass constraints can be achieved when the observed $R$ is larger. Lastly, we show the process of constraining lens mass using the observed values of $\Delta t$ and $R$ of two double-peaked FRB sources, i.e. FRB 121002 and FRB 130729, as references, although the double-peaked profiles are not necessarily caused by strong lensing.

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provides valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the ONC. We also re-measure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using HST ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 136 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of ($\sigma_{v_\alpha}$, $\sigma_{v_\delta}$, $\sigma_{v_r}$) = ($1.76\pm0.12$, $2.16^{+0.14}_{-0.15}$, $2.54^{+0.16}_{-0.17}$) km s$^{-1}$. Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction, but not in the line of sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a non-spherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north-south direction following the filament, which has been well studied in previous literature, and an elongation in the line of sight to tangential velocity direction. These 3-D kinematics, coupled with estimates of source masses, will allow future studies to determine the dominant formation mechanism, differentiating between models such as competitive accretion and turbulent fragmentation.

Aditya-L1 is India's first solar mission with Visible Emission Line Coronagraph (VELC) consisting of three spectral channels taking high-resolution spectroscopic observations of the inner corona up to 1.5 Rsun at 5303 \AA, 7892 \AA, and 10747 \AA. In this work, we present the strategy for the slit width optimization for the VELC using synthetic line profiles by taking into account the instrument characteristics and coronal conditions for log(T) varying from 6 to 6.5. The synthetic profiles are convolved with simulated instrumental scattered light and noise to estimate the signal-to-noise ratio (SNR), which will be crucial to design the future observation plans. We find that the optimum slit width for VELC turns out to be 50 microns providing sufficient SNR for observations in different solar conditions. We also analyzed the effect of plasma temperature on the SNR at different heights in the VELC field of view for the optimized slit width. We also studied the expected effect of the presence of a CME on the spectral channel observations. This analysis will help to plan the science observations of VELC in different solar conditions.

Exoplanet imaging missions utilizing an external occulter (starshade) for starlight suppression require precise alignment between the telescope and starshade, necessitating maintenance of the starshade orbit during observations. Differential lateral acceleration between the two spacecraft serves as a proxy for fuel use and number of required interruptions to the observation. Comparison against results from high fidelity simulations of stationkeeping validates the use of this easy-to-compute proxy. Among starshade positions constrained to the surface of a sphere centered about the telescope, minima of differential lateral acceleration lie on a great circle and its corresponding poles. We present a closed expression for telescope to star vectors requiring minimal stationkeeping for observation from a telescope at an arbitrary position

Numerical models have shown that disc dispersal via internal photoevaporation driven by the host star can successfully reproduce the observed pile-up of warm Jupiters near 1-2 au. However, since a range of different mechanisms have been proposed to cause the same feature, clear observational diagnostics of disc dispersal leaving an imprint in the observed distribution of giant planets could help constraining the dominant mechanisms. We aim to assess the impact of disc dispersal via X-ray driven photoevaporation (XPE) onto giant planet separations in order to provide theoretical constraints on the location and size of any possible features related to this process within their observed orbital distribution. For this purpose, we perform a set of 1D population syntheses with varying initial conditions and correlate the gas giants' final parking locations with the X-ray luminosities of their host stars in order to quantify observables of this process within the $a$-$L_\mathrm{x}$-plane of these systems. We find that XPE indeed creates an underdensity of gas giants near the gravitational radius, with corresponding pile-ups in- and/or outside of this location. However, the size and location of these features are strongly dependent on the choice of initial conditions in our model, such as the assumed formation location of the planets. XPE can strongly affect the migration process of giant planets and leave potentially observable signatures within the observed orbital separations of giant planets. However, due to the simplistic approach employed in our model, which lacks a self-consistent treatment of planet formation within an evolving disc, a quantitative analysis of the final planet population orbits is not possible. Our results however strongly motivate future studies to include realistic disc dispersal mechanisms into global planet population synthesis models.

We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev-Zel'dovich-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 dataset. With signal-to-noise of 62 (43) for galaxy (weak lensing) profiles over scales of about $0.2-20h^{-1}$ Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: 1. The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. 2. The full mass profile is also consistent with the simulations; hence it can constrain alternative dark matter models that modify the mass distribution of clusters. 3. The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. This can be used to constrain processes such as quenching and tidal disruption that alter the galaxy distribution inside the halo, and scale-dependent features in the transition regime outside the halo. We measure the dependence of the profile shapes on the galaxy sample, redshift and cluster mass. We extend the Diemer \& Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation and cosmology are discussed.

Six images of IRC+10216 taken by the Hubble Space Telescope at three epochs in 2001, 2011, and 2016 are compared in the rest frame of the central carbon star. An accurate astrometry has been achieved with the help of Gaia Data Release 2. The positions of the carbon star in the individual epochs are determined using its known proper motion, defining the rest frame of the star. In 2016, a local brightness peak with compact and red nature is detected at the stellar position. A comparison of the color maps between 2016 and 2011 epochs reveals that the reddest spot moved along with the star, suggesting a possibility of its being the dusty material surrounding the carbon star. Relatively red, ambient region is distributed in an $\Omega$ shape and well corresponds to the dusty disk previously suggested based on near-infrared polarization observations. In a larger scale, differential proper motion of multiple ring-like pattern in the rest frame of the star is used to derive the average expansion velocity of transverse wind components, resulting in $\sim$ 12.5 km s$^{-1}$ ($d$/123 pc), where $d$ is the distance to IRC+10216. Three dimensional geometry is implied from its comparison with the line-of-sight wind velocity determined from half-widths of submillimeter emission line profiles of abundant molecules. Uneven temporal variations in brightness for different searchlight beams and anisotropic distribution of extended halo are revisited in the context of the stellar light illumination through a porous envelope with postulated longer-term variations for a period of $\lesssim10$ years.

In this work we examine M dwarf rotation rates at a range of ages to establish benchmarks for Mdwarf gyrochronology. This work includes a sample of 713 spectroscopically-classified M0-M8 dwarfs with new rotation rates measured from K2 light curves. We analyzed data and recover rotation rates for 179 of these objects. We add these to rotation rates for members of clusters with known ages (5-700 Myr), as well as objects assumed to have field ages ($>$1 Gyr). We use Gaia DR2 parallax and (G-GRP) photometry to create color-magnitude diagrams to compare objects across samples. We use color-period plots to analyze the period distributions across age, as well as incorporate Halpha equivalent width and tangential velocity where possible to further comment on age dependence. We find that the age of transition from rapid to slow rotation in clusters, which we define as an elbow in the period-color plots, depends on spectral type. Later spectral types transition at older ages: M4 for Praesepe at 700 Myr, one of the oldest clusters for which M dwarf rotation rates have been measured. The transition from active to inactive Halpha equivalent width also occurs at this elbow, as objects transition from rapid rotation to the slowly rotating sequence. Redder or smaller stars remain active at older ages. Finally, using Gaia kinematics we find evidence for rotation stalling for late Ms in the field sample, suggesting the transition happens much later for mid to late-type M dwarfs.

We identified 8 additional stars as members of the Helmi stream (HStr) in the combined GALAH+ DR3 and $Gaia$ EDR3 catalog. By consistently reevaluating claimed members from the literature, we consolidate a sample of 22 HStr stars with parameters determined from high-resolution spectroscopy and spanning a considerably wider (by $\sim$0.5 dex) metallicity interval ($-2.5 \lesssim \rm[Fe/H] < -1.0$) than previously reported. Our study focuses on $\alpha$ (Mg and Ca) and neutron-capture (Ba and Eu) elements. We find that the chemistry of HStr is typical of dwarf spheroidal (dSph) galaxies, in good agreement with previous $N$-body simulations of this merging event. Stars of HStr constitute a clear declining sequence in $\rm[\alpha/Fe]$ for increasing metallicity up to $\rm[Fe/H] \sim -1.0$. Moreover, stars of HStr show a median value of $+$0.5 dex for $\rm[Eu/Fe]$ with a small dispersion ($\pm$0.1 dex). Every star analyzed with $\rm[Fe/H] < -1.2$ belong to the $r$-process enhanced ($\rm[Eu/Fe] > +0.3$ and $\rm[Ba/Eu] < 0.0$) metal-poor category, providing remarkable evidence that, at such low-metallicity regime, stars of HStr experienced enrichment in neutron-capture elements predominantly via $r$-process nucleosynthesis. Finally, the extended metallicity range also suggests an increase in $\rm[Ba/Eu]$ for higher $\rm[Fe/H]$, in conformity with other surviving dwarf satellite galaxies of the Milky Way.

The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) different ML techniques, also on equal footing, to optimize the immense Big Data parameter space created by these many predictors; third, the establishment of a robust, three-pronged communication effort oriented toward policy makers, space-weather stakeholders and the wider public. FLARECAST pledged to make all its data, codes and infrastructure openly available worldwide. The combined use of 170+ properties (a total of 209 predictors are now available) in multiple ML algorithms, some of which were designed exclusively for the project, gave rise to changing sets of best-performing predictors for the forecasting of different flaring levels. At the same time, FLARECAST reaffirmed the importance of rigorous training and testing practices to avoid overly optimistic pre-operational prediction performance. In addition, the project has (a) tested new and revisited physically intuitive flare predictors and (b) provided meaningful clues toward the transition from flares to eruptive flares, namely, events associated with coronal mass ejections (CMEs). These leads, along with the FLARECAST data, algorithms and infrastructure, could help facilitate integrated space-weather forecasting efforts that take steps to avoid effort duplication. In spite of being one of the most intensive and systematic flare forecasting efforts to-date, FLARECAST has not managed to convincingly lift the barrier of stochasticity in solar flare occurrence and forecasting: solar flare prediction thus remains inherently probabilistic.

HD 106906 is a young, binary stellar system, located in the Lower Centaurus Crux (LCC) group. This system is unique among discovered systems in that it contains an asymmetrical debris disk, as well as an 11 M$_{Jup}$ planet companion, at a separation of $\sim$735 AU. Only a handful of other systems are known to contain both a disk and directly imaged planet, where HD 106906 is the only one in which the planet has apparently been scattered. The debris disk is nearly edge on, and extends roughly to $>$500 AU, where previous studies with HST have shown the outer regions to have high asymmetry. To better understand the structure and composition of the disk, we have performed a deep polarimetric study of HD 106906's asymmetrical debris disk using newly obtained $H$-, $J$-, and $K1$-band polarimetric data from the Gemini Planet Imager (GPI). An empirical analysis of our data supports a disk that is asymmetrical in surface brightness and structure, where fitting an inclined ring model to the disk spine suggests that the disk may be highly eccentric ($e\gtrsim0.16$). A comparison of the disk flux with the stellar flux in each band suggests a blue color that also does not significantly vary across the disk. We discuss these results in terms of possible sources of asymmetry, where we find that dynamical interaction with the planet companion, HD 106906b, is a likely candidate.

The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. 23 meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as a HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of about 156 cm diameter with high bulk density about 2.85 g/cm3, a relatively low albedo pV about 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of about 0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth-impacting orbit via the nu_6 resonance. The impact that ejected 2018 LA in an orbit towards Earth occurred 22.8 +/- 3.8 Ma ago. Zircons record a concordant U-Pb age of 4563 +/- 11 Ma and a consistent 207Pb/206Pb age of 4563 +/- 6 Ma. A much younger Pb-Pb phosphate resetting age of 4234 +/- 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin.

Black hole low mass X-ray binaries in their hard spectral state are found to display two different correlations between the radio emission from the compact jets and the X-ray emission from the inner accretion flow. Here, we present a large data set of quasi-simultaneous radio and X-ray observations of the recently discovered accreting black hole MAXI J1348-630 during its 2019/2020 outburst. Our results span almost six orders of magnitude in X-ray luminosity, allowing us to probe the accretion-ejection coupling from the brightest to the faintest phases of the outburst. We find that MAXI J1348-630 belongs to the growing population of outliers at the highest observed luminosities. Interestingly, MAXI J1348-630 deviates from the outlier track at $L_{\rm X} \lesssim 7 \times 10^{35} (D / 2.2 \ {\rm kpc})^2$ erg s$^{-1}$ and ultimately rejoins the standard track at $L_{\rm X} \simeq 10^{33} (D / 2.2 \ {\rm kpc})^2$ erg s$^{-1}$, displaying a hybrid radio/X-ray correlation, observed only in a handful of sources. However, for MAXI J1348-630 these transitions happen at luminosities much lower than what observed for similar sources (at least an order of magnitude). We discuss the behaviour of MAXI J1348-630 in light of the currently proposed scenarios and we highlight the importance of future deep monitorings of hybrid correlation sources, especially close to the transitions and in the low luminosity regime.

We present results from neutral atomic hydrogen (HI) observations of Hydra I, the first cluster observed by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. For the first time we show that WALLABY can reach its final survey sensitivity. Leveraging the sensitivity, spatial resolution and wide field of view of WALLABY, we identify a galaxy, ESO 501-G075, that lies near the virial radius of Hydra I and displays an HI tail. ESO 501-G075 shows a similar level of morphological asymmetry as another cluster member, which lies near the cluster centre and shows signs of experiencing ram pressure. We investigate possible environmental processes that could be responsible for producing the observed disturbance in the HI morphology of ESO 501-G075. We rule out tidal interactions, as ESO 501-G075 has no nearby neighbours within $\sim0.34$Mpc. We use a simple model to determine that ram pressure can remove gas from the disc at radii $r\gtrsim25$kpc. We conclude that, as ESO 501-G075 has a typical HI mass compared to similar galaxies in the field and its morphology is compatible with a ram pressure scenario, ESO 501-G075 is likely recently infalling into the cluster and in the early stages of experiencing ram pressure.

The streaming instability (SI) is a mechanism to aerodynamically concentrate solids in protoplanetary disks and trigger the formation of planetesimals. The SI produces strong particle clumping if the ratio of solid to gas surface density -- an effective metallicity -- exceeds a critical value. This critical value depends on particle sizes and disk conditions such as radial drift-inducing pressure gradients and levels of turbulence. To quantify these thresholds, we perform a suite of vertically-stratified SI simulations over a range of dust sizes and metallicities. We find a critical metallicity as low as 0.4% for the optimum particle sizes and standard radial pressure gradients (normalized value of $\Pi = 0.05$). This sub-Solar metallicity is lower than previous results due to improved numerical methods and computational effort. We discover a sharp increase in the critical metallicity for small solids, when the dimensionless stopping time (Stokes number) is $\leq 0.01$. We provide simple fits to the size-dependent SI clumping threshold, including generalizations to different disk models and levels of turbulence. We also find that linear, unstratified SI growth rates are a surprisingly poor predictor of particle clumping in non-linear, stratified simulations, especially when the finite resolution of simulations is considered. Our results widen the parameter space for the SI to trigger planetesimal formation.

We studied the absorption features of CO lines against the continuum originating from the heated dust in the obscuring tori around active galactic nuclei (AGNs). We investigated the formation of absorption lines corresponding to the CO rotational transitions using three-dimensional non-LTE line transfer simulations considering the dust thermal emission. As in Papers I--III of this series, we performed post-processed radiative transfer calculations using the "radiation-driven fountain model" (wada2016}, which yields a geometrically thick obscuring structure around the nucleus. This model is consistent with the spectral energy distribution of the nearest type-2 Seyfert galaxy, the Circinus galaxy. We found that the continuum-subtracted channel maps of $J = 4-3$ and higher transitions show absorption regions along the disk mid-plane for an edge-on viewing angle. The spectra consist of multiple absorption and emission features, reflecting the internal inhomogeneous and turbulent structure of the torus. The deepest absorption feature is caused by the gas on the near-side of the torus between $r =10$ and 15 pc, which is located in front of the AGN-heated dust inside $r \simeq 5$ pc. We also found that a spatial resolution of 0.5--1.0 pc is necessary to resolve the absorption features. Moreover, the inclination angle must be close to the edge-on angle (i.e., $\sim 85^\circ$) to observe the absorption features. The findings of the present study imply that combining our radiation-hydrodynamic model with high-resolution observations of CO (7-6) by ALMA can provide new information about the internal structure of the molecular tori in nearby AGNs.

Chondrules are often surrounded by fine-grained rims or igneous rims. The properties of these rims reflect their formation histories. While the formation of fine-grained rims is modeled by the accretion of dust grains onto chondrules, the accretion should be followed by the growth of dust grains due to the shorter growth timescale than the accretion. In this paper, we investigate the formation of rims, taking into account the growth of porous dust aggregates. We estimate the rim thickness as a function of the chondrule fraction at a time when dust aggregate accretion onto chondrules is switched to collisions between these chondrules. Our estimations are consistent with the measured thicknesses of fine-grained rims in ordinary chondrites. However, those of igneous rims are thicker than our estimations. The thickness of igneous rims would be enlarged in remelting events.

The Galileo mission to Jupiter discovered magnetic signatures associated with hidden sub-surface oceans at the moons Europa and Callisto using the phenomenon of magnetic induction. These induced magnetic fields originate from electrically conductive layers within the moons and are driven by Jupiter's strong time-varying magnetic field. The ice giants and their moons are also ideal laboratories for magnetic induction studies. Both Uranus and Neptune have a strongly tilted magnetic axis with respect to their spin axis, creating a dynamic and strongly variable magnetic field environment at the orbits of their major moons. Although Voyager 2 visited the ice giants in the 1980s, it did not pass close enough to any of the moons to detect magnetic induction signatures. However, Voyager 2 revealed that some of these moons exhibit surface features that hint at recent geologically activity, possibly associated with sub-surface oceans. Future missions to the ice giants may therefore be capable of discovering sub-surface oceans, thereby adding to the family of known ocean worlds in our solar system. Here, we assess magnetic induction as a technique for investigating sub-surface oceans within the major moons of Uranus. Furthermore, we establish the ability to distinguish induction responses created by different interior characteristics that tie into the induction response: ocean thickness, conductivity, and depth, and ionospheric conductance. The results reported here demonstrate the possibility of single-pass ocean detection and constrained characterization within the moons of Miranda, Ariel, and Umbriel, and provide guidance for magnetometer selection and trajectory design for future missions to Uranus.

Interstellar gas heating is a powerful cosmology-independent observable for exploring the parameter space of primordial black holes (PBHs) formed in the early Universe that could constitute part of the dark matter (DM). We provide a detailed analysis of the various aspects for this observable, such as PBH emission mechanisms. Using observational data from the Leo T dwarf galaxy, we constrain the PBH abundance over a broad mass-range, $M_{\rm PBH} \sim \mathcal{O}(1) M_{\odot}-10^7 M_{\odot}$, relevant for the recently detected gravitational wave signals from intermediate-mass BHs. We also consider PBH gas heating of systems with bulk relative velocity with respect to the DM, such as Galactic clouds.

We utilize the robust membership determination algorithm, ML-MOC, on the precise astrometric and deep photometric data from Gaia Early Data Release 3 within a region of radius 5$^{\circ}$ around the center of the intermediate-age galactic open cluster NGC 752 to identify its member stars. We report the discovery of the tidal tails of NGC 752, extending out to $\sim$35 pc on either side of its denser central region and following the cluster orbit. From comparison with PARSEC stellar isochrones, we obtain the mass function of the cluster with a slope, $\chi=-1.26\pm0.07$. The high negative value of $\chi$ is indicative of a disintegrating cluster undergoing mass-segregation. $\chi$ is more negative in the intra-tidal regions as compared to the outskirts of NGC 752. We estimate a present day mass of the cluster, M$\rm_{C}=297\pm10$ M$_{\odot}$. Through mass-loss due to stellar evolution and tidal interactions, we further estimate that NGC 752 has lost nearly 95.2-98.5 % of its initial mass, $\rm M_{i}~=~0.64~-2~\times~10^{4}~M_{\odot}$.

The statistical characterization of the distribution of visible matter in the universe is a central problem in modern cosmology. In this respect, a crucial question still lacking a definitive answer concerns how large are the greatest structures in the universe. This point is closely related to whether or not such a distribution can be approximated as being homogeneous on large enough scales. Here we assess this problem by considering the size distribution of superclusters of galaxies and by leveraging on the properties of Zipf-Mandelbrot law, providing a novel approach which complements standard analysis based on the correlation functions. We find that galaxy superclusters are well described by a pure Zipf's law with no deviations and this implies that all the catalogs currently available are not sufficiently large to spot a truncation in the power-law behavior. This finding provides evidence that structures larger than the greatest superclusters already observed are expected to be found when deeper redshift surveys will be completed. As a consequence the scale beyond which galaxy distribution crossovers toward homogeneity, if any, should increase accordingly

Radio interferometer arrays with non-homogeneous element patterns are more difficult to calibrate compared to the more common homogeneous array. In particular, the non-homogeneity of the patterns has significant implications on the computational tractability of evaluating the calibration solutions. We apply the A-stacking technique to this problem and explore the trade-off to be made between the calibration accuracy and computational complexity. Through simulations, we show that this technique can be favourably applied in the context of an SKA-Low station. We show that the minimum accuracy requirements can be met at a significantly reduced computational cost, and this cost can be reduced even further if the station calibration timescale is relaxed from 10 minutes to several hours. We demonstrate the impact antenna designs with differing levels of non-homogeneity have on the overall computational complexity in addition to some cases where calibration performs poorly.

High-energy cosmic rays are observed indirectly by detecting the extensive air showers initiated in Earth's atmosphere. The interpretation of these observations relies on accurate models of air shower physics, which is a challenge and an opportunity to test QCD under extreme conditions. Air showers are hadronic cascades, which eventually decay into muons. The muon number is a key observable to infer the mass composition of cosmic rays. Air shower simulations with state-of-the-art QCD models show a significant muon deficit with respect to measurements; this is called the Muon Puzzle. The origin of this discrepancy has been traced to the composition of secondary particles in hadronic interactions. The muon discrepancy starts at the TeV scale, which suggests that this change in hadron composition is observable at the Large Hadron Collider. An effect that can potentially explain the puzzle has been discovered at the LHC, but needs to be confirmed for forward produced hadrons with LHCb, and with future data on oxygen beams.

Context. It is now well-established that small, rocky planets are common around low-mass stars. However, the detection of such planets is challenged by the short-term activity of the host stars. Aims. The HArps-N red Dwarf Exoplanet Survey (HADES) program is a long-term project at the Telescopio Nazionale Galileo aimed at the monitoring of nearby, early-type, M dwarfs, using the HARPS-N spectrograph to search for small, rocky planets. Methods. A total of 174 HARPS-N spectroscopic observations of the M0.5V-type star GJ 9689 taken over the past seven years have been analysed. We combined these data with photometric measurements to disentangle signals related to the stellar activity of the star from possible Keplerian signals in the radial velocity data. We run an MCMC analysis, applying Gaussian Process regression techniques to model the signals present in the data. Results. We identify two periodic signals in the radial velocity time series, with periods of 18.27 d, and 39.31 d. The analysis of the activity indexes, photometric data, and wavelength dependency of the signals reveals that the 39.31 d signal corresponds to the stellar rotation period. On the other hand, the 18.27 d signal shows no relation to any activity proxy or the first harmonic of the rotation period. We, therefore, identify it as a genuine Keplerian signal. The best-fit model describing the newly found planet, GJ 9689 b, corresponds to an period P$_{\rm b}$ = 18.27 $\pm$ 0.01 d, and a minimum mass M$_{\rm P}\sin i$ = 9.65 $\pm$ 1.41 M$_{\oplus}$.

Finding and characterising extrasolar Earth analogues will rely on interpretation of the planetary system's environmental context. The total budget and fractionation between C-H-O species sensitively affect the climatic and geodynamic state of terrestrial worlds, but their main delivery channels are poorly constrained. We connect numerical models of volatile chemistry and pebble coagulation in the circumstellar disk with the internal compositional evolution of planetesimals during the primary accretion phase. Our simulations demonstrate that disk chemistry and degassing from planetesimals operate on comparable timescales and can fractionate the relative abundances of major water and carbon carriers by orders of magnitude. As a result, individual planetary systems with significant planetesimal processing display increased correlation in the volatile budget of planetary building blocks relative to no internal heating. Planetesimal processing in a subset of systems increases the variance of volatile contents across planetary systems. Our simulations thus suggest that exoplanetary atmospheric compositions may provide constraints on $when$ a specific planet formed.

Tensions in cosmological parameters measurement motivate a revisit of the effects of instrumental systematics. In this article, we focus on the Pearson's correlation coefficient of the cosmic microwave background temperature and polarization E modes $\mathcal{R}_\ell^{\rm TE}$ which has the property of not being biased by multiplicative instrumental systematics. We build a $\mathcal{R}_\ell^{\rm TE}$-based likelihood for the Planck data, and present the first constraints on $\Lambda$CDM parameters from the correlation coefficient. Our results are compatible with parameters derived from a power spectra based likelihood. In particular the value of the Hubble parameter $H_0$ characterizing the expansion of the Universe today, 67.5 $\pm$ 1.3 km/s/Mpc, is consistent with the ones inferred from standard CMB analysis. We also discuss the consistency of the Planck correlation coefficient with the one computed from the most recent ACTPol power spectra.

It is generally believed that high energy radiation (power-law components) can be mostly produced by a hot corona gas in the accreting black holes. There is a very popular hybrid disk radial coupling model that the inner part of cool Keplerian disk (or Shakura-Sunyaev disk) can produce advection-dominated accretion flow or corona-like structure, which can also generate outflows/jets. Here we argue that this simple coupling model cannot explain the whole hardness-intensity diagram of the spectral states and their transitions, and associated jets of a $X-$ray binary. Based on recent theoretical works on advective disk structures around a black hole, as well as many observational behaviors of a source, we conclude that there should be a third component of hot accretion flow with the radial coupling model, which can successfully explain all the spectral states and transitions. Interestingly, this model also provides a new scenario for the jet generation, launching, and evolution during the states with help of created barrier of the inner flow. We have also find out the jet kinetic power expression with our new jet generation scenario.

Future surveys such as the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will observe an order of magnitude more astrophysical transient events than any previous survey before. With this deluge of photometric data, it will be impossible for all such events to be classified by humans alone. Recent efforts have sought to leverage machine learning methods to tackle the challenge of astronomical transient classification, with ever improving success. Transformers are a recently developed deep learning architecture, first proposed for natural language processing, that have shown a great deal of recent success. In this work we develop a new transformer architecture, which uses multi-head self attention at its core, for general multi-variate time-series data. Furthermore, the proposed time-series transformer architecture supports the inclusion of an arbitrary number of additional features, while also offering interpretability. We apply the time-series transformer to the task of photometric classification, minimising the reliance of expert domain knowledge for feature selection, while achieving results comparable to state-of-the-art photometric classification methods. We achieve a weighted logarithmic-loss of 0.507 on imbalanced data in a representative setting using data from the Photometric LSST Astronomical Time-Series Classification Challenge (PLAsTiCC). Moreover, we achieve a micro-averaged receiver operating characteristic area under curve of 0.98 and micro-averaged precision-recall area under curve of 0.87.

Spectroscopic observations of 32 HII regions in the spiral galaxy NGC 3963 and the barred irregular galaxy NGC 7292 were carried out with the 2.5-m telescope of the Caucasus Mountain Observatory of the Sternberg Astronomical Institute using the Transient Double-beam Spectrograph with a dispersion of 1A/pixel and a spectral resolution of 3A. These observations were used to estimate the oxygen and nitrogen abundances and the electron temperatures in HII regions through modern strong-line methods. In general, the galaxies have oxygen and nitrogen abundances typical of stellar systems with similar luminosities, sizes, and morphology. However, we have found some peculiarities in chemical abundance distributions in both galaxies. The distorted outer segment of the southern arm of NGC 3963 shows an excess oxygen and nitrogen abundances. Chemical elements abundances in NGC 7292 are constant and do not depend on the galactocentric distance. These peculiarities can be explained in terms of external gas accretion in the case of NGC 3963 and major merging for NGC 7292.

New results from long-term optical photometric observations of the pre-main sequence star GM Cep from UX Orionis type are reported. During ongoing photometric monitoring of the GM Cep four deep minimums in brightness are observed. The collected multicolour photometric data shows the typical of UXor variables colour reversal during the minimums in brightness. Recent $BVRI$ photometric observations of GM Cep have been collected from November 2014 to October 2020.

We present optical and near-infrared (NIR, $YJH$-band) observations of 42 Type Ia supernovae (SNe Ia) discovered by the untargeted intermediate Palomar Transient Factory (iPTF) survey. This new data-set covers a broad range of redshifts and host galaxy stellar masses, compared to previous SN Ia efforts in the NIR. We construct a sample, using also literature data at optical and NIR wavelengths, to examine claimed correlations between the host stellar masses and the Hubble diagram residuals. The SN magnitudes are corrected for host galaxy extinction using either a global total-to-selective extinction ratio, $R_V$=2.0 for all SNe, or a best-fit $R_V$ for each SN individually. Unlike previous studies which were based on a narrower range in host stellar mass, we do not find evidence for a "mass-step", between the color- and stretch-corrected peak $J$ and $H$ magnitudes for galaxies below and above $\log(M_{*}/M_{\odot}) = 10$. However, the mass-step remains significant ($3\sigma$) at optical wavelengths ($g,r,i$) when using a global $R_V$, but vanishes when each SN is corrected using their individual best-fit $R_V$. Our study confirms the benefits of the NIR SN Ia distance estimates, as these are largely exempted from the empirical corrections dominating the systematic uncertainties in the optical.

The turbulent dynamics of nearby and extragalactic gas structures can be studied with the column density power spectrum, which is often described by a broken power-law.In an extragalactic context, the breaks in the power spectra have been interpreted to constrain the disc scale height, which marks a transition from 2D disc-like to 3D motion. However, this interpretation has recently been questioned when accounting for instrumental effects. We use numerical simulations to study the spatial power spectra of isolated galaxies and investigate the origins of the break scale. We split the gas into various phases and analyze the time evolution of the power spectrum characteristics, such as the slope(s) and the break scale. We find that the break scale is phase dependent. The physics traced by the break scale also differ: in the warm gas it marks the transition from 2D (disk-like) to 3D (isotropic) turbulence. In the cold gas, the break scale traces the typical size of molecular clouds. We further show that the break scale almost never traces the disc scale height. We study turbulent properties of the ISM to show that, in the case where the break scale traces a transition to isotropic turbulence, the fraction of required accretion energy to sustain turbulent motions in the ISM increases significantly. Lastly, we demonstrate through simulated observations that it is crucial to account for observational effects, such as the beam and instrumental noise, in order to accurately recover the break scale in real observations.

The motion data of the S-stars around the Galactic center gathered in the last 28 yr imply that Sgr A* hosts a supermassive compact object of about $4\times 10^6$ $M\odot$, a result awarded with the Nobel Prize in Physics 2020. A non-rotating black hole (BH) nature of Sgr A* has been uncritically adopted since the S-star orbits agree with Schwarzschild geometry geodesics. The orbit of S2 has served as a test of General Relativity predictions such as the gravitational redshift and the relativistic precession. The central BH model is, however, challenged by the G2 post-peripassage motion and by the lack of observations on event-horizon-scale distances robustly pointing to its univocal presence. We have recently shown that the S2 and G2 astrometry data are better fitted by geodesics in the spacetime of a self-gravitating dark matter (DM) core - halo distribution of 56 keV-fermions, "darkinos", which also explains the outer halo Galactic rotation curves. This Letter confirms and extends this conclusion using the astrometry data of the $17$ best-resolved S-stars, thereby strengthening the alternative nature of Sgr A* as a dense core of darkinos.

The high-redshift quasar PMN J0909+0354 ($z=3.288$) is known to have a pc-scale compact jet structure, based on global 5-GHz very long baseline interferometry (VLBI) observations performed in 1992. Its kpc-scale structure was studied with the Karl G. Jansky Very Large Array (VLA) in the radio and the Chandra space telescope in X-rays. Apart from the north-northwestern jet component seen in both the VLA and Chandra images at $2.3''$ separation from the core, there is another X-ray feature at $6.48''$ in the northeastern (NE) direction. To uncover more details and possibly structural changes in the inner jet, we conducted new observations at 5 GHz using the European VLBI Network (EVN) in 2019. These data confirm the northward direction of the one-sided inner jet already suspected from the 1992 observations. A compact core and multiple jet components were identified that can be traced up to $\sim0.25$ kpc projected distance towards the north, while the structure becomes more and more diffuse. A comparison with arcsec-resolution imaging with the VLA shows that the radio jet bends by $\sim30^\circ$ between the two scales. The direction of the pc-scale jet as well as the faint optical counterpart found for the newly-detected X-ray point source (NE) favors the nature of the latter as a background or foreground object in the field of view. However, the extended ($\sim160$ kpc) emission around the positions of the quasar core and NE detected by the Wide-field Infrared Survey Explorer (WISE) in the mid-infrared might suggest physical interaction of the two objects.

Contamination by polarized foregrounds is one of the biggest challenges for future polarized Cosmic Microwave Background (CMB) surveys and the potential detection of primordial $B$-modes. Future experiments, such as Simons Observatory (SO) and CMB-S4, will aim at very deep observations in relatively small ($f_{\rm sky} \sim 0.1$) areas of the sky. In this work, we investigate the forecasted performance, as a function of the survey field location on the sky, for regions over the full sky, balancing between polarized foreground avoidance and foreground component separation modeling needs. To do this, we simulate observations by a SO-like experiment, and measure the error bar on the detection of the tensor-to-scalar ratio, $\sigma(r)$, with a pipeline that includes a parametric component separation method, Correlated Component Analysis (CCA), and the use of the Fisher information matrix. We forecast the performance over 192 survey areas covering the full sky and also for optimized low-foreground regions. We find that modeling the Spectral Energy Distribution (SED) of foregrounds is the most important factor, and any mismatch will result in residuals and bias in the primordial $B$-modes. At these noise levels, $\sigma(r)$ is not especially sensitive to the level of foreground contamination, provided the survey targets the least contaminated regions of the sky close to the Galactic Poles.

Thermal phase variations of short period planets indicate that they are not spherical cows: day-to-night temperature contrasts range from hundreds to thousands of degrees, rivaling their vertical temperature contrasts. Nonetheless, the emergent spectra of short-period planets have typically been fit using one-dimensional (1D) spectral retrieval codes that only account for vertical temperature gradients. The popularity of 1D spectral retrieval codes is easy to understand: they are robust and have a rich legacy in Solar System atmospheric studies. Exoplanet researchers have recently introduced multi-dimensional retrieval schemes for interpreting the spectra of short-period planets, but these codes are necessarily more complex and computationally expensive than their 1D counterparts. In this paper we present an alternative: phase-dependent spectral observations are inverted to produce longitudinally resolved spectra that can then be fitted using standard 1D spectral retrieval codes. We test this scheme on the iconic phase-resolved spectra of WASP-43b and on simulated JWST observations using the open-source pyratbay 1D spectral retrieval framework. Notably, we take the model complexity of the simulations one step further over previous studies by allowing for longitudinal variations in composition in addition to temperature. We show that performing 1D spectral retrieval on longitudinally resolved spectra is more accurate than applying 1D spectral retrieval codes to disk-integrated emission spectra, despite being identical in terms of computational load. We find that for the extant Hubble and Spitzer observations of WASP-43b the difference between the two approaches is negligible but that JWST phase measurements should be treated with longitudinally \textbf{re}solved \textbf{spect}ral retrieval (ReSpect).

Uranus provides a unique laboratory to test our understanding of planetary atmospheres under extreme conditions. Multi-spectral observations from Voyager, ground-based observatories, and space telescopes have revealed a delicately banded atmosphere punctuated by storms, waves, and dark vortices, evolving slowly under the seasonal influence of Uranus' extreme axial tilt. Condensables like methane and hydrogen sulphide play a crucial role in shaping circulation, clouds, and storm phenomena via latent heat release through condensation, strong equator-to-pole gradients suggestive of equatorial upwelling and polar subsidence, and through forming stabilising layers that may decouple different circulation and convective regimes as a function of depth. Weak vertical mixing and low atmospheric temperatures associated with Uranus' negligible internal heat means that stratospheric methane photochemistry occurs in a unique high-pressure regime, decoupled from the influx of external oxygen. The low homopause also allows for the formation of an extensive ionosphere. Finally, the atmosphere provides a window on the bulk composition of Uranus - the ice-to-rock ratio, supersolar elemental and isotopic enrichments inferred from remote sensing and future \textit{in situ} measurements - providing key insights into its formation and subsequent migration. This review reveals the state of our knowledge of the time-variable circulation, composition, meteorology, chemistry, and clouds on this enigmatic `Ice Giant,' summarising insights from more than three decades of observations, and highlighting key questions for the next generation of planetary missions. As a hydrogen-dominated, intermediate-sized, and chemically-enriched world, Uranus could be our closest and best example of atmospheric processes on a class of worlds that may dominate the census of planets beyond our own Solar System.

In this letter, we present an update of a method for analysing arrival directions of ultra-high-energy cosmic rays (UHECRs) above the Greisen--Zatsepin--Kuz'min cut-off with a deep convolutional neural network developed originally in Kalashev, Pshirkov, Zotov (2020). Namely, we introduce energy as another variable employed in the analysis. This allows us to take into account the intrinsic uncertainties in energy of primary cosmic rays present in any experiment, which were not taken into account in the previous study, without any loss of quality of the classifier. We present the architecture of the new neural network, results of its application to mock maps of UHECR arrival directions and outline possible directions of a further improvement of the method.

I explore a triple-star scenario where a tight neutron star (NS) - NS binary system enters the envelope of a red supergiant (RSG) star and spirals-in towards its core. The two NSs accrete mass through accretion disks and launch jets that power a very luminous and long transient event, a common envelope jets supernova (CEJSN) event. Dynamical friction brings the two NSs to merge either in the RSG envelope or core. The total energy of the event, radiation and kinetic, is >10^{52}erg. The light curve stays luminous for months to years and a signal of gravitational waves might be detected. The ejecta contains freshly synthesized r-process elements not only from the NS-NS merger as in kilonova events, but possibly also from the pre-merger jets that the NSs launch inside the core, as in the r-process CEJSN scenario. This scenario shortens the time to NS-NS merger compared with that of kilonovae, and might somewhat ease the problem of the NS-NS r-process scenario to account for r-process nucleosynthesis in the early Universe. I estimate the ratio of NS-NS merger in CEJSN events to core collapse supernova (CCSN) events to be <10^{-6}-2x10^{-5}. However, because they are much more luminous I expect their detection fraction to that of CCSNe to be much larger than this number. This study calls for considering this and similar CEJSN scenarios in binary and in triple star systems when explaining peculiar and puzzling super luminous supernovae.

According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest stream in each halo as a function of varying limiting magnitude. We find that none of the models shows signatures of stellar streams at $\mu_{r}^{lim} \leq 25$ mag arcsec$^{-2}$. The stream detection increases significantly between 27 and 28 mag arcsec$^{-2}$. Nevertheless, even at 30 mag arcsec$^{-2}$, 13 percent of our models show no detectable streams. We study the properties of the brightest streams progenitors (BSPs). We find that BSPs are accreted within a broad range of infall times, from 1.6 to 10 Gyr ago, with only 25 percent accreted within the last 5 Gyrs; thus most BSPs correspond to relatively early accretion events. We also find that 37 percent of the BSPs survive to the present day. The median infall times for surviving and disrupted BSPs are 5.6 and 6.7 Gyr, respectively. We find a clear relation between infall time and infall mass of the BSPs, such that more massive progenitors tend to be accreted at later times. However, we find that the BSPs are not, in most cases, the dominant contributor to the accreted stellar halo of each galaxy.

Multiple scalar fields appear in vast modern particle physics and gravity models. When they couple to gravity non-minimally, conformal transformation is utilized to bring the theory into Einstein frame. However, the kinetic terms of scalar fields are usually not canonical, which makes analytic treatment difficult. Here we investigate under what conditions the theories can be transformed to the quasi-canonical form, in which case the effective metric tensor in field space is conformally flat. We solve the relevant nonlinear partial differential equations for arbitrary number of scalar fields and present several solutions that may be useful for future phenomenological model building, including the $\sigma$-model with a particular non-minimal coupling. We also find conformal flatness can always be achieved in some modified gravity theories, for example, Starobinsky model.

The quantum to classical transition of fluctuations in the early universe is still not completely understood. Some headway has been made incorporating the effects of decoherence and the squeezing of states, though the methods and procedures continue to be challenged. But new developments in the analysis of the most recent Planck data suggest that the primordial power spectrum has a cutoff associated with the very first quantum fluctuation to have emerged into the semi-classical universe from the Planck domain at about the Planck time. In this paper, we examine the implications of this result on the question of classicalization, and demonstrate that the birth of quantum fluctuations at the Planck scale would have been a `process' supplanting the need for a `measurement' in quantum mechanics. Emerging with a single wavenumber, these fluctuations would have avoided the interference between different degrees of freedom in a superposed state. Moreover, the implied scalar-field potential had an equation-of-state consistent with the zero active mass condition in general relativity, allowing the quantum fluctuations to emerge in their ground state with a time-independent frequency. They were therefore effectively quantum harmonic oscillators with classical correlations in phase space from the very beginning.

Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform ongoing efforts to construct an accurate model of the general-relativistic binary dynamics via an effective-one-body description and other semi-analytical approaches. Existing self-force methods are to a large extent specialised to bound, inspiral orbits. Here we develop a technique for (numerical) self-force calculations that can efficiently tackle scatter orbits. The method is based on a time-domain reconstruction of the metric perturbation from a scalar-like Hertz potential that satisfies the Teukolsky equation, an idea pursued so far only for bound orbits. The crucial ingredient in this formulation are certain jump conditions that (each multipole mode of) the Hertz potential must satisfy along the orbit, in a 1+1-dimensional multipole reduction of the problem. We obtain a closed-form expression for these jumps, for an arbitrary geodesic orbit in Schwarzschild spacetime, and present a full numerical implementation for a scatter orbit. In this paper we focus on method development, and go only as far as calculating the Hertz potential; a calculation of the self-force and its physical effects on the scatter orbit will be the subject of forthcoming work.

We examine an intriguing possibility that a single field is responsible for both inflation and dark matter, focussing on the minimal set--up where inflation is driven by a scalar coupling to curvature. We study in detail the reheating process in this framework, which amounts mainly to particle production in a quartic potential, and distinguish thermal and non--thermal dark matter options. In the non--thermal case, the reheating is impeded by backreaction and rescattering, making this possibility unrealistic. On the other hand, thermalized dark matter is viable, yet the unitarity bound forces the inflaton mass into a narrow window close to half the Higgs mass.

We study primordial cosmology with two scalar fields that participate in inflation at the same time, by coupling quantum gravity (i.e., the theory $R+R^{2}+C^{2}$ with the fakeon prescription/projection for $C^{2}$) to a scalar field with a quadratic potential. We show that there exists a perturbative regime that can be described by an asymptotically de Sitter, cosmic RG flow in two couplings. Since the two scalar degrees of freedom mix in nontrivial ways, the adiabatic and isocurvature perturbations are not RG invariant on superhorizon scales. It is possible to identify the correct perturbations by using RG invariance as a guiding principle. We work out the resulting power spectra of the tensor and scalar perturbations to the NNLL and NLL orders, respectively. An unexpected consequence of RG invariance is that the theory remains predictive. Indeed, the scalar mixing affects only the subleading corrections, so the predictions of quantum gravity with single-field inflation are confirmed to the leading order.

In classical general relativity astrophysical black holes can be affected by the superradiant instability when gravity is minimally coupled to a light bosonic field. The majority of phenomenological studies have focused on the idealized case in which the black hole is initially surrounded by a single mode superradiant seed. By studying the evolution of a scalar field with multiple modes initial data in a quasiadiabatic approximation, we show that the dynamics is more involved and depend on the initial seed energy and the amplitude ratio between the modes. We also present preliminary results of the dynamical evolution of a massive scalar field around a Newtonian and a fully relativistic, black hole binary.

We explore the tail of various waiting time datasets of processes that follow a nonstationary Poisson distribution with a sinusoidal driver. Analytically, we find that the distribution of large waiting times of such processes can be described using a power law slope of -2.5. We show that this result applies more broadly to any nonstationary Poisson process driven periodically. Examples of such processes include solar flares, coronal mass ejections, geomagnetic storms, and substorms. We also discuss how the power law specifically relates to the behavior of driver near its minima.

Direct observations of gravitational waves at frequencies around deci-Hertz will play a crucial role in fully exploiting the potential of multi-messenger astronomy. In this chapter, we discuss the detection landscape for the next several decades of the deci-Hertz gravitational-wave spectrum. We provide an overview of the experimental frontiers being considered to probe this challenging regime and the astrophysics and fundamental goals accessible towards them. This includes interferometric observatories in space with heliocentric and geocentric satellites, cubesats, lunar-based experiments and atom intereferometry. A major focus of this chapter is towards the technology behind DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) and its scientific pathfinder mission concept B-DECIGO.

We revisit a supersymmetric string model for space-time foam, in which bosonic open-string states, such as photons, can possess quantum-gravity-induced velocity fluctuations in vacuum. We argue that the suggestion of the light speed variation recently observed from gamma-ray burst photon time delays can serve as a crucial support for this string-inspired framework, through connecting the experimental finding with model predictions. We also derive the value of the effective quantum-gravity mass in this framework, and give a qualitative study on the model-dependent coefficients. Constraints from birefringent effects and/or photon decays, including the novel $\gamma$-decay constraint obtained here from the latest Tibet AS$\gamma$ near-PeV photon, are also found to be consistent with predictions in such a quantum-gravity scheme. Future observation that can testify further the theory is suggested.

Due to the absence of any definite signals of new physics at colliders and from precision measurements, it has gradually become more and more popular in the community to utilize the effective field theory (EFT) framework in searching for new physics in a model-independent manner. In this letter, working in the EFT framework and focusing on neutrino non-standard interactions (NSIs), we report our most recent results on these NSIs from considering terrestrial neutrino oscillation experiments Daya Bay, Double Chooz, RENO, T2K and NO$\nu$A, and precision measurements of $N_{\rm eff}$ from Planck and CMB-S4.

Dark matter Axion search with riNg Cavity Experiment (DANCE) was proposed. To search for axion-like particles, we aim to detect the rotation and oscillation of optical linear polarization caused by axion-photon coupling with a bow-tie cavity. DANCE can improve the sensitivity to axion-photon coupling constant $g_{a \gamma}$ for axion mass $m_a < 10^{-10} \, \rm{eV}$ by several orders of magnitude compared to the best upper limits at present. A prototype experiment DANCE Act-1 is in progress to demonstrate the feasibility of the method and to investigate technical noises. We assembled the optics, evaluated the performance of the cavity, and estimated the current sensitivity. If we observe for a year, we can reach $g_{a \gamma} \simeq 9 \times 10^{-7} \, \rm{GeV^{-1}}$ at $m_a \simeq 10^{-13} \, \rm{eV}$. The current sensitivity was believed to be limited by laser intensity noise at low frequencies and by mechanical vibration at high frequencies.

SU(2) gauge fields coupled to an axion field can acquire an isotropic background solution during inflation. We study homogeneous but anisotropic inflationary solutions in the presence of such (massless) gauge fields. A gauge field in the cosmological background may pose a threat to spatial isotropy. We show, however, that such models $\textit{generally}$ isotropize in Bianchi type-I geometry, and the isotropic solution is the attractor. Restricting the setup by adding an axial symmetry, we revisited the numerical analysis presented in Wolfson et.al (2020). We find that the reported numerical breakdown in the previous analysis is an artifact of parametrization singularity. We use a new parametrization that is well-defined all over the phase space. We show that the system respects the cosmic no-hair conjecture and the anisotropies always dilute away within a few e-folds.

We propose Genesis, a one-armed simplified Convolutional Neural Network (CNN)for exoplanet detection, and compare it to the more complex, two-armed CNN called Astronet. Furthermore, we examine how Monte Carlo cross-validation affects the estimation of the exoplanet detection performance. Finally, we increase the input resolution twofold to assess its effect on performance. The experiments reveal that (i)the reduced complexity of Genesis, i.e., a more than 95% reduction in the number of free parameters, incurs a small performance cost of about 0.5% compared to Astronet, (ii) Monte Carlo cross-validation provides a more realistic performance estimate that is almost 0.7% below the original estimate, and (iii) the twofold increase in input resolution decreases the average performance by about 0.5%. We conclude by arguing that further exploration of shallower CNN architectures may be beneficial in order to improve the generalizability of CNN-based exoplanet detection across surveys.

Dark photon not only provides a portal linking dark sector particles and ordinary matter but also is a well-motivated dark matter candidate. We propose to detect the dark photon dark matter through the inverse Compton-like scattering process $p+\gamma^\prime \to p+\gamma$. Thanks to the ultra-high energy primary cosmic rays, we find that such a method is able to probe the dark photon mass from $10^{-2}$ eV down to $10^{-19}$ eV with the expected sensitivity of eROSITA $X$-ray telescope, which can extend the current lower limit of dark photon mass from Jupiter's magnetic fields experiment by about three orders of magnitude.

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF.

We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced LIGO and Advanced Virgo during O3a, the first half of their third observing run. We study: 1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; 2) how the interpretation of individual high-mass events would change if they were found to be lensed; 3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and 4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses.

Analyses of near-surface air temperature T in Poland for 1781-2016 and in Tbilisi (Georgia) for 1881-2016 have been carried out. We show that the centenary warming effect in Poland and in Tbilisi has almost the same peculiarities. An average centenary warming effect deltaT = (1.08+/-0.29) C is observed in Poland and in Tbilisi for 1881-2016. A warming effect is larger in winter season (deltaT = ~1.15 C) than in other seasons (average warming effect for these seasons deltaT = ~0.95 C). We show that a centenary warming is mainly related to the change of solar activity (estimated by sunspot numbers (SSN) and total solar irradiance (TSI)); particularly, a time interval about ~70 years (1890-1960), when a correlation coefficients between 11 years smoothed SSN and T, and TSI and T are high, r = 0.66+/-0.07 and r = 0.73+/-0.07 for Poland and r = 0.82+/-0.05 and r = 0.90+/-0.05 for Tbilisi, respectively; in this period solar activity contributes decisively in the global warming. We show that a global warming effect equals zero based on the temperature T data in Poland for period 1781-1880, when human activities were relatively less than in 1881-2016. We recognize a few feeble ~20+/-3 years disturbances in the temperature changes for period 1885-1980, most likely related with the fluctuations of solar magnetic cycles. We distinguish the fluctuations of ~7-8 years in Poland's T data, possibly connected with local effects of the North Atlantic Oscillation.

Within Einstein's theory of gravity, any compact object heavier than a few solar masses must be a black hole. Any observation showing otherwise would imply either new physics beyond General Relativity or new exotic matter fields beyond the Standard Model, and might provide a portal to understand some puzzling properties of a black hole. We give a short overview on tests of the nature of dark compact objects with present and future gravitational-wave observations, including inspiral tests of the multipolar structure of compact objects and of their tidal deformability, ringdown tests, and searches for near-horizon structures with gravitational-wave echoes.