Papers for Wednesday, Jun 15 2022

We present a novel method to estimate galaxy physical properties from spectral energy distributions (SEDs), alternate to template fitting techniques and based on self-organizing maps (SOM) to learn the high-dimensional manifold of a photometric galaxy catalog. The method has been previously tested with hydrodynamical simulations in Davidzon et al. (2019) while here is applied to real data for the first time. It is crucial for its implementation to build the SOM with a high quality, panchromatic data set, which we elect to be the "COSMOS2020" galaxy catalog. After the training and calibration steps with COSMOS2020, other galaxies can be processed through SOM to obtain an estimate of their stellar mass and star formation rate (SFR). Both quantities result to be in good agreement with independent measurements derived from more extended photometric baseline, and also their combination (i.e., the SFR vs. stellar mass diagram) shows a main sequence of star forming galaxies consistent with previous studies. We discuss the advantages of this method compared to traditional SED fitting, highlighting the impact of having, instead of the usual synthetic templates, a collection of empirical SEDs built by the SOM in a "data-driven" way. Such an approach also allows, even for extremely large data sets, an efficient visual inspection to identify photometric errors or peculiar galaxy types. Considering in addition the computational speed of this new estimator, we argue that it will play a valuable role in the analysis of oncoming large-area surveys like Euclid or the Legacy Survey of Space and Time at the Vera Cooper Rubin Telescope.

The Ly$\alpha$ emission line is one of the main observables of galaxies at high redshift, but its output depends strongly on the neutral gas distribution and kinematics around the star-forming regions where UV photons are produced. We present observations of Ly$\alpha$ and 21-cm HI emission at comparable scales with the goal to qualitatively investigate how the neutral interstellar medium (ISM) properties impact Ly$\alpha$ transfer in galaxies. We have observed 21-cm HI at the highest angular resolution possible (~ 3" beam) with the VLA in two local galaxies from the Lyman Alpha Reference Sample. We contrast this data with HST Ly$\alpha$ imaging and spectroscopy, and MUSE and PMAS ionized gas observations. In LARS08, high intensity Ly$\alpha$ emission is co-spatial with high column density HI where dust content is the lowest. The Ly$\alpha$ line is strongly redshifted, consistent with velocity redistribution which allows Ly$\alpha$ escape from high column density neutral medium with low dust content. In eLARS01, high intensity Ly$\alpha$ emission is located in regions of low column density HI, below the HI data sensitivity limit ($<2\times10^{20}\,$cm$^{-2}$). The perturbed ISM distribution with low column density gas in front of the Ly$\alpha$ emission region plays an important role in the escape. In both galaxies, the faint Ly$\alpha$ emission ($\sim 1\times10^{-16}$erg.s$^{-1}$cm$^{-2}$arcsec$^{-2}$) traces intermediate H$\alpha$ emission regions where HI is found, regardless of the dust content. Dust seems to modulate, but not prevent, the formation of a faint Ly$\alpha$ halo. This study suggests the existence of scaling relations between dust, H$\alpha$, HI, and Ly$\alpha$ emission in galaxies.

We report new observations of "Hanny's Voorwerp" (hereafter HV) taken from the second data release of the LOFAR Two-metre Sky Survey (LoTSS). HV is a highly-ionised region in the environs of the galaxy IC2497, first discovered by the Galaxy Zoo project. The new 150MHz observations are considered in the context of existing multi-frequency radio data and archival narrow-band imaging from the Hubble Space Telescope, centred on the [Oiii] emission line. The combined sensitivity and spatial resolution of the LoTSS data -- which far exceed what was previously available at radio frequencies -- reveal clear evidence for large-scale extended emission emanating from the nucleus of IC2497. The radio jet appears to have punched a hole in the neutral gas halo, in a region co-located with HV. The new 150MHz data, alongside newly-processed archival 1.64GHz eVLA data, reveal that the extended emission has a steep spectrum, implying an age $>10^8$yr. The jet supplying the extended 150MHz structure must have "turned off" long before the change in X-ray luminosity reported in recent works. In this picture, a combination of jet activity and the influence of the radiatively efficient active galactic nucleus are responsible for the unusual appearance of HV.

The physical properties of galaxies are encoded within their spectral energy distribution and require comparison with models to be extracted. These models must contain a synthetic stellar population and, where infrared data is to be used, also consider prescriptions for energy reprocessing and re-emission by dust. While many such models have been constructed, there are few analyses of the impact of stellar population model choice on derived dust parameters, or vice versa. Here we apply a simple framework to compare the impact of these choices, combining three commonly-used stellar population synthesis models and three dust emission models. We compare fits to the ultraviolet to far-infrared spectral energy distributions of a validation sample of infrared-luminous galaxies. We find that including different physics, such as binary stellar evolution, in the stellar synthesis model can introduce biases and uncertainties in the derived parameters of the dust and stellar emission models, largely due to differences in the far-ultraviolet emission available for reprocessing. This may help to reconcile the discrepancy between the cosmic star formation rate and stellar mass density histories. Notably the inclusion of a dusty stellar birth cloud component in the dust emission model provides more flexibility in accommodating the stellar population model, as its reemission is highly sensitive to the ultraviolet radiation field spectrum and density. Binary populations favour a longer birth cloud dissipation timescale than is found when assuming only single star population synthesis.

We present 3-D hydrodynamical models of the evolution of superbubbles powered by stellar winds and supernovae from young coeval massive star clusters within low metallicity ($Z = 0.02$Z$_{\odot}$), clumpy molecular clouds. We explore the initial stages of the superbubble evolution, including the occurrence of pair-instability and core-collapse supernovae. Our aim is to study the occurrence of dust grain growth within orbiting dusty clumps, and in the superbubble's swept-up supershell. We also aim to address the survival of dust grains produced by sequential supernovae. The model accounts for the star cluster gravitational potential and self-gravity of the parent cloud. It also considers radiative cooling (including that induced by dust) and a state-of-the-art population synthesis model for the coeval cluster. As shown before, a superbubble embedded into a clumpy medium becomes highly distorted, expanding mostly due to the hot gas streaming through low density channels. Our results indicate that in the case of massive ($\sim10^7$M$_{\odot}$) molecular clouds, hosting a super star cluster ($\sim5.6\times10^5$M$_{\odot}$), grain growth increments the dust mass at a rate $\sim4.8\times10^{-5}$M$_{\odot}$ yr$^{-1}$ during the first $2.5$Myr of the superbubble's evolution, while the net contribution of pair-instability and core-collapse supernovae to the superbubble's dust budget is $\sim1200$M$_{\odot} (M_{SC}/5.6\times10^{5}$M$_{\odot})$, where $M_{SC}$ is the stellar mass of the starburst. Therefore, dust grain growth and dust injection by supernovae lead to create, without invoking a top-heavy initial mass function, massive amounts of dust within low-metallicity star-forming molecular clouds, in accordance with the large dust mass present in galaxies soon after the onset of cosmic reionization.

The 21cm line emitted by neutral hydrogen (HI) during the Dark Ages carries imprints of pristine primordial correlations. In models of inflation driven by a single, canonical scalar field, we show that a phase of ultra-slow-roll can lead to a null in all the primordial correlations at a specific wavenumber $k_\textrm{dip}$. We consider scenarios wherein the null in the correlations occurs over wave numbers $1 \lesssim k_\textrm{dip} \lesssim 10\,\mathrm{Mpc}^{-1}$, and examine the prospects of detecting such a damping in the HI signal due to the nulls at the level of power and bi-spectra in future observational missions.

The cosmic formation rate of long Gamma Ray Bursts (LGRBs) encodes the evolution, across cosmic times, of their progenitors' properties and of their environment. The LGRB formation rate and the luminosity function, with its redshift evolution, are derived by reproducing the largest set of observations collected in the last four decades, namely the observer-frame prompt emission properties of GRB samples detected by the Fermi and Compton Gamma Ray Observatory (CGRO) satellites and the redshift, luminosity and energy distributions of flux-limited, redshift complete, samples of GRBs detected by Swift. The model that best reproduces all these constraints consists of a GRB formation rate increasing with redshift $\propto (1+z)^{3.2}$, i.e. steeper than the star formation rate, up to $z\sim3$ followed by a decrease $\propto(1+z)^{-3}$. On top of this, our model predicts also a moderate evolution of the characteristic luminosity function break $\propto(1+z)^{0.6}$. Models with only luminosity or rate evolution are excluded at $>5\sigma$ significance. The cosmic rate evolution of LGRBs is interpreted as their preference to occur in environments with metallicity $12+\log(\rm O/H)<8.6$, consistently with theoretical models and host galaxy observations. The LGRB rate at $z=0$, accounting for their collimation, is $\rho_0=79^{+57}_{-33}$ Gpc$^{-3}$ yr$^{-1}$ (68% confidence interval). This corresponds to $\sim$1\% of broad-line Ibc supernovae producing a successful jet in the local Universe. This fraction increases up to $\sim$7% at $z\ge3$. Finally, we estimate that at least $\approx0.2-0.7$ yr$^{-1}$ of Swift and Fermi detected bursts at $z<0.5$ are jets observed slightly off-axis.

We present a novel method for including the effects of early (pre-supernova) feedback in simulations of galaxy evolution. Rather than building a model which attempts to match idealized, small-scale simulations or analytic approximations, we rely on direct observational measurements of the time-scales over which star-forming molecular clouds are disrupted by early feedback. We combine observations of the spatial de-correlation between molecular gas and star formation tracers on $\sim100$~pc scales with an analytic framework for the expansion of feedback fronts driven by arbitrary sources or mechanisms, and use these to constrain the time-scale and momentum injection rate by early feedback. This allows us to directly inform our model for feedback from these observations, sidestepping the complexity of multiple feedback mechanisms and their interaction below the resolution scale. We demonstrate that this new model has significant effects on the spatial clustering of star formation, the structure of the ISM, and the driving of outflows from the galactic plane, while preserving the overall regulation of the galaxy-integrated star formation rate. We find that this new feedback model results in galaxies that regulate star formation through the rapid disruption of star-forming clouds, rather than by highly efficient, global galactic outflows. We also demonstrate that these results are robust to stochasticity, degraded numerical resolution, changes in the star formation model parameters, and variations in the single free model parameter that is unconstrained by observations.

We analyse synthetic $^{12}$CO, $^{13}$CO, and [CII] emission maps of simulated molecular clouds of the SILCC-Zoom project, which include an on-the-fly evolution of H$_2$, CO, and C$^+$. We use simulations of hydrodynamical and magnetohydrodynamical clouds, both with and without stellar feedback. We introduce a novel post-processing of the C$^+$ abundance using CLOUDY, to account for further ionization states of carbon due to stellar radiation. We report the first self-consistent synthetic emission maps of [CII] in feedback bubbles, largely devoid of emission inside them, as recently found in observations. The C$^+$ mass is only poorly affected by stellar feedback but the [CII] luminosity increases by $50 - 85$ per cent compared to runs without feedback. Furthermore, we investigate the capability of the CO/[CII] line ratio as a tracer of the amount of H$_2$ in the clouds and their evolutionary stage. We obtain, for both $^{12}$CO and $^{13}$CO, no clear trend of the luminosity ratio, $L_\mathrm{CO}/L_\mathrm{[CII]}$. It can therefore \textit{not} be used as a reliable measure of the H$_2$ mass fraction. We note a monotonic relation between $L_\mathrm{CO}/L_\mathrm{[CII]}$ and the H$_2$ fraction when considering the ratio for individual pixels of our synthetic maps, but with large scatter. Moreover, we show that assuming chemical equilibrium results in an overestimation of H$_2$ and CO masses of up to 110 and 30 per cent, respectively, and in an underestimation of H and C$^+$ masses of 65 and 7 per cent, respectively. In consequence, $L_\mathrm{CO}$ would be overestimated by up to 50 per cent, and $L_\mathrm{C[II]}$ be underestimated by up to 35 per cent. Hence, the assumption of chemical equilibrium in molecular cloud simulations introduces intrinsic errors of a factor of up to $\sim2$ in chemical abundances, luminosities and luminosity ratios.

We conducted a spectral and temporal analysis of X-ray data from the Be X-ray binary pulsar SXP 15.6 located in the Small Magellanic Cloud based on NuSTAR, NICER and Swift observations during the 2021 outburst. We present for the first time the broadband X-ray spectra of the system based on simultaneous NuSTAR and NICER observations. Moreover we use monitoring data to study the spectral and temporal properties of the system during the outburst. Comparison of the evolution of the 2021 outburst with archival data reveals a consistent pattern of variability with multiple peaks occurring at time intervals similar to the orbital period of the system (~36 d). Our spectral analysis indicates that most of the energy is released at high energies above 10 keV while we found no cyclotron absorption line in the spectrum. Analysis of the spectral evolution during the outburst, we find that the spectrum is softer-when-brighter, which in turn reveals that the system is probably in the super-critical regime where the accretion column is formed. This places an upper limit to the magnetic field of the system of the order of 7$\times$10$^{11}$ G. The spin-evolution of the neutron star (NS) during the outburst is consistent with an NS with a low magnetic field ($\sim{5}\times$10$^{11}$ G), while there is evident orbital modulation which we modelled and derived the orbital parameters. We found the orbit to have a moderate eccentricity of ~0.3. Our estimates of the magnetic field are consistent with the lack of an electron cyclotron resonance scattering feature in the broadband X-ray spectrum.

We aim to determine the intrinsic far-Infrared (far-IR) emission of X-ray-luminous quasars over cosmic time. Using a 16 deg^2 region of the Stripe 82 field surveyed by XMM-Newton and Herschel Space Observatory, we identify 2905 X-ray luminous (LX > 10^42 erg/s) Active Galactic Nuclei (AGN) in the range z ~ 0-3. The IR is necessary to constrain host galaxy properties such as star formation rate (SFR) and gas mass. However, only 10% of our AGN are detected both in the X-ray and IR. Because 90% of the sample is undetected in the far-IR by Herschel, we explore the mean IR emission of these undetected sources by stacking their Herschel/SPIRE images in bins of X-ray luminosity and redshift. We create stacked spectral energy distributions from the optical to the far-IR, and estimate the median star formation rate, dust mass, stellar mass, and infrared luminosity using a fitting routine. We find that the stacked sources on average have similar SFR/L_bol ratios as IR detected sources. The majority of our sources fall on or above the main sequence line suggesting that X-ray selection alone does not predict the location of a galaxy on the main sequence. We also find that the gas depletion timescales of our AGN are similar to those of dusty star forming galaxies. This suggests that X-ray selected AGN host high star formation and that there are no signs of declining star formation.

We implement a black hole spin evolution and jet feedback model into SWIFT, a smoothed particle hydrodynamics code. The jet power is determined self-consistently assuming Bondi accretion, using a realistic, spin-dependant efficiency. The jets are launched along the spin axis of the black hole, resulting in natural reorientation and precession. We apply the model to idealised simulations of galaxy groups and clusters, finding that jet feedback successfully quenches gas cooling and star formation in all systems. Our group-size halo ($M_\mathrm{200}=10^{13}$ $\mathrm{M}_\odot$) is quenched by a strong jet episode triggered by a cooling flow, and it is kept quenched by a low-power jet fed from hot halo accretion. In more massive systems ($M_\mathrm{200}\geq 10^{14}$ $\mathrm{M}_\odot$), hot halo accretion is insufficient to quench the galaxies, or to keep them quenched after the first cooling episode. These galaxies experience multiple episodes of gas cooling, star formation and jet feedback. In the most massive galaxy cluster that we simulate ($M_\mathrm{200}=10^{15}$ $\mathrm{M}_\odot$), we find peak cold gas masses of $10^{10}$ $\mathrm{M}_\odot$ and peak star formation rates of a few times $100$ $\mathrm{M}_\odot\mathrm{yr}^{-1}$. These values are achieved during strong cooling flows, which also trigger the strongest jets with peak powers of $10^{47}$ $\mathrm{erg}\hspace{0.3mm}\mathrm{s}^{-1}$. These jets subsequently shut off the cooling flows and any associated star formation. Jet-inflated bubbles draw out low-entropy gas that subsequently forms dense cooling filaments in their wakes, as seen in observations.

The role played by recombination during a common-envelope event has long been a subject of debate. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30 per cent of the total recombination energy, and occurs much deeper in the stellar envelope. We investigate the distinct roles played by hydrogen and helium recombination in a common-envelope interaction experienced by a 12 solar mass red supergiant donor. We perform adiabatic, 3D hydrodynamical simulations that (i) include hydrogen, helium, and molecular hydrogen recombination, (ii) include hydrogen and helium recombination, (iii) include only helium recombination, and (iv) do not include recombination energy. By comparing these simulations, we find that the addition of helium recombination energy alone ejects 30 per cent more envelope mass, and leads to a 16 per cent larger post-plunge-in separation. Under the adiabatic assumption, adding hydrogen recombination energy increases the amount of ejected mass by a further 40 per cent, possibly unbinding the entire envelope, but does not affect the post-plunge separation. Most of the ejecta becomes unbound at relatively high (>70 per cent) degrees of hydrogen ionisation, where the hydrogen recombination energy is likely to expand the envelope instead of being radiated away.

Life in the clouds of Venus, if present in sufficiently high abundance, must be affecting the atmospheric chemistry. It has been proposed that abundant Venusian life could obtain energy from its environment using three possible sulfur energy-metabolisms. These metabolisms raise the possibility of Venus's enigmatic cloud-layer SO$_2$-depletion being caused by life. We here couple each proposed energy-metabolism to a photochemical-kinetics code and self-consistently predict the composition of Venus's atmosphere under the scenario that life produces the observed SO$_2$-depletion. Using this photo-bio-chemical kinetics code, we show that all three metabolisms can produce SO$_2$-depletions, but do so by violating other observational constraints on Venus's atmospheric chemistry. We calculate the maximum possible biomass density of sulfur-metabolising life in the clouds, before violating observational constraints, to be $\sim10^{-5}\,-\,10^{-3}\,{\rm mg\,m^{-3}}$. The methods employed are equally applicable to aerial biospheres on Venus-like exoplanets, planets that are optimally poised for atmospheric characterisation in the near future.

Context. Gaia has been in operations since 2014. The third Gaia data release expands from the early data release (EDR3) in 2020 by providing 34 months of multi-epoch observations that allowed us to probe, characterise and classify systematically celestial variable phenomena. Aims. We present a summary of the variability processing and analysis of the photometric and spectroscopic time series of 1.8 billion sources done for Gaia DR3. Methods. We used statistical and Machine Learning methods to characterise and classify the variable sources. Training sets were built from a global revision of major published variable star catalogues. For a subset of classes, specific detailed studies were conducted to confirm their class membership and to derive parameters that are adapted to the peculiarity of the considered class. Results. In total, 10.5 million objects are identified as variable in Gaia DR3 and have associated time series in G, GBP, and GRP and, in some cases, radial velocity time series. The DR3 variable sources subdivide into 9.5 million variable stars and 1 million Active Galactic Nuclei/Quasars. In addition, supervised classification identified 2.5 million galaxies thanks to spurious variability induced by the extent of these objects. The variability analysis output in the DR3 archive amounts to 17 tables containing a total of 365 parameters. We publish 35 types and sub-types of variable objects. For 11 variable types, additional specific object parameters are published. An overview of the estimated completeness and contamination of most variability classes is provided. Conclusions. Thanks to Gaia we present the largest whole-sky variability analysis based on coherent photometric, astrometric, and spectroscopic data. Later Gaia data releases will more than double the span of time series and the number of observations, thus allowing for an even richer catalogue in the future.

To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from astrophysical models of compact objects as seen by distant observers. We incorporate two different schemes, namely Monte Carlo radiative transfer, integrating geodesics from the astrophysical region to distant observers, and camera techniques with backwards integration from the observer to the emission region. The code is validated by successfully passing several test cases, among them: thin accretion disks and neutron star hot spot emission.

W-CDF-S, ELAIS-S1, and XMM-LSS will be three Deep-Drilling Fields (DDFs) of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), but their extensive multi-wavelength data have not been fully utilized as done in the COSMOS field, another LSST DDF. To prepare for future science, we fit source spectral energy distributions (SEDs) from X-ray to far-infrared in these three fields mainly to derive galaxy stellar masses and star-formation rates. We use CIGALE v2022.0, a code that has been regularly developed and evaluated, for the SED fitting. Our catalog includes 0.8 million sources covering $4.9~\mathrm{deg^2}$ in W-CDF-S, 0.8 million sources covering $3.4~\mathrm{deg^2}$ in ELAIS-S1, and 1.2 million sources covering $4.9~\mathrm{deg^2}$ in XMM-LSS. Besides fitting normal galaxies, we also select candidates that may host active galactic nuclei (AGNs) or are experiencing recent star-formation variations and use models specifically designed for these sources to fit their SEDs; this increases the utility of our catalog for various projects in the future. We calibrate our measurements by comparison with those in well-studied smaller regions and briefly discuss the implications of our results. We also perform detailed tests of the completeness and purity of SED-selected AGNs. Our data can be retrieved from a public website.

Now that detection of gravitational wave signals from the coalescence of extra-galactic compact binary star mergers has become nearly routine, it is intriguing to consider other potential gravitational wave signatures. Here we examine the prospects for discovery of continuous gravitational waves from fast-spinning neutron stars in our own galaxy and from more exotic sources. Potential continuous-wave sources are reviewed, search methodologies and results presented and prospects for imminent discovery discussed.

There is a wealth of data on live, undecayed 60Fe ($t_{1/2} = 2.6 \ \rm Myr$) in deep-sea deposits, the lunar regolith, cosmic rays, and Antarctic snow, which is interpreted as originating from the recent explosions of at least two near-Earth supernovae. We use the 60Fe profiles in deep-sea sediments to estimate the timescale of supernova debris deposition beginning $\sim 3$ Myr ago. The available data admits a variety of different profile functions, but in all cases the best-fit 60Fe pulse durations are $> 1.6$ Myr when all of the data is combined. This timescale far exceeds the $\lesssim 0.1$ Myr pulse that would be expected if 60Fe was entrained in supernova blast wave plasma. We interpret the long signal duration as evidence that 60Fe arrives in the form of supernova dust, whose dynamics are separate from but coupled to the evolution of the blast plasma. In this framework, the $> 1.6$ Myr timescale is that for dust stopping due to drag forces. This scenario is consistent with the simulations in Fry et. al (2020), where is magnetically trapped in supernova remnants and thereby confined in and near regions of the remnant dominated by supernova ejecta, where magnetic fields are low. This picture fits naturally with models of cosmic-ray injection of refractory elements as sputtered supernova dust grains, and implies that the recent 60Fe detections in cosmic rays complement the fragments of grains that survived to arrive on the Earth and Moon. Finally, we present possible tests for this scenario.

Rich observational phenomenology associated with Pulsar B in PSR J0737$-$3039A/B system resembles in many respects phenomena observed in the Earth and Jupiter magnetospheres, originating due to the wind-magnetosphere interaction. We consider particle dynamics in the fast corotating magnetosphere of Pulsar B, when the spin period is shorter than the third adiabatic period. We demonstrate that trapped particles occasionally experience large radial variations of the L-parameter (effective radial distance) due to the parametric interaction of the gyration motion with the large scale electric fields induced by the deformations of the magnetosphere, in what could be called a betatron-induced diffusion. The dynamics of particles from the wind of Pulsar A trapped inside Pulsar B magnetosphere is governed by Mathieu's equation, so that the parametrically unstable orbits are occasionally activated; particle dynamics is not diffusive per se. The model explains the high plasma density on the closed field lines of Pulsar B, and the fact that the observed eclipsing region is several times smaller than predicted by the hydrodynamic models.

The revolutionary power of future Rubin-LSST observations will allow us to significantly improve the physics of pulsating stars, including RR Lyrae. In this context, an updated theoretical scenario predicting all the relevant pulsation observables in the corresponding photometric filters is mandatory. The bolometric light curves based on a recently computed extensive set of nonlinear convective pulsation models for RR Lyrae stars, covering a broad range in metal content and transformed into the Rubin-LSST photometric system. Predicted Rubin-LSST mean magnitudes and pulsation amplitudes have been adopted to built the Bailey diagrams (luminosity amplitude vs period) and the color-color diagrams in these bands. The current findings indicate that the gLSST-rLSST, rLSST-iLSST colors obey to a well defined linear relation with the metal content. Moreover, the Period Luminosity relations display in the reddest filters (rLSST,iLSST,zLSST,yLSST) a significant dependence on the assumed metal abundance. In particular, more metal-rich RR Lyrae are predicted to be fainter at fixed period. Metal-dependent Period-Wesenheit relations for different combinations of optical and NIR filters are also provided. These represent powerful tools to infer individual distances independently of reddening uncertainties, once the metal abundance is known and no relevant deviations from the adopted extinction law occur. Finally, we also derived new linear and quadratic absolute magnitude metallicity relations (gLSST vs [Fe/H]) and the metallicity coefficient is consistent with previous findings concerning the B and the V band.

Interchange reconnection is thought to play an important role in determining the dynamics and material composition of the slow solar wind that originates from near coronal hole boundaries. To explore the implications of this process we simulate the dynamic evolution of a solar wind stream along a newly-opened magnetic flux tube. The initial condition is composed of a piecewise continuous dynamic equilibrium in which the regions above and below the reconnection site are extracted from steady-state solutions along open and closed field lines. The initial discontinuity at the reconnection site is highly unstable and evolves as a Riemann problem, decomposing into an outward-propagating shock and inward-propagating rarefaction that eventually develop into a classic N-wave configuration. This configuration ultimately propagates into the heliosphere as a coherent structure and the entire system eventually settles to a quasi-steady wind solution. In addition to simulating the fluid evolution we also calculate the time-dependent non-equilibrium ionization of oxygen in real time in order to construct in situ diagnostics of the conditions near the reconnection site. This idealized description of the plasma dynamics along a newly-opened magnetic field line provides a baseline for predicting and interpreting the implications of interchange reconnection for the slow solar wind. Notably, the density and velocity within the expanding N-wave are generally enhanced over the ambient wind, as is the O7+/O6+ ionization ratio, which exhibits a discontinuity across the reconnection site that is transported by the flow and arrives later than the propagating N-wave.

We present photometric and spectroscopic observations of the Type IIn supernova SN 2019zrk (also known as ZTF20aacbyec). The SN shows a $\gtrsim$ 100 day precursor, with a slow rise, followed by a rapid rise to M $\sim -19.2$ in the $r$ and $g$ bands. The post-peak light-curve decline is well fit with an exponential decay with a timescale of $\sim 39$ days, but it shows prominent undulations, with an amplitude of $\sim 1$ mag. Both the light curve and spectra are dominated by an interaction with a dense circumstellar medium (CSM), probably from previous mass ejections. The spectra evolve from a scattering-dominated Type IIn spectrum to a spectrum with strong P-Cygni absorptions. The expansion velocity is high, $\sim 16,000$ km s$^{-1}$, even in the last spectra. The last spectrum $\sim 110$ days after the main eruption reveals no evidence for advanced nucleosynthesis. From analysis of the spectra and light curves, we estimate the mass-loss rate to be $\sim 4 \times 10^{-2}$ M$_\odot$ yr$^{-1}$ for a CSM velocity of 100 km s$^{-1}$, and a CSM mass of $\gtrsim 1$ M$_\odot$. We find strong similarities for both the precursor, general light curve, and spectral evolution with SN 2009ip and similar SNe, although SN 2019zrk displays a brighter peak magnitude. Different scenarios for the nature of the 09ip-class of SNe, based on pulsational pair instability eruptions, wave heating, and mergers, are discussed. }

We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with previous studies, we find that the central region of the cloud has higher line widths at fixed size relative to the rest of the molecular cloud and to other LMC clouds, indicating an enhanced level of turbulent motions. However, there is no clear trend in gravitational boundedness (as measured by the virial parameter) with distance from R136. Structures observed in $^{13}$CO are spatially coincident with filaments and are close to a state of virial equilibrium. In contrast, CO structures vary greatly in virialization, with low CO surface brightness structures outside of the main filamentary network being predominantly unbound. The low surface brightness structures constitute ~10% of the measured CO luminosity; they may be shredded remnants of previously star-forming gas clumps, or alternatively the CO-emitting parts of more massive, CO-dark structures.

We present surface photometry and stellar kinematics of IC 3167, a dwarf galaxy hosting a lopsided weak bar and infalling into the Virgo cluster. We measured the bar radius and strength from broad-band imaging and bar pattern speed by applying the Tremaine-Weinberg method to stellar-absorption integral-field spectroscopy. We derived the ratio of the corotation radius to bar radius (R = 1.7 + 0.5 - 0.3) from stellar kinematics and bar pattern speed. The probability that the bar is rotating slowly is more than twice as likely as that the bar is fast. This allows us to infer that the formation of this bar was triggered by the ongoing interaction rather than to internal processes.

SMC X-1 has exhibited three super-orbital period excursions since the onset of X-ray monitoring beginning with RXTE's launch in 1995. NICER has recently probed a fourth observed excursion beginning in 2021 with our program Monitoring Observations of SMC X-1's Excursions (MOOSE). These sensitive new MOOSE data probe different super-orbital periods and phases within them. Spectral fits to the high-state continuum during April 2021 to January 2022 show that the intrinsic spectral shapes are characterised by a soft (kT~0.19 keV) disc component and a hard (Gamma~0.7) power-law tail. When the 2021-2022 NICER observations, taken during an excursion, are compared to 2016 XMM-Newton observations (outside of an excursion), we find little evidence for intrinsic spectral variability across the high-states, but find evidence for a >3 sigma change in the absorption, although we caution that there may be calibration differences between the two instruments. Thus, over different lengths of super-orbital periods, we see little evidence for intrinsic spectral changes in the high-state. Upcoming studies of the pulse profiles may shed light on the mechanism behind the excursions.

In a rigidly-rotating magnetohydrodynamic (MHD) system with convective turbulence, a large-scale dynamo, categorized as the $\alpha^2$-type, can be excited when the spin rate is large enough. In this paper, the rotational dependence of the $\alpha^2$-type dynamo and the cause of it are explored by mean-field (MF) dynamo models coupled with direct numerical simulations (DNSs) of MHD convections. Bearing the application to the solar/stellar dynamo in mind, we adopt a strongly-stratified polytrope as a model of the convective atmosphere. Our DNS models show that the $\alpha^2$-type dynamo is excited when ${\rm Ro} \lesssim 0.1$ where ${\rm Ro}$ is the Rossby number defined with the volume-averaged mean convective velocity. From the corresponding MF models, we demonstrate that the rotational dependence of the $\alpha^2$-type dynamo is mainly due to the change in the magnitude of the turbulent magnetic diffusion. With increasing the spin rate, the turbulent magnetic diffusion weakens while the $\alpha$-effect remains essentially unchanged over the convection zone, providing the critical point for the excitation of the large-scale dynamo. The ${\rm Ro}$-dependence of the stellar magnetic activity observable in the cool star is also discussed from the viewpoint of the rotational dependence of the turbulent electro-motive force. Overall our results suggest that, to get a better grasp of the stellar dynamo activity and its ${\rm Ro}$-dependence, it should be quantified how the convection velocity changes with the stellar spin rate with taking account of the rotational quenching and the Lorentz force feedback from the magnetic field on the convective turbulence.

We have investigated the variable stars in the field surrounding NGC 2355 based on the time-series photometric observation data. More than 3000 CCD frames were obtained in the V band spread over 13 nights with the Nanshan One-meter Wide-field Telescope. We have detected 88 variable stars, containing 72 new variable stars and 16 known variable stars. By analyzing these light curves, we classified the variable stars as follows: 26 eclipsing binaries, 52 pulsating stars, 4 rotating variables, and 6 unclear type variable stars for which their periods are much longer than the time baseline chosen. Employing Gaia DR2 parallax, kinematics, and photometry, the cluster membership of these variable stars were also analyzed for NGC 2355. In addition to the 11 variable members reported by Cantat-Gaudin et al. (2018), we identify 4 more variable member candidates located at the outer region of NGC 2355 and showed homogeneity in space positions and kinematic properties with the cluster members. The main physical parameters of NGC 2355 estimated from the two-color and color-magnitude diagrams are log(age/yr) = 8.9, E(B - V) = 0.24 mag, and [Fe/H] = - 0.07 dex.

We present sensitive $\rm NH_3$ (1,1)--(7,7) line images from the Karl G. Jansky Very Large Array toward successive shocks, which are associated with the blueshifted outflow lobe driven by the compact protobinary system L1157. Within a projection distance of 0.1 pc, our observations not only trace the quiescent and cold gas in the flattened envelope but also illustrate the complex physical and chemical processes that take place where the high-velocity jet impinges on its surrounding medium. Specifically, the $\rm NH_3$ ortho-to-para ratio is enhanced by a factor of 2--2.5 along the jet path, where the velocity offset between the line peak and the blueshifted wing reaches values as high as $\rm 10\,km\,s^{-1}$; it also shows a strong spatial correlation with the $\rm NH_3$ column density, which is enhanced to $\rm >10^{16}\,cm^{-2}$ toward the shock cavities. At a linear resolution of 1500\,au, our refined temperature map from the seven $\rm NH_3$ lines shows a gradient from the warm B0 eastern cavity wall ($\rm >120\,K$) to the cool cavity B1 and the earlier shock B2 ($\rm <80\,K$), indicating shock heating.

Chemical diversity of low-mass protostellar sources has so far been recognized, and environmental effects are invoked as its origin. In this context, observations of isolated protostellar sources without influences of nearby objects are of particular importance. Here, we report chemical and physical structures of the low-mass Class 0 protostellar source IRAS 16544$-$1604 in the Bok globule CB68, based on 1.3 mm ALMA observations at a spatial resolution of $\sim$70~au that were conducted as part of the large program FAUST. Three interstellar saturated complex organic molecules (iCOMs), CH$_3$OH, HCOOCH$_3$, and CH$_3$OCH$_3$, are detected toward the protostar. The rotation temperature and the emitting region size for CH$_3$OH are derived to be $131\pm11$~K and $\sim$10~au, respectively. The detection of iCOMs in close proximity to the protostar indicates that CB68 harbors a hot corino. The kinematic structure of the C$^{18}$O, CH$_3$OH, and OCS lines is explained by an infalling-rotating envelope model, and the protostellar mass and the radius of the centrifugal barrier are estimated to be $0.08-0.30$~$M_\odot$ and $< 30$ au, respectively. The small radius of the centrifugal barrier seems to be related to the small emitting region of iCOMs. In addition, we detect emission lines of c-C$_3$H$_2$ and CCH associated with the protostar, revealing a warm carbon chain chemistry (WCCC) on a 1000~au scale. We therefore find that the chemical structure of CB68 is described by a hybrid chemistry. The molecular abundances are discussed in comparison with those in other hot corino sources and reported chemical models.

We examine the redshift evolution of density environments around 2,163 radio galaxies with the stellar masses of $\sim10^{9}-10^{12} M_\odot$ between redshifts of $z=0.3-1.4$, based on the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and Faint Images of the Radio Sky at Twenty-cm (FIRST). We use the $k$-nearest neighbor method to measure the local galaxy number density around our radio galaxy sample. We find that the overdensities of the radio galaxies are weakly but significantly anti-correlated with redshift. This is consistent with the known result that the relative abundance of less-massive radio galaxies increases with redshift, because less-massive radio galaxies reside in relatively low density regions. Massive radio galaxies with stellar mass of $M_* >10^{11}M_\odot$ are found in high density environments compared with the control sample galaxies with radio-non-detection and matched-stellar-mass. Less-massive radio galaxies with $M_* <10^{11}M_\odot$ reside in average density environments. The fraction of the radio galaxies associated with the neighbors within a typical major merger scale, $<70$ kpc, is higher than (comparable to) that of the control galaxies at $M_* >10^{11}M_\odot$ ($M_* <10^{11}M_\odot$). We also find that the local densities around the radio galaxies are anti-correlated with the radio luminosities and black hole mass accretion rates at fixed stellar mass. These findings suggest that massive radio galaxies have matured through galaxy mergers in the past, and have supermassive black holes whose mass accretion almost ceased at $z>1.4$, while less-massive radio galaxies undergo active accretion just at this epoch, as they have avoided such merger events.

The detection of the first exoplanet paved the way into the era of transit photometry space missions with a revolutionary photometric precision that aim at discovering new exoplanetary systems around different types of stars. With this high precision, it is possible to derive very accurately the radii of exoplanets which is crucial for constraining their type and composition. However, it requires an accurate description of host stars, especially their center-to-limb variation of intensities (so called limb darkening) as it affects the planet-to-star radius ratio determination. We aim at improving the accuracy of limb darkening calculations for stars with a wide range of fundamental parameters. We used the recently developed 1D MPS-ATLAS code to compute model atmosphere structures and to synthesize stellar limb darkening on a very fine grid of stellar parameters. For the computations we utilized the most accurate information on chemical element abundances and mixing length parameters including convective overshoot. The stellar limb darkening was fitted using the two most accurate limb darkening laws: the power-2 and 4-parameters non-linear laws. We present a new extensive library of stellar model atmospheric structures, the synthesized stellar limb darkening curves, and the coefficients of parameterized limb-darkening laws on a very fine grid of stellar parameters in the Kepler, TESS, CHEOPS, and PLATO passbands. The fine grid allows overcoming the sizable errors introduced by the need to interpolate. Our computations of solar limb darkening are in a good agreement with available solar measurements at different view angles and wavelengths. Our computations of stellar limb darkening agree well with available measurements of Kepler stars. A new grid of stellar model structures, limb darkening and their fitted coefficients in different broad filters is provided in CDS.

We investigate the regularity of galaxy cluster gas density profiles and the link to the relation between core-excised luminosity, LXc, and mass from the Yx proxy, MYx, for 93 SZE-selected objects. The sample spans masses M500=[0.5 - 20] x 10e14 Msun, and lies at redshifts 0.05<z<1.13. Using XMM-Newton observations, we derive an average ICM density profile for the SZE-selected systems and determine its scaling with mass and redshift. This average profile evolves slightly stronger than self-similar (a_z = 2.09+/-0.02), and has significant dependence on mass (a_M = 0.22 +/- 0.01). Deviations from the average scaling with radius indicate different evolution for the core regions and the bulk. We measure the radial variation of the intrinsic scatter, finding a slight evolution with redshift. The average profile of the SZE-selected systems describes that of X-ray-selected systems at low redshift. The scaled core properties are positively skewed at later times, suggesting an increased incidence of centrally peaked objects at lower redshifts. The relation between LXc and MYx has an intrinsic scatter of 13%. Using simulations, we investigate the impact of selection effects, intrinsic scatter, and covariance on this relation. The slope is insensitive to selection and intrinsic scatter between quantities; however, the scatter is very dependent on the covariance between LXc and Yx. Accounting for our use of the Yx proxy to determine the mass, we estimate an upper limit to the intrinsic scatter with respect to the true mass of 22%. We probe the connection between the scatter in density profiles and that in the LXc-M relation. Our results suggest that the ICM bulk evolves approximately self-similarly, with the core regions evolving separately; indicate a variation of the gas content with mass; and show that LXc has a tight relation to the underlying mass.

Context. As part of the third Gaia data release, we present the contributions of the non-stellar and classification modules from the eighth coordination unit (CU8) of the Data Processing and Analysis Consortium, which is responsible for the determination of source astrophysical parameters using Gaia data. This is the third in a series of three papers describing the work done within CU8 for this release. Aims. For each of the five relevant modules from CU8, we summarise their objectives, the methods they employ, their performance, and the results they produce for Gaia DR3. We further advise how to use these data products and highlight some limitations. Methods. The Discrete Source Classifier (DSC) module provides classification probabilities associated with five types of sources: quasars, galaxies, stars, white dwarfs, and physical binary stars. A subset of these sources are processed by the Outlier Analysis (OA) module, which performs an unsupervised clustering analysis, and then associates labels with the clusters to complement the DSC classification. The Quasi Stellar Object Classifier (QSOC) and the Unresolved Galaxy Classifier (UGC) determine the redshifts of the sources classified as quasar and galaxy by the DSC module. Finally, the Total Galactic Extinction (TGE) module uses the extinctions of individual stars determined by another CU8 module to determine the asymptotic extinction along all lines of sight for Galactic latitudes |b| > 5 deg. Results. Gaia DR3 includes 1591 million sources with DSC classifications; 56 million sources to which the OA clustering is applied; 1.4 million sources with redshift estimates from UGC; 6.4 million sources with QSOC redshift; and 3.1 million level 9 HEALPixes of size 0.013 squared degree, where the extinction is evaluated by TGE.

We report on a 5$\sigma$ detection of polarized 3-6 keV X-ray emission from the supernova remnant Cassiopeia A with the Imaging X-ray Polarimetry Explorer (IXPE) mission. The overall polarization degree of 1.8 $\pm$ 0.3% is detected by summing over a large region assuming circular symmetry for the polarization vectors. Correcting for thermal X-ray emission implies a 2.4% polarization degree for the synchrotron component. For the X-ray synchrotron dominated forward shock region the polarization degree of the synchrotron component is close to 5%. A pixel-by-pixel search for polarization provides a few tentative detections from discrete regions at the 3$\sigma$ confidence level. Given the number of pixels, the significance is insufficient to claim a detection for individual pixels, but implies considerable turbulence on small scales. Cas A's X-ray continuum emission is dominated by synchrotron radiation from regions within $10^{17}$ cm of the forward and reverse shocks. We find that i) the measured polarization angle corresponds to a radially-oriented magnetic field similar to what has been inferred from radio observations; ii) the X-ray polarization degree is lower than in the radio band (around 5%). Since shock compression should impose a tangential magnetic field structure, the IXPE results imply that magnetic-fields are reoriented within approximately $10^{17}$ cm of the shock. If the magnetic-field alignment is due to locally enhanced acceleration near quasi-parallel shocks, the preferred X-ray polarization angle suggests a size of $3 \times 10^{16}$ cm for cells with radial magnetic fields.

Magnetic fields play an important role in the evolution of molecular clouds and star formation. Using the Velocity Gradient Technique (VGT) model, we measured the magnetic field in Orion A using the 12CO, 13CO, and C18O (1-0) emission lines at a scale of 0.07 pc. The measured B-field shows an east-west orientation that is perpendicular to the integral shaped filament of Orion A at large scale. The VGT magnetic fields obtained from 13CO and C18O are in agreement with the B-field that is measured from the Planck 353 GHz dust polarization at a scale of 0.55 pc. Removal of density effects by using a Velocity Decomposition Algorithm can significantly improve the accuracy of the VGT in tracing magnetic fields with the 12CO (1-0) line. The magnetic field strength of seven sub-clouds, OMC-1, OMC-2, OMC-3, OMC-4, OMC-5, L 1641-N, and NGC 1999 has also been estimated with the Davis-Chandrasekhar-Fermi (DCF) and MM2 technique, and these are found to be in agreement with previous results obtained from dust polarization at far-infrared and sub-millimeter wavelengths. At smaller scales, the VGT proves a good method to measure magnetic fields.

A recent work (arXiv:2104.14481) has found a statistically significant transition in the Baryonic Tully-Fisher relation (BTFR) using low redshift data ($z<0.1$), with the transitions occurring at about 9 and 17 Mpc. Motivated by this finding, we carry out a variant of this analysis by fitting the data to an augmented BTFR, where both the exponent as well as normalization constant vary as a function of distance. We find that both the exponent and normalization constant show only a marginal variation with distance, and are consistent with a constant value, to within $2\sigma$. We also checked to see if there is a statistically significant difference between the BTFR results after bifurcating the dataset at distances of 9 and 17 Mpc. We find that almost all the sets of subsamples obey the BTFR with $\chi^2$/dof close to 1 and the best-fit parameters consistent across the subsamples. Only the subsample with $D<17$ Mpc shows a marginal discrepancy (at $1.75\sigma$) with respect to the BTFR. Therefore, we do not find any evidence for statistically significant differences in the BTFR at distances of 9 and 17 Mpc.

Galaxy morphological classification is a fundamental aspect of galaxy formation and evolution studies. Various machine learning tools have been developed for automated pipeline analysis of large-scale surveys, enabling a fast search for objects of interest. However, crowded regions in the image may pose a challenge as they can lead to bias in the learning algorithm. In this Research Note, we present galmask, an open-source package for unsupervised galaxy masking to isolate the central object of interest in the image. galmask is written in Python and can be installed from PyPI via the pip command.

We present qrpca, a fast and scalable QR-based principal component analysis package. The software, written in both R and python languages, makes use of torch for internal matrix computations, and enables GPU acceleration, when available. qrpca provides similar functionalities to prcomp (R) and sklearn (python) packages respectively. A benchmark test shows that qrpca can achieve computational speeds 10-20 $\times$ faster for large dimensional matrices than default implementations, and is at least twice as fast for a standard decomposition of spectral data cubes. The qrpca source code is made freely available to the community.

Correlated noise could impact the search for the gravitational wave background at future Earth-based gravitational-wave detectors. Due to the small distance ($\sim$ 400 m) between the different interferometers of the Einstein Telescope, correlated seismic noise could have a significant effect. To this extent, we study the seismic correlations at the Earth's surface, as well as underground, between seismometers and geophones separated by several hundreds of meters, in the frequency range 0.05 Hz - 50 Hz. Based on these correlated seismic fields we predict the levels of correlated Newtonian noise (NN). We construct upper limits on the allowed seismic coupling function such that correlated seismic noise does not affect the search for an isotropic gravitational wave background. Assuming a facility located 300 m below the surface, the impact on the search for a gravitational wave background of correlated NN from Rayleigh waves are found to be problematic up to $\sim$ 5 Hz. The NN from body waves, however, constitutes a serious threat to the search of a gravitational wave background. Correlated NN from body waves could be up to five to seven orders of magnitude above the planned sensitivity at $\sim$ 3 Hz and it could impede any search for a gravitational wave background below 40 Hz. With a factor 10 of NN reduction via NN cancellation in each interferometer, the effects of the NN on the stochastic search could be eliminated above 30 Hz.

The X-ray spectra of active galactic nuclei (AGNs) often exhibit an excess of emission above the primary power-law at energies <~ $2$ keV. Two models for the origin of this `soft excess' are ionized relativistic reflection from the inner accretion disc and Comptonization of thermal emission in a warm corona. Here, we introduce reXcor, a new AGN X-ray ($0.3$-$100$ keV) spectral fitting model that self-consistently combines the effects of both ionized relativistic reflection and the emission from a warm corona. In this model, the accretion energy liberated in the inner disc is distributed between a warm corona, a lamppost X-ray source, and the accretion disc. The emission and ionized reflection spectrum from the inner $400$ $r_g$ of the disc is computed, incorporating the effects of relativistic light-bending and blurring. The resulting spectra predict a variety of soft excess shapes and sizes that depend on the fraction of energy dissipated in the warm corona and lamppost. We illustrate the use of reXcor by fitting to the joint XMM-Newton and NuSTAR observations of the Seyfert 1 galaxies HE 1143-1820 and NGC 4593, and find that both objects require a warm corona contribution to the soft excess. Eight reXcor table models, covering different values of accretion rate, lamppost height and black hole spin, are publicly available through the XSPEC website. Systematic use of reXcor will provide insight into the distribution of energy in AGN accretion flows.

This is the second paper in a series using data from tens of thousands S0 galaxies of the local Universe ($z\lesssim 0.1$) retrieved from the NASA-Sloan Atlas. It builds on the outcomes of the previous work, which introduced a new classification scheme for these objects based on the principal component analysis (PCA) of their optical spectrum and its projections on to the first two eigenvectors or principal components (the PC1$\unicode{x2013}$PC2 diagram). We provide a comprehensive characterization of the activity of present-day S0s throughout both the broad-band PC1$\unicode{x2013}$PC2 spectral classifier and the conventional narrow-line BPT/WHAN ones, contrasting the different types of activity classes they define, and present an alternative diagram that exploits the concordance between WHAN and PCA demarcations. The analysis is extended to the mid-infrared, radio and X-ray wavelengths by crossmatching our core sample with data from the WISE, FIRST, XMM$\unicode{x2013}$Newton, and Chandra surveys. This has allowed us to carry out a thorough comparison of the most important activity diagnostics in the literature over different wavebands, discuss their similarities and differences, and explore the connections between them and with parameters related to star formation and black hole accretion. In particular, we find evidence that the bulk of nebular emission from radio and X-ray detected S0$\unicode{x2013}$Seyfert and LINER systems is not driven by star birth, while the dominant ionising radiation for a number of LINERs might come from post-AGB stars. These and other outcomes from the present work should be transferable to other morphologies.

Post-common envelope binary systems evolve when matter is transferred from the primary star at a rate that cannot be accommodated by its secondary companion. A common envelope forms which is subsequently ejected resulting in a system with a binary period frequently between 2 and 3 hours. Where circumbinary companions are predicted, it remains unclear whether they form before or after the common envelope ejection. From observations of eclipse time variations (ETVs), exoplanet databases e.g. NASA Exoplanet Archive, list typically a dozen systems with confirmed circumbinary planets. Here we examine seven of these systems, discuss other possible causes and consider whether, for these dynamic systems, the ETV methodology is a reliable indicator of planetary companions. The systems selected were those where we could determine precise eclipse timings, free from significant extraneous effects such as pulsations, and present 163 new times of minima permitting us to test existing models. Over thirty circumbinary models have been proposed for these seven systems and note all, other than the latest model for NY Vir which remains to be fully tested, fail within a year to accurately predict eclipse times. In examining alternative mechanisms we find that magnetic effects could contribute significantly in two of the seven systems studied. We conclude that the structure of these dynamic systems, with the extreme temperature differences and small binary separations, are not fully understood and that many factors may contribute to the observed ETVs.

Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between color and luminosity ($\beta$) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here we take advantage of state-of-the-art modeling of galaxy properties to characterize dust parameters (dust attenuation $A_V$, and a parameter describing the dust law slope $R_V$) for the Dark Energy Survey (DES) SN Ia host galaxies using the publicly available \texttt{BAGPIPES} code. Utilizing optical and infrared data of the hosts alone, we find three key aspects of host dust that impact SN Ia cosmology: 1) there exists a large range ($\sim1-6$) of host $R_V$ 2) high stellar mass hosts have $R_V$ on average $\sim0.7$ lower than that of low-mass hosts 3) there is a significant ($>3\sigma$) correlation between the Hubble diagram residuals of red SNe Ia that when corrected for reduces scatter by $\sim13\%$ and the significance of the ``mass step'' to $\sim1\sigma$. These represent independent confirmations of recent predictions based on dust that attempted to explain the puzzling ``mass step'' and intrinsic scatter ($\sigma_{\rm int}$) in SN Ia analyses. We also find that red-sequence galaxies have both lower and more peaked dust law slope distributions on average in comparison to non red-sequence galaxies. We find that the SN Ia $\beta$ and $\sigma_{\rm int}$ both differ by $>3\sigma$ when determined separately for red-sequence galaxy and all other galaxy hosts. The agreement between fitted host-$R_V$ and SN Ia $\beta$ \& $\sigma_{\rm int}$ suggests that host dust properties play a major role in SN Ia color-luminosity standardization and supports the claim that SN Ia intrinsic scatter is driven by $R_V$ variation.

We present spatially-resolved morphological properties of [CII] 158 $\mu$m, [OIII] 88 $\mu$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $\mu$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $\beta = 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5\sigma$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies.

The third data release (DR3) of the European Space Agency satellite Gaia provides coordinates, parallaxes, proper motions, and radial velocities for a sample of $\sim 34$ million stars. We use the combined 6-dimensional phase space information to search for hypervelocity stars (HVSs), unbound stars accelerated by dynamical processes happening in the Galactic Centre. By looking at the kinematics of Gaia DR3 stars in Galactocentric coordinates and by integrating their orbits in the Galactic potential, we do not identify any HVS candidates with a velocity higher than $700$ km s$^{-1}$ and robustly observed kinematics. Assuming a scenario wherein the interaction between a stellar binary and the massive black hole Sgr A$^*$ is responsible for HVS ejections from the Galactic Centre, we derive degenerate limits on the ejection rate of HVSs and the slope of the initial mass function of the primary star among binaries in the Galactic Centre. Our results indicate that the HVS ejection rate is $\lesssim 8\times10^{-5}$ yr$^{-1}$ assuming a Salpeter mass function, and this upper limit becomes progressively smaller for an increasingly top-heavy mass distribution. A fiducial HVS ejection rate of $10^{-4}$ yr$^{-1}$ prefers a mass function slope $\lesssim -2.35$, disfavouring previously claimed top-heavy initial mass functions among stars in the Galactic Centre.

The Fermi paradox has given rise to various attempts to explain why no evidence of extraterrestrial civilisations was found so far on Earth and in our Solar System. Here, we present a dynamical model for the development of such civilisations, which accounts for self-destruction, colonisation and astrophysical destruction mechanisms of civilisations including gamma-ray bursts, type Ia and type II supernovae as well as radiation from the supermassive black hole. We adopt conservative estimates regarding the efficiency of such processes and find that astrophysical effects can influence the development of intelligent civilisations and change the number of systems with such civilisations by roughly a factor of 2; potentially more if the feedback is enhanced. Our results show that non-equilibrium evolution allows for solutions in-between extreme cases such as "rare Earth" or extreme colonisation, including scenarios with civilisation fractions between 10^{-2} and 10^{-7}. These would imply still potentially large distances to the next such civilisations, particularly when persistence phenomena are being considered. As previous studies, we confirm that the main uncertainties are due to the lifetime of civilisations as well as the assumed rate of colonisation. For SETI-like studies, we believe that unbiased searches are needed considering both the possibilities that the next civilisations are nearby or potentially very far away.

We discuss baryogenesis in scenarios where the Universe is reheated to temperatures $\lesssim 100\,$GeV by the decay of long-lived massive particles into energetic SM particles. Before its thermalization, the center-of-mass energy in collisions between such a particle and a particle from the ambient plasma can be higher than the typical sphaleron mass, even if the temperature of the plasma itself is much lower. Optimistic estimates for the high energy enhancement of the sphaleron cross section suggest that successful baryogenesis is possible for reheating temperatures as low as $0.1\text{-}1\,$GeV. With a simple extension of the SM, sufficient baryon production can be achieved even if more pessimistic results for the sphaleron rate are correct. In both cases this scenario can be probed in collider and cosmic-ray experiments. We briefly discuss the possible origin of the required CP violation.

Cosmic strings that couple to neutrinos may account for a portion of the high-energy astrophysical neutrino (HEAN) flux seen by IceCube. Here, we calculate the observed spectrum of neutrinos emitted from a population of cosmic string loops that contain quasi-cusps, -kinks, or kink-kink collisions. We consider two broad neutrino emission models: one where these string features emit a neutrino directly, and one where they emit a scalar particle which then eventually decays to a neutrino. In either case, the spectrum of cosmic string neutrinos does not match that of the observed HEAN spectrum. We thus find that the maximum contribution of cosmic string neutrinos, through these two scenarios, to be at most $\sim 45$ % of the observed flux. However, we also find that the presence of cosmic string neutrinos can lead to bumps in the observed neutrino spectrum. Finally, for each of the models presented, we present the viable parameter space for neutrino emission.

In this Letter we demonstrate that a future accelerator-based neutrino experiment such as DUNE can greatly increase its sensitivity to a variety of new physics scenarios by operating in a mode where the proton beam impinges on a beam dump. We consider two new physics scenarios, namely light dark matter (LDM) and axion-like particles (ALPs) and show that by utilizing a dump mode at DUNE, unexplored new regions of parameter space can be probed with an exposure of only 3 months with half of its expected initial beam power. Specifically, target-less configuration of future high intensity neutrino experiments will probe the parameter space for thermal relic DM as well as the QCD axion (DFSZ and KSVZ). The strength of such configuration in the context of new physics searches stems from the fact that the neutrino flux is significantly reduced compared to that of the target, resulting in much smaller backgrounds from neutrino interactions. We have verified this in detail by explicitly computing neutrino fluxes which we make publicly available in order to facilitate further studies with a target-less configuration.

We present a new, testable picture of the early universe in finite nonlocal quantum gravity, which is conformally invariant at the classical and quantum levels. The high-energy regime of the theory consists of two phases, a conformally invariant trans-Planckian phase and a sub-Planckian or Higgs phase described by an action quadratic in the Ricci tensor and where the cosmos evolves according to the standard radiation-dominated model. In the first phase, all the issues of the hot big bang such as the singularity, flatness, and horizon problems find a universal and simple non-inflationary solution by means of Weyl conformal invariance, regardless of the microscopic details of the theory. In the second phase, once conformal symmetry is spontaneously broken, logarithmic quantum corrections to the action make both the primordial tensor spectrum (from graviton fluctuations) and the scalar spectrum (from thermal fluctuations) quasi scale invariant. Although nonlocal quantum gravity is an explicit realization of this scenario, any finite conformal theory with quadratic limit would yield the same results. In particular, the value of the scalar spectral index only depends on a beta function and a scale ratio, and is consistent with observations. The theory also predicts a positive tensor spectral index $n_{\rm t} =1-n_{\rm s}$, a fixed tensor-to-scalar ratio $r_{0.05}\gtrsim 0.01$, and a blue-tilted stochastic gravitational-wave background, all universal predictions testable in the immediate or near future.

It is widely accepted that curvature singularity resolution should be a feature of quantum gravity. We present a class of time-dependent asymptotically flat spherically symmetric metrics that model gravitational collapse in quantum gravity. The metrics capture intuitions associated with the dynamics of singularity resolution, and horizon formation and evaporation following a matter bounce. A parameter in the metric associated with the speed of the bounce determines black hole lifetime as a power of its mass; this includes the Hawking evaporation time $M^{3}$.

A key goal in cosmology in the upcoming decade will be to form a better understanding of the accelerated expansion of the Universe. Upcoming surveys, such as the Vera C. Rubin Observatory's 10-year Legacy Survey of Space and Time (LSST), Euclid and the Square Killometer Array (SKA) will deliver key datasets required to tackle this and other puzzles in contemporary cosmology. With this data, constraints of unprecedented power will be put on different models of dark energy and modified gravity. In this context it is crucial to understand how screening mechanisms, which hide the deviations of these theories from the predictions of general relativity in local experiments, affect structure formation. In this work we approach this problem by using a combination of analytic and numerical methods to describe chameleon screening in the context of cosmic voids. We apply a finite element method code, SELCIE, to solve the chameleon equation of motion for a number of void profiles derived from observational data and simulations. The obtained results indicate a complex relationship between the properties of cosmic voids and the size of the chameleon acceleration of a test particle. We find that the fifth force on a test particle in a void is primarily related to the depth and the inner density gradient of the void. For realistic void profiles, the obtained chameleon-to-Newtonian acceleration ratios range between $a_{\phi}/a_{\rm Newt} \approx 10^{-6} - 10^{-5}$. However, it should be noted that in unusually deep voids with large inner density gradients, the acceleration ratios can be significantly higher. We also discuss the optimal density profiles for detecting the fifth force in the upcoming observational surveys.

The detection of kilohertz-band gravitational waves promises discoveries in astrophysics, exotic matter, and cosmology. To improve the kilohertz quantum noise-limited sensitivity of interferometric gravitational-wave detectors, we investigate nondegenerate internal squeezing: optical parametric oscillation inside the signal-recycling cavity with distinct signal and idler frequencies. We use an analytic Hamiltonian model to show that this stable, all-optical technique is tolerant to decoherence from optical detection loss and that it, with its optimal readout scheme, is feasible for broadband sensitivity enhancement.

In the present article, we show that a simple modification to the Einstein-Hilbert action can explain the possibility of mutual interaction between the cosmic fluids. That is achieved considering the Weyl Integrable Spacetime in the background of a nonflat Friedmann-Lemna\^{i}tre-Robertson-Walker geometry for the universe. We show that widely-known phenomenological interacting cosmological scenarios can naturally appear in this context. We then performed the dynamical system analysis of the underlying cosmological scenario and explored many possibilities extracted from this gravitational theory.

Cosmological Collider Physics provides information on high-energy elementary particles by observing the spacetime fluctuations generated by inflation through the cosmic microwave background radiation. In other words, it is a method to investigate physics on energy scales that cannot be reached by terrestrial accelerators by means of precise observations of the universe. In this paper, we focus on the case where the GUT scale is close to the energy scale of inflation, and calculate the non-Gaussianity due to scalar fluctuations produced by the Higgs boson in GUT. The results are found to be consistent with the current observed restrictions on non-Gaussianity without a drastic fine tuning of parameters, which suggests the existence of Higgs boson in GUT.

WIMPs can be captured in compact stars such as white dwarves (WDs) leading to an increase in the star luminosity through their annihilation process. We show that when the WIMP interacts with the nuclear targets within the WD through inelastic scattering and its mass exceeds a few tens GeV the data on low-temperature large-mass WDs in the Messier 4 globular cluster can probe values of the mass splitting $\delta\lesssim$ 40 MeV. Such value largely exceeds those ensuing from direct detection and from solar neutrino searches. We apply such improved constraint to the specific DM scenario of a self-conjugate bi-doublet in the Left-Right Symmetric Model (LRSM), where the standard $SU(2)_L$ group with coupling $g_L$ is extended by an additional $SU(2)_R$ with coupling $g_R$. We show that bounds from WDs significantly reduce the cosmologically viable parameter space of such scenario, in particular requiring $g_R>g_L$. For instance, for $g_R/g_L$ = 1.8 we find the two viable mass ranges 1.2 TeV $\lesssim m_\chi\lesssim$ 3 TeV and 5 TeV $\lesssim m_\chi\lesssim$ 10 TeV, when the charged $SU(2)_R$ gauge boson mass $M_{W_2}$ is lighter than $\simeq$ 12 TeV. We also discuss the ultraviolet completion of the LRSM model, when the latter is embedded in a Grand Unified Theory. We show that such low-energy parameter space and compatibility to proton-decay bounds require a non-trivial extension of the particle content of the minimal model. We provide a specific example where $M_{W_2}\lesssim$ 10 TeV is achieved by extending the LRSM at high energy with color triplets that are singlets under all other groups, and $g_R/g_L>$1 is obtained by introducing $SU(2)_L$ triplets with no $SU(2)_R$ counterparts, i.e. by breaking the symmetry between the multiplets of $SU(2)_L$ and $SU(2)_R$.

We discuss the discovery potential of JUNO experiment for neutrino lines from MeV dark matter (DM) annihilation and decay in a model independent way. We find that JUNO will be able to give severe constraints on the cross section of DM annihilating into neutrinos and on the lifetime of DM decaying into neutrinos. More concretely, with $20$ years of data-taking in the fiducial volume $17$ kton, the cross section will be constrained as strong as $4\times 10^{-26}\,{\rm cm^{3}\,sec^{-1}}$ for the mass of a DM particle $m_{\chi} \simeq 15\,{\rm MeV}$ at $90\,\%$ C.L., which might be strong enough to test thermal production mechanism of DM particles for such range of DM mass. The lifetime will be constrained as strong as $1\times 10^{24}\,{\rm sec}$ for the mass of a DM particle $m_{\chi} \simeq 100\,{\rm MeV}$ at $90\,\%$ C.L..

Over the last two decades, motivations for modified gravity have emerged from both theoretical and observational levels. $f(R)$ and Chern-Simons gravity have received more attention as they are the simplest generalization. However, $f(R)$ and Chern-Simons gravity contain only the scalar degree of freedom and, as a result, do not include other modes of modified theories of gravity. In contrast, quadratic gravity (also referred to as Stelle gravity) is the most general second-order modification to 4-D general relativity and contains massive tensor modes that are not present in $f(R)$ and Chern-Simons gravity. Using two different physical settings $-$ the gravitational wave energy-flux measured by the detectors and the backreaction of the emitted gravitational radiation on the spacetime of the remnant black hole $-$ we demonstrate that massive tensor modes carry more energy than scalar modes. Our analysis shows that the effects are pronounced for intermediate-mass black holes, which are prime targets for LISA.

The discrepancies between data on rare $b$-hadron decays, controlled by the underlying neutral-current transitions $b\to s\ell^+\ell^- (\ell = e, \mu)$, and the corresponding Standard Model predictions constitute one of the most intriguing hints for new physics. Leptoquarks are prime candidates to solve these anomalies and, in particular, the scalar leptoquark, $S_3$, triplet under $SU(2)_L$ with hypercharge $Y=-1/3$, provides a very good fit to data. Here, for the first time, we entertain the possibility that the same scalar leptoquark, responsible for the LFU anomalies, is the portal to a dark sector consisting of two additional vector-like fermions, one of which is a candidate for the cosmological dark matter. We study two scenarios, where the dark matter candidate belongs to an $SU(2)_L$ singlet and triplet respectively, and discuss the theory parameter space in the context of the dark matter candidate's relic density and prospects for direct and indirect dark matter searches. Direct detection rates are highly suppressed, and generically below the neutrino floor. Current observations with, and future prospects for, high-energy gamma-ray telescopes such as HESS and the Cherenkov Telescope Array are much more promising, as they already provide powerful constraints on the models under consideration, and will potentially probe the full parameter space in the future.

Non-standard neutrino interactions with a massive boson can produce the bosons in the core of core-collapse supernovae (SNe). After the emission of the bosons from the SN core, their subsequent decays into neutrinos can modify the SN neutrino flux. We show future observations of neutrinos from a next galactic SN in Super-Kamiokande (SK) and Hyper-Kamiokande (HK) can probe flavor-universal non-standard neutrino couplings to a light boson, improving the previous limit from the SN 1987A neutrino burst by several orders of magnitude. We also discuss sensitivity of the flavor-universal non-standard neutrino interactions in future observations of diffuse neutrinos from all the past SNe, known as the diffuse supernova neutrino background (DSNB). According to our analysis, observations of the DSNB in HK, JUNO and DUNE experiments can probe such couplings by a factor of $\sim 2$ beyond the SN 1987A constraint. However, our result is also subject to a large uncertainty concerning the precise estimation of the DSNB.

We propose a novel idea using time-of-flight measurements to probe the interaction of dark matter (DM) with Standard Model particles including neutrinos ($\nu$). When a supernova (SN) explosion occurs in a galaxy, DM with mass $m_\chi \lesssim \mathcal{O}$(MeV) in the halo can be boosted by SN$\nu$ to relativistic speed. The majority of the SN$\nu$ boosted DM (BDM) arrives the Earth with a delayed time determined solely by $m_\chi$ and independent of the size of the cross section. These BDM can interact with targets in neutrino and DM experiments and produce afterglow events post the detection of SN$\nu$ therein. The BDM events contain unique time-dependent features allowing for the measurement of $m_\chi$. New cross section constraints on $\sqrt{\sigma_{\chi e} \sigma_{\chi\nu}}$ are derived from SN1987a in Large Magellanic Cloud with data from the Kamiokande and Super-Kamiokande experiments. Potential sensitivities for a next Galactic SN with Hyper-Kamiokande are projected. Our results improve the existing bounds on $\sigma_{\chi e}$ by 1-3 orders of magnitude for $m_\chi\lesssim \mathcal{O}(100\,{\rm keV})$ for cases where $\sigma_{\chi\nu}$ and $\sigma_{\chi e}$ are of similar strength.

We study the magnetic moment of leptons in extremely hot universe and superdense media of stars at high temperatures. Anomalous magnetic moment of charged leptons is inversely proportional to its mass, whereas, the induced dipole moment of neutral leptons is directly proportional to their mass. Neutral massive point particles exert nonzero magnetic moment as a higher order effect, which is smaller than the anomalous magnetic moment of a charged particle of their same flavor partner. All leptons acquire some extra mass due to their interaction with the medium and affect the magnetic moment accordingly. We compare the contribution to the magnetic moment of various leptons due to their temperature and chemical potential dependent masses. It is shown that the magnetic moment contributions are non-ignorable for lighter leptons and heavy neutrinos. These calculations are very important to study particle propagation in the early universe and in superdense stellar media.

The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. In \cite{BretJPP2018}, we devised a model for the kinetic history of the plasma through the shock front, allowing to self-consistently compute the downstream anisotropy, hence the density jump, in terms of the upstream parameters. This model dealt with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only 1 switch-on solution within a limited range of Alfv\'{e}n Mach numbers, our model offers 2 kinds of solutions within a slightly different range of Alfv\'{e}n Mach numbers. These 2 solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold.