Papers for Wednesday, Nov 10 2021

In this work we introduce group-equivariant self-attention models to address the problem of explainable radio galaxy classification in astronomy. We evaluate various orders of both cyclic and dihedral equivariance, and show that including equivariance as a prior both reduces the number of epochs required to fit the data and results in improved performance. We highlight the benefits of equivariance when using self-attention as an explainable model and illustrate how equivariant models statistically attend the same features in their classifications as human astronomers.

We study the evolution of baryonic gas in cosmologically growing dark matter halos. To accurately model both the inner and outer regions of the halos, we use a dark matter density profile that transitions smoothly from the NFW profile within the virial radius to a more realistic flat profile far beyond the halo. We construct a dark matter gravitational potential consistent with this density profile, and we use a "cosmological" potential that accounts for gas evolution consistent with Hubble expansion at large radii. Gas is initialized with a density $\approx$ 0.2 times the dark matter density, consistent with the universal baryon fraction $\rho_{\rm g}/(\rho_{\rm g}+\rho_{\rm DM}) \approx 0.17$. We study the formation of the virial shock and evolution of the baryon fraction, including the effects of radiative cooling and AGN jet feedback. The feedback is powered by the accretion of cold gas onto a central supermassive black hole (SMBH). The cores of the halo exhibit heating and cooling cycles, whose strength and duration depend on the feedback efficiency and the halo mass. The central SMBH initially grows exponentially with time in the early quasar phase, but the growth slows down at later times. The baryon fraction in the core decreases with increasing feedback efficiency and decreasing halo mass. While the halo outskirts evolve self-similarly, the core density is non-evolving, in agreement with cluster observations. We analyze the correlations between the properties of the gas and the central SMBH, and explore the existence of a fundamental plane.

One signature of an expanding universe is the time-variation of the cosmological abundances of its different components. For example, a radiation-dominated universe inevitably gives way to a matter-dominated universe, and critical moments such as matter-radiation equality are fleeting. In this paper, we point out that this lore is not always correct, and that it is possible to obtain a form of "stasis" in which the relative cosmological abundances $\Omega_i$ of the different components remain unchanged over extended cosmological epochs, even as the universe expands. Moreover, we demonstrate that such situations are not fine-tuned, but are actually global attractors within certain cosmological frameworks, with the universe naturally evolving towards such long-lasting periods of stasis for a wide variety of initial conditions. The existence of this kind of stasis therefore gives rise to a host of new theoretical possibilities across the entire cosmological timeline, ranging from potential implications for primordial density perturbations, dark-matter production, and structure formation all the way to early reheating, early matter-dominated eras, and even the age of the universe.

Hydrodynamical simulations show that the ram-pressure stripping in galaxy clusters fosters a strong interaction between stripped interstellar medium (ISM) and the surrounding medium, with the possibility of intracluster medium (ICM) cooling into cold gas clouds. Exploiting the MUSE observation of three jellyfish galaxies from the GAs Stripping Phenomena in galaxies with MUSE (GASP) survey, we explore the gas metallicity of star-forming clumps in their gas tails. We find that the oxygen abundance of the stripped gas decreases as a function of the distance from the parent galaxy disk; the observed metallicity profiles indicate that more than 40% of the most metal-poor stripped clouds are constituted by cooled ICM, in qualitative agreement with simulations that predict mixing between the metal-rich ISM and the metal-poor ICM.

The anisotropies of the Stochastic Gravitational Wave Background (SGWB) produced by merging compact binaries constitute a possible new probe of the Large-Scale Structure (LSS) of the Universe. However, a large shot noise contribution, caused by the discreteness of the GW emitters, is expected to cover the anisotropies induced by the LSS. In this work, we investigate the potential of cross-correlating forthcoming high precision measurements of the SGWB energy density and the Cosmic Microwave Background (CMB) lensing convergence in order to mitigate the effect of shot noise. Combining a detailed modeling of stellar and galactic astrophysics with a novel framework to distribute the GW emitters in the sky, we compute the auto- and cross-correlation power spectra for the two cosmic fields, we evaluate the shot noise contribution and we make a prediction for the signal-to-noise ratio. The results of our analysis show that the SGWB energy density is significantly correlated with the CMB lensing convergence and that the cross-correlation between these two cosmic fields is a powerful tool to reduce the impact of shot noise, paving the way to the detection of the intrinsic SGWB anisotropies.

We revisit the possibility that the stochastic common-spectrum process recently detected by the NANOGrav pulsar timing array experiment could be due to primordial gravitational waves (GWs). A na\"{i}ve extrapolation down to interferometer scales of the blue GW spectrum required to explain NANOGrav consistently with Cosmic Microwave Background (CMB) observations would strongly violate upper limits on the stochastic GW background (SGWB) amplitude from LIGO/Virgo. In combination with the fact that there are over 19 decades in frequency between CMB and interferometer scales, this motivates us to move beyond the commonly adopted approximation of a pure power-law GW spectrum. We consider a broken power-law parametrization for the SGWB spectrum, which turns from blue to red above the break frequency: while phenomenological, this choice maps to various well-motivated early-Universe models, including scenarios featuring non-instantaneous reheating or a non-standard background expansion following reheating. After a detailed discussion of the contribution of the resulting SGWB to the early-Universe radiation energy density, we constrain the broken power-law model against a wide variety of multi-frequency cosmological and GW observations. We find that this phenomenological model is able to explain the NANOGrav signal while remaining in agreement with upper limits on the tensor-to-scalar ratio on CMB scales, Big Bang Nucleosynthesis constraints on the early-Universe radiation energy density, and upper limits on the SGWB amplitude on interferometer scales. We briefly discuss the very bright prospects for testing this model with next-generation probes across the GW frequency landscape, which motivate further exploring connections to specific well-motivated early-Universe models.

We investigate shock acceleration in a realistic astrophysical environment with density inhomogeneities. The turbulence induced by the interaction of the shock precursor with upstream density fluctuations amplifies both upstream and downstream magnetic fields via the turbulent dynamo. The dynamo-amplified turbulent magnetic fields (a) introduce variations of shock obliquities along the shock face, (b) enable energy gain through a combination of shock drift and diffusive processes, (c) give rise to various spectral indices of accelerated particles, (d) regulate the diffusion of particles both parallel and perpendicular to the magnetic field, and (e) increase the shock acceleration efficiency. Our results demonstrate that upstream density inhomogeneities and dynamo amplification of magnetic fields play an important role in shock acceleration, and thus shock acceleration depends on the condition of the ambient interstellar environment. The implications on understanding radio spectra of supernova remnants are also discussed.

Analyses of absorption from disk winds and atmospheres in accreting compact objects typically treat the central emitting regions in these systems as point sources relative to the absorber. This assumption breaks down if the absorbing gas is located within $few \times 1000\cdot GM/{c}^{2}$, in which case a small component of the absorber's Keplerian motion contributes to the velocity-width of absorption lines. Here, we demonstrate how this velocity-broadening effect can be used to constrain the sizes of central engines in accreting compact objects via a simple geometric relationship, and develop a method for modeling this effect. We apply this method on the Chandra/HETG spectra of three ultra-compact and short period neutron star X-ray binaries in which evidence of gravitationally redshifted absorption, owing to an inner-disk atmosphere, has recently been reported. The significance of the redshift is above $5\sigma$ for XTE J1710$-$281 (this work) and 4U 1916$-$053, and is inconsistent with various estimates of the relative radial velocity of each binary. For our most sensitive spectrum (XTE J1710$-$281), we obtain a 1$\sigma$ upper bound of 310 $\text{km}$ $\text{s}^{-1}$ on the magnitude of this geometric effect and a central engine of size ${R}_{CE} < 60 ~ GM/{c}^{2}$ (or, $< 90 ~ GM/{c}^{2}$ at the $3\sigma$ level). These initial constraints compare favorably to those obtained via microlensing in quasars and approach the sensitivity of constraints via relativistic reflection in neutron stars. This sensitivity will increase with further exposures, as well as the launch of future microcalorimeter and grating missions.

We measure the central kinematics for the dwarf spheroidal galaxy Leo I using integrated-light measurements and previously published data. We find a steady rise in the velocity dispersion from $300^{\prime\prime}$ into the center. The integrated-light kinematics provide a velocity dispersion of $11.76\pm0.66$ km/s inside $75^{\prime\prime}$. After applying appropriate corrections to crowding in the central regions, we achieve consistent velocity dispersion values using velocities from individual stars. Crowding corrections need to be applied when targeting individual stars in high density stellar environments. From integrated light, we measure the surface brightness profile and find a shallow cusp towards the center. Axisymmetric, orbit-based models measure the stellar mass-to-light ratio, black hole mass and parameters for a dark matter halo. At large radii it is important to consider possible tidal effects from the Milky Way so we include a variety of assumptions regarding the tidal radius. For every set of assumptions, models require a central black hole consistent with a mass $(3.3 \pm 2) \times 10^6\, M_\odot$. The no-black-hole case for any of our assumptions is excluded at over 95% significance, with $6.4<\Delta\chi^2<14$. A black hole of this mass would have significant effect on dwarf galaxy formation and evolution. The dark halo parameters are heavily affected by the assumptions for the tidal radii, with the circular velocity only constrained to be above 30 km/s. Reasonable assumptions for the tidal radius result in stellar orbits consistent with an isotropic distribution in the velocities. These more realistic models only show strong constraints for the mass of the central black hole.

Active galactic nuclei (AGN) at the centres of galaxies can cycle between periods of activity and of quiescence. Characterising the duty-cycle of AGN is crucial for understanding their impact on the evolution of the host galaxy. For radio AGN, their evolutionary stage can be identified from a combination of morphological and spectral properties. We summarise the results we have obtained in the last few years by studying radio galaxies in various crucial phases of their lives, such as remnant and restarted sources. We used morphological information derived from LOw Frequency ARray (LOFAR) images at 150 MHz, combined with resolved spectral indices maps, obtained using recently released images at 1400 MHz from the APERture Tile In Focus (Apertif) phased-array feed system installed on the Westerbork Synthesis Radio Telescope. Our study, limited so far to the Lockman Hole region, has identified radio galaxies in the dying and restarted phases. We found large varieties in their properties, relevant for understanding their evolutionary stage. We started by quantifying their occurrences, the duration of the 'on' (active) and 'off' (dying) phase, and we compared the results with models of the evolution of radio galaxies. In addition to these extreme phases, the resolved spectral index images can also reveal interesting secrets about the evolution of apparently normal radio galaxies. The spectral information can be connected with, and used to improve, the Fanaroff--Riley classification, and we present one example of this, illustrating what the combination of the LOFAR and Apertif surveys now allow us to do routinely.

We search for the signature of shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the recently released Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass <M200m> = 1e14.9 msol and redshift <z> = 0.55. We analyze in parallel a set of zoom-in hydrodynamical simulations from The Three Hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = $1.08 \pm 0.09$, measured at $3.1\sigma$ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = $4.58 \pm 1.24$ at $2.0\sigma$ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyze clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (Rsh_acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, Rsh_acc/Rsp > $2.16 \pm 0.59$.

The X-ray emission of neutron stars enables a probe of their temperatures, geometries, and magnetospheric properties. The current number of X-ray emitting pulsars is insufficient to rule out observational biases that may arise from poorly known distance, age, or location of the neutron stars. One approach to overcome such biases is to create a distance-limited sample with sufficiently deep observations. With the aim of better sampling of the nearby ($\lesssim$ 2kpc) neutron stars population, we started an XMM-Newton survey of pulsars to measure their X-ray fluxes or derive respective constraining upper limits. We investigated 14 nearby pulsars for potential X-ray counterparts using different energy bands and detectors. In addition to our new XMM-Newton data, we also considered archival data and catalogs. We discuss source properties and also check for alternative counterparts to the X-ray sources. In our new XMM-Newton data, we found two pulsar counterpart candidates with significance above over $4\sigma$ and one candidate with $3.5\sigma$ by combining EPIC camera detection likelihoods. We also report the detection of potential X-ray counterparts to 8 radio pulsars in the 4XMM-DR10 catalog which have not been reported in the literature.

We present resolved HI synthesis maps from the Very Large Array (VLA) of three interacting dwarf systems: the NGC 3664 dwarf pair, the NGC 3264 dwarf pair, and the UGC 4638 dwarf triplet. All three dwarf systems are captured at various stages of interaction and span a range of environments. We detect clear hallmarks of tidal interactions through the presence of HI bridges, and diffuse HI extensions that surround the dwarfs. We overlay the HI data on Pan-STARRS r-band images and find further evidence of tidal interactions through coincident distorted HI and tidal stellar features in NGC 3264 and UGC 4638, and an unwound spiral arm pointing towards its smaller companion in NGC 3264. In UGC 4638, both the gas and diffuse stars are extended to similar radii east of the primary, which could indicate that the smaller dwarf in the system has already completed one pass through the primary. We additionally find that our three systems, and those from the Local Volume TiNy Titans survey, are not HI deficient and thus the interaction has not resulted in a loss of gas from the systems. A comparison with non-interacting dwarf galaxies shows that the interactions have a significant impact on the kinematics of the systems. Our new resolved HI kinematics, combined with detailed stellar and HI morphologies, provide crucial constraints for future dynamical modelling of hierarchical mergers and the baryon cycle at the low-mass scale.

Local measurements of the Hubble parameter obtained from the distance ladder at low redshift are in tension with global values inferred from cosmological standard rulers. A key role in the tension is played by the assumptions on the cosmological history, in particular on the origin of dark energy. Here we consider a scenario where dark energy originates from the amplification of quantum fluctuations of a light field in inflation. We show that spatial correlations inherited from inflationary quantum fluctuations can reduce the Hubble tension down to one standard deviation, thus relieving the problem with respect to the standard cosmological model. Upcoming missions, like Euclid, will be able to test the predictions of models in this class.

The Geneva 7-colour photometric system is successfully applied to the study of various astrophysical objects. It measures the slope of the Paschen continuum, the Balmer discontinuity, and blocking absorption due to hydrogen or metallic lines. One of its greatest strengths is its intrinsic homogeneity. A new catalogue of the available measurements was generated, 30 years after the last publication. The identifications for the individual stars were cross-checked on the basis of the $Gaia$ and 2MASS catalogues. The high precision coordinates together with proper motions (if available) are included, for the first time, in the catalogue. Special caution was exercised with binaries and high-proper motion stars. The catalogue includes 42911 entries of highly accurate photometry. The data of this catalogue can be used for various applications such as new calibrations of astrophysical parameters, the standardization of new observations, and as additional information for ongoing and forthcoming all-sky surveys, such as the Transiting Exoplanet Survey Satellite (TESS).

We present details of the design, simulation, and initial test results of prototype detectors for the fourth-generation receiver of the South Pole Telescope (SPT). Optimized for the detection of key secondary anisotropies of the cosmic microwave background (CMB), SPT-3G+ will measure the temperature and polarization of the mm/sub-mm sky at 220, 285, and 345 GHz, beyond the peak of the CMB blackbody spectrum. The SPT-3G+ focal plane will be populated with microwave kinetic inductance detectors (MKIDs), allowing for significantly increased detector density with reduced cryogenic complexity. We present simulation-backed designs for single-color dual-polarization MKID pixels at each SPT-3G+ observation frequency. We further describe design choices made to promote resonator quality and uniformity, enabling us to maximize the available readout bandwidth. We also discuss aspects of the fabrication process that enable rapid production of these devices and present an initial dark characterization of a series of prototype devices.

Although radio observations have been historically seen as less valuable than optical observations, today's broadband radio spectra of peaked spectrum sources reveal detailed physics from within the inner region of the galaxy, on spatial scales beyond what an optical telescope can resolve. Peaked radio spectra are thought to be evolving into large scale radio galaxies, although an over-abundance of the most compact sources reveals that a significant fraction are confined within their host galaxies. Furthermore, at the lowest luminosities, these sources are largely unknown, and may reveal the small scale precursors of FR-I galaxies. Here I summarise the previous work exploring the properties of low luminosity peaked radio sources, and the future work that extends on this within even deeper radio observations of well studied fields.

We report the radio continuum properties for several samples of peaked spectrum (PS) radio sources. Broadband spectra of the objects were analysed using the RATAN-600 six-frequency (1.2-22 GHz) observations and available literature data, obtained on a time scale of 20-30 years. We discuss statistical differences in radio properties for several AGN types with peaked spectra and focus on PS quasars at high redshifts ($z > 3$). We confirm that a relatively small fraction (1-2%) of bright PS sources can be considered as genuine GPSs when they have been monitored densely and for a long time. The contamination of GPS source samples by blazars is getting stronger as the redshift increases, and we confirm that it is underestimated due to lack of systematic multifrequency observations.

We present a simple but powerful technique for the analysis of polarized emission from radio galaxies and other objects. It is based on the fact that images of Stokes parameters often contain considerably more information than is available in polarized intensity and angle maps. In general, however, the orientation of the Stokes parameters will not be matched to the position angles of structures in the source. Polarization tomography, the technique presented in this paper, consists of making a series of single linear Stokes parameter images, $S(\rho)$, where each image is rotated by an angle $\rho$ from the initial orientation of Q and U. Examination of these images, in a series of still frames or a movie, reveals often hidden patterns of polarization angles, as well as structures that were obscured by the presence of overlapping polarized emission. We provide both cartoon examples and a quick look at the complex polarized structure in Cygnus A.

A new silicon monoxide ($^{28}$Si$^{16}$O) line list covering infrared, visible and ultraviolet regions called SiOUVenIR is presented. This line list extends the infrared EBJT ExoMol line list by including vibronic transitions to the $A\,{}^{1}\Pi$ and $E\,{}^{1}\Sigma^{+}$ electronic states. Strong perturbations to the $A\,{}^{1}\Pi$ band system are accurately modelled through the treatment of 6 dark electronic states: $C\,{}^{1}\Sigma^{-}$, $D\,{}^{1}\Delta$, $a\,{}^{3}\Sigma^{+}$, $b\,{}^{3}\Pi$, $e\,{}^{3}\Sigma^{-}$ and $d\,{}^{3}\Delta$. Along with the $X\,{}^{1}\Sigma^{+}$ ground state, these 9 electronic states were used to build a comprehensive spectroscopic model of SiO using a combination of empirical and ab initio curves, including the potential energy (PE), spin-orbit (SO), electronic angular momentum (EAM) and (transition) dipole moment curves. The ab initio PE and coupling curves, computed at the multireference configuration interaction (MRCI) level of theory, were refined by fitting their analytical representations to 2617 experimentally derived SiO energy levels determined from 97 vibronic bands belonging to the $X$-$X$, $E$-$X$ and $A$-$X$ electronic systems through the MARVEL procedure. 112 observed forbidden transitions from the $C$-$X$, $D$-$X$, $e$-$X$, and $d$-$X$ bands were assigned using our predictions, and these could be fed back into the MARVEL procedure. The SiOUVenIR line list was computed using published ab initio transition dipole moments for the $E$-$X$ and $A$-$X$ bands; the line list is suitable for temperatures up to 10,000 K and for wavelengths longer than 140 nm. SiOUVenIR is available from www.exomol.com and the CDS database.

The light curve of the microlensing event KMT-2021-BLG-0912 exhibits a very short anomaly relative to a single-lens single-source form. We investigate the light curve for the purpose of identifying the origin of the anomaly. We model the light curve under various interpretations. From this, we find four solutions, in which three solutions are found under the assumption that the lens is composed of two masses (2L1S models), and the other solution is found under the assumption that the source is comprised of a binary-star system (1L2S model). The 1L2S model is ruled out based on the contradiction that the faint source companion is bigger than its primary, and one of the 2L1S solutions is excluded from the combination of the relatively worse fit, blending constraint, and lower overall probability, leaving two surviving solutions with the planet/host mass ratios of $q\sim 2.8\times 10^{-5}$ and $\sim 1.1\times 10^{-5}$. A subtle central deviation supports the possibility of a tertiary lens component, either a binary companion to the host with a very large or small separation or a second planet lying near the Einstein ring, but it is difficult to claim a secure detection due to the marginal fit improvement, lack of consistency among different data sets, and difficulty in uniquely specifying the nature of the tertiary component. With the observables of the event, it is estimated that the masses of the planet and host are $\sim (6.9~M_\oplus, 0.75~M_\odot)$ according to one solution and $\sim (2.8~M_\oplus, 0.80~M_\odot)$ according to the other solution, indicating that the planet is a super Earth around a K-type star, regardless of the solution.

Despite their ubiquity throughout the Universe, the formation of S0 galaxies remains uncertain. Recent observations have revealed that S0 galaxies make up a diverse population which is difficult to explain with a single formation pathway, suggesting that the picture of how these galaxies form is more complicated than originally envisioned. Here we take advantage of the latest hydrodynamical cosmological simulations and follow up these studies with an investigation into the formation histories of S0s in IllustrisTNG. We first classify IllustrisTNG galaxies in a way which is fully consistent with the observations, and reproduce the observed photometric and environmental distributions seen for the S0 population. We then trace the formation histories of S0 galaxies back through time, identifying two main distinct pathways; those which experienced gas stripping via group infalls (37 percent of S0s) or significant merger events (57 percent). We find that those forming via mergers feature a transient star-forming ring, whose present-day occurrence rate matches observations. We find that these formation pathways together can reproduce the range in rotational support in observed S0s, concluding that there are two main formation pathways for S0 galaxies.

We examine the cosmological constraining power from two cross-correlation probes between galaxy and CMB surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langle\delta_g\kappa_{\rm CMB}\rangle$ and source galaxy weak lensing shear with CMB lensing convergence $\langle\gamma\kappa_{\rm CMB}\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to systematic effects that could be present in galaxy-only or CMB-only analyses. We study how the constraining power of $\langle\delta_g\kappa_{\rm CMB}\rangle+\langle\gamma\kappa_{\rm CMB}\rangle$ changes as the galaxy samples and CMB dataset qualitatively change from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, cosmological constraints from $\langle\delta_g\kappa_{\rm CMB}\rangle$ and $\langle\gamma\kappa_{\rm CMB}\rangle$ are comparable in Stage-III; but as we move to Stage-IV, shot noise from the galaxy density becomes subdominant to cosmic variance, preventing the contribution from $\langle\delta_g\kappa_{\rm CMB}\rangle$ to further improve. This implies that to maximize the cosmological constraints from future $\langle\delta_g\kappa_{\rm CMB}\rangle+\langle \gamma\kappa_{\rm CMB}\rangle$ analyses, we should focus more on the requirements on $\langle\gamma \kappa_{\rm CMB}\rangle$ instead of $\langle\delta_g\kappa_{\rm CMB}\rangle$. In addition, the selection of the lens galaxy sample in $\langle\delta_g\kappa_{\rm CMB}\rangle$ should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

We present high resolution observations of CS(J=1-0), H13CO+ (J=1-0), and SiO(v=0:J=1-0) lines, together with the 49GHz and 86GHz continuum emissions, toward W49N carried out with Nobeyama Millimeter Array. We identified 11 CS, 8 H13CO+, and 6 SiO clumps with radii of 0.1-0.5pc. The CS and H13CO+ clumps are mainly divided into two velocity components, one at 4kms-1 and the other at 12kms-1, while the SiO clumps have velocities between the two components. The SiO emission is distributed toward the UCHII ring, where the 4kms-1 component clumps of CS and H13CO+ also exist. The 12kms-1 component clumps of CS are detected at the east and west of the UCHII ring with an apparent hole toward the ring. The clump masses vary from 4.4x10^2 M_SUN to 4.9x10^4 M_SUN with the mean values of 0.94x10^4M_SUN, 0.88x10^4M_SUN, and 2.2x10^4M_SUN for the CS, H13CO+, and SiO clumps, respectively. The total masses derived from CS, H13CO+, and SiO clumps are 1.0x10^5M_SUN, 0.70x10^5M_SUN, and 1.3x10^5 M_SUN, respectively, which agree well with the corresponding virial masses of 0.71x10^5M_SUN, 1.3x10^5M_SUN, and 0.88x10^5M_SUN, respectively. The average molecular hydrogen densities of the clumps are 0.90x10^6 cm-3, 1.4x10^6cm-3, and 7.6x10^6 cm-3 for the CS, H13CO+ and SiO clumps, respectively. The density derived from the SiO clumps seems significantly higher than those from the others, probably because the SiO emission is produced in high density shocked regions. The free fall time scale of the clumps is estimated to be ~3x10^4 yr, which gives an accretion rate of 3x10-3-1M_SUN yr-1 onto a stellar core. The observed clumps are, if they are undergoing free fall, capable of producing dozens of massive stars in the next 10^5 yr. We propose a view that pre-existing two clouds collided with each other almost face-on to produce the observed clumps and triggered the burst of massive star formation in W49N.

In this paper, we report white dwarfs identified in the 5th Data Release of the Large Area Multi-Object fibre Spectroscopic Telescope, including spectral types of DA, DB, DC, DZ, and so on. There are 2 625 DA spectra of 2 281 DA stars, 182 DB spectra of 166 DB stars, 62 DC spectra of 58 DC stars, 36 DZ spectra of 33 DZ stars and many other types identified, in addition to our previous paper (Data Release 2). Among those sources, 393 DA stars and 46 DB stars are new identifications after cross-matching with the literature. In order to select DA candidates, we use the classification result from the LAMOST pipeline, colour-colour cut method and a random forest machine learning method. For DBs, since there is no template for DB in the pipeline model, a random forest machine learning method is chosen to select candidates. All the WD candidates have been visually checked individually. The parameters of effective temperature, surface gravity, mass, and cooling age have been estimated for relatively high signal-to-noise ratio DAs and DBs. The peaks of the DA and DB mass distributions are found to be around 0.62Msun and 0.65Msun, respectively. Finally, the data and method we used to select white dwarf candidates for the second phase of LAMOST survey are also addressed in this paper.

Planets with large moon(s) or those in the habitable zone of low-mass stars may experience much stronger tidal force and tide-induced ocean mixing than that on Earth. Thus, the vertical diffusivity (or, more precisely, diapycnal diffusivity) on such planets, which represents the strength of vertical mixing in the ocean, would be greater than that on Earth. In this study, we explore the effects of extremely high diffusivity on the ocean circulation and surface climate of Earth-like planets in one asynchronous rotation orbit. The response of planetary climate to 10 and 100 times greater vertical diffusivity than that found on Earth is investigated using a fully coupled atmosphere-ocean general circulation model. In order to perform a clear comparison with the climate of modern Earth, Earth's orbit, land-sea configuration, and present levels of greenhouse gases are included in the simulations. We find that a larger vertical diffusivity intensifies the meridional overturning circulation (MOC) in the ocean, which transports more heat to polar regions and melts sea ice there. Feedback associated with sea ice, clouds, and water vapor act to further amplify surface warming. When the vertical diffusivity is 10 (100) times the present-day value, the magnitude of MOC increases by $\approx$3 (18) times, and the global-mean surface temperature increases by $\approx$4$^{\circ}$C (10$^{\circ}$C). This study quantifies the climatic effect of an extremely strong vertical diffusivity and confirms an indirect link between planetary orbit, tidal mixing, ocean circulation, and surface climate. Our results suggest a moderate effect of varying vertical ocean mixing on planetary climate.

Solar eruptions are explosive release of coronal magnetic field energy as manifested in solar flares and coronal mass ejection. Observations have shown that the core of eruption-productive regions are often a sheared magnetic arcade, i.e., a single bipolar configuration, and, particularly, the corresponding magnetic polarities at the photosphere are elongated along a strong-gradient polarity inversion line (PIL). It remains unclear what mechanism triggers the eruption in a single bipolar field and why the one with a strong PIL is eruption-productive. Recently, using high accuracy simulations, we have established a fundamental mechanism of solar eruption initiation that a bipolar field as driven by quasi-static shearing motion at the photosphere can form an internal current sheet, and then fast magnetic reconnection triggers and drives the eruption. Here we investigate the behavior of the fundamental mechanism with different photospheric magnetic flux distributions, i.e., magnetograms, by combining theoretical analysis and numerical simulation. Our study shows that the bipolar fields of different magnetograms, as sheared continually, all exhibit similar evolutions from the slow storage to fast release of magnetic energy in accordance with the fundamental mechanism, which demonstrates the robustness of the mechanism. We further found that the magnetograms with stronger PIL produce larger eruptions, and the key reason is that the sheared bipolar fields with stronger PIL can achieve more non-potentiality, and their internal current sheet can form at a lower height and with a larger current density, by which the reconnection can be more efficient. This also provides a viable trigger mechanism for the observed eruptions in active region with strong PIL.

A simple analytical argument is proposed for a possible explanation of the characteristics of the lateral shower age ($s$) of proton ($p$)/nuclei-initiated showers. The analytical argument states that lateral density distribution (LDD) of electrons of a $p$-initiated shower is due to superposition of several electromagnetic (EM) sub-showers developed at a very early stage in the atmosphere from the decay of neutral pions ($\pi^{0}$s). Thanks to the superposition property of the electron LDD in a $p$ shower, a plausible analytical parametrization has been worked out by giving well represented analytic function for the electron LDDs of $p$- and $\pi^{0}$-initiated showers. Based on cosmic ray extensive air shower simulations, we have validated how the various characteristics of $s$ can be understood in the context of the present analytical argument. The $s$ parameter of a $p$ shower and its correlations with the shower ages of electron- and $\pi^{0}$-initiated showers supports the idea that the result of superposition of several EM sub-showers initiated by $\pi^{0}$s with varied energies at a very early stage might produce the LDD of electrons of a $p$ shower. It is also noticed with the simulated data that the stated feature still persists concerning the notion of the local shower age parameter.

In this work we investigate the possibility of transporting material to the NEO region via the 8:3 MMR with Jupiter, potentially even material released from the dwarf planet Ceres. By applying the FLI map method to the 8:3 MMR region in the orbital plane of Ceres, we were able to distinguish between stable and unstable orbits. Subsequently, based on the FLI maps (for mean anomaly $M=60^\circ$ and also $M=30^\circ$), 500 of the most stable and 500 of the most unstable particles were integrated for $15\,Myr$ for each map. Long-term integration in the case of $M=60^\circ$ showed that most of the stable particles evolved, in general, in uneventful ways with only 0.8\% of particles reaching the limit of q $\leq$ 1.3 $AU$. However, in the case of $M=30^\circ$, a stable evolution was not confirmed. Over 40\% of particles reached orbits with q $\leq$ 1.3 $AU$ and numerous particles were ejected to hyperbolic orbits or orbits with a > 100 $AU$. The results for stable particles indicate that short-term FLI maps are more suitable for finding chaotic orbits, than for detecting the stable ones. A rough estimate shows that it is possible for material released from Ceres to get to the region of 8:3 MMR with Jupiter. A long-term integration of unstable particles in both cases showed that transportation of material via 8:3 MMR close to the Earth is possible.

The spatial distribution between dark matter and baryonic matter of the Universe is biased or deviates from each other. In this work, by comparing the results derived from IllustrisTNG and WIGEON simulations, we find that many results obtained from TNG are similar to those from WIGEON data, but differences between the two simulations do exist. For the ratio of density power spectrum between dark matter and baryonic matter, as scales become smaller and smaller, the power spectra for baryons are increasingly suppressed for WIGEON simulations; while for TNG simulations, the suppression stops at $k=15-20h{\rm Mpc}^{-1}$, and the power spectrum ratios increase when $k>20h{\rm Mpc}^{-1}$. The suppression of power ratio for WIGEON is also redshift-dependent. From $z=1$ to $z=0$, the power ratio decreases from about 70% to less than 50% at $k=8h{\rm Mpc}^{-1}$. For TNG simulation, the suppression of power ratio is enhanced with decreasing redshifts in the scale range $k>4h{\rm Mpc}^{-1}$, but is nearly unchanged with redshifts in $k<4h{\rm Mpc}^{-1}$ These results indicate that turbulent heating can also have the consequence to suppress the power ratio between baryons and dark matter. Regarding the power suppression for TNG simulations as the norm, the power suppression by turbulence for WIGEON simulations is roughly estimated to be 45% at $k=2h{\rm Mpc}^{-1}$, and gradually increases to 69% at $k=8h{\rm Mpc}^{-1}$, indicating the impact of turbulence on the cosmic baryons are more significant on small scales.

Context. Most of the modeling of interstellar dust infrared emission spectrum is done by assuming some variations around a single temperature greybody approximation. For example, the foreground modeling of Planck mission maps involves a unique dust temperature T along a given line-of-sight with a unique emissivity index \b{eta}. The two parameters are then fitted and therefore variable from one line-of-sight to the other. Aims. Our aim is to go beyond that modeling in an economical way. Methods. We model the dust spectrum with a temperature distribution around the mean value and show that only the second temperature moment matters. Following Pitrou and Stebbins (2014), we advocate the use of the temperature logarithm as the proper variable. Results. If the interstellar medium is not too heterogeneous, there is a universal analytical spectrum, which is derived here, that goes beyond the greybody assumption. We show how the Cosmic Microwave Background radiatively interacts with the dust spectrum (a non-negligible corrective term at millimeter wavelengths). Finally, we construct a universal ladder of discrete temperatures which gives a minimal and fast description of dust emission spectra as measured by photometric mapping instruments that lends itself to an almost linear fitting. This data modeling can include contributions from the Cosmic Infrared Background fluctuations

Spherical asymmetries are prevalent within the outflows of AGB stars. Since binary interaction with a stellar or planetary companion is thought to be the underlying mechanism behind large-scale structures, we included the effects of UV radiation originating from a stellar companion in our chemical kinetics model. The one-dimensional model provides a first approximation of its effects on the chemistry throughout the outflow. The presence of a close-by stellar companion can strongly influence the chemistry within the entire outflow. Its impact depends on the intensity of the radiation (set by the stellar radius and blackbody temperature) and on the extinction the UV radiation experiences (set by the outflow density, density structure, and assumed radius of dust formation). Parent species can be photodissociated by the companion, initiating a rich photon-driven chemistry in the inner parts of the outflow. The outcome depends on the balance between two-body reactions and photoreactions. If two-body reactions dominate, chemical complexity within the outflow increases. This can make the abundance profiles of daughters appear like those of parents, with a larger inner abundance and a gaussian decline. If photoreactions dominate, the outflow can appear molecule-poor. We model three stellar companions. The impact of a red dwarf companion is limited. Solar-like companions show the largest effect, followed by a white dwarf. A stellar companion can also lead to the formation of unexpected species. The outflow's molecular content, especially combined with abundance profiles, can indicate a stellar companion's presence. Our results pave the way for further outflow-specific (three-dimensional) model development.

The Earth Trojans are co-orbitals librating around the Lagrange points $L_4$ or $L_5$ of the Sun-Earth system. Although many numerical studies suggest that they can maintain their dynamical status and be stable on timescales up to a few tens of thousands of years or even longer, they remain an elusive population. Thus far only one transient member (2010 TK$_7$) has been discovered serendipitously. Here, we present a dynamical study of asteroid 2020 XL$_5$. With our meticulous followup astrometric observations of the object, we confirmed that it is a new Earth Trojan. However, its eccentric orbit brings it close encounters with Venus on a frequent basis. Based on our N-body integration, we found that the asteroid was captured into the current Earth Trojan status in the 15th century, and then it has a likelihood of 99.5% to leave the $L_4$ region within the next $\sim$10 kyr. Therefore, it is most likely that 2020 XL$_5$ is dynamically unstable over this timescale.

In this work we explore and test new formulations of cosmological scenarios in $f(Q)$ theories. In these settings, the non-metricity scalar ($Q$) is the main source of gravity and Friedmann equations are modified to account for the associated degrees of freedom. This work focuses first on the derivation, and then theoretical and observational analysis of two such (new) exact cosmological models; they both display a non-standard behaviour in which an additional parameter encoding non-metricity effects acts in the fashion of a screened cosmological constant. One of the new settings has the same background evolution as the well know DGP cosmological model, while the other resembles the former considerably, although its origin is purely phenomenological. We use the Markov Chain Montecarlo method combined with standard statistical techniques to perform observational astrophysical tests relying upon background data, specifically these are Type Ia Supernovae luminosities and direct Hubble data (from cosmic clocks), along with Cosmic Microwave Background shift and Baryon Acoustic Oscillations data. In addition, we compute some of the cosmographic parameters and other discriminators with the purpose of refining our knowledge about these models in the light of their theoretical and observational signatures, and this allows for a better comparison with the (concordance) $\Lambda$CDM setup. We conclude that these scenarios do not show signatures indicating a departure from the $\Lambda$CDM behaviour.

We present CoLoRe, a public software package to efficiently generate synthetic realisations of multiple cosmological surveys. CoLoRe can simulate the growth of structure with different degrees of accuracy, with the current implementation supporting lognormal fields, first, and second order Lagrangian perturbation theory. CoLoRe simulates the density field on an all-sky light-cone up to a desired maximum redshift, and uses it to generate multiple 2D and 3D maps: galaxy positions and velocities, lensing (shear, magnification, convergence), integrated Sachs-Wolfe effect, line intensity mapping, and line of sight skewers for simulations of the Lyman-$\alpha$ forest. We test the accuracy of the simulated maps against analytical theoretical predictions, and showcase its performance with a multi-survey simulation including DESI galaxies and quasars, LSST galaxies and lensing, and SKA intensity mapping and radio galaxies. We expect CoLoRe to be particularly useful in studies aiming to characterise the impact of systematics in multi-experiment analyses, quantify the covariance between different datasets, and test cross-correlation pipelines for near-future surveys.

We show that upon applying Palatini $f(R)$, characterized by an $\alpha R^2$-term, within a scenario motivated by a temporal variation of strong coupling constant, then one gets a quadratic kinetic energy. We do not drop this term, but rather study two extreme cases: $\alpha <<1$ and $\alpha >>1$. In both cases one can generate an inflationary paradigm, legitimately called thus a k-inflation, which fits the Planck 2018 data.

PSZ2 G091.83+26.11 is a massive galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile [1]. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at higher redshifts (where we expect the fraction of merging objects to be higher), it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of PSZ2 G091.83+26.11 by the NIKA2 camera to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption.

We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relativity. We first perform a series of numerical tests to show that the gravitational wave (GW) signal obtained is in excellent agreement with the one expected from theory. This benchmark is done for well studied scenarios (such as binary stellar systems). We then apply our code to calculate the GW signals from tidal disruption events (TDEs), finding that our results are consistent with the theoretical estimates obtained in previous studies for selected parameters. We illustrate interesting results from the catalogue, where we stress how the gravitational signal is affected by variations of some parameters (like black hole spin, stellar orbital eccentricity and inclination). The full catalogue is available online. It is intended to be a living catalogue.

We estimate the mass accretion rate (MAR) of the 346 clusters of galaxies in the HectoMAP Cluster Survey. The clusters span the redshift range $0.17-0.42$ and the $M_{200}$ mass range $\approx (0.5 - 3.5)\cdot 10^{14}$M$_\odot$. The MAR estimate is based on the caustic technique along with a spherical infall model. Our analysis extends the measurement of MARs for 129 clusters at $z<0.3$ from the Cluster Infall Regions in the Sloan Digital Sky Survey (CIRS) and the Hectospec Cluster Survey (HeCS) to redshift $z \sim 0.42$. Averaging over redshift, low-mass clusters with $M_{200}\sim 0.7\cdot 10^{14}$M$_\odot$ accrete $\sim 3\cdot 10^4$M$_\odot$yr$^{-1}$; more massive clusters with $M_{200}\sim 2.8\cdot 10^{14}$M$_\odot$ accrete $\sim 1\cdot 10^5$M$_\odot$yr$^{-1}$. Low- and high-mass clusters increase their MAR by $\sim 46\%$ and $\sim 84\%$, respectively, as the redshift increases from $z\sim 0.17-0.29$ to $z\sim 0.34-0.42$. The MARs at fixed redshift increase with mass and MARs at fixed mass increase with redshift in agreement with $\Lambda$CDM cosmological model for hierarchical structure formation. We consider the extension of MAR measurements to $z \sim 1$.

The geology of Earth and super-Earth sized planets will, in many cases, only be observable via their atmospheres. Here, we use the creation of volcanic atmospheres as a key window into planetary geochemistry. We couple volcanic outgassing with atmospheric chemistry models to simulate the growth of C-O-H-S-N atmospheres in thermochemical equilibrium, aiming to establish what information about the planet's mantle fO2 and bulk silicate H/C ratio can be determined by atmospheric observation. Warm (800 K) volcanic atmospheres develop distinct compositional groups as the mantle fO2 is varied, which can be identified using sets of (often minor) indicator species: Class O, representing an oxidised mantle and containing SO2 and sulfur allotropes; Class I, formed by intermediate mantle fO2's and containing CO2, CH4, CO and COS; and Class R, produced by reduced mantles, containing H2, NH3 and CH4. These atmospheric classes are largely independent of the bulk silicate H/C ratio. However, the H/C ratio does affect the dominant atmospheric constituent, which can vary between H2, H2O, CO2 and CH4 once the chemical composition has stabilised to a point where it no longer changes substantially with time. This final state is dependent on the mantle fO2, the H/C ratio, and time since the onset of volcanism. Superchondritic H/C enrichment to the level of Earth (H/C = 0.99 +/- 0.42) and higher can only be inferred for planets with reduced mantles producing Class R atmospheres. On warm, volcanically active planets, mantle fO2 could be identifiable from atmospheric observations using JWST.

The current studies of cosmic rays are focused on most energetic particles entering the atmosphere and producing a single Extensive Air Shower (EAS). There are, however, models predicting that interactions of high energy particles may result in Cosmic-Ray Ensembles (CRE) created far from the Earth. They could be observed as some number of correlated air showers of relatively low energies spread over a large area. The objective of the Cosmic Ray Extremely Distributed Observatory (CREDO) is to search for CRE using all available data from different detectors and observatories including even small but numerous detectors spread over large areas. Interpretation of such measurements require precise information on properties of EAS in a very wide energy spectrum. Low energy EAS are analysed using events from CORSIKA, the program performing air shower simulations. The primary cosmic ray particle energy range extends from 1TeV up to 4000TeV. The secondary particles at the ground level are studied in order to obtain their density fluctuations and correlations in location. Although the fluctuations observed in multiplicity distributions are consistent with random the more detailed analysis reveals that near a selected particle the density of other particles is enhanced over that expected in the absence of correlations. The results of this analysis may be useful in further calculations, for example to obtain probability of detection of an EAS without special simulations.

The potential detection of ppb levels phosphine (PH3) in the clouds of Venus through millimeter-wavelength astronomical observations is extremely surprising as PH3 is an unexpected component of an oxidized environment of Venus. A thorough analysis of potential sources suggests that no known process in the consensus model of Venus' atmosphere or geology could produce PH3 at anywhere near the observed abundance. Therefore, if the presence of PH3 in Venus' atmosphere is confirmed, it is highly likely to be the result of a process not previously considered plausible for Venusian conditions. The source of atmospheric PH3 could be unknown geo- or photochemistry, which would imply that the consensus on Venus' chemistry is significantly incomplete. An even more extreme possibility is that strictly aerial microbial biosphere produces PH3. This paper summarizes the Venusian PH3 discovery and the scientific debate that arose since the original candidate detection one year ago.

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, the binary tidal parameter at the threshold ($\Lambda_{th}$), the maximum mass of nonrotating NSs, and the radii of reference mass NSs. We correct the systematic errors in previously reported fitting coefficients that were obtained with approximate general-relativity simulations. We combine the EOS-insensitive relations, phenomenological constraints on NS properties and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and 1.6 M$_\odot$ NS of 9.81 km and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability ($\Lambda$) of a 1.4 $M_\odot$ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters -- fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We introduce new methods to constrain the upper as well as lower limit of NS maximum mass using future BNS detections and their identification as prompt or delayed collapse. With future observations it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

We present a new methodology to analyse in a comprehensive way large-scale and supernovae (or any other standard candle) surveys. Our approach combines galaxy and supernova position and redshift data with supernova peculiar velocities, obtained through their magnitude scatter, to construct a 6x2pt analysis which includes six power spectra. The 3x3 correlation matrix of these spectra expresses exhaustively the information content of the surveys at the linear level. We proceed then to forecast the performance of future surveys like LSST and 4MOST with a Fisher Matrix analysis, adopting both a model-dependent and a model-independent approach. We compare the performance of the 6x2pt approach to the traditional one using only galaxy clustering and some recently proposed combinations of galaxy and supernovae data and quantify the possible gains by optimally extracting the linear information. We show that the 6x2pt method shrinks the uncertainty area in the $\sigma_8, \gamma$ plane by more than half when compared to the traditional method. The combined clustering and velocity data on the growth of structures has uncertainties at similar levels to those of the CMB but exhibit orthogonal degeneracies, and the combined constraints yield improvements of factors of 5 in each of the five cosmological parameters here considered. Concerning the model-independent results, we find that our method can improve the constraints on $H(z)/H_0$ in all redshift bins by more than 70% with respect to the galaxy clustering alone, reaching a precision of 3--4% at high redshifts.

Following optical pulses ($\lambda=405~\text{nm}$) on titanium nitride (TiN) Microwave Kinetic Inductance Detectors (MKIDs) cooled down at temperatures $T \le T_c / 20$ ($T_c \simeq 4.6~\text{K}$), we observe a large phase-response highlighting two different modes simultaneously that are nevertheless related. The first corresponds to the well-known transition of cooper-pair breaking into quasi-particles which produces a known phase response. This is immediately followed by a large inverse response lasting several hundreds of microseconds to several milliseconds depending on the temperature. We propose to model this inverse pulse as the thermal perturbation of the superconductor and interaction with two level system (TLS) that reduces the dielectric constant which in turns modify the capacitance and therefore the resonance frequency. The ratio of the TLS responding to the illumination is on the order of that of the area of the inductor to the whole resonator

Our knowledge about the photometric properties of bulges in late-type galaxies (LTGs) is founded upon image decomposition into a S\'ersic model for the central luminosity excess of the bulge and an exponential model for the underlying disk. We argue that the standard practice of adopting an exponential model for the disk all the way to its center is inadequate because it implicitly neglects the fact of star formation (SF) quenching (SFQ) in the centers of LTGs. Extrapolating the fit for the observable star-forming zone of the disk (outside the bulge) inwardly overestimates the true surface brightness of the disk in its SF-quenched central zone. We refer to this effect as Dio. The primary consequence of the neglect of Dio in bulge-disk decomposition studies is the oversubtraction of the disk underneath the bulge, leading to a systematic underestimation of the luminosity of the latter. Framed in the picture of galaxy downsizing and inside-out SFQ, Dio is expected to differentially impact galaxies across redshift and stellar mass M*, thus leading to systematic and complex biases in the scatter and slope of various galaxy scaling relations. We conjecture that correction for Dio will lead to a downbending of the bulge vs. super-massive black hole (SMBH) relation for galaxies below log(M*/Msolar)~10.7. A decreasing M(SMBH)/M* ratio with decreasing M* would help consistently explain the scarcity and weakness of accretion-powered nuclear activity in low-mass spiral galaxies. A well detectable Dio (~2 r mag) can emerge early on through inward migration of SF clumps from the disk in combination with a strong contrast of emission-line equivalent widths between the quenched proto-bulge and its SF periphery. Spatially resolved studies with the JWST, ELT, and Euclid could therefore offer key insights into the chronology and physical drivers of SFQ in the early phase of galaxy assembly. (abridged)

We report the results of new transit observations for the three hot Jupiter-like planets HATP-36b, HATP-56b and WASP-52b respectively. Transit timing variations (TTVs) are presented for these systems based on observations that span the period 2016 - 2020. The data were collected with the 0.6 m telescope at Adiyaman University (ADYU60, Turkey) and the 1.0 m telescope at T\"UB\.ITAK National Observatory (TUG, Turkey). Global fits were performed to the combined light curves for each system along with the corresponding radial velocity (RV) data taken from the literature. The extracted parameters (for all three systems) are found to be consistent with the values from previous studies. Through fits to the combined mid-transit times data from our observations and the data available in the literature, an updated linear ephemeris is obtained for each system. Although a number of potential outliers are noted in the respective O-C diagrams, the majority of the data are consistent within the 3$\sigma$ confidence level implying a lack of convincing evidence for the existence of additional objects in the systems studied.

We present a feasibility study for the high-redshift galaxy part of the Science Verification Campaign with the 220-440 GHz DESHIMA 2.0 integrated superconducting spectrometer on the ASTE telescope. The first version of the DESHIMA 2.0 chip has been recently manufactured and tested in the lab. Based on these realistic performance measurements, we evaluate potential target samples and prospects for detecting the [CII] and CO emission lines. The planned observations comprise two distinct, but complementary objectives: (1) acquiring spectroscopic redshifts for dusty galaxies selected in far-infrared/mm-wave continuum surveys; (2) multi-line observations to infer physical conditions in dusty galaxies.

Linear acceleration emission is one of the mechanisms proposed to explain the intense pulsar radio emissions. This mechanism is however not well understood due to a lack of its proper mathematical analyses, e.g., of the collective plasma response and the resulting emission power. We utilize 1D relativistic particle-in-cell simulations to derive the emission properties of two instabilities in neutron star magnetospheres, relativistic beam instability and interactions of plasma bunches/clouds. We found that the emission power by plasma bunch interactions exceeds emission due to streaming instability by seven orders of magnitude. The wave power generated by a plasma bunch interaction can be obtained as large as $3.4\times10^{19}$ W. It alone can account for the total radio power emitted by typical pulsars ($10^{18}-10^{22}$ W). The emission of the plasma bunch has a number of features of the observed pulsar radiation. Its spectrum is characterized by an almost flat profile for lower frequencies and a power-law with an index $\approx-2.5$ for higher frequencies. The angular width of the radiation decreases with increasing frequency. The generated wave power depends on the pulsar rotation angle. It can cause fine structures in the observed intensity as it fluctuates between positive and negative wave interference as a function of the emission angle.

B, V, i, and z bandpass observations were collected in late 2020 for three RRab type stars: UU Ceti, UW Gruis, and W Tucanae. The period-luminosity (PL) relationships of Catalen et al. (2004) and Caceres \& Catelan (2008) were applied to derive distances. These were found to be in reasonable agreement with the Gaia Early DR3 distances, lending confidence to use of the PL relationships. Fourier decompositions were applied to data from the TESS space telescope to derive, using stepwise linear regression, an empirical relationship between terms of the decomposition and the pulsation period with metallicity [Fe/H]. TESS data were available for UU Cet and W Tuc out of the three studied studies. The derived equation gave metallicities in line with the literature for both stars, lending confidence to their usage in the PL-derived distances.

The distribution of member stars in the surroundings of an Open Cluster (OC) can shed light on the process of its formation, evolution and dissolution. The analysis of structural parameters of OCs as a function of their age and position in the Galaxy brings constraints on theoretical models of cluster evolution. The Gaia catalogue is very appropriate to find members of OCs at large distance from their centers. We aim at revisiting the membership lists of OCs from the solar vicinity, in particular by extending these membership lists to the peripheral areas thanks to Gaia EDR3. We used the clustering algorithm HDBSCAN on Gaia parallaxes and proper motions to systematically look for members up to 50 pc from the cluster centers. We fitted a King's function on the radial density profile of these clusters and a Gaussian Mixture Model on their two dimensional distribution of members. We also evaluated the degree of mass segregation of the clusters. Our methodology performs well on 389 clusters out of the 467 selected ones. We report the detection of vast coronae around almost all the clusters and the detection of 71 OCs with tidal tails, multiplying by more than four the number of such structures identified. We find the size of the cores to be on average smaller for old clusters than for young ones. Also, the overall size of the clusters seems to slightly increase with age while the fraction of stars in the halo seems to decrease. As expected the mass segregation is more pronounced in the oldest clusters but a clear trend with age is not seen. OCs are more extended than previously expected, regardless of their age. The decrease in the proportion of stars populating the clusters halos highlights the different cluster evaporation processes and the short timescales they need to affect the clusters. Reported parameters all depend on cluster ages but can not be described as single functions of time.

It was found that satellites of nearby galaxies can form flattened co-rotating structures called disks of satellites or planes of satellites. Their existence is not expected by the current galaxy formation simulations in the standard dark-matter-based cosmology. On the contrary, modified gravity offers a promising alternative: the objects in the disks of satellites are tidal dwarf galaxies, that is small galaxies that form from tidal tails of interacting galaxies. After introducing the topic, we review here our work on simulating the formation of the disks of satellites of the Milky Way and Andromeda galaxies. The initial conditions of the simulation were tuned to reproduce the observed positions, velocities and disk orientations of the galaxies. The simulation showed that the galaxies had a close flyby 6.8Gyr ago. One of the tidal tails produced by the Milky Way was captured by Andromeda. It formed a cloud of particles resembling the disk of satellites at Andromeda by its size, orientation, rotation and mass. A hint of a disk of satellites was formed at the Milky Way too. In addition, the encounter induced a warp in the disk of the simulated Milky Way that resembles the real warp by its magnitude and orientation. We present here, for the first time, the proper motions of the members of the disk of satellites of Andromeda predicted by our simulation. Finally, we point out some of the remaining open questions which this hypothesis for the formation of disks of satellites brings up.

Precise measurement of the fundamental parameters of stellar systems, including mass and radius, depends critically on how well the stellar distances are known. Astrometry from space provides parallax measurements of unprecented accuracy, from which distances can be derived, initially from the Hipparcos mission, with further refinement of that analysis provided by van Leeuwen in 2007. The publication of the Gaia DR2 catalogue promises a dramatic improvement in the available data. We have recalculated the dynamical masses of a sample of 1 700 close visual binary stars using Gaia DR2 and compared the results with masses derived from both the original and enhanced Hipparcos data. We show the van Leeuwen analysis yields results close to those of Gaia DR2, but the latter are significantly more accurate. We consider the impact of the Gaia DR2 parallaxes on our understanding of the sample of visual binaries.

Modified Newtonian dynamics by Milgrom is a paradigm for explaining the rotation curves of spiral galaxies and various other large scale structures. This paradigm includes several different theories. Here we present Milgrom's modified inertia (MI) theory in terms of a simple and tractable non-conservative Newtonian dynamics, which is useful in obtaining observable predictions of MI. It is found that: 1) Modified inertia theory is equivalent to a Newtonian theory, with a non-conservative gravitational field, and dark matter density; 2) The tidal force in the equivalent Newtonian dynamics is non-conservative, and its effect on a binary system in free fall in the gravitational field of a spheroid is addressed. We also discuss attempts to restore conservation in MI.

In a collisionless plasma, the energy distribution function of plasma particles can be strongly affected by turbulence, in particular, it can develop a non-thermal power-law tail at large energies. We argue that turbulence with initially relativistically strong magnetic perturbations (magnetization parameter $\sigma \gg 1$) quickly evolves into a state with ultra-relativistic plasma temperature but mildly relativistic turbulent fluctuations. We present a phenomenological and numerical study suggesting that in this case, the exponent $\alpha$ in the power-law particle energy distribution function, $f(\gamma)d\gamma\propto \gamma^{-\alpha}d\gamma$, depends on magnetic compressibility of turbulence. Our analytic prediction for the scaling exponent $\alpha$ is in good agreement with the numerical results.

Gravitational waves emitted by distorted black holes, such as those arising from the coalescence of binary black holes, or a falling compact star into a supermassive black hole, carry not only information about the corresponding spacetime but also the information about the environment surrounding the black holes. In this paper, we study the effects of the dark matter halos with three different density profiles on the gravitational axial perturbations of a Schwarzschild-like black hole. For this purpose, we first consider modified Schwarzschild black holes with three different dark matter profiles and derive the equation of motion of the axial perturbations of the modified Schwarzschild metric. It is shown that by ignoring the dark matter perturbations, a Regge-Wheeler-like master equation with a modified potential for the axial perturbation can be obtained explicitly. Then we calculate the complex frequencies of the quasi-normal modes of the Schwarzschild-like black hole in the dark matter halos by applying the sixth-order WKB method. The corresponding gravitational wave spectra with the effects of the dark matter halos have also been discussed.

We consider neutrino evolution master equations in dense moving and polarized matter consisted of electrons, neutrons, protons and neutrinos. We also take into account the neutrino magnetic moment interaction with a magnetic field. We point out the mechanisms responsible for the neutrino spin precession and provide the expressions for the corresponding interaction Hamiltonians that should be taken into account in theoretical treatments of collective neutrino oscillations.

We show how a modified Friedmann equation, originating from a model of the universe built from a certain $W_3$ algebra, has the potential to explain the difference between the Hubble constants extracted from CMB data and from supernova data.

We study the efficiency of radiation in driving jets around black holes. Including general relativity for the radiation driving, we also show that the radiation field is affected by strong gravitational field in non linear manner, making general relativistic analysis more significant. We obtained internal shocks in jets close to the base, as a result of non-conical cross section and nature of radiation field on jet dynamics. Theoretical evidence of internal shocks is significant, as these are required to explain high energy tail of the spectra of radio sources. Under thermal and radiative driving, jets with electron-positron composition are obtained to be achieving relativistic speeds up to Lorentz factors $\gamma \sim 10$ while for electron-proton composition it is $\gamma \sim 2$ for luminous discs. We also showed that extragalactic jets around AGNs are faster than those around microquasars. We obtain that through Compton scattering, relativistic jets generated with non relativistic temperatures, are efficiently heated by temperatures up to $T \sim 10^{11-12} \rm K$. Further, we obtained transonic solutions with relativistic terminal speeds, for bound state of the jet at the base (that is, generalized Bernoulli parameter $E < 1$), where gravity is dominant over thermal driving. In Thomson scattering regime no transonic solutions could be obtained. This shows that radiation, in fact, has greater role as it not only accelerates the unbound matter of the jets, but the acceleration is effective even in absence of any other accelerating agent such as thermal acceleration. A detailed analysis of dependence of jet variables upon various parameters like plasma composition, magnetic pressure in the disc, luminosity, accretion rate, jet geometry and jet launching energy etc, is carried out.

We show that, if the principle of causality is violated, a system that is dissipative in a reference frame may look anti-dissipative in another one. By contrast, if causality holds, perturbations which are damped in one reference frame, are damped in all reference frames, whereas systems that are unstable in one reference frame, are unstable in all reference frames. These statements are mathematical theorems, valid in the non-linear regime, for all deterministic field theories, independently from the dynamical details. The present results solve and clarify the conceptual difficulties related with modelling dissipation in special and general relativity, showing that sub-luminality sets the ultimate boundary between exploitable and non-exploitable theories.

The Einstein and Maxwell equations are both systems of hyperbolic equations which need to satisfy a set of elliptic constraints throughout evolution. However, while electrodynamics (EM) and magnetohydrodynamics (MHD) have benefited from a large number of evolution schemes that are able to enforce these constraints and are easily applicable to curvilinear coordinates, unstructured meshes, or N-body simulations, many of these techniques cannot be straightforwardly applied to existing formulations of the Einstein equations. We develop a 3+1 a formulation of the Einstein equations which shows a striking resemblance to the equations of relativistic MHD and to EM in material media. The fundamental variables of this formulation are the frame fields, their exterior derivatives, and the Nester-Witten and Sparling forms. These mirror the roles of the electromagnetic 4-potential, the electromagnetic field strengths, the field excitations and the electric current. The role of the lapse function and shift vector, corresponds exactly to that of the scalar electric potential. The formulation, that we name dGREM (for differential forms, General Relativity and Electro-Magnetism), is manifestly first order and flux-conservative, which makes it suitable for high-resolution shock capturing schemes and finite-element methods. Being derived as a system of equations in exterior derivatives, it is directly applicable to any coordinate system and to unstructured meshes, and leads to a natural discretization potentially suitable for the use of machine-precision constraint propagation techniques such as the Yee algorithm and constrained transport. Due to these properties, we expect this new formulation to be beneficial in simulations of many astrophysical systems, such as binary compact objects and core-collapse supernovae as well as cosmological simulations of the early universe.

Scientific analysis for the gravitational-wave detector LISA will require theoretical waveforms from extreme-mass-ratio inspirals (EMRIs) that extensively cover all possible orbital and spin configurations around astrophysical Kerr black holes. However, on-the-fly calculations of these waveforms have not yet overcome the high dimensionality of the parameter space. To confront this challenge, we present a user-ready EMRI waveform model for generic (eccentric and inclined) orbits in Kerr spacetime, using an analytical self-force approach. Our model accurately covers all EMRIs with arbitrary inclination and black hole spin, up to modest eccentricity ($\lesssim 0.3$) and separation ($\gtrsim2$--$10M$ from the last stable orbit). In that regime, our waveforms are accurate at the leading `adiabatic' order, and they approximately capture transient self-force resonances that significantly impact the gravitational-wave phase. The model fills an urgent need for extensive waveforms in ongoing data-analysis studies, and its individual components will continue to be useful in future science-adequate waveforms.