Papers for Thursday, Sep 29 2022

Simulating time-domain observations of gravitational wave (GW) detector environments will allow for a better understanding of GW sources, augment datasets for GW signal detection and help in characterizing the noise of the detectors, leading to better physics. This paper presents a novel approach to simulating fixed-length time-domain signals using a three-player Wasserstein Generative Adversarial Network (WGAN), called DVGAN, that includes an auxiliary discriminator that discriminates on the derivatives of input signals. An ablation study is used to compare the effects of including adversarial feedback from an auxiliary derivative discriminator with a vanilla two-player WGAN. We show that discriminating on derivatives can stabilize the learning of GAN components on 1D continuous signals during their training phase. This results in smoother generated signals that are less distinguishable from real samples and better capture the distributions of the training data. DVGAN is also used to simulate real transient noise events captured in the advanced LIGO GW detector.

Simulations show that the orbits of planets are readily disrupted in dense star-forming regions; planets can be exchanged between stars, become free-floating and then be captured by other stars. However, dense star-forming regions also tend to be populous, containing massive stars that emit photoionising radiation, which can evaporate the gas in protoplanetary discs. We analyse N-body simulations of star-forming regions containing Jovian-mass planets and determine the times when their orbits are altered, when they become free-floating, and when they are stolen or captured. Simultaneously, we perform calculations of the evolution of protoplanetary discs when exposed to FUV radiation fields from massive stars in the same star-forming regions. In almost half (44 per cent) of the planetary systems that are disrupted - either altered, captured, stolen or become free-floating, we find that the radius of the protoplanetary disc evolves inwards, or the gas in the disc is completely evaporated, before the planets' orbits are disrupted. This implies that planets that are disrupted in dense, populous star-forming regions are more likely to be super Earths or mini Neptunes, as Jovian mass planets would not be able to form due to mass loss from photoevaporation. Furthermore, the recent discoveries of distant Jovian mass planets around tightly-packed terrestrial planets argue against their formation in populous star-forming regions, as photoevaporation would preclude gas giant planet formation at distances of more than a few au.

Galaxy clusters grow by accreting galaxies as individual objects, or as members of a galaxy group. These groups can strongly impact galaxy evolution, stripping the gas from galaxies, and enhancing the rate of galaxy mergers. However, it is not clear how the dynamics and structure of groups are affected when they interact with a large cluster, or whether all group members necessarily experience the same evolutionary processes. Using data from TheThreeHundred project, a suite of 324 hydrodynamical resimulations of large galaxy clusters, we study the properties of 1340 groups passing through a cluster. We find that half of group galaxies become gravitationally unbound from the group by the first pericentre, typically just 0.5-1 Gyr after cluster entry. Most groups quickly mix with the cluster satellite population; only 8% of infalling group haloes later leave the cluster, although for nearly half of these, all of their galaxies have become unbound, tidally disrupted or merged into the central by this stage. The position of galaxies in group-centric phase space is also important -- only galaxies near the centre of a group ($r\lesssim0.7R_{200}$) remain bound once a group is inside a cluster, and slow-moving galaxies in the group centre are likely to be tidally disrupted, or merge with another galaxy. This work will help future observational studies to constrain the environmental histories of group galaxies. For instance, groups observed inside or nearby to clusters have likely approached very recently, meaning that their galaxies will not have experienced a cluster environment before.

Gravitational-wave (GW) detections are starting to reveal features in the mass distribution of double compact objects. The lower end of the black hole (BH) mass distribution is especially interesting as few formation channels contribute here and because it is more robust against variations in the cosmic star formation than the high mass end. In this work we explore the stable mass transfer channel for the formation of GW sources with a focus on the low-mass end of the mass distribution. We conduct an extensive exploration of the uncertain physical processes that impact this channel. We note that, for fiducial assumptions, this channel reproduces the peak of the GW-observed binary BH merger rate remarkably well and predicts a cutoff mass that coincides with the upper edge of the purported neutron star BH mass gap. The peak and cutoff mass are a consequence of unique properties of this channel, namely (1) the requirement of stability during the mass transfer phases, and (2) the complex way in which the final compact object masses scale with the initial mass. We provide an analytical expression for the cutoff in the primary component mass and show that this adequately matches our numerical results. Our results imply that selection effects resulting from the formation channel alone can provide an explanation for the purported NS-BH mass gap in GW detections. This provides an alternative to the commonly adopted view that the gap is results from supernova fallback during BH formation.

We present a nonlinear study of the inflationary epoch based on numerical lattice simulations. Lattice simulations are a well-known tool in primordial cosmology, and they have been extensively used to study the reheating epoch after inflation. We generalize this known machinery to the inflationary epoch. Being this the first simulation of the inflationary epoch much before the end of inflation, the first part of the thesis focuses on the minimal single-field model of inflation. We discuss the conceptual and technical ingredients needed to simulate inflation on a lattice. The simulation is used to reproduce the nearly scale-invariant spectrum of scalar perturbations, as well as the oscillations in the power spectrum caused by a step in the potential. In the second part, we focus on the more complicated axion-U(1) model of inflation and present the first lattice simulation of this model during the deep inflationary epoch. We use the simulation to discover new properties of primordial scalar perturbations from this model. In the linear regime of the theory, we find high-order non-Gaussianity (beyond trispectrum) to be key to describing the statistical properties of scalar perturbations. Conversely, we find perturbations to be nearly Gaussian in the nonlinear regime of the theory. This relaxes existing constraints from the overproduction of primordial black holes, allowing for a gravitational waves signal in the observable range of upcoming experiments such as LISA. Our results show that lattice simulations can be a powerful tool to study the inflationary epoch and its observational signatures.

A fundamental assumption adopted in nearly every extragalactic emission-line study is that the attenuation of different emission lines can be described by a single attenuation curve. Here we show this assumption fails in many cases with important implications for derived results. We developed a new method to measure the differential nebular attenuation among three kinds of transitions: the Balmer lines of hydrogen, high-ionization transitions, and low-ionization transitions. This method bins the observed data in a multidimensional space spanned by attenuation-insensitive line ratios. Within each small bin, the variations in line ratios are mainly driven by the variations in the nebular attenuation. This allows us to measure the nebular attenuation using both forbidden lines and Balmer lines. We applied this method to a sample of 2.4 million H II region spaxels from SDSS-IV MaNGA. We found that the attenuation of high ionization lines and Balmer lines can be well described by a single Fitzpatrick (1999) extinction curve with $R_V=3.1$. However, no single attenuation curve can simultaneously account for all three transitions. This strongly suggests that different lines have different effective attenuations, likely because spectroscopy at kilo-parsecs resolution mixes multiple regions with different intrinsic line ratios and different levels of attenuation. As a result, the assumption that different lines follow the same attenuation curve breaks down. Using a single attenuation curve determined by Balmer lines to correct attenuation-sensitive forbidden line ratios could bias the nebular parameters derived by 0.06--0.25 dex at $A_V = 1$, depending on the details of the dust attenuation model. Observations of a statistically large sample of H II regions with high spatial resolutions and large spectral coverage are vital for improved modeling and deriving accurate corrections for this effect.

Recent observations suggest galaxies may ubiquitously host a molecular component to their multiphase circumgalactic medium (CGM). However, the structure and kinematics of the molecular CGM remains understudied theoretically and largely unconstrained observationally. Recent work suggests molecular gas clouds with efficient cooling survive acceleration in hot winds similar to atomic clouds. Yet the pressure-driven fragmentation of molecular clouds when subjected to external shocks or undergoing cooling remains unstudied. We perform radiative, inviscid hydrodynamics simulations of clouds perturbed out of pressure equilibrium to explore the process of hydrodynamic fragmentation to molecular temperatures. We find molecular clouds larger than a critical size can shatter into a mist of tiny droplets, with the critical size deviating significantly from the atomic case. We find that cold clouds shatter only if the sound crossing time exceeds the local maximum of the cooling time ~8000 K. Moreover, we find evidence for a universal mechanism to 'shatter' cold clouds into a 'mist' of tiny droplets as a result of rotational fragmentation -- a process we dub 'splintering.' Our results have implications for resolving the molecular phase of the CGM in observations and cosmological simulations.

We present a measurement of the intrinsic space density of intermediate redshift ($z\sim0.5$), massive ($M_{*} \sim 10^{11} \ \text{M}_{\odot}$), compact ($R_{e} \sim 100$ pc) starburst ($\Sigma_{SFR} \sim 1000 \ \text{M}_{\odot} \ \text{yr}^{-1} \text{kpc}^{-1}$) galaxies with tidal features indicative of them having undergone recent major mergers. A subset of them host kiloparsec scale, $>1000 \ \text{km}\ \text{s}^{-1}$ outflows and have little indication of AGN activity, suggesting that extreme star formation can be a primary driver of large-scale feedback. The aim for this paper is to calculate their space density so we can place them in a better cosmological context. We do this by empirically modeling the stellar populations of massive, compact starburst galaxies. We determine the average timescale for which galaxies that have recently undergone an extreme nuclear starburst would be targeted and included in our spectroscopically selected sample. We find that massive, compact starburst galaxies targeted by our criteria would be selectable for $\sim 148 ^{+27}_{-24}$ Myr and have an intrinsic space density $n_{\text{CS}} \sim (1.1^{+0.5}_{-0.3}) \times 10^{-6} \ \ \text{Mpc}^{-3}$. This space density is broadly consistent with our $z\sim0.5$ compact starbursts being the most extremely compact and star forming low redshift analogs of the compact star forming galaxies in the early Universe as well as them being the progenitors to a fraction of intermediate redshift post starburst and compact quiescent galaxies.

The recent observation of the anomalous X-ray pulsar 4U 0142+61 by the Imaging X-ray Polarimetry Explorer (IXPE) opened up a new avenue to study magnetars, neutron stars endowed with superstrong magnetic fields ($B\gtrsim 10^{14}$ G). The detected polarized X-rays exhibit an intriguing 90$^\circ$ linear polarization swing from low photon energies ($E\lesssim 4$ keV) to high photon energies ($E\gtrsim 5.5$ keV). We show that this swing can be naturally explained by photon polarization mode conversion at the vacuum resonance in the magnetar atmosphere; the resonance arises from the combined effects of plasma-induced birefringence and QED-induced vacuum birefringence in strong magnetic fields, and the mode conversion is analogous to the MSW effect for neutrino oscillations. This explanation suggests that the atmosphere of 4U 0142 be composed of partially ionized heavy elements (such as Fe), and the surface magnetic field be comparable or less than $10^{14}$~G, consistent with the dipole field inferred from the measured spindown. It also implies that the spin axis of 4U 0142+61 is aligned with its velocity direction.

A recently-discovered class of outflows, extremely high-velocity outflows (EHVOs), may be key to understanding feedback processes as it is likely the most powerful in terms of mass-energy. These EHVOs have been observed at redshifts 1.052 < z_em < 7.641, but the potential connection with outflows in emission had not been studied. We find that EHVOs, albeit their small numbers at the moment, appear to show distinct CIV and HeII properties. In particular, EHVOs are more predominant in quasars with large blueshifts of the CIV emission line, suggesting a connection between emission and absorption outflowing signatures for these extreme outflows. We also find incipient trends with the maximum velocity of the outflows, which is similar to what has been previously found in BALQSOs, but now extending previous studies to speeds up to ~0.2c. We find the bolometric luminosities, Eddington ratios, and black hole masses of our sample are overall very similar from the general quasar population upon considering their CIV emission properties. This is close to the case for HeII EW as we observe a tentative upper limit to the HeII strength for a quasar to host an EHVO. This study shows that extreme outflows such as EHVOs appear in quasars that are clearly a distinct class from the overall BALQSO population, and solidify the relation between outflows observed in emission and in absorption.

Dwarf galaxies are ubiquitous throughout the universe and are extremely sensitive to various forms of internal and external feedback. Over the last two decades, the census of dwarf galaxies in the Local Group and beyond has increased markedly. While hydrodynamic simulations (e.g. FIRE II, MINT Justice League) have reproduced the observed dwarf properties down to the ultra-faints, such simulations require extensive computational resources to run. In this work, we constrain the standard physical implementations in the semi-analytic model Galacticus to reproduce the observed properties of the Milky Way satellites down to the ultra-faint dwarfs found in SDSS. We run Galacticus on merger trees from our high-resolution N-body simulation of a Milky Way analog. We determine the best-fit parameters by matching the cumulative luminosity function and luminosity-metallicity relation from both observations and hydrodynamic simulations. With the correct parameters, the standard physics in Galacticus can reproduce the observed luminosity function and luminosity-metallicity relation of the Milky Way dwarfs. In addition, we find a multi-dimensional match with half-light radii, velocity dispersions and mass-to-light ratios at z = 0 down to M_V <= -6 (L >= 10^4 L_solar). In addition to successfully reproducing the properties of the z = 0 Milky Way satellite population, our modeled dwarfs have star formation histories which are consistent with those of the Local Group dwarfs.

In this paper we analyze Sloan g,r,i archival imaging data of M33 taken by Hartman et al. (2006) using Megacam at the Canada-France-Hawaii Telescope. To determine the distance to the M33 galaxy, we performed several analytical steps to identify its Cepheid population. We used the Lomb-Scargle algorithm to find periodicity and visually identified 1989 periodic variable stars. Since Cepheids occupy a specific region of the color-magnitude diagram, to differentiate Cepheids from other variables we used the expected position of the Cepheid instability strip to down-select Cepheids in M33 from other variables. This led to our sample of 1622 variables, the largest Cepheid sample known in M33 to date. We further classified these Cepheids into different subclasses, and used the fundamental mode Cepheids to estimate distance moduli for M33 in different filters: {\mu}= 25.044+/-0.083 mag in the g filter, {\mu}= 24.886+/-0.074 mag in the r filter, and {\mu}= 24.785+/-0.068 mag for the i filter. These results are in agreement with previous results.

The young and well-studied planetary nebula NGC 7027 harbors significant molecular gas that is irradiated by luminous, point-like UV (central star) and diffuse (shocked nebular) X-ray emission. This nebula represents an excellent subject to investigate the molecular chemistry and physical conditions within photon- and X-ray-dominated regions (PDRs and XDRs). As yet, the exact formation routes of CO$^+$ and HCO$^+$ in PN environments remain uncertain. Here, we present $\sim$2$"$ resolution maps of NGC 7027 in the irradiation tracers CO$^+$ and HCO$^+$, obtained with the IRAM NOEMA interferometer, along with SMA CO and HST 2.12~$\mu$m H$_2$ data for context. The CO$^+$ map constitutes the first interferometric map of this molecular ion in any PN. Comparison of CO$^+$ and HCO$^+$ maps reveal strikingly different emission morphologies, as well as a systematic spatial displacement between the two molecules; the regions of brightest HCO$^+$, found along the central waist of the nebula, are radially offset by $\sim$1$"$ ($\sim$900 au) outside the corresponding CO$^+$ emission peaks. The CO$^+$ emission furthermore precisely traces the inner boundaries of the nebula's PDR (as delineated by near-IR H$_2$ emission), suggesting that central star UV emission drives CO$^+$ formation. The displacement of HCO$^+$ radially outward with respect to CO$^+$ is indicative that dust-penetrating soft X-rays are responsible for enhancing the HCO$^+$ abundance in the surrounding molecular envelope, forming an XDR. These interferometric CO$^+$ and HCO$^+$ observations of NGC 7027 thus clearly establish the spatial distinction between the PDR and XDR formed (respectively) by intense UV and X-ray irradiation of molecular gas.

We present 3.6 and 4.5 um photometry for a set of 61 standard stars observed by Spitzer's Infrared Spectrograph (IRS). The photometry was obtained with the Infrared Array Camera (IRAC) on Spitzer in order to help tie the calibration of IRAC and the IRS, which had been anchored to the calibration of the Multiband Infrared Photometer for Spitzer (MIPS). The wavelength range of the IRS data only slightly overlaps with the IRAC 4.5 um band and not at all with the 3.6 um band. Therefore, we generated synthetic spectra from spectral templates of stars with the same spectral types and luminosity classes as our sample stars, normalized to the IRS data at 6-7 um, and compared those to the observed photometry. The new IRAC observations of IRS standard stars demonstrate that the two instruments are calibrated to within 1% of each other.

To understand the formation of stars from clouds of molecular gas, one essentially needs to know two things: What gas collapses, and how long it takes to do so. We address these questions by embedding pseudo-Lagrangian tracer particles in three simulations of self-gravitating turbulence. We identify prestellar cores at the end of the collapse, and use the tracer particles to rewind the simulations to identify the preimage gas for each core at the beginning of each simulation. This is the first of a series of papers, wherein we present the technique and examine the first question: What gas collapses? For the preimage gas at the t=0, we examine a number of quantities; the probability distribution function (PDF) for several quantities, the structure function for velocity, several length scales, the volume filling fraction, the overlap between different preimages, and fractal dimension of the preimage gas. Analytic descriptions are found for the PDFs of density and velocity for the preimage gas. We find that the preimage of a core is large and sparse, and we show that gas for one core comes from many turbulent density fluctuations and a few velocity fluctuations. We find that binary systems have preimages that overlap in a fractal manner. Finally, we use the density distribution to derive a novel prediction of the star formation rate.

We present polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3, G43.3-0.2, G41.1-0.3, and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. We detect polarization from G46.8-0.3, G43.3-0.2 and G39.2-0.3 but find upper limits at the 1% level of Stokes I for G41.1-0.3. For G46.8-0.3 and G39.2-0.3 the fractional polarization varies on small scales from 1% to ~6%. G43.3-0.2 is less polarized with fractional polarization <~3%. We find upper limits at the 1% level for the brighter regions in each SNR with no evidence for associated enhanced Faraday depolarization. We observe significant variation in Faraday depth and fractional polarization on angular scales down to the resolution limit of 16". Approximately 6% of our polarization detections from G46.8-0.3 and G39.2-0.3 exhibit two-component Faraday rotation and 14% of polarization detections in G43.3-0.2 are multi-component. For G39.2-0.3 we find a bimodal Faraday depth distribution with a narrow peak and a broad peak for all polarization detections as well as for the subset with two-component Faraday rotation. We identify the narrow peak with the front side of the SNR and the broad peak with the back side. Similarly, we interpret the observed Faraday depth distribution of G46.8-0.3 as a superposition of the distributions from the front side and the back side. We interpret our results as evidence for a partially filled shell with small-scale magnetic field structure and internal Faraday rotation.

We present GTC MEGARA high-dispersion integral field spectroscopic observations of the nova remnant QU\,Vul, which provide a comprehensive 3D view of this nova shell. The tomographic analysis of the H$\alpha$ emission reveals a complex physical structure characterized by an inhomogeneous and clumpy distribution of the material within this shell. The overall structure can be described as a prolate ellipsoid with an axial ratio of 1.4$\pm$0.2, a major axis inclination with the line of sight of $12^{\circ}\pm6^{\circ}$, and polar and equatorial expansion velocities $\approx$560 km~s$^{-1}$ and 400$\pm$60 km s$^{-1}$, respectively. The comparison of the expansion velocity on the plane of the sky with the angular expansion implies a distance of 1.43$\pm$0.23 kpc. The ionized mass is found to be $\approx 2\times 10^{-4}$ M$_\odot$, noting that the information on the 3D distribution of material within the nova shell has allowed us to reduce the uncertainty on its filling factor. The nova shell is still in its free expansion phase, which can be expected as the ejecta mass is much larger than the swept-up circumstellar medium mass. The 3D distribution and radial velocity of material within the nova shell provide an interpretation of the so-called "castellated" line profiles observed in early optical spectra of nova shells, which can be attributed to knots and clumps moving radially along different directions.

We here study the level of albedo variegation on the nucleus of Comet 67P/Churyumov-Gerasimenko. This is done by fitting the parameters of a standard photometric phase function model to disk-average radiance factor data in images acquired by the Rosetta/OSIRIS Narrow Angle Camera in the orange filter. Local discrepancies between the observed radiance factor and the disk-average solution are interpreted as a proxy W of the local single-scattering albedo. We find a wide range 0.02<W<0.09 around an average of W=0.055. The observed albedo variegation is strongly correlated with nucleus morphology - smooth terrain is brighter, and consolidated terrain is darker, than average. Furthermore, we find that smooth terrain darken prior to morphological changes, and that stratigraphically low terrain (with respect to the centre of each nucleus lobe) is brighter than stratigraphically high terrain. We propose that the observed albedo variegation is due to differences in porosity and the coherent effect: compaction causes small brighter particles to act collectively as larger optically effective particles, that are darker. Accordingly, we consider the dark consolidated terrain materials more compacted than smooth terrain materials, and darkening of the latter is due to subsidence.

This paper describes an adaptation of the Optimal Localized Averaging (OLA) inversion technique, originally developed for geo- and helioseismological applications, to the interpretation of solar spectroscopic data. It focuses on inverting the thermodynamical properties of the solar atmosphere assuming that the atmosphere and radiation field are in Local Thermodynamic Equilibrium (LTE). We leave inversions for magnetic field and non-LTE inversions for future work. The advantage with the OLA method is that it computes solutions that are optimally resolved (in depth) with minimal cross-talk error between variables. Additionally, the method allows for direct assessment of the vertical resolution of the inverted solutions. The primary challenges faced when adapting the method to spectroscopic inversions originate with the possible large amplitude differences between the atmospheric model used to initiate the inversion and the underlying atmosphere it aims to recover, necessitating the development of an iterative scheme. Here we describe the iterative OLA method we have developed for both single and multi-variable inversions and demonstrate its performance on simulated data and synthesized spectra. We note that when carrying out multi-variable inversions, employing response function amplification factors can address the inherent spectral-sensitivity bias that makes it hard to invert for less spectrally-sensitive variables. The OLA method can, in most cases, reliably invert as well as or better than the frequently employed Stokes Inversion based on Response functions (SIR) scheme, however some difficulties remain. In particular, the method struggles to recover large-scale offsets in the atmospheric stratification. We propose future strategies to improve this aspect.

The JWST early release data show unexpected high stellar mass densities of massive galaxies at 7<z<11, a high star formation efficiency is probably needed to explain this. However, such a high star formation efficiency would greatly increase the number of ionizing photons, which would be in serious conflict with current cosmic microwave background (CMB) measurements of cosmic reionization history. To solve this problem, we explore the fuzzy dark matter (FDM), which is composed of ultra-light scalar particles, e.g. ultra-light axions, and calculate its halo mass function and stellar mass density for different axion masses. We find that the FDM model can effectively suppress the formation of small halos and galaxies, so that with higher star formation efficiency, both the JWST data at $z\sim8$ and the optical depth of CMB scattering can be simultaneously matched. We find that an axion mass $m_a\simeq 5\times10^{-23}$ and $10^{-21}\ \rm eV$ are needed to fit the Planck and WMAP results, respectively. We also find that the JWST data at $z\sim10$ are still too high to fit in this scenario. We note that the estimated mean redshift of the sample may have large uncertainty, that it can be as low as $z\sim9$ depending on adopted spectral energy distribution (SED) templates and photometric-redshift code. Besides, the warm dark matter with $\sim$keV mass can also be an alternative choice, since it should have similar effects on halo formation as the FDM.

We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. As a result of these observations % that the observations provided, the state-of-the-art theoretical understanding of the formation and evolution of young stellar objects is described. In particular, we emphasize the importance of the coupling, decoupling, and re-coupling between weakly ionized gas and the magnetic field on appropriate scales. This highlights the complex and intimate relationship between gravitational collapse and magnetic fields in young protostars.

Solar flares, especially the M- and X-class flares, are often associated with coronal mass ejections (CMEs). They are the most important sources of space weather effects, that can severely impact the near-Earth environment. Thus it is essential to forecast flares (especially the M-and X-class ones) to mitigate their destructive and hazardous consequences. Here, we introduce several statistical and Machine Learning approaches to the prediction of the AR's Flare Index (FI) that quantifies the flare productivity of an AR by taking into account the numbers of different class flares within a certain time interval. Specifically, our sample includes 563 ARs appeared on solar disk from May 2010 to Dec 2017. The 25 magnetic parameters, provided by the Space-weather HMI Active Region Patches (SHARP) from Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO), characterize coronal magnetic energy stored in ARs by proxy and are used as the predictors. We investigate the relationship between these SHARP parameters and the FI of ARs with a machine-learning algorithm (spline regression) and the resampling method (Synthetic Minority Over-Sampling Technique for Regression with Gaussian Noise, short by SMOGN). Based on the established relationship, we are able to predict the value of FIs for a given AR within the next 1-day period. Compared with other 4 popular machine learning algorithms, our methods improve the accuracy of FI prediction, especially for large FI. In addition, we sort the importance of SHARP parameters by Borda Count method calculated from the ranks that are rendered by 9 different machine learning methods.

The anomalous heating of the solar upper atmosphere is one of the eight key problems in modern astronomy. Moreover, the stratification of the solar atmosphere is an outstanding key-problem in solar physics. In this study, a hot butterfly-like pattern is found to run through the chromosphere to the corona lying right on top of the magnetic butterfly pattern of sunspots in the photosphere. We thus propose to introduce the term butterfly body to describe the butterfly diagram in the 3-dimensional atmosphere. Besides, we discuss the so-called polar brightening in different layers. It is found to be statistically in anti-phase with the solar cycle in the photosphere and the chromosphere, while in phase with the solar cycle in the corona. Accordingly, we describe the role and relationship of solar magnetic elements of different magnetic flux strengths to explain the statistical structuring of the solar atmosphere with the butterfly body over the solar cycle.

We develop a framework for using clustering-based redshift inference (cluster-$z$) to measure the evolving galaxy luminosity function (GLF) and galaxy stellar mass function (GSMF) using WISE W1 ($3.4\mu m$) mid-infrared photometry and positions. We use multiple reference sets from the Galaxy And Mass Assembly (GAMA) survey, Sloan Digital Sky Survey (SDSS) and Baryon Oscillation Spectroscopic Survey (BOSS). Combining the resulting cluster-$z$s allows us to enlarge the study area, and by accounting for the specific properties of each reference set, making best use of each reference set to produce the best overall result. Thus we are able to measure the GLF and GSMF over $\sim 7500\, \mathrm{deg}^2 $ of the Northern Galactic Cap (NGC) up to $z<0.6$. Our method can easily be adapted for new studies with fainter magnitudes, which pose difficulties for the derivation of photo-$z$s. The measurement of the GSMF is currently limited by the models for k-corrections and mass-to-light ratios, rather than more complicated effects tied to the evolution of the differential galaxy bias. With better statistics in future surveys this technique is a strong candidate for studies with new emerging data from, e.g. the Vera C. Rubin Observatory, the Euclid mission or the Nancy Grace Roman Space Telescope.

IGR J16320-4751 is a highly obscured HMXB source containing a very slow neutron star ($P_{spin}\sim1300$ sec) orbiting its supergiant companion star with a period of $\sim$9 days. It shows high column density ($N_{H}\sim2-5\times10^{23}$ $cm^{-2}$) in the spectrum, and a large variation in flux along the orbit despite not being an eclipsing source. We report on some peculiar timing and spectral features from archival XMM-Newton observation of this source including 8 observations taken during a single orbit. The pulsar shows large timing variability in terms of average count rate from different observations, flaring activity, sudden changes in count rate, cessation of pulsation, and variable pulse profile even from observations taken a few days apart. We note that IGR J16320-4751 is among a small number of sources for which this temporary cessation of pulsation in the light curve has been observed. A time-resolved spectral analysis around the segment of missing pulse shows that variable absorption is deriving such behavior in this source. Energy resolved pulse profiles in the 6.2-6.6 keV band which has a partial contribution from Fe K$_\alpha$ photons, show strong pulsation. However, a more systematic analysis reveals a flat pulse profile from the contribution of Fe K$_\alpha$ photons in this band implying a symmetric distribution for the material responsible for this emission. Soft excess emission below 3 keV is seen in 6 out of 11 spectra of XMM-Newton observations.

We present LAMOST~J041920.07+072545.4 (hereafter J0419), a close binary consisting of a bloated extremely low mass pre-white dwarf (pre-ELM WD) and a compact object with an orbital period of 0.607189~days. The large-amplitude ellipsoidal variations and the evident Balmer and He~I emission lines suggest a filled Roche lobe and ongoing mass transfer. No outburst events were detected in the 15 years of monitoring of J0419, indicating a very low mass transfer rate. The temperature of the pre-ELM, $T_\mathrm{eff} = 5793_{-133}^{+124}\,\rm K$, is cooler than the known ELMs, but hotter than most CV donors. Combining the mean density within the Roche lobe and the radius constrained from our SED fitting, we obtain the mass of the pre-ELM, $M_1 = 0.176\pm 0.014\,M_\odot$. The joint fitting of light and radial velocity curves yields an inclination angle of $i = 66.5_{-1.7}^{+1.4}$ degrees, corresponding to the compact object mass of $M_2 = 1.09\pm 0.05\,M_\odot$. The very bloated pre-ELM has a smaller surface gravity ($\log g = 3.9\pm 0.01$, $R_1 = 0.78 \pm 0.02\,R_\odot$) than the known ELMs or pre-ELMs. The temperature and the luminosity ($L_\mathrm{bol} = 0.62_{-0.10}^{+0.11}\,L_\odot$) of J0419 are close to the main sequence, which makes the selection of such systems through the HR diagram inefficient. Based on the evolutionary model, the relatively long period and small $\log g$ indicate that J0419 could be close to the "bifurcation period" in the orbit evolution, which makes J0419 to be a unique source to connect ELM/pre-ELM WD systems, wide binaries and cataclysmic variables.

Primordial black holes (PBHs) whose masses are in $\sim[10^{-15}M_\odot,10^{-11}M_{\odot}]$ have been extensively studied as a candidate of whole dark matter (DM). One of the probes to test such a PBH-DM scenario is scalar-induced stochastic gravitational waves (GWs) accompanied with the enhanced primordial fluctuations to form the PBHs with frequency peaked in the mHz band being targeted by the LISA mission. In order to utilize the stochastic GWs for checking the PBH-DM scenario, it needs to exactly relate the PBH abundance and the amplitude of the GWs spectrum. Recently in Kitajima et al., the impact of the non-Gaussianity of the enhanced primordial curvature perturbations on the PBH abundance has been investigated based on the peak theory, and they found that a specific non-Gaussian feature called the exponential tail significantly increases the PBH abundance compared with the Gaussian case. In this work, we investigate the spectrum of the induced stochastic GWs associated with PBH DM in the exponential-tail case. In order to take into account the non-Gaussianity properly, we employ the diagrammatic approach for the calculation of the spectrum. We find that the amplitude of the stochastic GW spectrum is slightly lower than the one for the Gaussian case, but it can still be detectable with the LISA sensitivity. We also find that the non-Gaussian contribution can appear on the high-frequency side through their complicated momentum configurations. Although this feature emerges under the LISA sensitivity, it might be possible to obtain information about the non-Gaussianity from GW observation with a deeper sensitivity such as the DECIGO mission.

Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mission. Potential flares are identified in two steps: first, we search the light curves for at least two subsequent data points exceeding a 5{\sigma} threshold above the running mean. For these two cadences, we subtract the "quiet" stellar flux from the Kepler pixel data to obtain new images where the potential flare is the main light source. In the second step, we use a Bayesian approach to fit the point spread function of the instrument to determine the most likely location of the flux excess on the detector. We applied our method to 5862 main-sequence stars with near-solar effective temperatures. We found 2274 events exceeding the 5-sigma in at least two consecutive points in the light curves. Applying the second step reduced this number to 342 superflares. Of these, 283 flares happened on 178 target stars, 47 events are associated with fainter background stars, and in 10 cases, the flare location cannot be distinguished between the target and a background star. We also present cases where flares have been reported previously but our technique could not attribute them to the target star. We conclude that 1) identifying outliers in the light curves alone is insufficient to attribute them to stellar flares and 2) flares can only be uniquely attributed to a certain star when the instrument pixel-level data together with the point spread function are taken into account. As a consequence, previous flare statistics are likely contaminated by instrumental effects and unresolved astrophysical sources.

We derive the local stellar formation history from the Gaia-defined 40 pc white dwarf sample. This is currently the largest volume-complete sample of white dwarfs for which spectroscopy is available, allowing for classification of the chemical abundances at the photosphere, and subsequently accurate determination of the atmospheric parameters. We create a population synthesis model and show that a uniform stellar formation history for the last ~10.5 Gyr provides a satisfactory fit to the observed distribution of absolute Gaia G magnitudes. To test the robustness of our derivation, we vary various assumptions in the population synthesis model, including the initial mass function, initial-to-final mass relation, kinematic evolution, binary fraction and white dwarf cooling timescales. From these tests, we conclude that the assumptions in our model have an insignificant effect on the derived relative stellar formation rate as a function of look-back time. However, the onset of stellar formation (age of Galactic disc) is sensitive to a variety of input parameters including the white dwarf cooling models. Our derived stellar formation history gives a much better fit to the absolute Gaia G magnitudes than most previous studies.

We examine the particle acceleration mechanism in the Mpc-scale bridge between Abell 399 and Abell 401 and assess in particular if the synchrotron emission originates from first-order or second-order Fermi re-acceleration. We use deep (~40 hours) LOw-Frequency ARray (LOFAR) observations from Abell 399 and Abell 401 and apply improved direction-dependent calibration to produce deep radio images at three different resolutions at 144 MHz. With a point-to-point analysis we find in the bridge trends between the radio emission from our new maps and X-ray emission from an XMM Newton observation. By analyzing our observations and results, we argue that second-order Fermi re-acceleration is currently the most favoured process to explain the emission from the radio bridge, where past AGN activity may be responsible for the supply of fossil plasma needed for in-situ re-acceleration. The radio halos from Abell 401 and Abell 399 are also consistent with a second-order Fermi re-acceleration model.

We measure the expansion rate of the Universe at redshift $z=2.3$ from the anisotropy of Lyman-$\alpha$ (Ly$\alpha$) forest correlations measured by the Sloan Digital Sky Survey (SDSS). Our result is the most precise from large-scale structure at $z>1$. In flat $\Lambda$CDM we determine the matter density to be $\Omega_\mathrm{m}=0.36^{+0.03}_{-0.04}$ from Ly$\alpha$ alone, a factor of two tighter than baryon acoustic oscillation results from the same data. Using a nucleosynthesis prior, we measure the Hubble constant to be $H_0=63.2\pm2.5$ km/s/Mpc. In combination with other SDSS tracers, we find $H_0=67.2\pm0.9$ km/s/Mpc and measure the dark energy equation-of-state parameter to be $w=-0.90\pm0.12$. Our work opens a new avenue for constraining cosmology at high redshift.

According to the number of detected bursts, fast radio bursts (FRBs) can be classified into two categories, i.e., one-off FRBs and repeating ones. We make a statistical comparison of these two categories based on the first FRB catalog of the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project. Using the Anderson-Darling, Kolmogrov-Smirnov, and Energy statistic tests, we find significant statistical differences ($p$-value $<$ 0.001) of the burst properties between the one-off FRBs and the repeating ones. More specifically, after controlling for distance, we find that the peak luminosities of one-off FRBs are, on average, higher than the repeating ones; the pulse temporal widths of repeating FRBs are, on average, longer than the one-off ones. The differences indicate that these two categories could have distinct physical origins. Moreover, we discuss the sub-populations of FRBs and provide statistical evidence to support the existence of sub-populations in one-off FRBs and in repeating ones.

In this work we make use of available Integral Field Unit (IFU) spectroscopy and slit spectra of several nearby galaxies. The pre-existing empirical R and S calibrations for abundance determinations are constructed using a sample of HII regions with high quality slit spectra. In this paper, we test the applicability of those calibrations to the IFU spectra. We estimate the calibration-based abundances obtained using both the IFU and the slit spectroscopy for eight nearby galaxies. The median values of the slit and IFU spectra-based abundances in bins of 0.1 in fractional radius Rg (normalized to the optical radius) of a galaxy are determined and compared. We find that the IFU and the slit spectra-based abundances obtained through the R calibration are close to each other, the mean value of the differences of abundances is 0.005 dex and the scatter in the differences is 0.037 dex for 38 datapoints. The S calibration can produce systematically underestimated values of the IFU spectra-based abundances at high metallicities, the mean value of the differences is -0.059 dex for 21 datapoints, while at lower metallicities the mean value of the differences is -0.018 dex and the scatter is 0.045 dex for 36 data points. This evidences that the R calibration produces more consistent abundance estimations between the slit and the IFU spectra than the S calibration. We find that the same calibration can produce close estimations of the abundances using IFU spectra obtained with different spatial resolution and different spatial samplings. This is in line with the recent finding that the contribution of the diffuse ionized gas to the large aperture spectra of HII regions has a secondary effect.

Asteroseismic observations of internal stellar rotation have indicated a substantial lack of angular momentum transport in theoretical models of subgiant and red-giant stars. Accurate core and surface rotation rate measurements are therefore needed to constrain internal transport processes included in the models. We eliminate substantial systematic errors of asteroseismic surface rotation rates found in previous studies. We propose a new objective function for the Optimally Localized Averages method of rotational inversions for red-giant stars, which results in more accurate envelope rotation rate estimates obtained from the same data. We use synthetic observations from stellar models across a range of evolutionary stages and masses to demonstrate the improvement. We find that our new inversion technique allows us to obtain estimates of the surface rotation rate that are independent of the core rotation. For a star at the base of the red-giant branch, we reduce the systematic error from about 20% to a value close to 0, assuming constant envelope rotation. We also show the equivalence between this method and the method of linearised rotational splittings. Our new rotational inversion method substantially reduces the systematic errors of red-giant surface rotation rates. In combination with independent measures of the surface rotation rate, this will allow better constraints to be set on the internal rotation profile. This will be a very important probe to further constrain the internal angular momentum transport along the lower part of the red-giant branch.

The ESA global astrometry space mission Gaia has been monitoring the position of a billion stars since 2014. The analysis of such a massive dataset is challenging in terms of the data processing involved. In particular, the blind detection and characterization of single or multiple companions to stars (planets, brown dwarfs, or stars) using Gaia astrometry requires highly efficient algorithms. In this article, we present a set of analytical methods to detect and characterize companions in scanning space astrometric time series as well as via a combination of astrometric and radial velocity time series. We propose a general linear periodogram framework and we derive analytical formulas for the false alarm probability (FAP) of periodogram peaks. Once a significant peak has been identified, we provide analytical estimates of all the orbital elements of the companion based on the Fourier decomposition of the signal. The periodogram, FAP, and orbital elements estimates can be computed for the astrometric and radial velocity time series separately or in tandem. These methods are complementary with more accurate and more computationally intensive numerical algorithms (e.g., least-squares minimization, Markov chain Monte Carlo, genetic algorithms). In particular, our analytical approximations can be used as an initial condition to accelerate the convergence of numerical algorithms. Our formalism has been partially implemented in the Gaia exoplanet pipeline for the third Gaia data release. Since the Gaia astrometric time series are not yet publicly available, we illustrate our methods on the basis of Hipparcos data, together with on-ground CORALIE radial velocities, for three targets known to host a companion: HD 223636 (HIP 117622), HD 17289 (HIP 12726), and HD 3277 (HIP 2790).

We investigate the properties of the mixed-mode (RRd) RR Lyrae (RRL) variables in the Fornax dwarf spheroidal (dSph) galaxy by using $B$- and $V$-band time series collected over twenty-four years. We compare the properties of the RRds in Fornax with those in the Magellanic Clouds and in nearby dSphs, with special focus on Sculptor. We found that the ratio of RRds over the total number of RRLs decreases with metallicity. Typically, dSphs have very few RRds with 0.49$\ltsim P_0 \ltsim $0.53 days, but Fornax fills this period gap in the Petersen diagram (ratio between first overtone over fundamental period versus fundamental period). We also found that the distribution in the Petersen diagram of Fornax RRds is similar to SMC RRds, thus suggesting that their old stars have a similar metallicity distribution. We introduce the Period-Amplitude RatioS (PARS) diagram, a new pulsation diagnostics independent of distance and reddening. We found that LMC RRds in this plane are distributed along a short- and a long-period sequence that we identified as the metal-rich and the metal-poor component. These two groups are also clearly separated in the Petersen and Bailey (luminosity amplitude versus logarithmic period) diagrams. These circumstantial evidence indicates that the two groups have different evolutionary properties. All the pulsation diagnostics adopted in this investigation suggest that old stellar populations in Fornax and Sculptor dSphs underwent different chemical enrichment histories. Fornax RRds are similar to SMC RRds, while Sculptor RRds are more similar to the metal-rich component of the LMC RRds.

Searching for compact objects (black holes, neutron stars, or white dwarfs) in the Milky Way is essential for understanding the stellar evolution history, the physics of compact objects, and the structure of our Galaxy. Compact objects in binaries with a luminous stellar companion are perfect targets for optical observations. Candidate compact objects can be achieved by monitoring the radial velocities of the companion star. However, most of the spectroscopic telescopes usually obtain stellar spectra at a relatively low efficiency, which makes a sky survey for millions of stars practically impossible. The efficiency of a large-scale spectroscopic survey, the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST), presents a specific opportunity to search for compact object candidates, i.e., simply from the spectroscopic observations. Late-type K/M stars are the most abundant populations in our Galaxy. Owing to the relatively large Keplerian velocities in the close binaries with a K/M-dwarf companion, a hidden compact object could be discovered and followed-up more easily. In this study, compact object candidates with K/M-dwarf companions are investigated with the LAMOST low-resolution stellar spectra. Based on the LAMOST Data Release 5, we obtained a sample of $56$ binaries, each containing a K/M-dwarf with a large radial velocity variation $\Delta V_{\rm R} > 150~{\rm km~s}^{-1}$. Complemented with the photometric information from the Transiting Exoplanet Survey Satellite, we derived a sample of $35$ compact object candidates, among which, the orbital periods of $16$ sources were revealed by the light curves. Considering two sources as examples, we confirmed that a compact object existed in the two systems by fitting the radial velocity curve. This study demonstrates the principle and the power of searching for compact objects through LAMOST.

We aim to quantify the effect of chemistry on the infall velocity in the prestellar core L1544. Previous observational studies have found evidence for double-peaked line profiles for the rotational transitions of several molecules, which cannot be accounted for with presently available models for the physical structure of the source, without ad hoc up-scaling of the infall velocity. We ran one-dimensional hydrodynamical simulations of the collapse of a core with L1544-like properties (in terms of mass and outer radius) using a state-of-the-art chemical model with a very large chemical network combined with an extensive description of molecular line cooling, determined via radiative transfer simulations, with the aim of determining whether these expansions of the simulation setup (as compared to previous models) can lead to a higher infall velocity. After running a series of simulations where the simulation is sequentially simplified, we found that the infall velocity is almost independent of the size of the chemical network or the approach to line cooling. We conclude that chemical evolution does not have a large impact on the infall velocity, and that the higher infall velocities that are implied by observations may be the result of the core being more dynamically evolved than what is now thought, or alternatively the average density in the simulated core is too low. However, chemistry does have a large influence on the lifetime of the core, which varies by about a factor of two across the simulations and grows longer when the chemical network is simplified. Therefore, although the model is subject to several sources of uncertainties, the present results clearly indicate that the use of a small chemical network leads to an incorrect estimate of the core lifetime, which is naturally a critical parameter for the development of chemical complexity in the pre-collapse phase.

Ocean-vaporizing impacts of chemically reduced planetesimals onto the early Earth have been suggested to catalyse atmospheric production of reduced nitrogen compounds and trigger prebiotic synthesis despite an oxidized lithosphere. While geochemical evidence supports a dry, highly reduced late veneer on Earth, the composition of late-impacting debris around lower-mass stars is subject to variable volatile loss as a result of their hosts' extended pre-main sequence phase. We perform simulations of late-stage planet formation across the M-dwarf mass spectrum to derive upper limits on reducing bombardment epochs in Hadean analog environments. We contrast the Solar System scenario with varying initial volatile distributions due to extended primordial runaway greenhouse phases on protoplanets and desiccation of smaller planetesimals by internal radiogenic heating. We find a decreasing rate of late-accreting reducing impacts with decreasing stellar mass. Young planets around stars $\leq$ 0.4 $M_\odot$ experience no impacts of sufficient mass to generate prebiotically relevant concentrations of reduced atmospheric compounds once their stars have reached the main sequence. For M-dwarf planets to not exceed Earth-like concentrations of volatiles, both planetesimals and larger protoplanets must undergo extensive devolatilization processes and can typically emerge from long-lived magma ocean phases with sufficient atmophile content to outgas secondary atmospheres. Our results suggest that transiently reducing surface conditions on young rocky exoplanets are favoured around FGK- stellar types relative to M-dwarfs.

We study the global agreement between the most recent observations of the Cosmic Microwave Background temperature and polarization anisotropies angular power spectra released by the Atacama Cosmology Telescope and the Planck satellite in various cosmological models that differ by the inclusion of different combinations of additional parameters. By using the Suspiciousness statistic, we show that the global "CMB tension" between the two experiments, quantified at the Gaussian equivalent level of $\sim 2.5\,\sigma$ within the baseline $\Lambda$CDM, is reduced at the level of $1.8\sigma$ when the effective number of relativistic particles ($N_{\rm eff}$) is significantly less than the standard value, while it ranges between $2.3\,\sigma$ and $3.5\,\sigma$ in all the other extended models.

High sensitivity interstellar scintillation and polarization observations of PSR~B0656+14 made at three epochs over a year using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) show that the scattering is dominated by two different compact regions. We identify the one nearer to the pulsar with the shell of the Monogem Ring, thereby confirming the association. The other is probably associated with the Local Bubble. We find that the observed position angles of the pulsar spin axis and the space velocity are significantly different, with a separation of $19\fdg3\pm$0\fdg8, inconsistent with a previously published near-perfect alignment of $1\degr\pm 2\degr$. The two independent scattering regions are clearly defined in the secondary spectra which show two strong forward parabolic arcs. The arc curvatures imply that the scattering screens corresponding to the outer and inner arcs are located approximately 28~pc from PSR B0656+14 and 185~pc from the Earth, respectively. Comparison of the observed Doppler profiles with electromagnetic simulations shows that both scattering regions are mildly anisotropic. For the outer arc, we estimate the anisotropy $A_R$ to be approximately 1.3, with the scattering irregularities aligned parallel to the pulsar velocity. For the outer arc, we compare the observed delay profiles with delay profiles computed from a theoretical strong-scattering model. Our results suggest that the spatial spectrum of the scattering irregularities in the Monogem Ring is flatter than Kolmogorov, but further observations are required to confirm this.

The purely gravitational evidence supporting the need for dark matter (DM) particles is compelling and based on Galactic to cosmological scale observations. Thus far, the promising WIMP scenarios particles have eluded detection, motivating alternative models for DM. We consider the scenarios involving the super-heavy dark matter (SHDM) that potentially can be emitted by primordial black holes (PBH) and can decay or annihilate into ultra-high energy (UHE) neutrinos and photons. The observation of a population of photons with energies $E\ge 100$ EeV would imply the existence of completely new physical phenomena, or shed some light on DM models. Only the ultra-high energy cosmic ray observatories have the capabilities to detect such UHE decay products via the measurements of UHE photon induced extensive air showers. Using the upper bound on the flux of UHE cosmic rays beyond $10^{11.3}$ GeV implying $J(>10^{11.3}~{\rm{GeV}})< 3.6\times 10^{-5}$ km$^{-2}$sr$^{-1}$y$^{-1}$, at the $90\%$ C.L. reported by the Pierre Auger Observatory, we obtain global limits on the lifetime of the DM particles with masses $10^{15}\le M_{X} \le 10^{17}$ GeV. The constraints derived here are new and cover a region of the parameter space not yet explored. We compare our results with the projected constraints from future POEMMA and JEM-EUSO experiments, in order to quantify the improvement that will be obtained by these missions. Moreover, assuming that an epoch of early PBHs domination introduces a unique spectral break,$f_{\ast}$, in the gravitational wave (GW) spectrum, which frequency is related to the SHDM mass, we map potential probes and limits of the DM particles masses on the $f_{\ast}-M_{X}$ parameter space.

Extending the results of a previous paper \citep{DelPopolo2005}, by taking into account the role of dynamical friction, we recovered the luminosity-temperature relation (LTR). While by assuming self-similarity, a scaling law in which $L\propto T^2$ is obtained, observations show that the relation between luminosity and temperature is steeper, ${L \propto T^ {\simeq 3}}$. This difference can be explained in terms of energy input by non-gravitational processes, like pre-heating, supernovae feedback, and heating from AGN. In this paper, we studied the LTR by means of a modified version of the punctuated equilibria model \citep{Cavaliere1999}, taking into account in addition dynamical friction, thus extending the approach found in \citep{DelPopolo2005}. The result is a non-self-similar LTR with a bend at $\simeq 2$ keV, with a slope $2.76 \pm 0.18$ at larger energies and $3.4 \pm 0.18$ at energies smaller than 2 keV. This result is in agreement with the XXL survey \citep{Giles2016}. Moreover the steeper slopes at smaller energies is in agreement with some studies claiming a further steepening of the LTR at the low mass end. We also compared the results of our model with the 400d groups sample, finding that in groups the slope is slightly steeper than in clusters, %and namely $3.35 \pm 0.3$, in agreement with the \citep{Zou2016} study for the 400d groups sample, that gives a slope $3.29 \pm 0.33$.

The DAMA/LIBRA observation of an annual modulation in the detection rate compatible with that expected for dark matter particles from the galactic halo has accumulated evidence for more than twenty years. It is the only hint of a direct detection of the elusive dark matter, but it is in strong tension with the negative results of other very sensitive experiments, requiring ad-hoc scenarios to reconcile all the present experimental results. Testing the DAMA/LIBRA result using the same target material, NaI(Tl), removes the dependence on the particle and halo models and is the goal of the ANAIS-112 experiment, taking data at the Canfranc Underground Laboratory in Spain since August 2017 with 112.5 kg of NaI(Tl). At very low energies, the detection rate is dominated by non-bulk scintillation events and careful event selection is mandatory. This article summarizes the efforts devoted to better characterize and filter this contribution in ANAIS-112 data using a boosted decision tree (BDT), trained for this goal with high efficiency. We report on the selection of the training populations, the procedure to determine the optimal cut on the BDT parameter, the estimate of the efficiencies for the selection of bulk scintillation in the region of interest (ROI), and the evaluation of the performance of this analysis with respect to the previous filtering. The improvement achieved in background rejection in the ROI, but moreover, the increase in detection efficiency, push the ANAIS-112 sensitivity to test the DAMA/LIBRA annual modulation result beyond 3$\sigma$ with three-year exposure, being possible to reach 5$\sigma$ by extending the data taking for a few more years than the scheduled 5 years which were due in August 2022.

It has been shown some years ago that dark matter haloes outskirts are characterized by very steep density profiles in a very small radial range. This feature has been interpreted as a pile up of at a similar location of different particle orbits, namely splashback material at half an orbit after collapse. Adhikari et al. (2014), obtained the location of the splashback radius through a very simple model, namely calculating a dark matter shell trajectory in the secondary infall model while it crosses a growing, NFW profile shaped, dark matter halo. Since they imposed a halo profile instead of calculating it from the trajectories of the shells of dark matter, they were not able to find the dark matter profile around the splashback radius. In the present paper, we use an improved spherical infall model taking into shell crossing, and several physical effects like ordered, and random angular momentum, dynamical friction, adiabatic contraction, etc. This allow us to determine the density profile from the inner to outer region, and study the behavior of the outer density profile. We will compare the density profiles, and the logarithmic slope of the density profile with the results of Diemer \& Kravtsov (2014) simulations, finding a good agreement between the prediction of the model and the simulations.

We use spatially-resolved spectroscopy of a distant giant gravitational arc to test orientation effects on MgII absorption equivalent width (EW) and covering fraction (kappa) in the circumgalactic medium of a foreground star-forming galaxy (G1) at z~0.77. Forty-two spatially-binned arc positions uniformly sample impact parameters (D) to G1 between 10 and 30 kpc and azimuthal angles alpha between 30 and 90 degrees (minor axis). We find an EW-D anti-correlation, akin to that observed statistically in quasar absorber studies, and an apparent correlation of both EW and kappa with alpha, revealing a non-isotropic gas distribution. In line with our previous results on MgII kinematics suggesting the presence of outflows in G1, at minimum a simple 3-D static double-cone model (to represent the trace of bipolar outflows) is required to recreate the EW spatial distribution. The D and alpha values probed by the arc cannot confirm the presence of a disc, but the data highly disfavor a disc alone. Our results support the interpretation that the EW-alpha correlation observed statistically using other extant probes is partly shaped by bipolar metal-rich winds.

This review aims at proposing to the field an overview of the Cusp-core problem, including a discussion of its advocated solutions, assessing how each can satisfactorily provide a description of central densities. Whether the Cusp-core problem reflects our insufficient grasp on the nature of dark matter, of gravity, on the impact of baryonic interactions with dark matter at those scales, as included in semi-analytical models or fully numerical codes, the solutions to it can point either to the need for a paradigm change in cosmology, or to to our lack of success in ironing out the finer details of the $\Lambda$CDM paradigm.

Radial colour gradients within galaxies arise from gradients of stellar age, metallicity and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy datasets and can be used to constrain galaxy formation models. Here we measure colour gradients for low-redshift galaxies (z<0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ~93,000 galaxies with spectroscopic redshifts and ~574,000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5 Re to Re, which corresponds to the inner disk of typical late type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average $g-r$ and $r-z$ colour, $M_r$, $M_\star$, and sSFR. Trends of gradient strength with $M_r$ ($M_\star$) show an inflection around $M_r\sim-21$ ($\log_{10} \, M_\star/\mathrm{M_\odot}\sim10.5$). Below this mass, colour gradients become steeper with increasing $M_\star$, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star}\sim10^{8}\,\mathrm{M_\odot}$. We compare our results to age and metallicity gradients in two datasets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA datasets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

We have presented NuSTAR and Swift observations of the newly discovered Be/X-ray pulsar eRASSU J052914.9-662446. This is the first detailed study of the temporal and spectral properties of the pulsar using 2020 observations. A coherent pulsation of 1411.5$\pm$0.5 s was detected from the source. The pulse profile was found to resemble a simple single peaked feature which may be due to emission from the surface of the neutron star only. Pulse profiles are highly energy dependent. The variation of the pulse fraction of the pulse profiles are found to be non-monotonic with energy. The 0.5-20 keV Swift and NuSTAR simultaneous can be fitted well with power-law modified by high energy cutoff of $\sim$ 5.7 keV. The NuSTAR luminosity in the 0.5-79 keV energy range was $\sim$ 7.9x10 35 erg/s. The spectral flux in 3-79 keV shows modulation with the pulse phase.

We study the power spectrum of the comoving curvature perturbation $\cal R$ in the model that glues two linear potentials of different slopes, originally proposed by Starobinsky. We find that the enhanced power spectrum reaches its maximum at the wavenumber which is $\pi$ times the junction scale. The peak is $\sim2.61$ times larger than the ultraviolet plateau. We also show that its near-peak behavior can be well approximated by a constant-roll model, once we define the effective ultra-slow-roll $e$-folding number appropriately by considering the contribution from non-single-clock phase only. Such an abrupt transition to non-attractor phase can leave some interesting characteristic features in the energy spectrum of the scalar-induced gravitational waves, which are detectable in the space-borne interferometers if the primordial black holes generated at such a high peak are all the dark matter.

The major interactions are known to trigger star formation in galaxies and alter their colour. We study the major interactions in filaments and sheets using the SDSS data to understand the influence of large-scale environments on the galaxy interactions. We identify the galaxies in filaments and sheets using the local dimension and also find the major pairs residing in these environments. The star formation rate and colour of the interacting galaxies as a function of pair separation are separately analyzed in filaments and sheets. The analysis is repeated for three volume limited samples covering different magnitude ranges. The major pairs residing in the filaments show a significantly higher star formation rate (SFR) and bluer colour than those residing in the sheets up to the projected pair separation of $\sim 50$ kpc. We observe a complete reversal of this behaviour for both the SFR and colour of the galaxy pairs having a projected separation larger than 50 kpc. Some earlier studies report that the galaxy pairs align with the filament axis. Such alignment inside filaments indicates anisotropic accretion that may cause these differences. We do not observe these trends in the brighter galaxy samples. The pairs in filaments and sheets from the brighter galaxy samples trace relatively denser regions in these environments. The absence of these trends in the brighter samples may be explained by the dominant effect of the local density over the effects of the large-scale environment.

We propose that the interaction between the axion-like particles (ALPs) and photons can be a possible origin of the optical circular polarization (CP) in blazars. The non-detection of the optical CP at the level of $0.1\%$ can be used to place a constraint on the ALP-photon coupling $g_{a\gamma}\cdot B_\mathrm{T0}\lesssim7.9\times10^{-12}~\mathrm{G\cdot GeV}^{-1}$ for $m_{a}\lesssim 10^{-13}~\mathrm{eV}$, which depends on the magnetic field model of the blazar jet. This constraint would be stringent for the blazar models with a large magnetic field strength, such as hadronic radiation models. We also perform an analysis for the possible observations of the optical CP in two blazars, and find that they could be explained by the ALP-photon coupling of $\mathcal{O}(10^{-12})~\mathrm{GeV}^{-1}$. As an outlook, our analysis can be improved by further researches on the radiation models of blazars and high-precision simultaneous measurements of the optical linear polarization and CP.

The forest of Lyman-$\alpha$ absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable -- the 1D flux power spectrum (FPS) -- should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard $\Lambda$CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-accelerated code {\sc Cholla} to study the evolution of the Lyman-$\alpha$ forest under a wide range of physically-motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts $4.0\lesssim z\lesssim 5.2$ (Boera et al. 2019) yields a lower limit $m_{\rm WDM}>3.1$ keV (95 percent CL) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a best-fit thermal relic mass of $m_{\rm WDM}=4.5_{-1.4}^{+45}$ keV (95 percent CL). We find that the suppression of the FPS from free-streaming saturates at $k\gtrsim 0.1\,$s km$^{-1}$ because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) at $144$ MHz, is a daunting, time-consuming process, even with extensive manpower. Using a machine learning approach, we aim to automate the radio component association of large ($> 15$ arcsec) radio components. We turned the association problem into a classification problem and trained an adapted Fast region-based convolutional neural network to mimic the expert annotations from the first LoTSS data release. We implemented a rotation data augmentation to reduce overfitting and simplify the component association by removing unresolved radio sources that are likely unrelated to the large and bright radio components that we consider using predictions from an existing gradient boosting classifier. For large ($> 15$ arcsec) and bright ($> 10$ mJy) radio components in the LoTSS first data release, our model provides the same associations for $85.3\%\pm0.6$ of the cases as those derived when astronomers perform the association manually. When the association is done through public crowd-sourced efforts, a result similar to that of our model is attained. Our method is able to efficiently carry out manual radio-component association for huge radio surveys and can serve as a basis for either automated radio morphology classification or automated optical host identification. This opens up an avenue to study the completeness and reliability of samples of radio sources with extended, complex morphologies.

This paper continues a series reporting different aspects of the behaviour of disc galaxy simulations that support spiral instabilities. The focus in this paper is to demonstrate how linear spiral instabilities saturate and decay, and how the properties of the disc affect the limiting amplitude of the spirals. Once again, we employ idealized models that each possess a single instability that we follow until it has run its course. Remarkably, we find a tight correlation between the growth rate of the mode and its limiting amplitude, albeit from only six simulations. We show that non-linear orbit deflections near corotation cause the mode to saturate, and that the more time available in a slowly-growing mode creates the critical deflections at lower amplitude. We also find that scattering at the inner Lindblad resonance is insignificant until after the mode has saturated. Our objective in this series of papers, which we believe we have now achieved, has been to develop a convincing and well-documented account of the physical behaviour of the spiral patterns that have been observed in simulations by others, and by ourselves, for many decades. Understanding the simulations is an important step towards the greater objective, which is to find observational evidence from galaxies that could confront the identified mechanism.

Am stars are often components of short-period binary systems, where tidal interactions would result in low rotational velocities and help to develop the chemical peculiarities observed. However, the origin of single Am stars and Am stars that belong to wide binary systems is unclear. There is very recent evidence of an Am star hosting a hot-brown dwarf likely synchronized and other possible Am stars hosting hot-Jupiter planets. We wonder if these hot-low mass companions could play a role in the development of an Am star, that is to say, if they could help to mitigate the "single Am" problem. We studied a sample of 19 early-type stars, 7 of them hosting hot-brown dwarfs and 12 of them hosting hot-Jupiter planets. We detected 4 Am stars in our sample (KELT-19A, KELT-17, HATS-70 and TOI-503) and 2 possible Am stars (TOI-681 and HAT-P-69). In particular, we detected the new Am star HATS-70 which hosts a hot-brown dwarf, and rule out this class for the hot-Jupiter host WASP-189, both showing different composition than previously reported. We estimated the incidence of Am stars within stars hosting hot-brown dwarfs (50-75%) and within stars hosting hot-Jupiters (20-42%). The incidence of Am stars hosting hot-brown dwarfs resulted higher than the frequency of Am stars in general. This would imply that the presence of hot-brown dwarfs could play a role in the development of Am stars and possibly help to mitigate the "single Am" problem, different to the case of hot-Jupiter planets. Notably, these results would also indicate that the search for hot-brown dwarfs may be benefited by targeting single Am stars or Am stars in wide binary systems. We encourage the analysis off additional early-type stars hosting hot-companions in order to improve the significance of the initial trends found here. [abridged]

We have computed the deflection angles caused by 195 objects in the solar system, including 177 satellites and eight asteroids. Twenty-one satellites and six asteroids can bend light from distant compact extragalactic sources by more than 0.1 $\mu$as, and fourteen satellites and the asteroid Ceres can deflect light by more than 1.0 $\mu$as. We calculated the zones and durations of perturbations posed by the gravitational fields of five planets (excluding Earth, Jupiter, and Saturn), Pluto, and Ceres, where the perturbations would affect astrometry measured with the Squared Kilometre Array (SKA). Perturbed zones with deflection angles larger than 0.1 and 1.0 $\mu$as appear as ribbons. Their widths range from dozens of degrees for Uranus, Neptune, and Venus to several degrees or less for other objects at 0.1 $\mu$as, and from $\sim$ 16$^{\circ}$ for Venus to several degrees or less for other objects at 1.0 $\mu$as. From the calculated perturbation durations, the influence of the gravitational fields of selected objects can be divided into four levels: hardly affect SKA astrometry (I), may have little effect (II), may have a great effect (III) on single-epoch astrometry, and may greatly affect both single- and multi-epoch astrometry (IV). The objects corresponding to these levels are Ceres (I), Pluto (II), Mercury and Mars (III), and other objects (IV).

First order phase transitions in the early universe could produce a gravitational-wave background that might be detectable by the Laser Interferometer Space Antenna (LISA). Such an observation would provide evidence for physics beyond the Standard Model. We study the ability of LISA to observe a gravitational-wave background from phase transitions in the presence of an extragalactic foreground from binary black hole mergers throughout the universe, a galactic foreground from white dwarf binaries, and LISA noise. Modelling the phase transition gravitational wave background as a double broken power law, we use the deviance information criterion as a detection statistic, and Fisher matrix and Markov Chain Monte Carlo methods to assess the measurement accuracy of the parameters of the power spectrum. While estimating all the parameters associated with the gravitational-wave backgrounds, foregrounds, and LISA noise, we find that LISA could detect a gravitational-wave background from phase transitions with a peak frequency of 1 mHz and normalized energy density amplitude of $\OmPeak \simeq 3 \times 10^{-11}$. With $\OmPeak \simeq 10^{-10}$, the signal is detectable if the peak frequency is in the range $4 \times 10^{-4}$ to $9 \times 10^{-3}$ Hz, and the peak amplitude and frequency can be estimated to an accuracy of 10\% to 1\%.

The presence of acoustic waves in models of compressible flows can present complications for analytical and numerical analysis. Therefore, several methods have been developed to filter out these waves, leading to various "sound-proof" models, including the Boussinesq, anelastic and pseudo-incompressible models. We assess the validity of each of these approximate models for describing magnetic buoyancy in the context of the solar interior. A general sound-proof model is introduced and compared to the fully compressible system in a number of asymptotic regimes, including both non-rotating and rotating cases. We obtain specific constraints that must be satisfied in order that the model captures the leading-order behaviour of the fully compressible system. We then discuss which of the existing sound-proof models satisfy these constraints, and in what parameter regimes. We also present a variational derivation of the pseudo-incompressible MHD model, demonstrating its underlying Hamiltonian structure.

No existing spherical convolutional neural network (CNN) framework is both computationally scalable and rotationally equivariant. Continuous approaches capture rotational equivariance but are often prohibitively computationally demanding. Discrete approaches offer more favorable computational performance but at the cost of equivariance. We develop a hybrid discrete-continuous (DISCO) group convolution that is simultaneously equivariant and computationally scalable to high-resolution. While our framework can be applied to any compact group, we specialize to the sphere. Our DISCO spherical convolutions not only exhibit $\text{SO}(3)$ rotational equivariance but also a form of asymptotic $\text{SO}(3)/\text{SO}(2)$ rotational equivariance, which is more desirable for many applications (where $\text{SO}(n)$ is the special orthogonal group representing rotations in $n$-dimensions). Through a sparse tensor implementation we achieve linear scaling in number of pixels on the sphere for both computational cost and memory usage. For 4k spherical images we realize a saving of $10^9$ in computational cost and $10^4$ in memory usage when compared to the most efficient alternative equivariant spherical convolution. We apply the DISCO spherical CNN framework to a number of benchmark dense-prediction problems on the sphere, such as semantic segmentation and depth estimation, on all of which we achieve the state-of-the-art performance.

We consider the Hill four-body problem where three oblate, massive bodies form a relative equilibrium triangular configuration, and the fourth, infinitesimal body orbits in a neighborhood of the smallest of the three massive bodies. We regularize collisions between the infinitesimal body and the smallest massive body, via McGehee coordinate transformation. We describe the corresponding collision manifold and show that it undergoes a bifurcation when the oblateness coefficient of the small massive body passes through the zero value.

We consider an extension of the Standard Model that explains the neutrino masses and has a rich dark matter phenomenology. The model has two dark matter candidates, a vector WIMP and a fermion FIMP, and the sum of their relic densities matches the total dark matter abundance. We extensively study the dark matter production mechanisms and its connection with the neutrino sector, together with various bounds from present and future experiments. The extra scalar field in the model may induce a first-order phase transition in the early Universe. We study the production of stochastic gravitational waves associated with the first-order phase transition. We show that the phase transition can be strong, and thus the model may satisfy one of the necessary conditions for a successful electroweak baryogenesis. Detectability of the phase transition-associated gravitational waves is also discussed.

In a subclass of Horndeski theories with the speed of gravity equivalent to that of light, we study gravitational radiation emitted during the inspiral phase of compact binary systems. We compute the waveform of scalar perturbations under a post-Newtonian expansion of energy-momentum tensors of point-like particles that depend on a scalar field. This scalar mode not only gives rise to breathing and longitudinal polarizations of gravitational waves, but it is also responsible for scalar gravitational radiation in addition to energy loss associated with transverse and traceless tensor polarizations. We calculate the Fourier-transformed gravitational waveform of two tensor polarizations under a stationary phase approximation and show that the resulting waveform reduces to the one in a parametrized post-Einsteinian (ppE) formalism. The ppE parameters are directly related to a scalar charge in the Einstein frame, whose existence is crucial to allow the deviation from General Relativity (GR). We apply our general framework to several concrete theories and show that a new theory of spontaneous scalarization with a higher-order scalar kinetic term leaves interesting deviations from GR that can be probed by the observations of gravitational waves emitted from neutron star-black hole binaries. If the scalar mass exceeds the order of typical orbital frequencies $\omega \simeq 10^{-13}$ eV, which is the case for a recently proposed scalarized neutron star with a self-interacting potential, the gravitational waveform practically reduces to that in GR.

The chirality imbalance in QCD is spontaneously induced by a repulsive axial-vector interaction from the instanton anti-instanton pairing at high temperature above the chiral phase transition, and vanishes at low temperature. Phase transition of the chirality imbalance is always a first-order one in the early universe with $\mathcal {P}$ and $\mathcal {CP}$ violation. The spectra of gravitational waves and the formation of the primordial black holes from this first-order phase transition is investigated in this work, and the effect of a strong magnetic field is also analyzed. The gravitational waves produced by chirality imbalance can be detected by LISA, Taiji and DECIGO, with the peak energy density locating in the range of $10^{-11}$ to $10^{-9}$ and the peak frequencies lying in the range of $10^{-5}$ Hz to $10^{-2}$ Hz. The spectrum with larger axial vector coupling strength and stronger magnetic field has higher peak energy density and lower peak frequency. According to this trend, the gravitational waves spectra might also be able to be detected by SKA, IPTA and EPTA. Based on the mechanism of postponement of the false vacuum decay, it is scarcely possible to form PBHs in our model with typical parameters because the phase transition completes in an extreme short time and thus the false vacuum energy density decays sharply.

We calculate the energetics of compressible and sub-Alfv\'enic turbulence based on the dynamics of coherent cylindrical fluid parcels. We show that parallel and perpendicular magnetic fluctuations are generalized coordinates of the local perturbed Lagrangian of a magnetized fluid, and prove analytically that the bulk of the magnetic energy transferred to kinetic is the energy stored in the coupling between the initial and fluctuating magnetic field, $\vec{B}_{0} \cdot \delta \vec{B}/4\pi$. The analytical relations are consistent with numerical data up to second order terms.