Papers for Friday, Mar 26 2021

The efficiency of the transport of angular momentum and chemical elements inside intermediate-mass stars lacks proper calibration, thereby introducing uncertainties on a star's evolutionary pathway. Improvements require better estimation of stellar masses, evolutionary stages, and internal mixing properties. We aim to develop a neural network approach for asteroseismic modelling and test its capacity to provide stellar masses, ages, and overshooting parameter for a sample of 37 $\gamma$ Doradus stars. Here, our goal is to perform the parameter estimation from modelling of individual periods measured for dipole modes with consecutive radial order. We have trained neural networks to predict theoretical pulsation periods of high-order gravity modes as well as the luminosity, effective temperature, and surface gravity for a given mass, age, overshooting parameter, diffusive envelope mixing, metallicity, and near-core rotation frequency. We have applied our neural networks for Computing Pulsation Periods and Photospheric Observables, C-3PO, to our sample and compute grids of stellar pulsation models for the estimated parameters. We present the near-core rotation rates (from literature) as a function of the inferred stellar age and critical rotation rate. We assess the rotation rates of the sample near the start of the main sequence assuming rigid rotation. Furthermore, we measure the extent of the core overshoot region and find no correlation with mass, age, or rotation. The neural network approach developed in this study allows for the derivation of stellar properties dominant for stellar evolution -- such as mass, age, and extent of core-boundary mixing. It also opens a path for future estimation of mixing profiles throughout the radiative envelope, with the aim to infer those profiles for large samples of $\gamma$ Doradus stars.

The spatial decorrelation of dense molecular gas and young stars observed on $\lesssim 1$ kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC300, and we show that this agreement is due to different aspects of feedback, including H$_2$ dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affect the correlation on $\lesssim 100$ pc scales, while supernovae play a significant role on $\gtrsim 100$ pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per freefall time, $\epsilon_{\rm ff}$, which provides a strong observational constraint on $\epsilon_{\rm ff}$ when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.

The existence of supermassive black holes (SMBHs) with masses greater than $\sim 10^{9}M_{\odot}$ at high redshift ($z\gtrsim 7$) is difficult to be accommodated in standard astrophysical scenarios. We study the possibility that (nearly) totally dissipative self-interacting dark matter (tdSIDM)--in rare, high density dark matter fluctuations in the early Universe--produces SMBH seeds through catastrophic collapse. We use a semi-analytic model, tested and calibrated by a series of N-body simulations of isolated dark matter halos, to compute the collapse criteria and timescale of tdSIDM halos, where dark matter loses nearly all of its kinetic energy in a single collision in the center-of-momentum frame. Applying this model to halo merger trees, we empirically assign SMBH seeds to halos and trace the formation and evolution history of SMBHs. We make predictions for the quasar luminosity function, the $M_{\rm BH}$-$\sigma_{\rm v}^{\ast}$ relation, and cosmic SMBH mass density at high redshift and compare them to observations. We find that a dissipative dark matter interaction cross-section of $\sigma/m \sim 0.05~\rm cm^2/g$ is sufficient to produce the SMBHs observed in the early Universe while remaining consistent with ordinary SMBHs in the late Universe.

We apply for the first time a two-dimensional fitting statistic, $\tau^2$, to rotational evolution models (REMs) of stars (0.1 to 1.3 $M_{\odot}$) on the period-mass plane. The $\tau^2$ statistic simultaneously considers all cluster rotation data to return a goodness of fit, allowing for data-driven improvement of REMs. We construct data sets for Upper Sco, the Pleiades and Praesepe, to which we tune our REMs. We use consistently determined stellar masses (calculated by matching $K_\textrm{s}$ magnitudes to isochrones) and literature rotation periods. As a first demonstration of the $\tau^2$ statistic, we find the best-fitting gyrochronology age for Praesepe, which is in good agreement with the literature. We then systematically vary three parameters which determine the dependence of our stellar wind torque law on Rossby number in the saturated and unsaturated regimes, and the location of the transition between the two. By minimising $\tau^2$, we find best-fit values for each parameter. These values vary slightly between clusters, mass determinations and initial conditions, highlighting the precision of $\tau^2$ and its potential for constraining REMs, gyrochronology, and understanding of stellar physics. Our resulting REMs, which implement the best-possible fitting form of a broken power-law torque, are statistically improved on previous REMs using similar formulations, but still do not simultaneously describe the observed rotation distributions of the lowest masses, which have both slow and fast rotators by the Praesepe age, and the shape of the converged sequence for higher masses. Further complexity in the REMs is thus required to accurately describe the data.

The W-CDF-S and ELAIS-S1 fields will be two of the LSST Deep Drilling fields, but the availability of spectroscopic redshifts within these two fields is still limited on deg^2 scales. To prepare for future science, we use EAZY to estimate photometric redshifts (photo-zs) in these two fields based on forced-photometry catalogs. Our photo-z catalog consists of ~0.8 million sources covering 4.9 deg^2 in W-CDF-S and ~0.8 million sources covering 3.4 deg^2 in ELAIS-S1, among which there are ~0.6 (W-CDF-S) and ~0.4 (ELAIS-S1) million sources having signal-to-noise-ratio (SNR) > 5 detections in more than 5 bands. By comparing photo-zs and available spectroscopic redshifts, we demonstrate the general reliability of our photo-z measurements. Our photo-z catalog is publicly available at \doi{10.5281/zenodo.4603178}.

In this paper, we consider three phantom dark energy models, in the context of interaction between the dark components namely cold dark matter (CDM) and dark energy (DE). The first model, known as $w_{\textrm{d}}$CDM can induce a big rip singularity (BR) while the two remaining induce future abrupt events known as the Little Rip (LR) and Little Sibling of the Big Rip (LSBR). These phantom DE models can be distinguished by their equation of state. We invoke a new phenomenon such as the interaction between CDM and DE given that it could solve or alleviate some of the problems encountered in standard cosmology. We aim to find out the effect of such an interaction on the cosmological parameters of the studied models, as well as, the persistence or the disappearance of the singularity and the abrupt events induced by the models under study. We choose an interaction term proportional to DE density, i.e. $Q=\lambda H \rho_{\textrm{d}}$, since the case where $Q\propto \rho_{\textrm{m}}$ could lead to a large scale instability at early time. We also do not claim at all that $Q=\lambda H \rho_{\textrm{d}}$ is the ideal choice since it suffers from a negative CDM density in the future. By the use of a Markov Chain Monte Carlo (MCMC) approach, and by assuming a flat FLRW Universe, we constrain the cosmological parameters of each of the three phantom DE models studied. Furthermore, by the aid of the corrected Akaike Information Criterion ($\text{AIC}_{c}$) tool, we compare our phantom DE models. Finally, a perturbative analysis of phantom DE models under consideration is performed based on the best fit background parameters.

Large samples of galaxy clusters provide knowledge of both astrophysics in the most massive virialised environments and the properties of the cosmological model that defines our Universe. However, an important issue that affects the interpretation of galaxy cluster samples is the role played by the selection waveband and the potential for this to introduce a bias in the physical properties of clusters thus selected. We aim to investigate waveband-dependent selection effects in the identification of galaxy clusters by comparing the X-ray Multi-Mirror (XMM) Ultimate Extra-galactic Survey (XXL) and Subaru Hyper Suprime-Cam (HSC) CAMIRA cluster samples identified from a common 22.6 deg2 sky area. We compare 150 XXL and 270 CAMIRA clusters in a common parameter space defined by X-ray aperture brightness and optical richness. We find that 71/150 XXL clusters are matched to the location of a CAMIRA cluster, the majority of which (67/71) display richness values N>15 that exceed the CAMIRA catalogue richness threshold. We find that 67/270 CAMIRA clusters are matched to the location of an XXL cluster (defined within XXL as an extended X-ray source). Of the unmatched CAMIRA clusters, the majority display low X-ray fluxes consistent with the lack of an XXL counterpart. However, a significant fraction (64/107) CAMIRA clusters that display high X-ray fluxes are not asociated with an extended source in the XXL catalogue. We demonstrate that this disparity arises from a variety of effects including the morphological criteria employed to identify X-ray clusters and the properties of the XMM PSF.

We have studied the evolution of different timing and spectral properties of the X-ray pulsar 2S 1417$-$624 during the recent outburst in January 2021 based on the Neutron Star Interior Composition Explorer (NICER) and Swift observations. The spin period during the outburst was $P \sim17.3622$ s based on the NICER data and the period decreases slowly with time. The evolution of the spin period and pulsed flux is studied with Fermi/GBM during the outburst. The pulse profile shows strong energy dependence and variability. The pulse profile shows multiple peaks and dips which evolve significantly with energy. The pulsed fraction shows a positive correlation with energy. The evolution of the spectral state is also studied during the outburst. The NICER energy spectrum is well described with a model of $-$ broken power-law and a blackbody emission along with a photo-electric absorption component.

Both simulations and observations have shown that turbulence is a pervasive phenomenon in cosmic scenarios, yet it is particularly difficult to model numerically due to its intrinsically multiscale character which demands high resolutions. Additionally, turbulence is tightly connected to the dynamical state and the formation history of galaxies and galaxy clusters, producing a diverse phenomenlogy which requires large samples of such structures to attain robust conclusions. In this work, we use an adaptive mesh refinement (AMR) cosmological simulation to explore the generation and dissipation of turbulence in galaxy clusters, in connection to its assembly history. We find that major mergers, and more generally accretion of gas, is the main process driving turbulence in the ICM. We have especially focused on solenoidal turbulence, which can be quantified through enstrophy. Our results seem to confirm a scenario for its generation which involves baroclinicity and compression at the external (accretion) and internal (merger) shocks, followed by vortex stretching downstream of them. We have also looked at the infall of mass to the cluster beyond its virial boundary, finding that gas follows trajectories with some degree of helicity, as it has already developed some vorticity in the external shocks.

$\tau$ Sco, a well-studied magnetic B-type star in the Upper Sco association, has a number of surprising characteristics. It rotates very slowly and shows nitrogen excess. Its surface magnetic field is much more complex than a purely dipolar configuration which is unusual for a magnetic massive star. We employ the CMFGEN radiative transfer code to determine the fundamental parameters and surface CNO and helium abundances. Then, we employ MESA and GENEC stellar evolution models accounting for the effects of surface magnetic fields. To reconcile $\tau$ Sco's properties with single-star models, an increase is necessary in the efficiency of rotational mixing by a factor of 3 to 10 and in the efficiency of magnetic braking by a factor of 10. The spin down could be explained by assuming a magnetic field decay scenario. However, the simultaneous chemical enrichment challenges the single-star scenario. Previous works indeed suggested a stellar merger origin for $\tau$ Sco. However, the merger scenario also faces similar challenges as our magnetic single-star models to explain $\tau$ Sco's simultaneous slow rotation and nitrogen excess. In conclusion, the single-star channel seems less likely and versatile to explain these discrepancies, while the merger scenario and other potential binary-evolution channels still require further assessment as to whether they may self-consistently explain the observables of $\tau$ Sco.

Vortices are one of the most promising mechanisms to locally concentrate millimeter dust grains and allow the formation of planetesimals through gravitational collapse. The outer disk around the binary system HD 142527 is known for its large horseshoe structure with azimuthal contrasts of 3-5 in the gas surface density and of about 50 in the dust. Using 13CO and C18O J = 3-2 transition lines, we detect kinematic deviations to the Keplerian rotation, which are consistent with the presence of a large vortex around the dust crescent, as well as a few spirals in the outer regions of the disk. Comparisons with a vortex model suggest velocity deviations up to 350 m/s after deprojection compared to the background Keplerian rotation, as well as an extension of about 40 au radially on both sides of the vortex and 200 degrees azimuthally, yielding an azimuthal-to-radial aspect ratio of 5. Another alternative for explaining the vortex-like signal implies artificial velocity deviations generated by beam smearing in association with variations of the gas velocity due to gas pressure gradients at the inner and outer edges of the circumbinary disk. The two scenarios are not currently discernible and, for this purpose, would probably require the use of multiple lines at a higher spatial resolution. The beam smearing effect, due to the finite spatial resolution of the observations and gradients in the line emission, should be common in observations of protoplanetary disks and may lead to misinterpretations of the gas velocity, in particular around ring-like structures.

We report on our extensive photometry and imaging of Comet 46P/Wirtanen during its 2018/19 apparition and use these data to constrain modeling of Wirtanen's activity. Narrowband photometry was obtained on nine epochs from 2018 October through 2019 March as well as 10 epochs during the 1991, 1997, and 2008 apparitions. The ensemble photometry reveals a typical composition and a secular decrease in activity since 1991. Production rates were roughly symmetric around perihelion for the carbon-bearing species (CN, C$_3$, and C$_2$), but steeper for OH and NH outbound. Our imaging program emphasized CN, whose coma morphology and lightcurve yielded rotation periods reported in a companion paper (Farnham et al., PSJ, 2, 7). Here, we compare the gas and dust morphology on the 18 nights for which observations of additional species were obtained. The carbon-bearing species exhibited similar morphology that varied with rotation. OH and NH had broad, hemispheric brightness enhancements in the tailward direction that did not change significantly with rotation, which we attribute to their originating from a substantial icy grain component. We constructed a Monte Carlo model that replicates the shape, motion, and brightness distribution of the CN coma throughout the apparition with a single, self-consistent solution in principal axis rotation. Our model yields a pole having (R.A., Decl.) = 319$^\circ$, $-$5$^\circ$ (pole obliquity of 70$^\circ$) and two large sources (radii of 50$^\circ$ and 40$^\circ$) centered at near-equatorial latitudes and separated in longitude by $\sim$160$^\circ$. Applications of the model to explain observed behaviors are discussed.

The Ly$\alpha$ forest provides one of the best means of mapping large-scale structure at high redshift, including our tightest constraint on the distance-redshift relation before cosmic noon. We describe how the large-scale correlations in the Ly$\alpha$ forest can be understood as an expansion in cumulants of the optical depth field, which itself can be related to the density field by a bias expansion. This provides a direct connection between the observable and the statistics of the matter fluctuations which can be computed in a systematic manner. We discuss the way in which complex, small-scale physics enters the predictions, the origin of the much-discussed velocity bias and the `renormalization' of the large-scale bias coefficients. Our calculations are within the context of perturbation theory, but we also make contact with earlier work using the peak-background split. Using the structure of the equations of motion we demonstrate, to all orders in perturbation theory, that the large-scale flux power spectrum becomes the linear spectrum times the square of a quadratic in the cosine of the angle to the line of sight. Unlike the case of galaxies, both the isotropic and anisotropic pieces receive contributions from small-scale physics.

Parker (1972) first proposed that coronal heating was the necessary outcome of an energy flux caused by the tangling of coronal magnetic field lines by photospheric flows. In this paper we discuss how this model has been modified by subsequent numerical simulations outlining in particular the substantial differences between the "nanoflares" introduced by Parker and "elementary events", defined here as small-scale spatially and temporally isolated heating events resulting from the continuous formation and dissipation of field-aligned current sheets within a coronal loop. We present numerical simulations of the compressible 3-D MHD equations using the HYPERION code. We use two clustering algorithms to investigate the properties of the simulated elementary events: an IDL implementation of a Density-Based Spatial Clustering of Applications with Noise (DBSCAN) technique; and our own Physical Distance Clustering (PDC) algorithm. We identify and track elementary heating events in time, both in temperature and in Joule heating space. For every event we characterize properties such as: density, temperature, volume, aspect ratio, length, thickness, duration and energy. The energies of the events are in the range $10^{18}-10^{21}$ ergs, with durations shorter than 100 seconds. A few events last up to 200 seconds and release energies up to $10^{23}$ ergs. While high temperature are typically located at the flux tube apex, the currents extend all the way to the footpoints. Hence a single elementary event cannot at present be detected. The observed emission is due to the superposition of many elementary events distributed randomly in space and time within the loop.

We present the systematic spectral analyses of gamma-ray bursts (GRBs) detected by the Fermi Gamma-Ray Burst Monitor (GBM) during its first ten years of operation. This catalog contains two types of spectra; time-integrated spectral fits and spectral fits at the brightest time bin, from 2297 GRBs, resulting in a compendium of over 18000 spectra. The four different spectral models used for fitting the spectra were selected based on their empirical importance to the shape of many GRBs. We describe in detail our procedure and criteria for the analyses, and present the bulk results in the form of parameter distributions both in the observer frame and in the GRB rest frame. 941 GRBs from the first four years have been re-fitted using the same methodology as that of the 1356 GRBs in years five through ten. The data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC).

As more exoplanets are being discovered around ultracool dwarfs, understanding their magnetic activity -- and the implications for habitability -- is of prime importance. To find stellar flares and photometric signatures related to starspots, continuous monitoring is necessary, which can be achieved with spaceborn observatories like the Transiting Exoplanet Survey Satellite (TESS). We present an analysis of TRAPPIST-1 like ultracool dwarfs with TESS full-frame image photometry from the first two years of the primary mission. A volume-limited sample up to 50 pc is constructed consisting of 339 stars closer than 0.5 mag to TRAPPIST-1 on the Gaia colour-magnitude diagram. The 30-min cadence TESS light curves of 248 stars were analysed, searching for flares and rotational modulation caused by starspots. The composite flare frequency distribution of the 94 identified flares shows a power law index similar to TRAPPIST-1, and contains flares up to $E_\mathrm{TESS} = 3 \times 10^{33}$ erg. Rotational periods shorter than 5 days were determined for 42 stars, sampling the regime of fast rotators. The ages of 88 stars from the sample were estimated using kinematic information. A weak correlation between rotational period and age is observed, which is consistent with magnetic braking.

We present a photometric and kinematical analysis of two and poorly studied open clusters; Koposov 12 (FSR 802) and Koposov 43 (FSR 848) by using cross-matched data from PPMXL and Gaia DR2 catalog. We use astrometric parameters to identify 285 and 310 cluster members for Koposov 12 and Koposov 43, respectively. Using the extracted member candidates and isochrone fitting to near-infrared (J, H, Ks) and Gaia DR2 bands (G, GBP, GRP), and Color Magnitude Diagrams (CMDs), we have estimated ages: log (age/yr) = 9.00 +/- 0.20 and 9.50 +/- 0.20, and distances d = 1850 +/- 43 pc and 2500 +/- 50 pc for Koposov 12 and Koposov 43, respectively, assuming Solar metallicity (Z=0.019). The estimated masses of the cluster derived using initial mass function and synthetic CMD are 364 +/- 19 M_sun and 352 +/- 19 M_sun. We have also computed their velocity ellipsoid parameters based on (3x3) matrix elements (mu_ij).

The second supernova that forms double-neutron-star systems is expected to occur in a progenitor that is ultra-stripped due to binary interactions. Thus, the secondary neutron star's mass as well as the post-supernova binary's orbital parameters will depend on the nature of the collapsing progenitor core. Since neutron stars are in the strong-gravity regime, their binding energy makes up a significant fraction of their total mass-energy budget. The second neutron star's binding energy may thus provide a unique insight as to whether its progenitor was a low-mass iron core or an oxygen-neon-magnesium core. I obtain relations for the baryonic mass and binding energy incorporating both a hadronic equation-of-state catalog as well as recent multi-messenger neutron-star observations. With these relations, I obtain the first tight constraints on the baryonic mass and binding energy of three neutron stars that are thought to have formed from an ultra-stripped progenitor. With these constraints, I test if each neutron star is consistent with forming from an ONeMg core that undergoes an electron-capture supernova. From these tests, I find that this scenario can be ruled out for one of three neutron stars. Neutron-star properties and the dense-matter equation of state can thus help distinguish binary formation scenarios.

Due to the finite amount of observational data, the best-fit parameters corresponding to the reconstructed cluster mass have uncertainties. In turn, these uncertainties affect the inferences made from these mass models. Following our earlier work, we have studied the effect of such uncertainties on the singularity maps in simulated and actual galaxy clusters. The mass models for both simulated and real clusters have been constructed using grale. The final best-fit mass models created using grale give the simplest singularity maps and a lower limit on the number of point singularities that a lens has to offer. The simple nature of these singularity maps also puts a lower limit on the number of three image (tangential and radial) arcs that a cluster lens has. Hence, we estimate the number of galaxy sources giving rise to the three image arcs, which can be observed with the James Webb Space Telescope (JWST). We find that we expect to observe at least 20-30 tangential and 5-10 radial three-image arcs in the Hubble Frontier Fields cluster lenses with the JWST.

We present constraints on extensions to the $\Lambda$CDM cosmological model from measurements of the $E$-mode polarization auto-power spectrum and the temperature-$E$-mode cross-power spectrum of the cosmic microwave background (CMB) made using 2018 SPT-3G data. The extensions considered vary the primordial helium abundance, the effective number of relativistic degrees of freedom, the sum of neutrino masses, the relativistic energy density and mass of a sterile neutrino, and the mean spatial curvature. We do not find clear evidence for any of these extensions, from either the SPT-3G 2018 dataset alone or in combination with baryon acoustic oscillation and \textit{Planck} data. None of these model extensions significantly relax the tension between Hubble-constant, $H_0$, constraints from the CMB and from distance-ladder measurements using Cepheids and supernovae. The addition of the SPT-3G 2018 data to \textit{Planck} reduces the square-root of the determinants of the parameter covariance matrices by factors of $1.3 - 2.0$ across these models, signaling a substantial reduction in the allowed parameter volume. We also explore CMB-based constraints on $H_0$ from combined SPT, \textit{Planck}, and ACT DR4 datasets. While individual experiments see some indications of different $H_0$ values between the $TT$, $TE$, and $EE$ spectra, the combined $H_0$ constraints are consistent between the three spectra. For the full combined datasets, we report $H_0 = 67.49 \pm 0.53\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$, which is the tightest constraint on $H_0$ from CMB power spectra to date and in $4.1\,\sigma$ tension with the most precise distance-ladder-based measurement of $H_0$. The SPT-3G survey is planned to continue through at least 2023, with existing maps of combined 2019 and 2020 data already having $\sim3.5\times$ lower noise than the maps used in this analysis.

Shock waves associated with fast coronal mass ejections (CMEs) accelerate solar energetic particles (SEPs) in the long duration, gradual events that pose hazards to crewed spaceflight and near-Earth technological assets, but the source of the CME shock-accelerated plasma is still debated. Here, we use multi-messenger observations from the Heliophysics System Observatory to identify plasma confined at the footpoints of the hot, core loops of active region 11944 as the source of major gradual SEP events in January 2014. We show that the elemental composition signature detected spectroscopically at the footpoints explains the measurements made by particle counting techniques near Earth. Our results localize the elemental fractionation process to the top of the chromosphere. The plasma confined closest to that region, where the coronal magnetic field strength is high (a few hundred Gauss), develops the SEP composition signature. This source material is continually released from magnetic confinement and accelerated as SEPs following M-, C-, and X-class flares.

The search for the first stars formed from metal-free gas in the universe is one of the key issues in astronomy because it relates to many fields, such as the formation of stars and galaxies, the evolution of the universe, and the origin of elements. It is not still clear if metal-free first stars can be found in the present universe. These first stars are thought to exist among extremely metal-poor stars in the halo of our Galaxy. Here we propose a new scenario for the formation of low-mass first stars that have survived until today and observational counterparts in our Galaxy. The first stars in binary systems, consisting of massive- and low-mass stars, are examined using stellar evolution models, simulations of supernova ejecta colliding with low-mass companions, and comparisons with observed data. These first star survivors will be observed as metal-rich halo stars in our Galaxy. We may have identified a candidate star in the observational database where elemental abundances and kinematic data are available. Our models also account for the existence of several solar-metallicity stars in the literature having space velocities equivalent to the halo population. The proposed scenario demands a new channel of star formation in the early universe and is a supplementary scenario for the origin of the known metal-poor stars.

We carried a detailed spectral and temporal study of blazar, \mbox{PKS\,0903-57} using the \emph{Fermi}-LAT and \emph{Swift}-XRT/UVOT observations, during its brightest flaring period MJD\,58931--58970. During this period, the maximum daily averaged $\gamma$-ray flux ($\rm F_{0.1-500\,GeV}$) of $\rm9.42\times10^{-6}\,ph\,cm^{-2}\,s^{-1}$ is observed on MJD\,58951.5, the highest $\gamma$-ray flux detected from \mbox{PKS\,0903-57} till now. Several high-energy (HE) photons ($>\,10$\,GeV) consistent with the source location at high probability (>\,99\%) are detected, and the $\gamma$-ray lightcurve in the active state shows multiple substructures with asymmetric profile. In order to understand the possible physical scenario responsible for the flux enhancement, we carried a detailed broadband spectral study of \mbox{PKS\,0903-57} by choosing different flux states from its active period. Neglecting the multi-band variability in each of the selected time intervals, we could reproduce their averaged broadband SEDs with a one-zone leptonic model whose parameters were derived with a $\chi^2$-fit. We found that the broadband SED during different flux states can be reproduced by the synchrotron, synchrotron-self-Compton (SSC) and External-Compton (EC) processes. The seed photons for EC process from BLR or IR torus provide acceptable fits to the GeV spectrum in all the flux states; however, the detection of HE photons together with the equipartition condition suggest that the EC/IR process is a more likely scenario. Further, a detailed comparison between the fit parameters shows that the flux enhancement from quiescent-state to the flaring-state is mostly related to increase in the bulk Lorentz factor of the emission region and change in the break energy of the source spectrum.

Rate constants have been measured for the C(3P) + CH3CN reaction between 50 K and 296 K using a continuous-flow supersonic reactor. C(3P) atoms were created by the in-situ pulsed laser photolysis of CBr4 at 266 nm, while the kinetics of C(3P) atom loss were followed by direct vacuum ultra-violet laser induced fluorescence at 115.8 nm. Secondary measurements of product H(2S) atom formation were also made, allowing absolute H-atom yields to be obtained by comparison with those obtained for the C(3P) + C2H4 reference reaction. In parallel, quantum chemical calculations were performed to obtain the various complexes, adducts and transition states relevant to the title reaction over the triplet potential energy surface, allowing us to better understand the preferred reaction pathways. The reaction is seen to be very fast, with measured rate constants in the range (3-4) x 10-10 cm3 s-1 with little or no observed temperature dependence. As the C + CH3CN reaction is not considered in current astrochemical networks, we test its influence on interstellar methyl cyanide abundances using a gas-grain dense interstellar cloud model. Its inclusion leads to predicted CH3CN abundances that are significantly lower than the observed ones.

We investigate the clustering property of primordial black holes (PBHs) in a scenario where PBHs can explain the existence of supermassive black holes (SMBHs) at high redshifts. We analyze the angular correlation function of PBHs originating from fluctuations of a spectator field which can be regarded as a representative model to explain SMBHs without conflicting with the constraint from the spectral distortion of cosmic microwave background. We argue that the clustering property of PBHs can give a critical test for models with PBHs as the origin of SMBHs and indeed show that the spatial distribution of PBHs in such a scenario is highly clustered, which suggests that those models may be disfavored from observations of SMBHs although a careful comparison with observational data would be necessary.

The coronal mass ejections (CMEs) from the Sun are known for their space weather and geomagnetic consequences. Among all CMEs, so-called radio-loud (RL) and halo CMEs are considered the most energetic in the sense that they are usually faster and wider than the general population of CMEs. Hence the study of RL and halo CMEs has become important and the prediction of their occurrence rate in a future cycle will give a warning in advance. In the present paper, the occurrence rates of RL and halo CMEs in solar cycle (SC) 25 are predicted. For this, we obtained good correlations between the numbers of RL and halo CMEs in each year and the yearly mean sunspot numbers in the previous two cycles. The predicted values of sunspot numbers in SC 25 by NOAA/NASA were considered as representative indices and the corresponding numbers of RL and halo CMEs have been determined using linear relations. Our results show that the maximum number of RL and halo CMEs will be around 39 $/pm$ 3 and 45 $/pm$ 4, respectively. Removing backside events, a set of front-side events was also considered separately and the front-side events alone in SC 25 are predicted again. The peak values of front-side RL and halo events have been estimated to be around 31 $/pm$ 3 and 29 $/pm$ 3 respectively. These results are discussed in comparison with the predicted values of sunspots by different authors.

Orbital monitoring of exoplanetary and stellar systems is fundamental for analysing their architecture, dynamical stability and evolution, and mechanisms of formation. Current high-contrast extreme-adaptive optics imagers like SPHERE, GPI, and SCExAO+CHARIS explore the population of giant exoplanets and brown dwarf and stellar companions beyond typically 10 au, covering generally a small fraction of the orbit (<20%) leading to degeneracies and biases in the orbital parameters. Precise and robust measurements over time of the position of the companions are critical, which require good knowledge of the instrumental limitations and dedicated observing strategies. The homogeneous dedicated calibration strategy for astrometry implemented for SPHERE has facilitated high-precision studies by its users since its start of operation in 2014. As the precision of exoplanet imaging instruments is now reaching milliarcseconds and is expected to improve with the upcoming facilities, we initiated a community effort, triggered by the SPHERE experience, to share lessons learned for high-precision astrometry in direct imaging. A homogeneous strategy would strongly benefit the VLT community, in synergy with VLTI instruments like GRAVITY/GRAVITY+, future instruments like ERIS and MAVIS, and in preparation for the exploitation of the ELT's first instruments MICADO, HARMONI, and METIS.

We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose.

Thermal-infrared measurements of asteroids are crucial for deriving the objects' sizes, albedos, and also the thermophysical properties of the surface material. Depending on the available data, a range of simple to complex thermal models are applied to achieve specific science goals. However, testing these models is often a difficult process and the uncertainties of the derived parameters are not easy to estimate. Here, we make an attempt to verify a widely accepted thermophysical model (TPM) against unique thermal infrared (IR), full-disk, and well-calibrated measurements of the Moon. The data were obtained by the High-resolution InfraRed Sounder (HIRS) instruments on board a fleet of Earth weather satellites that serendipitously scan over the Moon. We found 22 Moon intrusions, taken in 19 channels between 3.75 micron and 15.0 micron, and over a wide phase angle range from -73.1 deg to +73.8 deg. The similarity between these Moon data and typical asteroid spectral-IR energy distributions allows us to benchmark the TPM concepts and to point out problematic aspects. The TPM predictions match the HIRS measurements within 5% (10% at the shortest wavelengths below 5 micron when using the Moon's known properties (size, shape, spin, albedo, thermal inertia, roughness) in combination with a newly established wavelength-dependent hemispherical emissivity. In the 5-7.5 micron and in the 9.5 to 11 micron ranges, the global emissivity model deviates considerably from the known lunar sample spectra. Our findings will influence radiometric studies of near-Earth and main-belt asteroids in cases where only short-wavelength data (from e.g., NEOWISE, the warm Spitzer mission, or ground-based M-band measurements) are available. The new, full-disk IR Moon model will also be used for the calibration of IR instrumentation on interplanetary missions (e.g., for Hayabusa-2) and weather satellites.

We performed a systematic search for broad-velocity-width molecular features (BVFs) in the disk part of our Galaxy by using the CO J = 1--0 survey data obtained with the Nobeyama Radio Observatory 45 m telescope. From this search, 58 BVFs were identified. In comparisons with the infrared and radio continuum images, 36 BVFs appeared to have both infrared and radio continuum counterparts, and 15 of them are described as molecular outflows from young stellar objects in the literature. In addition, 21 BVFs have infrared counterparts only, and eight of them are described as molecular outflows in the literature. One BVF (CO 16.134--0.553) does not have any luminous counterpart in the other wavelengths, which suggests that it may be an analog of high-velocity compact clouds in the Galactic center.

The abstract is available at: this http URL

The disruption of asteroids and comets produces cm-sized meteoroids that end up impacting the Earth's atmosphere and producing bright fireballs that might have associated shock waves or, in geometrically-favorable occasions excavate craters that put them into unexpected hazardous scenarios. The astrometric reduction of meteors and fireballs to infer their atmospheric trajectories and heliocentric orbits involves a complex and tedious process that generally requires many manual tasks. To streamline the process, we present a software package called SPMN 3D Fireball Trajectory and Orbit Calculator (3D-FireTOC), an automatic Python code for detection, trajectory reconstruction of meteors, and heliocentric orbit computation from video recordings. The automatic 3D-FireTOC package comprises of a user interface and a graphic engine that generates a realistic 3D representation model, which allows users to easily check the geometric consistency of the results and facilitates scientific content production for dissemination. The software automatically detects meteors from digital systems, completes the astrometric measurements, performs photometry, computes the meteor atmospheric trajectory, calculates the velocity curve, and obtains the radiant and the heliocentric orbit, all in all quantifying the error measurements in each step. The software applies corrections such as light aberration, refraction, zenith attraction, diurnal aberration and atmospheric extinction. It also characterizes the atmospheric flight and consequently determines fireball fates by using the $\alpha - \beta$ criterion that analyses the ability of a fireball to penetrate deep into the atmosphere and produce meteorites. We demonstrate the performance of the software by analyzing two bright fireballs recorded by the Spanish Fireball and Meteorite Network (SPMN).

We present high signal-to-noise, full polarization pulse profiles for 40 bright, 'slowly'-rotating (non-recycled) pulsars using the new Ultra-Wideband Low-frequency (UWL; 704-4032 MHz) receiver on the Parkes radio telescope. We obtain updated and accurate interstellar medium parameters towards these pulsars (dispersion measures and Faraday rotation measures), and reveal Faraday dispersion towards PSR J1721-3532 caused by interstellar scattering. We find general trends in the pulse profiles including decreasing fractional linear polarization and increasing degree of circular polarization with increasing frequency, consistent with previous studies, while also revealing new features and frequency evolution. This demonstrates results that can be obtained using UWL monitoring observations of slow pulsars, which are valuable for improving our understanding of pulsar emission and the intervening interstellar medium. The calibrated data products are publicly available.

Light curves for more than century optical photometric observations of the blazar OJ 287 reveals strong flares with a quasi-period of about 12 years. For a long time, this period has been interpreted by processes in a binary black hole system. We propose an alternative explanation for this period, which is based on Doppler factor periodic variations of the emitting region caused by jet helicity. Using multi-epoch very large baseline interferometry (VLBI) observations carried out in a framework of the MOJAVE (Monitoring Of Jets in Active galactic nuclei with VLBA Experiments) program and other VLBA (Very Long Baseline Array) archival experiments at the observing frequency of 15 GHz, we derived geometrical parameters of the jet helix. To reach an agreement between the VLBI and photometric optical observation data, the jet component motion at a small angle to the radial direction is necessary. Such non-radial motion is observed and, together with the jet helical shape, can be naturally explained by the development of the Kelvin-Helmholtz instability in the parsec-scale outflow. In this case, the true precession of the OJ 287 jet may manifest itself in differences between the peak flux values of the 12-year optical flares. A possibility to create this precession due to Lense-Thirring effect of a single supermassive black hole is also discussed.

This paper reports 209 O-type stars found with LAMOST. All 135 new O-type stars discovered so far with LAMOST so far are given. Among them, 94 stars are firstly presented in this sample. There are 1 Iafpe star, 5 Onfp stars, 12 Oe stars, 1 Ofc stars, 3 ON stars, 16 double-lined spectroscopic binaries, and 33 single-lined spectroscopic binaries. All O-type stars are determined based on LAMOST low-resolution spectra (R ~ 1800), with their LAMOST median-resolution spectra (R~7500) as supplements.

The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to be affected not only by the stellar intrinsic properties but also by the inclination angle of the stellar rotation axis. Until now such an inclination effect on S-index has remained largely unexplored. To fill this gap, we develop a physics-based model to calculate S-index, focusing on the Sun. Using the distributions of solar magnetic features derived from observations together with Ca II H & K spectra synthesized in non-local thermodynamic equilibrium, we validate our model by successfully reconstructing the observed variations of solar S-index over four activity cycles. Further, using the distribution of magnetic features over the visible solar disk obtained from surface flux transport simulations, we obtain S-index time series dating back to 1700 and investigate the effect of inclination on S-index variability, both on the magnetic activity cycle and the rotational timescales. We find that when going from an equatorial to a pole-on view, the amplitude of S-index variations decreases weakly on the activity cycle timescale and strongly on the rotational timescale (by about 22% and 81%, respectively, for a cycle of intermediate strength). The absolute value of S-index depends only weakly on the inclination. We provide analytical expressions that model such dependencies.

The hypothesis that the late Universe is isotropic and homogeneous is adopted by most cosmological studies. The expansion rate $H_0$ is thought to be spatially constant, while bulk flows are often presumed to be negligible compared to the Hubble expansion, even at local scales. Their effects on the redshift-distance conversion are hence usually ignored. Any deviation from this consensus can strongly bias the results of such studies and thus the importance of testing these assumptions cannot be understated. Scaling relations of galaxy clusters can be effectively used for that. In previous works, we observed strong anisotropies in cluster scaling relations, whose origins remain ambiguous. By measuring many different cluster properties, several scaling relations with different sensitivities can be built. Nearly independent tests of cosmic isotropy and bulk flows are then feasible. We make use of up to 570 clusters with measured properties at X-ray, microwave, and infrared wavelengths, to construct 10 different cluster scaling relations (five of them presented for the first time) and test the isotropy of the local Universe. Through rigorous tests, we ensure that our analysis is not prone to generally known systematic biases and X-ray absorption issues. By combining all available information, we detect an apparent $9\%$ spatial variation in the local $H_0$ between $(l,b)\sim ({280^{\circ}}^{+35^{\circ}}_{-35^{\circ}},{-15^{\circ}}^{+20^{\circ}}_{-20^{\circ}})$ and the rest of the sky. The observed anisotropy has a nearly dipole form. Using Monte Carlo simulations, we assess the statistical significance of the anisotropy to be $>5\sigma$. This result could also be attributed to a $\sim 900$ km/s bulk flow which seems to extend out to at least $\sim 500$ Mpc. These two effects are indistinguishable until more high$-z$ clusters are observed by future all-sky surveys, such as eROSITA.

We present Quasarnet, a novel research platform that enables deployment of data-driven modeling techniques for the investigation of the properties of super-massive black holes. Black hole data sets -- observations and simulations -- have grown rapidly in the last decade in both complexity and abundance. However, our computational environments and tool sets have not matured commensurately to exhaust opportunities for discovery with these data. Our pilot study presented here is motivated by one of the fundamental open questions in understanding black hole formation and assembly across cosmic time - the nature of the black hole host galaxy and parent dark matter halo connection. To explore this, we combine and co-locate large, observational data sets of quasars, the high-redshift luminous population of accreting black holes, at z > 3 alongside simulated data spanning the same cosmic epochs in Quasarnet. We demonstrate the extraction of the properties of observed quasars and their putative dark matter parent halos that permit studying their association and correspondence. In this paper, we describe the design, implementation, and operation of the publicly queryable Quasarnet database and provide examples of query types and visualizations that can be used to explore the data. Starting with data collated in Quasarnet, which will serve as training sets, we plan to utilize machine learning algorithms to predict properties of the as yet undetected, less luminous quasar population. To that ultimate goal, here we present the first key step in building the BH-galaxy-halo connection that underpins the formation and evolution of supermassive black holes. All our codes and compiled data are available on the public Google Kaggle Platform.

Context. Class II methanol masers at 6.7 and 12.2 GHz occur close to high-mass young stellar objects (HMYSOs). When they are observed simultaneously, such studies may contribute to refining the characterisation of local physical conditions. Aims. We aim to search for the 12.2 GHz methanol emission in 6.7 GHz methanol masers that might have gone undetected in previous surveys of northern sky HMYSOs, mainly due to their variability. Contemporaneous observations of both transitions are used to refine the flux density ratio and examine the physical parameters. Methods. We observed a sample of 153 sites of 6.7 GHz methanol maser emission in the 12.2 GHz methanol line with the Torun 32 m radio telescope, using the newly built X-band receiver. Results. The 12.2 GHz methanol maser emission was detected in 36 HMYSOs, with 4 of them detected for the first time. The 6.7 GHz to 12.2 GHz flux density ratio for spectral features of the contemporaneously observed sources has a median value of 5.1, which is in agreement with earlier reports. The ratio differs significantly among the sources and for the periodic source G107.298+5.639 specifically, the ratio is weakly recurrent from cycle to cycle, but it generally reaches a minimum around the flare peak. This is consistent with the stochastic maser process, where small variations in the physical parameters along the maser path can significantly affect the ratio. A comparison of our data with historical results (from about ten years ago) implies significant (> 50%) variability for about 47% and 14% at 12.2 GHz and 6.7 GHz, respectively. This difference can be explained via the standard model of methanol masers.

During solar flares, a large flux of energetic electrons propagate from the tops of reconnecting magnetic flux tubes toward the lower atmosphere. Over the course of the electrons' transport, a co-spatial counter-streaming return current is induced, thereby balancing the current density. In response to the return current electric field, a fraction of the ambient electrons will be accelerated into the runaway regime. However, models describing the accelerated electron beam/return-current system have generally failed to take these suprathermal runaway electrons into account self-consistently. We develop a model in which an accelerated electron beam drives a steady-state, sub-Dreicer co-spatial return-current electric field, which locally balances the direct beam current and freely accelerates a fraction of background (return-current) electrons. The model is self-consistent, i.e., the electric field induced by the co-evolution of the direct beam and the runaway current is considered. We find that (1) the return current electric field can return a significant number of suprathermal electrons to the acceleration region, where they can be further accelerated to higher energies, runaway electrons can be a few tens of percent of the return current flux returning to the nonthermal beam's acceleration region, (2) the energy gain of the suprathermal electrons can be up to $10-35$\,keV, (3) the heating rate in the corona can be reduced by an order of magnitude in comparison to models which neglect the runaway component. The results depend on the injected beam flux density, the temperature and density of the background plasma.

We derive all contributions to the dispersion measure (DM) of electromagnetic pulses to linear order in cosmological perturbations, including both density fluctuations and relativistic effects. We then use this result to calculate the power spectrum of DM-based cosmological observables to linear order in perturbations. In particular, we study two cases: maps of the dispersion measure from a set of localized sources (including the effects of source clustering), and fluctuations in the density of DM-selected sources. The impact of most relativistic effects is limited to large angular scales, and is negligible for all practical applications in the context of ongoing and envisaged observational programs targeting fast radio bursts. We compare the leading contributions to DM-space clustering, including the effects of gravitational lensing, and find that the signal is dominated by the fluctuations in the free electron column density, rather than the local source clustering or lensing contributions. To compensate for the disappointing irrelevance of relativistic effects, we re-derive them in terms of the geodesic equation for massive particles in a perturbed Friedmann-Robertson-Walker metric.

We investigate clustering properties of dark matter halos and galaxies to search for optimal statistics and scales where possible departures from general relativity (GR) could be found. We use large N-body cosmological simulations to perform measurements based on the two-point correlation function (2PCF) in GR and in selected modified gravity (MG) structure formation scenarios. As a test-bed, we employ two popular beyond-GR models: $f(R)$ gravity and the normal branch of the Dvali-Gabadadze-Porrati (nDGP) braneworld. We study a range of simulated halo and galaxy populations and reveal a noticeable MG signal in the monopole and quadrupole moments of the redshift-space 2PCF, and in the so-called clustering wedges. However, once expressed in terms of the linear distortion parameter, $\beta$, the statistical significance of these signals largely diminishes due to a strong degeneracy between MG-enhanced clustering and modified tracer bias. To circumvent this, we consider statistics less dependent on the bias: relative clustering ratios. We generalize the monopole ratio proposed in earlier work to multipole moments and clustering wedges, and introduce a new estimator of the $\beta$ parameter. The clustering ratios we extract foster noticeable differences between MG and GR models, reaching a maximum deviation of 10\% at 2$\sigma$ significance for specific variants of $f(R)$ and nDGP. We show that such departures could be measured for $\beta$ if non-linear effects at intermediate scales are correctly modeled. Our study indicates that the clustering ratios give great promise to search for signatures of MG in the large-scale structure. We also find that the selection of an optimal tracer sample depends on a particular statistics and gravity model to be considered.

Recently, a nitrogen iceberg was proposed as a possible origin for the first interstellar object, 1I/2017 U1, also known as `Oumuamua. Here, we show that the mass budget in exo-Pluto planets necessary to explain the detection of `Oumuamua as a nitrogen iceberg chipped off from a planetary surface is comparable to the total mass in stars, making the scenario exceedingly unlikely.

We revisit to investigate shadows cast by Kerr-like wormholes. The boundary of the shadow is determined by unstable circular photon orbits. We find that, in certain parameter regions, the orbit is located at the throat of the Kerr-like wormhole, which was not considered in the literatures. In these cases, the existence of the throat alters the shape of the shadow significantly, and it will be much easier to differentiate it from that of a Kerr black hole.

We discuss the unitarity bounds on vectors coupled to currents whose non-conservation is due to mass terms, such as $U(1)_{L_\mu - L_\tau}$. Due to emission of many final state longitudinally polarized gauge bosons, inclusive rates grow exponentially fast in energy, leading to constraints that are only logarithmically dependent on the symmetry breaking mass term. This exponential growth is unique to Stueckelberg theories and reverts back to polynomial growth at energies above the mass of the radial mode. As an example, we demonstrate how the total inelastic cross section of the LHC beats out cosmological bounds to place the strongest limit on Stueckelberg $U(1)_{L_\mu - L_\tau}$ models for most masses below a keV. We also present a stronger, but more uncertain, bound coming from the validity of perturbation theory at the LHC.

We reconsider the widely held view that the Mannheim--Kazanas (MK) vacuum solution for a static, spherically-symmetric system in conformal gravity (CG) predicts flat rotation curves, such as those observed in galaxies, without the need for dark matter. The conformal equivalence of the MK and Schwarzschild--de-Sitter (SdS) metrics, where the latter does not predict flat rotation curves, raises concerns that the prediction in the MK frame may be a gauge artefact. This ambiguity arises from assuming that, in each frame, test particles have fixed rest mass and follow timelike geodesics, which are not conformally invariant. The mass of such particles must instead be generated dynamically through interaction with a scalar field, the energy-momentum of which means that the spacetime outside a static, spherically-symmetric matter distribution in CG is, in general, not given by the MK vacuum solution. A unique solution does exist, however, for which the scalar field energy-momentum vanishes and the metric retains the MK form. Nonetheless, we show that in both the Einstein and MK frames of this solution, in which the scalar field is constant or radially-dependent, respectively, massive particles follow timelike geodesics of the SdS metric, thereby resolving the apparent frame dependence of physical predictions and unambiguously yielding rotation curves with no flat region. We further find that the general form of the conformal transformation linking the Einstein and MK frames is unique in preserving the structure of any diagonal static, spherically-symmetric metric with a radial coefficient that is (minus) the reciprocal of its temporal one. We also comment briefly on how our analysis resolves the long-standing uncertainty regarding gravitational lensing in the MK metric. (Abridged)

We study the quasi-periodic oscillations from the accretion disk around rotating traversable wormholes by means of the resonance models. We investigate the linear stability of the circular geodesic orbits in the equatorial plane for a general class of wormhole geometries deriving analytical expressions for the epicyclic frequencies. Since wormholes can often mimic black holes in the astrophysical observations, we analyze the properties of the quasi-circular oscillatory motion in comparison with the Kerr black hole. We demonstrate that wormholes possess distinctive features, which can be observationally significant. It is characteristic for the Kerr black hole that the orbital and the epicyclic frequencies obey a constant ordering in the whole range of the spin parameter. In contrast, for wormhole spacetimes we can have various types of orderings between the frequencies in the different regions of the parametric space. This enables the excitation of much more diverse types of resonances including parametric and forced resonances of lower order, which could lead to stronger observable signals. In addition, for co-rotating orbits the resonances can be excited in a very close neighbourhood of the wormhole throat for a wide range of values of the angular momentum, making wormholes a valuable laboratory for testing strong gravity.

We test the validity of the Generalized Heisenberg's Uncertainty principle in the presence of strong gravitational fields nearby rotating black holes; Heisenberg's principle is supposed to require additional correction terms when gravity is taken into account, leading to a more general formulation also known as the Generalized Uncertainty Principle. Using as probe electromagnetic waves acquiring orbital angular momentum when lensed by a rotating black hole, we find from numerical simulations a relationship between the spectrum of the orbital angular momentum of light and the corrections needed to formulate the Generalized Uncertainty Principle, here characterized by the rescaled parameter $\beta_0$, a function of the Planck's mass and the bare mass of the black hole. Then, from the analysis of the observed twisted light due to the gravitational field of the compact object observed in M87*, we find new limits for the parameter $\beta_0$. With this method, complementary to black hole shadow circularity analyses, we obtain more precise limits from the experimental data of M87*, confirming the validity of scenarios compatible with General Relativity, within the uncertainties due to the experimental errors present in EHT data and those due to the numerical simulations and analysis.

The mechanism of the generation of dark matter and dark radiation from the evaporation of primordial black holes is very interesting. We consider the case of Kerr black holes to generalize previous results obtained in the Schwarzschild case. For dark matter, the results do not change dramatically and the bounds on warm dark matter apply similarly: in particular, the Kerr case cannot save the scenario of black hole domination for light dark matter. For dark radiation, the expectations for $\Delta N_{eff}$ do not change significantly with respect to the Schwarzschild case, but for an enhancement in the case of spin 2 particles: in the massless case, however, the projected experimental sensitivity would be reached only for extremal black holes.

In both Newtonian gravity and Einstein gravity there is no force on a test particle located inside a spherical cavity cut out of a static, spherically symmetric mass distribution. Inside the cavity exterior matter is decoupled and there is no external field effect that could act on the test particle. However, for potentials other than the Newtonian potential or for geometries other than Ricci flat ones this is no longer the case, and there then is an external field effect. We explore this possibility in various alternate gravity scenarios, and suggest that such (Machian) external field effects can serve as a diagnostic for gravitational theory.

The cosmology of a standard model (SM) gauge singlet complex scalar dark matter (DM), stabilized by a reflection symmetry, is studied including all renormalizable interactions that preserve the reflection symmetry but can break the larger global U(1) symmetry of DM number. We find an interesting interplay of the ensuing DM self-scatterings and annihilations in generating the present DM density, and possible particle-antiparticle asymmetry in the DM sector. The role of DM self-scatterings in determining its present density and composition is a novel phenomenon. The simultaneous presence of the self-scatterings and annihilations is required to obtain a non-zero asymmetry, which otherwise vanishes due to unitarity sum rules.

A first-order Electroweak Phase Transition (EWPT) could explain the observed baryon-antibaryon asymmetry and have a detectable gravitational wave signature, but this is not possible in the Standard Model. We therefore study the EWPT in the Standard Model Effective Field Theory (SMEFT) including dimension-six operators. A first-order EWPT has previously been shown to be possible in the SMEFT in scenarios with a tree-level barrier between minima, which requires a negative Higgs quartic coupling and a new physics scale low enough to raise questions about the validity of the EFT approach. In this work we show that a first-order EWPT is possible in a novel scenario where the barrier between minima is generated radiatively, the quartic coupling is positive, the scale of new physics is higher, and there is good agreement with experimental bounds. Our calculation is done in a consistent, gauge-invariant way, and we carefully analyze the scaling of parameters necessary to generate a barrier in the potential. We perform a global fit in the relevant parameter space and explicitly find the points with a first-order transition that agree with experimental data. We also briefly discuss the prospects for probing the allowed parameter space using di-Higgs production in colliders.