Papers for Friday, Nov 04 2022

We present stellar rotation periods for late K- and early M-dwarf members of the 4 Gyr old open cluster M67 as calibrators for gyrochronology and tests of stellar spin-down models. Using Gaia EDR3 astrometry for cluster membership and Pan-STARRS (PS1) photometry for binary identification, we build this set of rotation periods from a campaign of monitoring M67 with the Canada-France-Hawaii Telescope's MegaPrime wide field imager. We identify 1807 members of M67, of which 294 are candidate single members with significant rotation period detections. Moreover, we fit a polynomial to the period versus color-derived effective temperature sequence observed in our data. We find that the rotation of very cool dwarfs can be explained by a simple solid-body spin-down between 2.7 and 4 Gyr. We compare this rotational sequence to the predictions of gyrochronological models and find that the best match is Skumanich-like spin-down, P_rot \propto t^0.62, applied to the sequence of Ruprecht 147. This suggests that, for spectral types K7-M0 with near-solar metallicity, once a star resumes spinning down, a simple Skumanich-like is sufficient to describe their rotation evolution, at least through the age of M67. Additionally, for stars in the range M1-M3, our data show that spin-down must have resumed prior to the age of M67, in conflict with predictions of the latest spin-down models.

We report a correlation between the presence of a Gaia-Sausage-Enceladus (GSE) analogue and dark matter halo spin in the ARTEMIS simulations of Milky Way-like galaxies. The haloes which contain a large population of accreted stars on highly radial orbits (like the GSE) have lower spin on average than their counterparts with more isotropic stellar velocity distributions. The median modified spin parameters $\lambda^\prime$ differ by a factor of $\sim1.7$ at the present-day, with a similar value when the haloes far from virial equilibrium are removed. We also show that accreted stars make up a smaller proportion of the stellar populations in haloes containing a GSE analogue, and are stripped from satellites with stellar masses typically $\sim4$ times smaller. Our findings suggest that the higher spin of DM haloes without a GSE-like feature is due to mergers with large satellites of stellar mass $\sim10^{10}M_\odot$, which do not result in prominent radially anisotropic features like the GSE.

Surveys with the James Webb Space Telescope (JWST) will have the sensitivity to explore the assembly history of the first nuclear black holes (BHs), as they grow in mass from the scale of the first heavy BH seeds ($\rm M_{BH}\sim 10^4 - 10^6 M_\odot$), to the supermassive black holes (SMBHs, $\rm M_{BH} > 10^8 M_\odot$) powering luminous quasars out to $z \sim 7.5$. In this paper we provide predictions for the number of accreting BHs that would be observable with planned JWST surveys at $5 \le z < 15$. We base our study on the recently developed Cosmic Archaeology Tool (CAT), which allows us to model BH seeds formation and growth, while being consistent with the general population of AGNs and galaxies observed at $4 \le z \le 7$. We find that JWST planned surveys will provide a complementary view on the active BH population at $z > 5$, with JADES-Medium/-Deep being capable of detecting the numerous BHs that populate the faint-end of the distribution, COSMOS-Web sampling a large enough area to detect the rarest brightest systems, and CEERS/PRIMER bridging the gap between these two regimes. The relatively small field of view of the above surveys preferentially selects BHs with masses $\rm 6 \leq Log(M_{BH}/M_\odot) < 8$ at $7 \le z < 10$, residing in relatively metal poor galaxies ($\rm Log(Z/Z_\odot) \ge -2$). At $z \ge 10$, only JADES-Deep will have the sensitivity to detect growing BHs with masses $\rm 4 \leq Log(M_{BH}/M_\odot) < 6$, hosted in even more metal poor environments ($\rm -3 \leq Log(Z/Z_\odot) < -2$). In our model, the latter population corresponds to heavy BH seeds formed by the direct collapse of super-massive stars in their earliest phases of mass growth. Detecting these systems would provide unvaluable insights on the nature and early growth of the first BH seeds.

The buckling process in stellar bars is full of unsolved issues. We analyze the origin of the buckling instability in stellar bars using high-resolution N-body simulations. Previous studies have promoted the nonresonant firehose instability to be responsible for the vertical buckling. We have analyzed the buckling process in terms of the resonant excitation of stellar orbits in the bar, which pumps energy into vertical oscillations. We find that (1) the buckling is associated with an abrupt increase in the central mass concentration and triggers velocities along the bar and along its rotation axis. The velocity field projected on one of the main axes forms circulation cells and increases vorticity, which are absent in firehose instability; (2) The bending amplitude is nonlinear when measured by isodensity contours or curvature of the Laplace plane, which has a substantial effect on the stellar motions; (3) In the linear description, the planar and vertical 2:1 resonances appear only with the buckling and quickly reach the overlapping phase, thus supporting the energy transfer; (4) Using nonlinear orbit analysis, we analyze the stellar oscillations along the bar and along the rotation axis and find that stars cross the vertical 2:1 resonance simultaneously with the buckling. The overlapping planar and vertical 2:1 resonances trapping more than 25% of the bar particles provide the 'smoking gun' pointing to a close relationship between the bending of stellar orbits and the resonant action -- these particles provide the necessary ingredient assuring the cohesive response in the growing vertical asymmetry. We conclude that resonant excitation is important in triggering the buckling instability, and the contribution from the firehose instability should be reevaluated. Finally, we discuss some observational implications of buckling.

Diffuse radio emission has been detected in a considerable number of galaxy clusters and groups, revealing the presence of pervasive cosmic magnetic fields, and of relativistic particles in the large-scale structure (LSS) of the Universe. Since cluster radio emission is faint and steep spectrum, its observations are largely limited by the instrument sensitivity and frequency of observation, leading to a dearth of information, more so for lower-mass systems. The unprecedented sensitivity of recently commissioned low-frequency radio telescope arrays, aided by the development of advanced calibration and imaging techniques, have helped in achieving unparalleled image quality. At the same time, the development of sophisticated numerical simulations and the availability of supercomputing facilities have paved the way for high-resolution numerical modeling of radio emission, and the structure of the cosmic magnetic fields in LSS, leading to predictions matching the capabilities of observational facilities. In view of these rapidly-evolving scenerio in modeling and observations, in this review, we summarise the role of the new telescope arrays and the development of advanced imaging techniques and discuss the detections of various kinds of cluster radio sources. In particular, we discuss observations of the cosmic web in the form of supercluster filaments, studies of emission in poor clusters and groups of galaxies, and of ultra-steep spectrum sources. We also review the current theoretical understanding of various diffuse cluster radio sources and the associated magnetic field and polarization. As the statistics of detections improve along with our theoretical understanding, we update the source classification schemes based on their intrinsic properties. We conclude by summarising the role of the upgraded GMRT and our expectations from the upcoming Square Kilometre Array (SKA) observatories.

Turbulent Radiative Mixing Layers (TRMLs) form at the interface of cold, dense gas and hot, diffuse gas in motion with each other. TRMLs are ubiquitous in and around galaxies on a variety of scales, including galactic winds and the circumgalactic medium. They host the intermediate temperature gases that are efficient in radiative cooling, thus play a crucial role in controlling the cold gas supply, phase structure, and spectral features of galaxies. In this work, we introduce a simple parameterization of the effective turbulent conductivity and viscosity that enables us to develop a simple and intuitive analytic 1.5 dimensional model for TRMLs. Our analytic model reproduces the mass flux, total cooling, and phase structure of 3D simulations of TRMLs at a fraction of the computational cost. It also reveals essential insights into the physics of TRMLs, particularly the importance of the viscous dissipation of relative kinetic energy in balancing radiative cooling. This dissipation takes place both in the intermediate temperature phase, which offsets the enthalpy flux from the hot phase, and in the cold phase, which enhances radiative cooling. Additionally, our model provides a fast and easy way of computing the column density and surface brightness of TRMLs, which can be directly linked to observations.

Gravitational-wave observations of neutron star mergers can probe the nuclear equation of state by measuring the imprint of the neutron star's tidal deformability on the signal. We investigate the ability of future gravitational-wave observations to produce a precise measurement of the equation of state from binary neutron star inspirals. Since measurability of the tidal effect depends on the equation of state, we explore several equations of state that span current observational constraints. We generate a population of binary neutron stars as seen by a simulated Advanced LIGO-Virgo network, as well as by a planned Cosmic Explorer observatory. We perform Bayesian inference to measure the parameters of each signal, and we combine measurements across each population to determine $R_{1.4}$, the radius of a $1.4M_{\odot}$ neutron star. We find that with 321 signals the LIGO-Virgo network is able to measure $R_{1.4}$ to better than 2% precision for all equations of state we consider, however we find that achieving this precision could take decades of observation, depending on the equation of state and the merger rate. On the other hand we find that with one year of observation, Cosmic Explorer will measure $R_{1.4}$ to better than 0.6% precision. In both cases we find that systematic biases, such as from an incorrect mass prior, can significantly impact measurement accuracy and efforts will be required to mitigate these effects.

We present a systematic study of fast neutrino-flavor conversion (FFC) with both small-scale and large-scale numerical simulations in spherical symmetry. We find that FFCs can, in general, reach in a quasi-steady state, and these features in the non-linear phase are not characterized by the growth rate of FFC instability but rather angular structures of electron neutrino lepton number (ELN) and heavy one (XLN). Our result suggests that neutrinos can almost reach a flavor equipartition even in cases with low growth rate of instability (e.g., shallow ELN crossing) and narrow angular regions (in momentum space) where flavor conversions occur vigorously. This exhibits that ELN and XLN angular distributions can not provide a sufficient information to determine total amount of flavor conversion in neutrinos and antineutrinos of all flavors. Based on the results of our numerical simulations, we provide a new approximate scheme of FFC that is designed so that one can easily incorporate effects of FFCs in existing classical neutrino transport codes for the study of core-collapse supernova (CCSN) and binary neutron star merger (BNSM). The scheme has an ability to capture key features of quasi-steady state of FFCs without solving quantum kinetic neutrino transport, which will serve to facilitate access to FFCs for CCSN and BNSM theorists.

Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the large number of parameters needed to model the PTA and GW signals. As the data span and number of pulsars increase, this problem will only worsen. An alternative Monte Carlo sampling method, Hamiltonian Monte Carlo (HMC), utilizes Hamiltonian dynamics to produce sample proposals informed by first-order gradients of the model likelihood. This in turn allows it to converge faster to high dimensional distributions. We implement HMC as an alternative sampling method in our search for an isotropic stochastic GW background, and show that this method produces equivalent statistical results to similar analyses run with standard MCMC techniques, while requiring 100-200 times fewer samples. We show that the speed of HMC sample generation scales as $\mathcal{O}(N_\mathrm{psr}^{5/4})$ where $N_\mathrm{psr}$ is the number of pulsars, compared to $\mathcal{O}(N_\mathrm{psr}^2)$ for MCMC methods. These factors offset the increased time required to generate a sample using HMC, demonstrating the value of adopting HMC techniques for PTAs.

We report the discovery of a transient seen in a strongly lensed arc at redshift $z_{\rm s}=1.2567$ in \emph{Hubble Space Telescope} imaging of the Abell 370 galaxy cluster. The transient is detected at $29.51\pm0.14$ AB mag in a WFC3/UVIS F200LP difference image made using observations from two different epochs, obtained in the framework of the \emph{Flashlights} program, and is also visible in the F350LP band ($m_{\rm F350LP}\sim30.53$ AB). The transient is observed on the negative-parity side of the critical curve at a distance of $\sim 0.6''$ from it, greater than previous examples of lensed stars. The large distance from the critical curve yields a significantly smaller macro-magnification, but our simulations show that bright, O/B-type supergiants can reach sufficiently high magnifications to be seen at the observed position and magnitude. In addition, the observed transient image is a trailing image with an observer-frame time delay of $\sim+0.8$ days from its expected counterpart, so that any transient lasting for longer than that should have also been seen on the minima side and is thus excluded. This, together with the blue color we measure for the transient ($m_{\rm F200LP} - m_{\rm F350LP} \sim [-0.7,-1]$ AB mag), rules out most other transient candidates such as (kilo)novae, for example, and makes a lensed star the prime candidate. Assuming the transient is indeed a lensed star as suggested, many more such events should be detected in the near future in cluster surveys with the \emph{Hubble Space Telescope} and \emph{James Webb Space Telescope}.

We present the first detection of the [N II] 122 $\mu$m and [O III] 52 $\mu$m lines for a reionisation-epoch galaxy. Based on these lines and previous [C II] 158 $\mu$m and [O III] 88 $\mu$m measurements, we estimate an electron density of $\lesssim$ 500 cm$^{-3}$ and a gas-phase metallicity $Z/Z_\odot \sim 1.1 \pm 0.2$ for A1689-zD1, a gravitationally-lensed, dusty galaxy at $z$ = 7.133. Other measurements or indicators of metallicity so far in galaxy ISMs at $z \gtrsim$ 6 are typically an order of magnitude lower than this. The unusually high metallicity makes A1689-zD1 inconsistent with the fundamental metallicity relation, although there is likely significant dust obscuration of the stellar mass, which may partly resolve the inconsistency. Given a solar metallicity, the dust-to-metals ratio is a factor of several lower than expected, hinting that galaxies beyond $z \sim$ 7 may have lower dust formation efficiency. Finally, the inferred nitrogen enrichment compared to oxygen, on which the metallicity measurement depends, indicates that star-formation in the system is older than about 250 Myr, pushing the beginnings of this galaxy to $z >$ 10.

Main Belt Comets (MBCs) exhibit sublimation-driven activity while occupying asteroid-like orbits in the Main Asteroid Belt. MBCs and candidates show stronger clustering of their longitudes of perihelion around 15{\deg} than other objects from the Outer Main Belt (OMB). This potential property of MBCs could facilitate the discovery of new candidates by observing objects in similar orbits. We acquired deep r-band images of 534 targeted asteroids using the INT/WFC between 2018 and 2020. Our sample is comprised of OMB objects observed near perihelion, with longitudes of perihelion between 0{\deg} and 30{\deg} and orbital parameters similar to knowns MBCs. Our pipeline applied activity detection methods to 319 of these objects to look for tails or comae, and we visually inspected the remaining asteroids. Our activity detection pipeline highlighted a faint anti-solar tail-like feature around 2001 NL19 (279870) observed on 2018 November 07, six months after perihelion. This is consistent with cometary activity but additional observations of this object will be needed during its next perihelion to investigate its potential MBC status. If it is active our survey yields a detection rate of $\sim$1:300, which is higher than previous similar surveys, supporting the idea of dynamical clustering of MBCs. If not, it is consistent with previously estimated abundance rates of MBCs in the OMB (<1:500).

Over the entire 20th century, Eta Carinae (\ec) has displayed a unique spectrum, which recently has been evolving towards that of a typical LBV. The two competing scenarios to explain such evolution are: (1) a dissipating occulter in front of a stable star or (2) a decreasing mass loss rate of the star. The first mechanism simultaneously explains why the central star appears to be secularly increasing its apparent brightness while its luminosity does not change; why the Homunculus' apparent brightness remains almost constant; and why the spectrum seen in direct light is becoming more similar to that reflected from the Homunculus (and which resembles a typical LBV). The second scenario does not account for these facts and predicts an increase in the terminal speed of the wind, contrary to observations. In this work, we present new data showing that the P Cygni absorption lines are secularly strengthening, which is not the expected behaviour for a decreasing wind-density scenario. CMFGEN modelling of the primary's wind with a small occulter in front agrees with observations. One could argue that invoking a dissipating coronagraphic occulter makes this object even more peculiar than it already appears to be. However, on the contrary, it solves the apparent contradictions between many observations. Moreover, by assigning the long-term behaviour to circumstellar causes and the periodic variations due to binarity, a star more stable after the 1900s than previously thought is revealed, contrary to the earlier paradigm of an unpredictable object.

The lack of objects between $2\,M_{\odot}$ and $5\,M_{\odot}$ in the joint mass distribution of compact objects has been termed "mass gap", and attributed mainly to the characteristics of the supernova mechanism precluding their birth. However, recent observations show that a number of candidates reported to lie inside the "gap" may fill it, suggesting instead a paucity that may be real or largely a result of small number statistics. We quantify in this work the individual candidates and evaluate the joint probability of a mass gap. Our results show that an absolute mass gap is not present, to a very high confidence level. It remains to be seen if a relative paucity of objects stands in the future, and how this population can be related to the formation processes, which may include neutron star mergers, collapse of a neutron star to a black hole and others.

We present the compilation of the Fifth Catalogue of Nearby Stars (CNS5), based on astrometric and photometric data from Gaia EDR3 and Hipparcos, and supplemented with parallaxes from ground-based astrometric surveys carried out in the infrared. The aim of the CNS5 is to provide the most complete sample of objects in the solar neighbourhood. For all known stars and brown dwarfs in the 25 pc sphere around the Sun, basic astrometric and photometric parameters are given. Furthermore, we provide the colour-magnitude diagram and various luminosity functions of the stellar content in the solar neighbourhood, and characterise the completeness of the CNS5 catalogue. We compile a sample of stars and brown dwarfs which most likely are located within 25 pc of the Sun, taking space-based parallaxes from Gaia EDR3 and Hipparcos as well as ground-based parallaxes from Best et al. (2021), Kirkpatrick et al. (2021), and from the CNS4 into account. We develop a set of selection criteria to clean the sample from spurious sources. Furthermore, we show that effects of blending in the Gaia photometry, which affect mainly the faint and red sources in Gaia, can be mitigated, to reliably place those objects in a colour-magnitude diagram. We also assess the completeness of the CNS5 using a Kolmogorov-Smirnov test and derive observational optical and mid-infrared luminosity functions for the main-sequence stars and white dwarfs in the solar neighbourhood. The CNS5 contains 5931 objects, including 5230 stars (4946 main-sequence stars, 20 red giants and 264 white dwarfs) and 701 brown dwarfs. We find that the CNS5 catalogue is statistically complete down to 19.7 mag in G-band and 11.8 mag in W1-band absolute magnitudes, corresponding to a spectral type of L8.

We report on JWST commissioning observations of the transiting exoplanet HAT-P-14 b, obtained using the Bright Object Time Series (BOTS) mode of the NIRSpec instrument with the G395H/F290LP grating/filter combination ($3-5\mu$m). While the data were used primarily to verify that the NIRSpec BOTS mode is working as expected, and to enable it for general scientific use, they yield a precise transmission spectrum which we find is featureless down to the precision level of the instrument, consistent with expectations given HAT-P-14~b's small scale-height and hence expected atmospheric features. The exquisite quality and stability of the \emph{JWST/NIRSpec} transit spectrum -- almost devoid of any systematic effects -- allowed us to obtain median uncertainties of 50-60 ppm in this wavelength range at a resolution of $R=100$ in a single exposure, which is in excellent agreement with pre-flight expectations and close to the (or at the) photon-noise limit for a $J = 9.094$, F-type star like HAT-P-14. These observations showcase the ability of NIRSpec/BOTS to perform cutting-edge transiting exoplanet atmospheric science, setting the stage for observations and discoveries to be made in Cycle 1 and beyond.

In this study, we use archival data from HST, Chandra, XMM-Newton, and Swift-XRT, to probe the nature of 9 (X1-X9) candidate ULXs in NGC 1672. Our study focuses on using the precise source positions obtained via improved astrometry based on {\it Chandra} and HST observations to search for and identify optical counterparts for these ULXs.Unique optical counterparts are identified for X2 an{d X6; two potential counterparts were determined for X1, X5 and X7 within the respective error radii while no optical counterparts were found for the remaining four sources. Based on spectral energy distributions (SEDs), X-ray and optical temporal analyses, some evidences about the nature of X1 and X2 were obtained.

We present preliminary results from our long-baseline interferometry (LBI) survey to constrain the multiplicity properties of intermediate-mass A-type stars within 80pc. Previous multiplicity studies of nearby stars exhibit orbital separation distributions well-fitted with a log-normal with peaks > 15au, increasing with primary mass. The A-star multiplicity survey of De Rosa et al. (2014), sensitive beyond 30au but incomplete below 100 au, found a log-normal peak around 390au. Radial velocity surveys of slowly-rotating, chemically peculiar Am stars identified a significant number of very close companions with periods $\leq$ 5 days, ~ 0.1au, a result similar to surveys of O- and B-type primaries. With the improved performance of LBI techniques, we can probe these close separations for normal A-type stars where other surveys are incomplete. Our initial sample consists of 27 A-type primaries with estimated masses between 1.44-2.49M$_{\odot}$ and ages 10-790Myr, which we observed with the MIRC-X instrument at the CHARA Array. We use the open source software CANDID to detect five companions, three of which are new, and derive a companion frequency of 0.19$^{+0.11}_{-0.06}$ over mass ratios 0.25-1.0 and projected separations 0.288-5.481 au. We find a probability of 10$^{-6}$ that our results are consistent with extrapolations based on previous models of the A-star companion population, over mass ratios and separations sampled. Our results show the need to explore these very close separations to inform our understanding of stellar formation and evolution processes.

We aim to revisit the system AB Pic which has a known companion at the exoplanet/ brown-dwarf boundary. We based this study on a rich set of observations to investigate the companion's orbit and atmosphere. We composed a spectrum of AB Pic b merging archival VLT/SINFONI K-band data, with published spectra at J and H-band (SINFONI) and Lp-band (Magellan-AO), and photometric measurements (HST and Spitzer). We modeled the spectrum with ForMoSA, based on two atmospheric models: ExoREM and BT-SETTL13. We determined the orbital properties of b fitting the astrometric measurements from NaCo (2003 and 2004) and SPHERE (2015). The orbital solutions favor a semi-major axis of $\sim$190au viewed edge-on. With Exo-REM, we derive a T$_{eff}$ of 1700$\pm$50K and surface gravity of 4.5$\pm$0.3dex, consistent with previous works, and we report for the first time a C/O ratio of 0.58$\pm$0.08 ($\sim$solar). The posteriors are sensitive to the wavelength interval and the family of models used. Given the 2.1hr rotation period and our vsin(i) of $\sim$73km/s, we estimate for the first time the true obliquity to be $\sim$45 or $\sim$135deg, indicating a significant misalignment between the planet's spin and orbit orientations. Finally, the existence of a proper motion anomaly between the Hipparcos and Gaia eDR3 compared to our SPHERE detection limits and adapted radial velocity limits indicate the existence of a $\sim$6M$_{Jup}$ inner planet orbiting from 2 to 10au (40-200mas). The possible existence of an inner companion, together with the likely miss-alignment of the spin axis orientation, strongly favor a formation path by gravitational instability or core accretion within a disk closer inside followed by dynamical interactions. Confirmation and characterization of planet c and access to a broader wavelength coverage for planet b will be essential to probe the uncertainties associated with the parameters.

The hot plasma within merging galaxy clusters is predicted to be filled with shocks and turbulence that may convert part of their kinetic energy into relativistic electrons and magnetic fields generating synchrotron radiation. Analyzing Low Frequency Array (LOFAR) observations of the galaxy cluster Abell 2255, we show evidence of radio synchrotron emission distributed over very large scales of at least 5 megaparsec. The pervasive radio emission witnesses that shocks and turbulence efficiently transfer kinetic energy into relativistic particles and magnetic fields in a region that extends up to the cluster outskirts. The strength of the emission requires a magnetic field energy density at least 100 times higher than expected from a simple compression of primordial fields, presumably implying that dynamo operates efficiently also in the cluster periphery. It also suggests that nonthermal components may contribute substantially to the pressure of the intracluster medium in the cluster periphery.

The EFIELD instrument is part of the geophysics and meteorology sensor package DraGMet on the Dragonfly mission, which will explore the surface of Titan in the mid-2030s. EFIELD consists of two electrodes designed to passively record the AC electric field at each landing site. The exploration zone of Dragonfly will mostly consist of dune fields, covered with sand grains. Little is known on the properties of these grains, although Cassini-Huygens observations suggest they are mostly made of organic material produced by Titan's atmospheric photochemistry and evolved at the surface. Little is known also about dune formation and in general about the transport of sediments by winds. The latter much depends on inter-particle forces and therefore on how grains are charged by friction. We demonstrate here that the EFIELD experiment can bring new insights on these questions. We have developed a hydrodynamic-electrostatic model to simulate the trajectory of a wind-blown charged sand grain in the vicinity of an idealized EFIELD probe and to predict how such a grain flying close to the probe would affect its potential. We show that, in some conditions, the resulting perturbation will be strong enough to be detected by the EFIELD probe. More specifically, we find that the detection of typical charged wind-blown grains (200 microns) on Titan requires an instrument standard deviation noise inferior to 1mV, though occasional larger grains flying close to one electrode could be detected with a higher noise level. Furthermore, we propose a method to retrieve information on the charge and velocity of wind-blown charged grains detected by the EFIELD experiment. This method well applies for cases where the particle trajectory can be regarded as quasi-linear. We validate our inversion approach on both synthetic and experimental data obtained with a laboratory prototype of the EFIELD experiment.

In molecular cloud cores, the cosmic ray (CR) induced sputtering via CR ion-icy grain collision is one of the desorption processes for ice molecules from mantles around dust grains. The efficiency of this process depends on the incident CR ion properties as well as the physicochemical character of the ice mantle. Our main objective is the examination of the sputtering efficiency for H$_2$O and CO ices found in molecular cloud cores. In the calculation routine, we consider a multi-dimensional parameter space that consists of thirty CR ion types, five different CR ion energy flux distributions, two separate ice mantle components (pure H$_2$O and CO), three ice formation states, and two sputtering regimes (linear and quadratic). We find that the sputtering behavior of H$_2$O and CO ices is dominated by the quadratic regime rather than the linear regime, especially for CO sputtering. The sputtering rate coefficients for H$_2$O and CO ices show distinct variations with respect to the adopted CR ion energy flux as well as the grain size-dependent mantle depth. The maximum radius of the cylindrical latent region is quite sensitive to the effective electronic stopping power. The track radii for CO ice are much bigger than H$_2$O ice values. In contrast to the H$_2$O mantle, even relatively light CR ions ($Z \geq 4$) may lead to a track formation within the CO mantle, depending on ${\rm S_{\rm e,eff}}$. We suggest that the latent track formation threshold can be assumed as a separator between the linear and the quadratic regimes for sputtering.

As we further our studies on Gamma-ray bursts (GRBs), both on theoretical models and observational tools, more and more options begin to open for exploration of its physical properties. As transient events primarily dominated by synchrotron radiation, it is expected that the synchrotron photons emitted by GRBs should present some degree of polarization throughout the evolution of the burst. Whereas observing this polarization can still be challenging due to the constraints on observational tools, especially for short GRBs, it is paramount that the groundwork is laid for the day we have abundant data. In this work, we present a polarization model linked with an off-axis spreading top-hat jet synchrotron scenario in a stratified environment with a density profile $n(r)\propto r^ {-k}$. We present this model's expected temporal polarization evolution for a realistic set of afterglow parameters constrained within the values observed in the GRB literature for four degrees of stratification $k=0,1,1.5 {\rm \, and\,} 2$ and two magnetic field configurations with high extreme anisotropy. We apply this model and predict polarization from a set of GRBs exhibiting off-axis afterglow emission. In particular, for GRB 170817A, we use the available polarimetric upper limits to rule out the possibility of a extremely anisotropic configuration for the magnetic field.

We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model in comparison to cold dark matter (CDM) models from a weak lensing perspective using perturbative methods. We use Eulerian perturbation theory up to fourth order to compute the tree-level matter trispectrum and the one-loop matter spectrum and bispectrum from consistently chosen initial conditions. We go on to derive the respective lensing spectra, bispectra and trispectra in CDM and FDM in the context of a Euclid-like weak lensing survey. Finally, we compute the attainable cumulative signal-to-noise ratios and an estimate of the attainable $\chi^2$-functionals for distinguishing FDM from CDM at particle masses $m=10^{-21}$ eV, $m = 10^{-22}$ eV and $m = 10^{-23}$ eV. We find that a weak lensing survey can potentially be used to distinguish between the FDM and CDM cases up to a mass of $m = 10^{-22}$ eV.

Single flux density measurements at observed-frame sub-millimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive haloes from the Feedback in Realistic Environments (FIRE) project. Our starting point is four star-forming, central galaxies at seven redshifts between z=1.5 and z=4.5. We generate modified quiescent galaxies that have been quenched for 100Myr, 500Myr, or 1Gyr prior to each of the studied redshifts by re-assigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature, T_dust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumed T_dust~25K, we show that the ratio between dust mass and 1.3mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non star-forming galaxies. This can lead to an underestimation of dust mass (and, when sub-mm flux density is used as a proxy for molecular gas content, gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3mm flux density probes rest-frame wavelengths far from the Rayleigh-Jeans regime, and hence depends super-linearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses.

More than 5,000 exoplanets have been confirmed and among them almost 4,000 were discovered by the transit method. However, few transiting exoplanets have an orbital period greater than 100 days. Here we report a transit detection of Kepler-167 e, a "Jupiter analog" exoplanet orbiting a K4 star with a period of 1,071 days, using the Unistellar ground-based telescope network. From 2021 November 18 to 20, citizen astronomers located in nine different countries gathered 43 observations, covering the 16 hour long transit. Using a nested sampling approach to combine and fit the observations, we detected the mid-transit time to be UTC 2021 November 19 17:20:51 with a 1$\sigma$ uncertainty of 9.8 minutes, making it the longest-period planet to ever have its transit detected from the ground. This is the fourth transit detection of Kepler-167 e, but the first made from the ground. This timing measurement refines the orbit and keeps the ephemeris up to date without requiring space telescopes. Observations like this demonstrate the capabilities of coordinated networks of small telescopes to identify and characterize planets with long orbital periods.

The halo stars on highly radial orbits should inevitably pass the center regions of the Milky Way. Under the assumption that the stellar halo is in dynamical equilibrium and axisymmetric, we integrate the orbits of $\sim 10,000$ halo K-giants at $5\leq r \leq 50$ kpc cross-matched from LAMOST DR5 and $Gaia$ DR3. By carefully considering the selection function, we construct the stellar halo distribution at the entire regions of $r \leq 50$ kpc. We find that a double-broken power-law function well describes the stellar halo density distribution with shallower slopes in the inner regions and the two breaks at $r=10$ kpc and $r=25$ kpc, respectively. The stellar halo becomes flatter from outer to inner regions but has $q\sim 0.5$ at $r \lesssim 5$ kpc. The stellar halo becomes isotropic with a slight prograde rotation in the inner 5 kpc, and reaches velocity dispersions of $\sim 250\rm \ km\ s^{-1}$. We get a weak negative metallicity gradient of $-0.005$ dex kpc$^{-1}$ at $5\leq r \leq 50$ kpc, while there is an excess of relative metal-rich stars with [Fe/H]$>-1$ in the inner 10 kpc. The halo interlopers at $r \leq 5$ kpc from integration of our sample has a mass of $\sim1.2 \times 10^8\ M_{\odot}$ ($\sim 4.7 \times 10^7\ M_{\odot}$ at [Fe/H]$<-1.5$), which can explain 50-100% of the metal-poor stars with [Fe/H]$<-1.5$ directly observed in the Galactic central regions.

The magnetar SGR 1935+2154 is reported to have an anti-glitch, accompanied by fast radio bursts, and transient pulsed radio emission. In the wind braking model, this triplet event tell people that (1) SGR 1935+2154 does not have a strong particle wind and can be approximated by magnetic dipole braking in the persistent state; (2) Its anti-glitch is due to an enhanced particle wind, similar to the first anti-glitch in magnetars; (3) Its transient pulsed radio emission may due to enlarged and decreasing emission beam during the outburst. The enhanced particle acceleration potential and pulsar death line may not be the dominate factor.

Observations by the Interface Region Imaging Spectrograph (IRIS) in the Mg II h & k spectral lines have provided a new diagnostic window towards the knowledge of the complex physical conditions in the solar chromosphere. Theoretical efforts focused on understanding the behavior of these lines have allowed us to obtain a better and more accurate vision of the chromosphere. These efforts include forward modeling, numerical simulations, and inversions. In this paper, we focus our attention on the uncertainties associated with the thermodynamic model atmosphere obtained after the inversion of the Mg II h & k lines. We have used ~ 50;000 synthetic representative profiles of the IRIS2 database to characterize the most important source of uncertainties in the inversion process, viz.: the inherent noise of the observations, the random initialization of process, and the selection criteria in a high-dimensional space. We have applied a Monte Carlo approach to this problem. Thus, for a given synthetic representative profile, we have created five randomized noise realizations (representative of the most popular exposure times in the IRIS observations), and inverted these profiles five times with different inversion initializations. The resulting 25 inverted profiles, fits to noisy data, and model atmospheres are then used to determine the uncertainty in the model atmosphere, based on the standard deviation and empirical selection criteria for the goodness of fit. With this approach, the new uncertainties of the models available in the IRIS2 database are more reliable at the optical depths where the Mg II h & k lines are sensitive to changes in the thermodynamics.

Blue large-amplitude pulsators (BLAPs) are a newly discovered type of variable star. Their typical pulsation periods are on the order of a few tens of minutes, with relatively large amplitudes of 0.2-0.4 mag in optical bands, and their rates of period changes are on the order of $10^{-7} yr^{-1}$ (both positive and negative). They are extremely rare objects and attempts to explain their origins and internal structures have attracted a great deal of attention. Previous studies have proposed that BLAPs may be pre-white dwarfs, with masses around $0.3M_\odot$, or core-helium-burning stars in the range of $\sim 0.7-1.1M_\odot$. In this work, we use a number of MESA models to compute and explore whether BLAPs could be explained as shell helium-burning subdwarfs type B (SHeB sdBs). The models that best match existing observational constraints have helium core masses in the range of $\sim 0.45-0.5M_\odot$. Our model predicts that the positive rate of period change may evolve to negative. The formation channels for SHeB sdBs involve binary evolution and although the vast majority of BLAPs do not appear to be binaries (with the exception of HD 133729), the observational constraints are still very poor. Motivated by these findings, we explored the Roche lobe overflow channel. Of the 304 binary evolution models we computed, about half of them are able to produce SHeB sdBs in long-period binaries that evade detection from the limited observations that are currently available.

Removing optical and atmospheric blur from galaxy images significantly improves galaxy shape measurements for weak gravitational lensing and galaxy evolution studies. This ill-posed linear inverse problem is usually solved with deconvolution algorithms enhanced by regularisation priors or deep learning. We introduce a so-called "physics-based deep learning" approach to the Point Spread Function (PSF) deconvolution problem in galaxy surveys. We apply algorithm unrolling and the Plug-and-Play technique to the Alternating Direction Method of Multipliers (ADMM) with a Poisson noise model and use a neural network to learn appropriate priors from simulated galaxy images. We characterise the time-performance trade-off of several methods for galaxies of differing brightness levels, showing an improvement of 26% (SNR=20)/48% (SNR=100) compared to standard methods and 14% (SNR=20) compared to modern methods.

The azimuthal asymmetries of dust rings in protoplanetary disks such as a crescent around young stars are often interpreted as dust traps, and thus as ideal locations for planetesimal and planet formations. Whether such dust traps effectively promote planetesimal formation in disks around very-low-mass stars (VLM; a mass of $\lesssim$0.2~$M_\odot$) is debatable, as the dynamical and grain growth timescales in such systems are long. To investigate grain growth in such systems, we studied the dust ring with crescent around the VLM star ZZ~Tau~IRS using the Karl G. Jansky Very Large Array (JVLA) at centimeter wavelengths. Significant signals were detected around ZZ~Tau~IRS. To estimate the maximum grain size ($a_{\rm max}$) in the crescent, we compared the observed spectral energy distribution (SED) with SEDs for various $a_{\rm max}$ values predicted by radiative transfer calculations. We found $a_{\rm max} \gtrsim$~1~mm and $\lesssim$~60~$\mu$m in the crescent and ring, respectively, though our modeling efforts rely on uncertain dust properties. Our results suggest that grain growth occurred in the ZZ~Tau~IRS disk, relative to sub-micron-sized interstellar medium. Planet formation in crescent with mm-sized pebbles might proceed more efficiently than in other regions with sub-millimeter-sized pebbles via pebble accretion scenarios.

The studies of the complex molecular lines towards the hot molecular cores at millimeter and submillimeter wavelengths provide instructive knowledge about the chemical complexity in the interstellar medium (ISM). We presented the detection of the rotational emission lines of the complex nitrogen-bearing molecule ethyl cyanide (C$_{2}$H$_{5}$CN) towards the hot molecular core G10.47+0.03 using the Atacama Large Millimeter/Submillimeter Array (ALMA) band 4 observations. The estimated column density of C$_{2}$H$_{5}$CN towards G10.47+0.03 was (9.5$\pm$0.1)$\times$10$^{16}$ cm$^{-2}$ with the rotational temperature of 223.8$\pm$4.3 K. The estimated fractional abundance of C$_{2}$H$_{5}$CN with respect to H$_{2}$ towards G10.47+0.03 was 1.90$\times$10$^{-8}$. We observed that the estimated fractional abundance of C$_{2}$H$_{5}$CN is nearly similar to the simulated abundance of C$_{2}$H$_{5}$CN which was predicted by the three-phase warm-up model from Garrod (2013). We also discuss the possible formation mechanism of C$_{2}$H$_{5}$CN towards the hot molecular cores and we claimed the barrierless and exothermic radical-radical reaction between CH$_{2}$ and CH$_{2}$CN is responsible for the production of high abundant C$_{2}$H$_{5}$CN ($\sim$10$^{-8}$) towards G10.47+0.03 during the warm-up phases.

Cosmic Rays (CRs) interact with the diffuse gas, radiation, and magnetic fields in the interstellar medium (ISM) to produce electromagnetic emissions that are a significant component of the all-sky flux across a broad wavelength range. The Fermi Large Area Telescope (LAT) has measured these emissions at GeV $\gamma$-ray energies with high statistics. Meanwhile, the High-Energy Stereoscopic System (H.E.S.S.) telescope array has observed large-scale Galactic diffuse emission in the TeV $\gamma$-ray energy range. The emissions observed at GeV and TeV energies are connected by the common origin of the CR particles injected by the sources, but the energy dependence of the mixture from the general ISM (true `diffuse'), those emanating from the relatively nearby interstellar space about the sources, and the sources themselves, is not well understood. In this paper, we investigate predictions of the broadband emissions using the GALPROP code over a grid of steady-state 3D models that include variations over CR sources, and other ISM target distributions. We compare, in particular, the model predictions in the VHE ($\geq$100 GeV) $\gamma$-ray range with the H.E.S.S. Galactic plane survey (HGPS) after carefully subtracting emission from catalogued $\gamma$-ray sources. Accounting for the unresolved source contribution, and the systematic uncertainty of the HGPS, we find that the GALPROP model predictions agree with lower estimates for the HGPS source-subtracted diffuse flux. We discuss the implications of the modelling results for interpretation of data from the next generation Cherenkov Telescope Array (CTA).

The detection of benzonitrile (C6H5CN), 1- and 2-cyano-naphthalene (C10H7CN) in the cold, dark molecular cloud TMC-1 at centimetre (cm) wavelengths has opened up prospects for the detection of other N- and CN-containing polycyclic aromatic hydrocarbons (PAHs). In this light, the pure rotational spectra of N-pyrene (C15H9N), CN-pyrene (C15H9CN), N-coronene (C23H11N) and CN-coronene (C23H11CN) are reported here for the first time. The B3LYP/6-311+G(d,p) level of theory, in the Density Functional Theory (DFT) calculations, achieves the best performance for calculating the spectroscopic parameters and simulating the rotational spectra. The large permanent dipole moment of CN-PAHs makes them the most suitable PAH species for detection in the interstellar medium. Additionally, pyrene's smaller partition function makes CN-pyrene a prime candidate to be discovered in cold, dark molecular clouds such as the TMC-1. The present work sets a benchmark for theoretical rotational spectra of N- and CN-containing PAHs and may act as a guide for laboratory experiments and observational searches.

NGC 1266 is a lenticular galaxy (S0) hosting an active galactic nucleus (AGN), and known to contain a large amount of shocked gas. We compare the luminosity ratio of mid-\emph{J} CO lines to IR continuum with star-forming galaxies (SFGs), and then model the CO spectral line energy distribution (SLED). We confirm that in the mid- and high-\emph{J} regions ($J_{\rm up}$ = 4--13), the C-type shock ($v_{\rm s}$ = 25 km s$^{-1}$, $n_{\rm H}$ = $5\times10^{4}$ cm$^{-3}$) can reproduce the CO observations well. The galaxy spectral energy distribution (SED) is constructed and modeled by the code {\tt X-CIGALE} and obtains a set of physical parameters including the star formation rate (SFR, 1.17 $\pm$ 0.47 \emph{M$_{\odot}$}yr$^{-1}$). Also, our work provides SFR derivation of [C\,{\sc ii}] from the neutral hydrogen regions only (1.38 $\pm$ 0.14 $M_{\odot}$yr$^{-1}$). Previous studies have illusive conclusions on the AGN or starburst nature of the NGC 1266 nucleus. Our SED model shows that the hidden AGN in the system is intrinsically low-luminosity, consequently the infrared luminosity of the AGN does not reach the expected level. Archival data from \emph{NuSTAR} hard X-ray observations in the 3--79 keV band shows a marginal detection, disfavoring presence of an obscured luminous AGN and implying that a compact starburst is more likely dominant for the NGC 1266 nucleus.

The Kodaikanal Observatory has made synoptic observations of the Sun in white light since 1904, and these images are sketched on the Stonyhurst grids called sun charts. These continuous hand-drawn data sets are used for long-term studies of the Sun. This article investigates temporal and periodic variations of the monthly hemispheric sunspot number and sunspot group area for 1905--2016, covering solar cycles 14 to 24. We find that the temporal variations of the sunspot number and group area are different in each hemisphere and peak at different times of the solar cycle in the opposite hemisphere. For both the data sets, Cycle 19 shows maximum amplitude. For the sunspot number time series, Cycle 24 was the weakest, and Cycle 15 for the group area. The existence of double peaks and violation of the ``odd-even rule'' was found in both data sets. We have studied the periodic and quasi-periodic variations in both the time series by wavelet technique. We noticed that along with the fundamental mode of the $\sim$ 11~year cycle and polarity reversal period of 22~years, the sunspot activity data also exhibited several mid-term periodicities in the opposite hemispheres, specifically the Riger group and quasi-biennial periodicities. The temporal evolution of these detected quasi-periodicities also differs in the northern and southern hemispheres. We analyzed the data set statistically to understand the bulk properties and coupling between the opposite hemispheres. The study indicates that the two hemispheric data sets differ, but some dependency could be present.

TIRCAM2 is the facility near-infrared Imager at the Devasthal 3.6-m telescope in northern India, equipped with an Aladdin III InSb array detector. We have pioneered the use of TIRCAM2 for very fast photometry, with the aim of recording Lunar Occultations (LO). This mode is now operational and publicly offered. In this paper we describe the relevant instrumental details, we provide references to the LO method and the underlying data analysis procedures, and we list the LO events recorded so far. Among the results, we highlight a few which have led to the measurement of one small-separation binary star and of two stellar angular diameters. We conclude with a brief outlook on further possible instrumental developments and an estimate of the scientific return. In particular, we find that the LO technique can detect sources down to K~ 9 mag with SNR=1 on the DOT telescope. Angular diameters larger than ~ 1 milliarcsecond (mas) could be measured with SNR above 10, or K~6 mag. These numbers are only an indication and will depend strongly on observing conditions such as lunar phase and rate of lunar limb motion. Based on statistics alone, there are several thousands LO events observable in principle with the given telescope and instrument every year.

State-of-the-art techniques to identify H\alpha emission line sources in narrow-band photometric surveys consist of searching for H\alpha excess with reference to nearby objects in the sky (position-based selection). However, while this approach usually yields very few spurious detections, it may fail to select intrinsically faint and/or rare H\alpha-excess sources. In order to obtain a more complete representation of the heterogeneous emission line populations, we recently developed a technique to find outliers relative to nearby objects in the colour-magnitude diagram (CMD-based selection). By combining position-based and CMD-based selections, we built an updated catalogue of H\alpha-excess candidates in the northern Galactic Plane. Here we present spectroscopic follow-up observations and classification of 114 objects from this catalogue, that enable us to test our novel selection method. Out of the 70 spectroscopically confirmed H\alpha emitters in our sample, 15 were identified only by the CMD-based selection, and would have been thus missed by the classic position-based technique. In addition, we explore the distribution of our spectroscopically confirmed emitters in the Gaia CMD. This information can support the classification of emission line sources in large surveys, such as the upcoming WEAVE and 4MOST, especially if augmented with the introduction of other colours.

We report on the Imaging X-ray Polarimetry Explorer (IXPE) observation of the closest and X-ray brightest Compton-thick active galactic nucleus (AGN), the Circinus galaxy. We find the source to be significantly polarized in the 2--6 keV band. From previous studies, the X-ray spectrum is known to be dominated by reflection components, both neutral (torus) and ionized (ionization cones). Our analysis indicates that the polarization degree is $28 \pm 7$ per cent (at 68 per cent confidence level) for the neutral reflector, with a polarization angle of $18{\deg} \pm 5{\deg}$, roughly perpendicular to the radio jet. The polarization of the ionized reflection is unconstrained. A comparison with Monte Carlo simulations of the polarization expected from the torus shows that the neutral reflector is consistent with being an equatorial torus with a half-opening angle of 45{\deg}-55{\deg}. This is the first X-ray polarization detection in a Seyfert galaxy, demonstrating the power of X-ray polarimetry in probing the geometry of the circumnuclear regions of AGNs, and confirming the basic predictions of standard Unification Models.

We develop a simple analytic model for what happens when an O star (or compact cluster of OB stars) forms in a shock compressed layer and carves out an approximately circular hole in the layer, at the waist of a bipolar HII Region (HIIR). The model is characterised by three parameters: the half-thickness of the undisturbed layer, Zlay, the mean number-density of hydrogen molecules in the undisturbed layer, nlay, and the (collective) ionising output of the star(s), NdotLyC. The radius of the circular hole is given by WIF ~ 3.8 pc [Zlay/0.1pc]^{-1/6} [nlay/10^4cm^{-3}]^{-1/3} [NdotLyC/10^{49} s^{-1}]^{1/6} [t/Myr]^{2/3}. Similar power-law expressions are obtained for the rate at which ionised gas is fed into the bipolar lobes; the rate at which molecular gas is swept up into a dense ring by the shock front (SF) that precedes the ionisation front (IF); and the density in this dense ring. We suggest that our model might be a useful zeroth-order representation of many observed HIIRs. From viewing directions close to the midplane of the layer, the HIIR will appear bipolar. From viewing directions approximately normal to the layer it will appear to be a limb-brightened shell but too faint through the centre to be a spherically symmetric bubble. From intermediate viewing angles more complicated morphologies can be expected.

While a large fraction of the stars are in multiple systems, our understanding of the processes leading to the formation of these systems is still inadequate. Given the large theoretical uncertainties, observation plays a basic role. Here we discuss the contribution of high contrast imaging, and more specifically of the SPHERE instrument at the ESO Very Large Telescope, in this area. SPHERE nicely complements other techniques - in particular those exploiting Gaia and ALMA - in detecting and characterising systems near the peak of the distribution with separation and allows to capture snapshots of binary formation within disks that are invaluable for the understanding of disk fragmentation.

Here we present an initial look at the dynamics and stability of 178 multiplanetary systems which are already confirmed and listed in the NASA Exoplanet Archive. To distinguish between the chaotic and regular nature of a system, the value of the MEGNO indicator for each system was determined. Almost three-quarters of them could be labelled as long-term stable. Only 45 studied systems show chaotic behaviour. We consequently investigated the effects of the number of planets and their parameters on the system stability. A comparison of results obtained using the MEGNO indicator and machine-learning algorithm SPOCK suggests that the SPOCK could be used as an effective tool for reviewing the stability of multiplanetary systems. A similar study was already published by Laskar and Petit in 2017. We compared their analysis based on the AMD criterion with our results. The possible discrepancies are discussed.

We develop an exact formalism for the computation of the abundance of primordial black holes (PBHs) in the presence of local non-gaussianity (NG) in the curvature perturbation field. For the first time, we include NG going beyond the widely used quadratic and cubic approximations, and consider a completely generic functional form. Adopting threshold statistics of the compaction function, we address the computation of the abundance both for narrow and broad power spectra. While our formulas are generic, we discuss explicit examples of phenomenological relevance considering the physics case of the curvaton field. We carefully assess under which conditions the conventional perturbative approach can be trusted. In the case of a narrow power spectrum, this happens only if the perturbative expansion is pushed beyond the quadratic order (with the optimal order of truncation that depends on the width of the spectrum). Most importantly, we demonstrate that the perturbative approach is intrinsically flawed when considering broad spectra, in which case only the non-perturbative computation captures the correct result. Finally, we describe the phenomenological relevance of our results for the connection between the abundance of PBHs and the stochastic gravitational wave (GW) background related to their formation. As NGs modify the amplitude of perturbations necessary to produce a given PBHs abundance and boost PBHs production at large scales for broad spectra, modelling these effects is crucial to connect the PBH scenario to its signatures at current and future GWs experiments.

Neutrino astronomy saw its birth with the discovery by IceCube of a diffuse flux at energies above 60 TeV with intensity comparable to a predicted upper limit to the flux from extra-galactic sources of ultra-high energy cosmic rays (UHECRs). While such an upper limit corresponds to the case of calorimetric sources, in which cosmic rays lose all their energy into photo-pion production, the first statistically significant coincident observation between neutrinos and gamma rays was observed from a blazar of intriguing nature. A very-high-energy muon event, of most probable neutrino energy of 290 TeV for an $E^{-2.13}$ spectrum, alerted other observatories triggering a large number of investigations in many bands of the electromagnetic (EM) spectrum. A high gamma-ray state from the blazar was revealed soon after the event and in a follow-up to about 40 days. A posteriori observations also in the optical and radio bands indicated a rise of the flux from the TXS 0506+056 blazar. A previous excess of events of the duration of more than 100~d was observed by IceCube with higher significance than the alert itself. These observations triggered more complex modeling than simple one-zone proton synchrotron models for proton acceleration in jets of active galactic nuclei (AGNs) and more observations across the EM spectrum. A second piece of evidence was a steady excess of about 50 neutrino events with reconstructed soft spectrum in a sample of lower energy well-reconstructed muon events than the alert event. A hot spot was identified in a catalog of 110 gamma-ray intense emitters and starburst galaxies in a direction compatible with NGC 1068 with a significance of $2.9\sigma$. NGC 1068 hosts a mildly relativistic jet in a starburst galaxy, seen not from the jet direction but rather through the torus. This Seyfert II galaxy is at only 14.4~Mpc from the Earth. We discuss these observations.

We have carried out a comprehensive study of the temperature structure of plasmoids, which successively occurred in recurrent active region jets. The multithermal plasmoids were seen to be travelling along the multi-threaded spire as well as at the footpoint region in the EUV/UV images recorded by the Atmospheric Imaging Assembly (AIA). The Differential Emission Measure (DEM) analysis was performed using EUV AIA images, and the high-temperature part of the DEM was constrained by combining X-ray images from the X-ray telescope (XRT/Hinode). We observed a systematic rise and fall in brightness, electron number densities and the peak temperatures of the spire plasmoid during its propagation along the jet. The plasmoids at the footpoint (FPs) (1.0-2.5 MK) and plasmoids at the spire (SPs) (1.0-2.24 MK) were found to have similar peak temperatures, whereas the FPs have higher DEM weighted temperatures (2.2-5.7 MK) than the SPs (1.3-3.0 MK). A lower limit to the electron number densities of plasmoids - SPs (FPs) were obtained that ranged between 3.4-6.1$\times$10$^{8}$ (3.3-5.9$\times$10$^{8}$) cm$^{-3}$ whereas for the spire, it ranged from 2.6-3.2$\times$10$^{8}$ cm$^{-3}$. Our analysis shows that the emission of these plasmoids starts close to the base of the jet(s), where we believe that a strong current interface is formed. This suggests that the blobs are plasmoids induced by a tearing-mode instability.

It is now well-established that $\Lambda$CDM predictions overestimate gravitational lensing measurements around massive galaxies by about 30%, the so-called lensing is low problem. Using a state-of-the-art hydrodynamical simulation, we show that this discrepancy reflects shortcomings in standard structure formation models rather than tensions within the $\Lambda$CDM paradigm itself. Specifically, this problem results from ignoring a variety of galaxy formation effects in simple models, including assembly bias, segregation of satellite galaxies relative to dark matter, and baryonic effects on the matter distribution. Each of these contributes towards overestimating gravitational lensing and, when combined, these explain the amplitude and scale dependence of the lensing is low problem. We conclude that simplistic structure formation models are inadequate to interpret lensing and clustering together, and that it is crucial to employ more sophisticated models for the upcoming generation of large-scale surveys.

The study of the polarization direction is crucial in the issue of restoring the spatial structure of the magnetic field in the active galaxy parsec-scale jets. But, due to relativistic effects, the magnetic field projected onto the celestial sphere in the source reference frame cannot be assumed to be orthogonal to the observed direction of the electric vector in the wave. Moreover, the local axis of the jet component may not coincide with its motion direction, which affects the observed polarization direction. In this article, we analyze the transverse to jet distributions of the electric vector in the wave, obtained as a result of modeling with different jet kinematic and geometrical parameters for a helical magnetic field with a different twist angle and for a toroidal magnetic field in the center, surrounded by a varying thickness sheath, penetrated by a poloidal field. We obtained: 1) the shape of the electric vector transverse distribution depends in a complex way on the angles of the jet axis and the velocity vector with the line of sight; 2) ambiguity in determining the twist direction of the helical magnetic field under using only the distributions of the electric vector in the wave; 3) both considered magnetic field topologies can reproduce both the ``spine-sheath'' polarization structure and individual bright details with the longitudinal to the jet axis polarization direction.

Rocky planet compositions regulate planetary evolution by affecting core sizes, mantle properties, and melting behaviours. Yet, quantitative treatments of this aspect of exoplanet studies remain generally under-explored. We attempt to constrain the range of potential bulk terrestrial exoplanet compositions in the solar neighbourhood (<200 pc). We circumscribe probable rocky exoplanet compositions based on a population analysis of stellar chemical abundances from the Hypatia and GALAH catalogues. We apply a devolatilization model to simulate compositions of hypothetical, terrestrial-type exoplanets in the habitable zones around Sun-like stars, considering elements O, S, Na, Si, Mg, Fe, Ni, Ca, and Al. We further apply core-mantle differentiation by assuming constant oxygen fugacity, and model the consequent mantle mineralogy with a Gibbs energy minimisation algorithm. We report statistics on several compositional parameters and propose a reference set of (21) representative planet compositions for using as end-member compositions in imminent modelling and experimental studies. We find a strong correlation between stellar Fe/Mg and metallic core sizes, which can vary from 18 to 35 wt%. Furthermore, stellar Mg/Si gives a first-order indication of mantle mineralogy, with high-Mg/Si stars leading to weaker, ferropericlase-rich mantles, and low-Mg/Si stars leading to mechanically stronger mantles. The element Na, which modulates crustal buoyancy and mantle clinopyroxene fraction, is affected by devolatilization the most. While we find that planetary mantles mostly consist of Fe/Mg-silicates, core sizes and relative abundances of common minerals can nevertheless vary significantly among exoplanets. These differences likely lead to different evolutionary pathways among rocky exoplanets in the solar neighbourhood.

Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends previous works to a dynamic range in stellar mass and star formation rate (10^-4<SFR<10 Mo yr^-1) never explored so far. The main sequence (MS) relation derived for all star forming galaxies within one virial radius of the Virgo cluster has a slope comparable to that observed in other nearby samples of isolated objects, but has a dispersion ~3 times larger. The dispersion is tightly connected to the available amount of HI gas, with gas-poor systems located far below objects of similar stellar mass but with a normal HI content. When measured on unperturbed galaxies with a normal HI gas content, the relation has a slope a=0.92, an intercept b=-1.57, and a scatter ~0.40. We compare these observational results to the prediction of models. The observed scatter in the MS relation can be reproduced only after a violent and active stripping process such as ram-pressure that removes gas from the disc and quenches star formation on short (<1 Gyr) timescales. This rules out milder processes such as starvation. This interpretation is also consistent with the position of galaxies of different star formation activity and gas content within the phase-space diagram. We also show that the star forming regions formed in the stripped material outside perturbed galaxies are located well above the MS relation drawn by unperturbed systems. These HII regions, which might be at the origin of compact sources typical in rich environments, are living a starburst phase lasting only <50 Myr, later becoming quiescent systems.

There are still some important unanswered questions about the unexpected very high energy $\gamma$-ray signatures. To help understand the mechanism, focusing on the linear and quadratic perturbation mode for subliminal case, the present paper revisited the expected signature for the Lorentz invariance violation effects on $\gamma-\gamma$ absorption in TeV spectra of Gamma-ray bursts (GRBs). We note that the existence of minimum photon energy threshold for the pair production process leads up to a break energy, which is sensitive to the quantum gravity energy scale. We suggest that energy spectral break in the few tens of TeV is a rough observational diagnostic for the LIV effects. The expected spectra characteristics are applied to a GRB 221009A. The results show that the cosmic opacity with Lorentz invariance violation effects considered here is able to roughly reproduce the observed $\gamma$-ray spectra for the source, which enabled us to constrain the lower limit of the linear values of energy scale at $E_{\rm QG,1}=3.35\times10^{20}$ GeV for the linear perturbation and $E_{\rm QG,2}=9.19\times10^{13}$ GeV for the quadratic perturbation. This value corresponds to a break energy $E_{\rm \gamma, break,1}\simeq 55.95~\rm TeV$ for the linear and $E_{\rm \gamma, break,2}\simeq 73.66~\rm TeV$ for the quadratic in the observed frame respectively.

Mass transfer stability is an essential issue in binary evolution. Ge et al. studied critical mass ratios for dynamically stable mass transfer by establishing adiabatic mass loss model and found that the donor stars on the giant branches tend to be more stable than that based on the composite polytropic stellar model. We would investigate the influence of mass transfer stability on the formation and properties of DWD populations. We performed a series of binary population synthesis, where the critical mass ratios from adiabatic mass loss model (Ge's model) and that from the composite polytropic model are adopted, respectively. For Ge's model, most of the DWDs are produced from the stable non-conservative Roche lobe overflow plus common envelope (CE) ejection channel (RL+CE channel) regardless of the CE ejection efficiency $\alpha_{CE}$. While the results of the polytropic model strongly depend on the adopted value of $\alpha_{ CE}$. We find DWDs produced from the RL+CE channel have comparable WD masses and the mass ratio distribution peaks at around 1. Based on the magnitude-limited sample of DWDs, the space densities for the detectable DWDs and those with extremely low-mass WD (ELM WD) companions in Ge's model is $1347$ and $473 kpc^{-3}$, respectively, close to observations. While the polytropic model overpredicts space density of DWDs by a factor of about $2-3$. We also find that the results of DWD merger rate distribution in Ge's model reproduce the observations better than that of the polytropic model, and the merger rate of DWDs with ELM WD companions in the Galaxy is about $1.8\times 10^{-3} yr^{-1}$ in Ge's model, which is comparable to the observation estimation. We confirm that the mass transfer stability plays important roles in the formation and properties of DWD populations, and then in the progenitors of SNe Ia and detectable GW sources.

We report the first detection of cesium (Z = 55) in the atmosphere of a white dwarf. Around a dozen absorption lines of Cs IV, Cs V, and Cs VI have been identified in the Far Ultraviolet Spectroscopic Explorer spectrum of the He-rich white dwarf HD 149499B (Teff = 49,500 K, log g = 7.97). The lines have equivalent widths ranging from 2.3 to 26.9 m\r{A}. We performed a spectral synthesis analysis to determine the cesium content in the atmosphere. Non-LTE atmosphere models were computed by considering cesium explicitly in the calculations. For this purpose we calculated oscillator strengths for the bound-bound transitions of Cs IV-Cs VI with both AUTOSTRUCTURE (multiconfiguration Breit-Pauli) and GRASP2K (multiconfiguration Dirac-Fock) atomic structure codes as neither measured nor theoretical values are reported in the literature. We determined a cesium abundance of log N(Cs)/N(He) = -5.45(0.35), which can also be expressed in terms of the mass fraction log X(Cs) = -3.95(0.35).

We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly during the collapse, large grain formation being controlled by the turbulent differential velocity. When turbulence is included, only ambipolar diffusion remains efficient at removing the small grains from the distribution, brownian motion is only efficient as a standalone process. The macroscopic gas-dust drift is negligible for grain growth and only dynamically significant near the first Larson core. At high density, we find that the coagulated distribution is unaffected by the initial choice of dust distribution. Strong magnetic fields are found to enhance the small grains depletion, causing an important increase of the ambipolar diffusion. This hints that the magnetic field strength could be regulated by the small grain population during the protostellar collapse. Fragmentation could be effective for bare silicates, but its modeling relies on the choice of ill-constrained parameters. It is also found to be negligible for icy grains. When fragmentation occurs, it strongly affects the magnetic resistivities profiles. Dust coagulation is a critical process that needs to be fully taken into account during the protostellar collapse. The onset and feedback of fragmentation remains uncertain and its modeling should be further investigated.

It is well known that blazars can show variability on a wide range of time scales. This behavior can include periodicity in their $\gamma$-ray emission, whose clear detection remains an ongoing challenge, partly due to the inherent stochasticity of the processes involved and also the lack of adequately-well sampled light curves. We report on a systematic search for periodicity in a sample of 24 blazars, using twelve well-established methods applied to Fermi-LAT data. The sample comprises the most promising candidates selected from a previous study, extending the light curves from nine to twelve years and broadening the energy range analyzed from $>$1 GeV to $>$0.1 GeV. These improvements allow us to build a sample of blazars that display a period detected with global significance $\gtrsim3\,\sigma$. Specifically, these sources are PKS 0454$-$234, S5 0716+714, OJ 014, PG 1553+113, and PKS 2155$-$304. Periodic $\gamma$-ray emission may be an indication of modulation of the jet power, particle composition, or geometry but most likely originates in the accretion disk, possibly indicating the presence of a supermassive black hole binary system.

We present updated results constraining multiplicity demographics for the stellar population of the Orion Nebula Cluster (ONC, a high-mass, high-density star-forming region), across primary masses 0.08-0.7M$_{\odot}$. Our study utilizes archival Hubble Space Telescope data obtained with the Advanced Camera for Surveys using multiple filters (GO-10246). Previous multiplicity surveys in low-mass, low-density associations like Taurus identify an excess of companions to low-mass stars roughly twice that of the Galactic field and find the mass ratio distribution consistent with the field. Previously, we found the companion frequency to low-mass stars in the ONC is consistent with the Galactic field over mass ratios=0.6-1.0 and projected separations=30-160au, without placing constraints on the mass ratio distribution. In this study, we investigate the companion population of the ONC with a double point-spread function (PSF) fitting algorithm sensitive to separations larger than 10au (0.025") using empirical PSF models. We identified 44 companions (14 new), and with a Bayesian analysis, estimate the companion frequency to low-mass stars in the ONC =0.13$^{+0.05}_{-0.03}$ and the power law fit index to the mass ratio distribution =2.08$^{+1.03}_{-0.85}$ over all mass ratios and projected separations of 10-200au. We find the companion frequency in the ONC is consistent with the Galactic field population, likely from high transient stellar density states, and a probability of 0.002 that it is consistent with that of Taurus. We also find the ONC mass ratio distribution is consistent with the field and Taurus, potentially indicative of its primordial nature, a direct outcome of the star formation process.

As a consequence of galaxy clustering, close galaxies observed on the plane of the sky should be spatially correlated with a probability inversely proportional to their angular separation. In principle, this information can be used to improve photometric redshift estimates when spectroscopic redshifts are available for some of the neighbouring objects. Depending on the depth of the survey, however, such angular correlation is reduced by chance projections. In this work, we implement a deep learning model to distinguish between apparent and real angular neighbours by solving a classification task. We adopt a graph neural network architecture to tie together the photometry, the spectroscopy and the spatial information between neighbouring galaxies. We train and validate the algorithm on the data of the VIPERS galaxy survey, for which SED-fitting based photometric redshifts are also available. The model yields a confidence level for a pair of galaxies to be real angular neighbours, enabling us to disentangle chance superpositions in a probabilistic way. When objects for which no physical companion can be identified are excluded, all photometric redshifts' quality metrics improve significantly, confirming that their estimates were of lower quality. For our typical test configuration, the algorithm identifies a subset containing ~75% of high-quality photometric redshifts, for which the dispersion is reduced by as much as 50% (from 0.08 to 0.04), while the fraction of outliers reduces from 3% to 0.8%. Moreover, we show that the spectroscopic redshift of the angular neighbour with the highest detection probability provides an excellent estimate of the redshift of the target galaxy, comparable or even better than the corresponding template fitting estimate.

Flares and coronal mass ejections are powered by magnetic energy stored in coronal electric currents. Here, we explore the nature of coronal currents in observed and model active region (ARs) by studying manifestations of these currents in photospheric vector magnetograms. We employ Gauss's separation method, recently introduced to the solar physics literature, to partition the photospheric field into three distinct components, each arising from a separate source: (i) currents passing through the photosphere, (ii) currents flowing below it, and (iii) currents flowing above it. We refer to component (iii) as the photospheric imprint of coronal currents. In both AR 10930 and AR 11158, photospheric imprints exhibit large-scale, spatially coherent structures along these regions' central, sheared polarity inversion lines (PILs) that are consistent with coronal currents flowing horizontally above these PILs, similar to recent findings in AR 12673 by Schuck et al. (2022). We find similar photospheric imprints in a simple model of a non-potential AR with known currents. We find that flare-associated changes in photospheric imprints in AR 11158 accord with earlier reports that near-PIL fields become more horizontal, consistent with the "implosion" scenario. We hypothesize that this evolution effectively shortens, in an overall sense, current-carrying coronal fields, leading to decreased inductive energy (DIE) in the coronal field. We further hypothesize that, in the hours prior to flares, parts of the coronal field slowly expand, in a process we deem coronal inflation (CI) -- essentially, the inverse of the implosion process. Both of these hypotheses are testable with non-potential coronal field extrapolations.

We investigate the connection of the regulation of star formation and the cycling of baryons within and in and out of galaxies. We use idealized numerical simulations of Milky Way-mass galaxies, in which we systemically vary the galaxy morphology (bulge-to-total mass ratio) and stellar feedback strength (total eight setups with 80 simulations). By following individual gas parcels through the disk, spiral arms, and massive star-forming clumps, we quantify how gas moves and oscillates through the different phases of the interstellar medium (ISM) and forms stars. We show that the residence time of gas in the dense ISM phase ($\tau_{\rm SF}$), the nature of spiral arms (strength, number), and the clump properties (number, mass function, and young star fraction) depend on both the galaxy morphology and stellar feedback. Based on these results, we quantify signatures of the baryon cycle within galaxies using the temporal and spatial power spectrum density (PSD) of the star formation history (SFH). Stronger stellar feedback leads to more bursty star formation while the correlation timescale of the SFH is longer, because stronger feedback dissolves the dense, star-forming ISM phase, leading to a more homogeneous ISM and a decrease in $\tau_{\rm SF}$. The bulge strength has a similar effect: the deep gravitational potential in a bulge-dominant galaxy imposes a strong shear force that effectively breaks apart gas clumps in the ISM; this subsequently inhibits the fragmentation of cool gas and therefore the star formation in the disk, leading to a decrease in the spatial power on scales of $\sim$ 1 kpc. We conclude that measurements of the temporal and spatial PSD of the SFH can provide constraints on the baryon cycle and the star formation process.

Smoothed Particle Hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron-star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches which can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

Stars are born in dense molecular filaments irrespective of their mass. Compression of the ISM by shocks cause filament formation in molecular clouds. Observations show that a massive star cluster formation occurs where the peak of gas column density in a cloud exceeds 10^23 cm^-2. In this study, we investigate the effect of the shock-compressed layer duration on filament/star formation and how the initial conditions of massive star formation are realized by performing three-dimensional (3D) isothermal magnetohydrodynamics (MHD) simulations with {gas inflow duration from the boundaries (i.e., shock wave duration)} as a controlling parameter. Filaments formed behind the shock expand after the duration time for short shock duration models, whereas long duration models lead to star formation by forming massive supercritical filaments. Moreover, when the shock duration is longer than two postshock free-fall times, the peak column density of the compressed layer exceeds 10^23 cm^-2, and {the gravitational collapse of the layer causes that} the number of OB stars expected to be formed in the shock-compressed layer reaches the order of ten (i.e., massive cluster formation).

A large number of studies, all using Bayesian parameter inference from Markov Chain Monte Carlo methods, have constrained the presence of a decaying dark matter component. All such studies find a strong preference for either very long-lived or very short-lived dark matter. However, in this letter, we demonstrate that this preference is due to parameter volume effects that drive the model towards the standard $\Lambda$CDM model, which is known to provide a good fit to most observational data. Using profile likelihoods, which are free from volume effects, we instead find that the best-fitting parameters are associated with an intermediate regime where around $3 \%$ of cold dark matter decays just prior to recombination. With two additional parameters, the model yields an overall preference over the $\Lambda$CDM model of $\Delta \chi^2 \approx -2.8$ with \textit{Planck} and BAO and $\Delta \chi^2 \approx -7.8$ with the SH0ES $H_0$ measurement, while only slightly alleviating the $H_0$ tension. Ultimately, our results reveal that decaying dark matter is more viable than previously assumed, and illustrate the dangers of relying exclusively on Bayesian parameter inference when analysing extensions to the $\Lambda$CDM model.

The observations of optical and near-infrared counterparts of binary neutron star mergers not only enrich our knowledge about the abundance of heavy elements in the Universe, or help reveal the remnant object just after the merger as generally known, but also can effectively constrain dense nuclear matter properties and the equation of state (EOS) in the interior of the merging stars. Following the relativistic mean-field description of nuclear matter, we perform the Bayesian inference of the EOS and the nuclear matter properties using the first multi-messenger event GW170817/AT2017gfo, together with the NICER mass-radius measurements of pulsars. The kilonova is described by a radiation-transfer model with the dynamical ejecta, and light curves connect with the EOS through the quasi-universal relations between the ejecta properties (the ejected mass, velocity, opacity or electron fraction) and binary parameters (the mass ratio and reduced tidal deformability). It is found that the posterior distributions of the reduced tidal deformability from the AT2017gfo analysis display a bimodal structure, with the first peak enhanced by the GW170817 data, leading to slightly softened posterior EOSs, while the second peak cannot be achieved by a nuclear EOS with saturation properties in their empirical ranges. The inclusion of NICER data in our analyses results in stiffened EOS posterior because of the massive pulsar PSR J0740+6620. We give results at nuclear saturation density for the nuclear incompressibility, the symmetry energy and its slope, as well as the nucleon effective mass, from our analysis of those observational data.

The extraordinarily bright gamma-ray burst GRB 221009A was observed by a large number of observatories, from radio frequencies to gamma-rays. Of particular interest are the reported observations of photon-like air showers of very high energy: an 18 TeV event in LHAASO and a 251 TeV event at Carpet-2. Gamma rays at these energies are expected to be absorbed by pair-production events on background photons when travelling intergalactic distances. Several works have sought to explain the observations of these events, assuming they originate from GRB 221009A, by invoking axion-like particles (ALPs). We reconsider this scenario and account for astrophysical uncertainties due to poorly known magnetic fields and background photon densities. We find that, robustly, the ALP scenario cannot simultaneously account for an 18 TeV and a 251 TeV photon from GRB 221009A.

The Scalar Field Dark Matter (SFDM) model, also called Fuzzy, Wave, Bose-Einstein, Ultra-light Dark Matter, has received a lot of attention because it has been able to provide simpler and more natural explanations for various features of galaxies, such as the number of satellite galaxies and the cusp-core problem. We recently showed that this model is able to explain the vast polar orbits of satellite galaxies around their host, the so-called VPO, and to explain the X-ray and gamma-ray emissions in the vacuum regions of our galaxy, that is, the Fermi Bubbles. In all these phenomena the quantum character of SFDM has been crucial. In this work we study the quantum effects of SFDM at the cosmological level, to see these effects not only at the galactic scale, but also at the cosmological scale. Using a convenient ansatz, we were able to integrate the perturbed equations to show that the shape of the SFDM halos resembling atoms is a generic result. The main conclusion of this work is that quantum mechanics, the successful microworld theory, could also explain the dark side of the cosmos.

Common envelope (CE) evolution, which is crucial in creating short period binaries and associated astrophysical events, can be constrained by reverse modeling of such binaries' formation histories. Through analysis of a sample of well-constrained white dwarf (WD) binaries with low-mass primaries (7 eclipsing double WDs, 2 non-eclipsing double WDs, 1 WD-brown dwarf), we estimate the CE energy efficiency $\alpha_{\rm{CE}}$ needed to unbind the hydrogen envelope. We use grids of He- and CO-core WD models to determine the masses and cooling ages that match each primary WD's radius and temperature. Assuming gravitational wave-driven orbital decay, we then calculate the associated ranges in post-CE orbital period. By mapping WD models to a grid of red giant progenitor stars, we determine the total envelope binding energies and possible orbital periods at the point CE evolution is initiated, thereby constraining $\alpha_{\rm CE}$. Assuming He-core WDs with progenitors of 0.9 - 2.0 $M_\odot$, we find $\alpha_{\rm CE} \! \sim \! 0.2-0.4$ is consistent with each system we model. Significantly higher values of $\alpha_{\rm{CE}}$ are required for higher mass progenitors and for CO-core WDs, so these scenarios are deemed unlikely. Our values are mostly consistent with previous studies of post-CE WD binaries, and they suggest a nearly constant and low envelope ejection efficiency for CE events that produce He-core WDs.

We utilise theoretical models of Population III stellar+nebular spectra to investigate the prospects of observing and accurately identifying Population III galaxies with JWST using both deep imaging and spectroscopy. We investigate a series of different colour cuts, finding that a combination of NIRCam and MIRI photometry through the F444W-F560W, F560W-F770W colours offers the most robust identifier of potential $z=8$ Pop III candidates. We calculate that NIRCam will have to reach $\sim$28.5-30.0 AB mag depths (1-20 h), and MIRI F560W must reach $\sim$27.5-29.0 AB mag depths (10-100 h) to achieve $5\sigma$ continuum detections of $M_* = 10^6~\mathrm{M}_\odot$ Pop III galaxies at $z=8$. We also discuss the prospects of identifying Pop III candidates through slitless and NIRSpec spectroscopic surveys that target Ly$\alpha$, H$\beta$ and/or He II $\lambda 1640$. We find small differences in the H$\beta$ rest-frame equivalent width (EW) between Pop III and non-Pop III galaxies, rendering this diagnostic likely impractical. We find that the detection of high EW He II $\lambda 1640$ emission will serve as the definitive Pop III identifier, requiring (ultra-)deep integrations (10-250 h) with NIRSpec/G140M for $M_*=10^6~\mathrm{M}_\odot$ Pop III galaxies at $z=8$. With moderate ($\mu=$2-3) lensing and/or moderately massive ($M_*= 2$-$3\times10^6~\mathrm{M}_\odot$) Pop III galaxies, such line detections can be achieved in medium-sized JWST GO programs. However, MIRI F770W detections of Pop III galaxies will require substantial gravitational lensing ($\mu=10$) and/or fortuitous imaging of exceptionally massive ($M_* = 10^7~\mathrm{M}_\odot$) Pop III galaxies. Thus, NIRCam medium-band imaging surveys that can search for high EW He II $\lambda 1640$ emitters in photometry may perhaps be a viable alternative for finding Pop III candidates.

The upcoming Square Kilometer Array (SKA) is expected to produce humongous amount of data for undertaking H{\sc i}~science. We have developed an MPI-based {\sc Python} pipeline to deal with the large data efficiently with the present computational resources. Our pipeline divides such large H{\sc i}~21-cm spectral cubes into several small cubelets, and then processes them in parallel using publicly available H{\sc i}~source finder {\sc SoFiA-$2$}. The pipeline also takes care of sources at the boundaries of the cubelets and also filters out false and redundant detections. By comapring with the true source catalog, we find that the detection efficiency depends on the {\sc SoFiA-$2$} parameters such as the smoothing kernel size, linking length and threshold values. We find the optimal kernel size for all flux bins to be between $3$ to $5$ pixels and $7$ to $15$ pixels, respectively in the spatial and frequency directions. Comparing the recovered source parameters with the original values, we find that the output of {\sc SoFiA-$2$} is highly dependent on kernel sizes and a single choice of kernel is not sufficient for all types of H{\sc i}~galaxies. We also propose use of alternative methods to {\sc SoFiA-$2$} which can be used in our pipeline to find sources more robustly.

We present the GPry algorithm for fast Bayesian inference of general (non-Gaussian) posteriors with a moderate number of parameters. GPry does not need any pre-training, special hardware such as GPUs, and is intended as a drop-in replacement for traditional Monte Carlo methods for Bayesian inference. Our algorithm is based on generating a Gaussian Process surrogate model of the log-posterior, aided by a Support Vector Machine classifier that excludes extreme or non-finite values. An active learning scheme allows us to reduce the number of required posterior evaluations by two orders of magnitude compared to traditional Monte Carlo inference. Our algorithm allows for parallel evaluations of the posterior at optimal locations, further reducing wall-clock times. We significantly improve performance using properties of the posterior in our active learning scheme and for the definition of the GP prior. In particular we account for the expected dynamical range of the posterior in different dimensionalities. We test our model against a number of synthetic and cosmological examples. GPry outperforms traditional Monte Carlo methods when the evaluation time of the likelihood (or the calculation of theoretical observables) is of the order of seconds; for evaluation times of over a minute it can perform inference in days that would take months using traditional methods. GPry is distributed as an open source Python package (pip install gpry) and can also be found at https://github.com/jonaselgammal/GPry.

We generalize to reduced Horndeski theories of gravity, where gravitational waves (GWs) travel at the speed of light, the expression of a statistically homogeneous and unpolarized stochastic gravitational wave background (SGWB) signal measured as the correlation between the individual signals detected by two interferometers in arbitrary configurations. We also discuss some results found in the literature regarding cosmological distances in modified theories, namely, the simultaneous validity of a duality distance relation for GW signals and of the coincidence between the gravitational wave luminosity distance, based on the energy flux, and the distance inferred from the wave amplitude. This discussion allows us to conclude that the spectral energy density per unit solid angle of an astrophysical SGWB signal has the same functional dependency with the luminosity of each emitting source as in General Relativity (GR). Using the generalized expression of the GW energy-momentum tensor and the modified propagation law for the tensor modes, we conclude that the energy density of a SGWB maintains the same functional relation with the scale factor as in GR, provided that the modified theory coincides with GR in a given hypersurface of constant time. However, the relation between the detected signal and the spectral energy density is changed by the global factor $G_4(\varphi(t_0))$, thus potentially serving as a probe for modified gravity theories.

Traditionally, gravitational waves are detected with techniques such as matched filtering or unmodeled searches based on wavelets. However, in the case of generic black hole binaries with non-aligned spins, if one wants to explore the whole parameter space, matched filtering can become impractical, which sets severe restrictions on the sensitivity and computational efficiency of gravitational-wave searches. Here, we use a novel combination of machine-learning algorithms and arrive at sensitive distances that surpass traditional techniques in a specific setting. Moreover, the computational cost is only a small fraction of the computational cost of matched filtering. The main ingredients are a 54-layer deep residual network (ResNet), a Deep Adaptive Input Normalization (DAIN), a dynamic dataset augmentation, and curriculum learning, based on an empirical relation for the signal-to-noise ratio. We compare the algorithm's sensitivity with two traditional algorithms on a dataset consisting of a large number of injected waveforms of non-aligned binary black hole mergers in real LIGO O3a noise samples. Our machine-learning algorithm can be used in upcoming rapid online searches of gravitational-wave events in a sizeable portion of the astrophysically interesting parameter space. We make our code, AResGW, and detailed results publicly available at https://github.com/vivinousi/gw-detection-deep-learning.

Superstring/M-theory is the theory of quantum gravity that can provide the UV-completion to viable inflation models. We modify the Starobinsky inflation model by adding the Bel-Robinson tensor $T^{\mu\nu\lambda\rho}$ squared term proposed as the leading quantum correction inspired by superstring theory. The $(R+\frac{1}{6m^2}R^2 -\frac{\beta}{8m^6}T^2)$ model under consideration has two parameters: the inflaton mass $m$ and the string-inspired positive parameter $\beta$. We derive the equations of motion in the Friedmann-Lemaitre-Robertson-Walker universe and investigate its solutions. We find the physical bounds on the value of the parameter $\beta$ by demanding the absence of ghosts and consistency of the derived inflationary observables with the measurements of the cosmic microwave background radiation.

We consider the propagation of gravitational waves in the late-time Universe, in the presence of structure, and we also consider the cosmic fluid to be viscous. In this work, we investigate the cumulative effect of inhomogeneities and viscosity of the cosmic-fluid, on the observables associated with the sources of the gravitational waves. Employing Buchert's averaging procedure in the backreaction framework, we consider a model of inhomogeneous spacetime. Using the modified redshift versus distance relation, through the averaging process in the context of the model, we study the variation of the redshift-dependent part of the observed gravitational wave amplitude for different combinations of our model parameters, while simultaneously considering damping of the gravitational wave amplitude due to viscosity of the cosmic-fluid. Then, we investigate the differences occurring in the variation of the redshift-dependent part of the observed gravitational wave amplitude, due to consideration of viscous attenuation. We show that there are significant deviations after the inclusion of viscous attenuation in our analysis, depending on the chosen value of the coefficient of viscosity. Our result signifies the importance of the effect of viscosity, within the model of an inhomogeneous Universe, on precision measurements of parameters of compact-binary sources of gravitational waves.

Cosmological constraints on the sum of the neutrino masses can be relaxed if the number density of active neutrinos is reduced compared to the standard scenario, while at the same time keeping the effective number of neutrino species $N_{\rm eff}\approx 3$ by introducing a new component of dark radiation. We discuss a UV complete model to realise this idea, which simultaneously provides neutrino masses via the seesaw mechanism. It is based on a $U(1)$ symmetry in the dark sector, which can be either gauged or global. In addition to heavy seesaw neutrinos, we need to introduce $\mathcal{O}(10)$ generations of massless sterile neutrinos providing the dark radiation. Then we can accommodate active neutrino masses with $\sum m_\nu \sim 1$ eV, in the sensitivity range of the KATRIN experiment. We discuss the phenomenology of the model and identify the allowed parameter space. We argue that the gauged version of the model is preferred, and in this case the typical energy scale of the model is in the 10 MeV to few GeV range.

Scalarization is a mechanism that endows strongly self-gravitating bodies, such as neutron stars and black holes, with a scalar field configuration. It resembles a phase transition in that the scalar configuration only appears when a certain quantity that characterizes the compact object, e.g., its compactness or spin, is beyond a threshold. We provide a critical and comprehensive review of scalarization, including the mechanism itself, theories that exhibit it, its manifestation in neutron stars, black holes, and their binaries, potential extension to other fields, and a thorough discussion of future perspectives.

The International Astronomical Youth Camp (IAYC) is a 50-year old summer camp, where participants work independently on astronomy projects. Due to the ongoing COVID-19 pandemic, the 2020 and 2021 instalments of the IAYC were cancelled, a first in the camp's history. An online format was established dubbed the eIAYC, consisting of three types of activities: (1) an engagement series with astronomical talks and workshops; (2) small independent research projects; and (3) a non-astronomical program involving a range of social activities. Here we present the experience of adapting an in-person camp into an online alternative in order to further the IAYC's mission. We discuss organisational challenges, experiences with online engagement, and how the 2020 eIAYC informed plans for this year's eIAYC.

Objects orbiting in the presence of a rotating massive body experience a gravitomagnetic frame-dragging effect, known as the Lense-Thirring effect, that has been experimentally confirmed in the weak-field limit. In the strong-field limit, near the horizon of a rotating black hole, frame dragging becomes so extreme that all objects must co-rotate with the black hole's angular momentum. In this work, we perform general relativistic numerical simulations to identify observable signatures of frame dragging in the strong-field limit that appear when infalling gas is forced to flip its direction of rotation as it is being accreted. In total intensity images, infalling streams exhibit "S"-shaped features due to the switch in the tangential velocity. In linear polarization, a flip in the handedness of spatially resolved polarization ticks as a function of radius encodes a transition in the magnetic field geometry that occurs due to magnetic flux freezing in the dragged plasma. Using a network of telescopes around the world, the Event Horizon Telescope collaboration has demonstrated that it is now possible to directly image black holes on event horizon scales. We show that the phenomena described in this work would be accessible to the next-generation Event Horizon Telescope (ngEHT) and extensions of the array into space, which would produce spatially resolved images on event horizon scales with higher spatial resolution and dynamic range.

Single-photon atom gradiometry is a powerful experimental technique that can be employed to search for the oscillation of atomic transition energies induced by ultralight scalar dark matter (ULDM). In the sub-Hz regime the background is expected to be dominated by gravity gradient noise (GGN), which arises as a result of mass fluctuations around the experiment. In this work we model the GGN as surface Rayleigh waves and construct a likelihood-based analysis that consistently folds GGN into the sensitivity estimates of vertical atom gradiometers in the frequency window between 1 mHz and 1 Hz. We show that in certain geological settings GGN can be significantly mitigated when operating a multi-gradiometer configuration, which consists of three or more atom interferometers in the same baseline. Multi-gradiometer experiments, such as future versions of AION and MAGIS-100, have the potential to probe regions of scalar ULDM parameter space in the sub-Hz regime that have not been excluded by existing experiments.

Since the first detection of gravitational waves by the LIGO/VIRGO team, the related research field has attracted more attention. The spinning compact binaries system, as one of the gravitational-wave sources for broadband laser interferometers, has been widely studied by related researchers. In order to analyze the gravitational wave signals using matched filtering techniques, reliable numerical algorithms are needed. Spinning compact binaries systems in Post-Newtonian (PN) celestial mechanics have an inseparable Hamiltonian. The extended phase-space algorithm is an effective solution for the problem of this system. We have developed correction maps for the extended phase-space method in our previous work, which significantly improves the accuracy and stability of the method with only a momentum scale factor. In this paper, we will add more scale factors to modify the numerical solution in order to minimize the errors in the constants of motion. However, we find that these correction maps will result in a large energy bias in the subterms of the Hamiltonian in chaotic orbits, whose potential and kinetic energy, etc. are calculated inaccurately. We develop new correction maps to reduce the energy bias of the subterms of the Hamiltonian, which can instead improve the accuracy of the numerical solution and also provides a new idea for the application of the manifold correction in other algorithms.

Ultra-light particles, such as axions, form a macroscopic condensate around a highly spinning black hole by the superradiant instability. Due to its macroscopic nature, the condensate opens the possibility of detecting the axion through gravitational wave observations. However, the precise evolution of the condensate must be known for the actual detection. For future observation, we numerically study the influence of the self-interaction, especially interaction between different modes, on the evolution of the condensate in detail. First, we focus on the case when condensate starts with the smallest possible angular quantum number. For this case, we perform the non-linear calculation and show that the dissipation induced by the mode interaction is strong enough to saturate the superradiant instability, even if the secondary cloud starts with quantum fluctuations. Our result indicates that explosive phenomena such as bosenova do not occur in this case. We also show that the condensate settles to a quasi-stationary state mainly composed of two modes, one with the smallest angular quantum number for which the superradiant instability occurs and the other with the adjacent higher angular quantum number. We also study the case when the condensate starts with the dominance of the higher angular quantum number. We show that the dissipation process induced by the mode coupling does not occur for small gravitational coupling. Therefore, bosenova might occur in this case.

In the future situation, aiming to seek more resources, human beings decided to march towards the mysterious and bright starry sky, which opened the era of great interstellar exploration. According to the Outer Space Treaty, any exploration of celestial bodies should be aimed at promoting global equality and for the benefit of all nations. Firstly, we defined global equity and set a Unified Equity Index (UEI) model to measure it. We merge the factors with greater correlation, and finally, get 6 elements, and then use the entropy method (TEM) to find the dispersion of these elements in different countries. Then use principal component analysis (PCA) to reduce the dimensionality of the dispersion, and then use the scandalized index to obtain the global equity. Secondly, we simulated a future with asteroid mining and evaluated its impact on Unified Equity Index (UEI). Then, we divided the mineable asteroids into three classes with different mining difficulties and values, identified 28 mining entities including private companies, national and international organizations. We considered changes in the asteroid classes, mining capabilities and mining scales to determine the changes in the value of minerals mined between 2025 and 2085. We convert mining output value into mineral transaction value through allocation matrix. Based on grey relational analysis (GRA). Finally, we presented three possible versions of the future of asteroid mining by changing the conditions. We propose two sets of corresponding policies for changes in future trends in global fairness with asteroid mining. We test the separate and combined effects of these policies and find that they are positive, strongly supporting the effectiveness of our model.

LHAASO collaboration detected photons with energy above 10 TeV from the most recent gamma-ray burst (GRB), GRB221009A. Given the redshift of this event, $z\sim 0.15$, photons of such energy are expected to interact with the diffuse extragalactic background light (EBL) well before reaching Earth. In this paper we provide the first neutrino-related explanation of the most energetic 18 TeV event reported by LHAASO. We find that the minimal viable scenario involves both mixing and transition magnetic moment portal between light and sterile neutrinos. The production of sterile neutrinos occurs efficiently via mixing while the transition magnetic moment portal governs the decay rate in the parameter space where tree-level decays via mixing to non-photon final states are suppressed. Our explanation of this event, while being consistent with the terrestrial constraints, points to the non-standard cosmology.