Papers for Tuesday, Mar 14 2023

WD 1856+534 b (or WD 1856 b for short) is the first known transiting planet candidate around a white dwarf star. WD 1856 b is about the size of Jupiter, has a mass less than about 12 Jupiter masses, and orbits at a distance of about 2% of an astronomical unit. The formation and migration history of this object is still a mystery. Here, we present constraints on the presence of long-period companions (where we explored eccentricity, inclination, mass, and period for the possible companion) in the WD 1856+534 planetary system from Transit Timing Variations (TTVs). We show that existing transit observations can rule out planets with orbital periods less than about 500 days. With additional transit observations over the next decade, it will be possible to test whether WD 1856 also hosts additional long-period planets that could have perturbed WD 1856 b into its current close-in orbit.

Dwarf galaxies provide a unique opportunity for studying the evolution of the Milky Way (MW) and the Local Group as a whole. Analysing the running solar apex based on the kinematics of the MW satellites, we discovered an unexpected behaviour of the dipole term of the radial velocity distribution as a function of the galactocentric distance. The nearby satellites (<100 kpc) have a bulk motion with an amplitude of 140-230 km/s while the more distant ones show an isotropic distribution of the radial velocities. Such strong solar apex variations can not be explained by the net rotation of the satellites, as it would require an enormously high rotation rate (~970 km/s). If we exclude the LMC and its most closely related members from our sample, this does not suppress the bulk motion of the nearby satellites strongly enough. Nevertheless, we have demonstrated that the observed peculiar kinematics of the MW satellites can be explained by a perturbation caused by the first infall of the LMC. First, we `undone' the effect of the perturbation by integrating the orbits of the MW satellites backwards (forwards) with (without) massive LMC. It appears that the present-day peculiar enhancement of the solar apex in the inner halo is diminished the most in the case of 2x10^{11} Msun LMC. Next, in self-consistent high-resolution N-body simulations of the MW-LMC interaction, we found that the solar apex shows the observed behaviour only for the halo particles with substantial angular momentum, comparable to that of the MW satellites.

Context: this paper describes the detection of wide binary and multiple central stars (CSs) of Galactic planetary nebulae (PNe) using the most up-to-date data available from the Gaia Data Release 3 (Gaia DR3). Aims: the objective of this study is to benefit from the Gaia DR3's reliable measurements of parallax and proper motion to reveal the binary, ternary, and quadruple CSs of planetary nebulae. Methods: in our search for the binary and multiple CSs in the Gaia archive, we adopted the criteria provided in the literature to ensure that the CS and its partner(s) have comparable parallaxes and proper motions. Results: we have detected a total of 45 binary and multiple stellar systems coupled with the CSs of PNe. Based on the standard error of the parallax, this sample was divided into four categories: highest probable, probable, possible, and uncertain systems, which comprise 18, 8, 10, and 9 objects, respectively. Except for A35, NGC246 and IC 3568, the radial velocities of the CSs companions are unknown for our PNe sample. The radial velocity measurements of these three companion stars are comparable to their linked CSs. The results show the detection of a highly probable quadruple system, and a probable quadruple system: NGC6853 and PHRJ1129-6012, respectively. In addition, we found one highly probable (Fr 2-42), one probable (M 1-58) and two possible (IC 2553 and PHRJ1123-6030) ternary CS systems. The results further show that the primary components of eight wide and very wide binary systems are actually close binary systems. Moreover, the masses of the components of nine PN binary and multiple wide systems were calculated.

Type Ia supernovae (SNe Ia) are important cosmological probes and contributors to galactic nucleosynthesis, particularly of the iron group elements. To improve both their reliability as cosmological probes and to understand galactic chemical evolution, it is vital to understand the binary progenitor system and explosion mechanism. The classification of SNe Ia into Branch groups has led to some understanding of the similarities and differences among the varieties of observed SNe Ia. However, partly due to small sample size, little work has been done on the broad-line or 02bo group. We perform direct spectral analysis on the pre-maximum spectra of the broad-line SN 2019ein and compare and contrast it to the core-normal SN~2011fe. Both SN 2019ein and SN 2011fe were first observed spectroscopically within two days of discovery, allowing us to follow the spectroscopic evolution of both supernovae in detail. We find that the optical depths of the primary features of both the CN and BL supernovae are very similar, except that there is a velocity shift between them. We further examine the SN 2002bo-like subclass and show that for nine objects with pre-maximum spectra in the range -6 -- -2 days with respect to B-maximum all the emission peaks of the Si II {\lambda}6355 line of BL are blueshifted pre-maximum, making this a simple classification criterion.

The NASA Exoplanet Archive was searched for planets with an equilibrium temperature below 600 K, mass uncertainty less than 27 percent, and radius uncertainty less than 8 percent. This search produced 93 planets with mass from 0.3 to 1680 ME; and 101 planets if the Solar System planets are included. The characteristics of the sample in this catalog are: (1) 94 percent of the Terrestrial planets have mass less than 2.9 ME and radius less than 1.4 RE, (3) The sample has a small drop in population consistent with the previously identified radius gap from 1.5 to 2.0 RE, (4) Planets in the radius range 1.50 to 2.25 RE are consistent with either a gas-rich Terrestrial composition or a rock-ice Terrestrial composition with a supercritical hydrosphere and water mass fraction less than 20 percent, (5) A super-Neptune radius desert is observed for the radius range 4.5 to 7.5 RE, (6) Saturn composition planets have masses from 15 to 170 ME and radii from 7.9 to 10.1 RE, (7) A nearly barren sub-Saturn mass-radius desert is found in the sample as indicated by a lack of planets with mass exceeding 20 ME and radii in the range 4.0 to 7.5 RE, (8) Most Jupiter composition planets have radii between 10.9 and 12.4 RE and mass exceeding 200 ME, (9) With few exceptions, planet radius can be used as a proxy for planet composition classification into Terrestrial, gas-rich Terrestrial or supercritical hydrosphere Terrestrial, Rock-Ice Giant, and Gas Giant composition classes for this sample of Teq less than 600 K planets. The characteristics of this sample are consistent with several predictions of the core accretion model including the predicted values for the critical core mass for gas accretion and runaway accretion, the pebble isolation mass, and the Saturn mass desert.

It is difficult to accurately identify galaxy mergers and it is an even larger challenge to classify them by their mass ratio or merger stage. In previous work we used a suite of simulated mergers to create a classification technique that uses linear discriminant analysis (LDA) to identify major and minor mergers. Here, we apply this technique to 1.3 million galaxies from the SDSS DR16 photometric catalog and present the probability that each galaxy is a major or minor merger, splitting the classifications by merger stages (early, late, post-coalescence). We present publicly-available imaging predictor values and all of the above classifications for one of the largest-yet samples of galaxies. We measure the major and minor merger fraction ($f_{\mathrm{merg}}$) and build a mass-complete sample of galaxies, which we bin as a function of stellar mass and redshift. For the major mergers, we find a positive slope of $f_{\mathrm{merg}}$ with stellar mass and negative slope of $f_{\mathrm{merg}}$ with redshift between stellar masses of $10.5 < M_* (log\ M_{\odot}) < 11.6$ and redshifts of $0.03 < z < 0.19$. We are able to reproduce an artificial positive slope of the major merger fraction with redshift when we do not bin for mass or craft a complete sample, demonstrating the importance of mass completeness and mass binning. We determine that the positive trend of the major merger fraction with stellar mass is consistent with a hierarchical assembly scenario. The negative trend with redshift requires that an additional assembly mechanism, such as baryonic feedback, dominates in the local Universe.

The Swift Burst Alert Telescope (BAT) is a coded aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consists of eighty-channel detector plane histograms that accumulate photon counts over time periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey dataset between $14$ and $195$ keV. Survey data comprises $> 90\%$ of all BAT data by volume and allows for the tracking of long term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey dataset has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis python package which provides a modern, open-source pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the python code and exemplify its usefulness and accuracy by analyzing survey data from the Crab Pulsar, NGC 2992, and a previously uncataloged MAXI Transient. The BatAnalysis package allows for $\sim$ 18 years of BAT survey to be used in a systematic way to study a large variety of astrophysical sources.

By implementing a model of primordial dust emission, we predict dust-continuum fluxes for massive galaxy sources similar to those recently detected by JWST at $z \gtrsim 7$. Current upper flux limits, obtained with ALMA for some of these sources, can constrain gas metallicity and dust fraction of the first galaxies. Encouragingly, if assuming expected properties for typical first galaxies (i.e., dust-to-metal mass ratio: $D/M = 5 \times 10^{-3}$, gas metallicity: $Z_{\rm g} = 5 \times 10^{-3}~Z_{\odot}$, star formation efficiency: $\eta = 0.01$), model far-infrared (FIR) fluxes are consistent with current upper flux limits inferred from ALMA bands 6 and 7 ($\lesssim 10^4$ nJy). Such low $D/M$ values and metallicities are in agreement with some scenarios proposed in the literature to explain the non-detection of the FIR dust continuum for high-$z$ JWST galaxy candidates. On the other hand, higher values of model parameters $D/M$ ($\gtrsim 0.06$) and $Z_{\rm g}$ ($\gtrsim 5 \times 10^{-2}~Z_{\odot}$) are ruled out by observational data, unless a higher $\eta$ is assumed. According to our findings, ALMA multi-band observations could constrain the dust chemistry and dust grain size distribution in the early universe. In this context, future observational challenges would involve not only reaching higher FIR sensitivities, but also increasing the wavelength coverage by exploring distinct ALMA bands.

The focus of this investigation is to quantify the conversion of magnetic to thermal energy initiated by a quiet Sun cancellation event and to explore the resulting dynamics from the interaction of the opposite polarity magnetic features. We used imaging spectroscopy in the H$\alpha$ line, along with spectropolarimetry in the \ion{Fe}{I} 6173~{\AA} and \ion{Ca}{II} 8542~{\AA} lines from the Swedish Solar Telescope (SST) to study a reconnection-related cancellation and the appearance of a quiet Sun Ellerman bomb (QSEB). We observed, for the first time, QSEB signature in both the wings and core of the \ion{Fe}{I} 6173~{\AA} line. We also found that, at times, the \ion{Fe}{I} line-core intensity reaches higher values than the quiet Sun continuum intensity. From FIRTEZ-dz inversions of the Stokes profiles in \ion{Fe}{I} and \ion{Ca}{II} lines, we found enhanced temperature, with respect to the quiet Sun values, at the photospheric ($\log\tau_c$ = -1.5; $\sim$1000 K) and lower chromospheric heights ($\log\tau_c$ = -4.5; $\sim$360 K). From the calculation of total magnetic energy and thermal energy within these two layers it was confirmed that the magnetic energy released during the flux cancellation can support heating in the aforesaid height range. Further, the temperature stratification maps enabled us to identify cumulative effects of successive reconnection on temperature pattern, including recurring temperature enhancements. Similarly, Doppler velocity stratification maps revealed impacts on plasma flow pattern, such as a sudden change in the flow direction.

A long-standing scientific puzzle has been to explain the origin of the heaviest elements in the Universe and, more particularly, the production of the elements heavier than iron up to uranium. The rapid neutron capture process (or r-process) is known to synthesize about 50\% of these heavy elements and the long-lived actinides observed in our solar system and so-called metal-poor r-process-enhanced stars. This thesis aims to study the r-process nucleosynthesis and some of the uncertainties that still govern our predictions, namely, nuclear yields and heating rates. Our focus will be on the r-process in neutron star (NS) mergers, which is in the spotlight after recent ``multi-messenger'' observations, including the combined detection of gravitational waves from a NS-NS merger event and its subsequent electromagnetic counterpart. In this work, we base our nucleosynthesis calculations on hydrodynamical simulations of NS merger systems, which estimate the amount of ejected mass and its properties during ejection. Through our r-process calculations, we estimate the composition of this gravitationally unbound material, which can contribute to the r-process enrichment of the Galaxy. This work is divided into two studies, each addressing a remaining open question regarding the r-process nucleosynthesis. First, a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive decay of nuclei synthesized in the dynamical ejecta of NS-NS mergers is presented. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations...

Compact stars (CS) are stellar remnants of massive stars. Inside CSs the density is so high that matter is in subatomic form composed of nucleons. With increase of density of matter towards the centre of the objects other degrees of freedom like hyperons, heavier non-strange baryons, meson condensates may appear. Not only that at higher densities, the nucleons may get decomposed into quarks and form deconfined strange quark matter (SQM). If it is so then CSs may contain SQM in the core surrounded by nucleonic matter forming hybrid stars (HSs). However, the nature and composition of matter inside CSs can only be inferred from the astrophysical observations of these CSs. Recent astrophysical observations in terms of CS mass-radius (M-R) relation and gravitational wave (GW) observation indicate that the matter should be soft in the intermediate density range and stiff enough at higher density range to attain the maximum possible mass above $2~M_\odot$ which is not compatible with pure hadronic equation of states (EOSs). Consequently, we study the HS properties with different models of SQM and find that within vector bag model considering density dependent bag parameter, the model goes well with the astrophysical observations so far.

Disc stars from the Gaia DR3 RVS sample are selected to explore the phase spiral as a function of position in the Galaxy. The data reveal a two-armed phase spiral pattern in the local $z-v_z$ plane inside the solar radius, which appears clearly when colour-coded by $\langle v_R \rangle (z,v_z)$: this is characteristic of a breathing mode that can in principle be produced by in-plane non-axisymmetric perturbations. We note that, on the contrary, the phase spiral pattern becomes single armed when outside the solar radius. When a realistic analytic model with an axisymmetric background potential plus a realistic steadily rotating central bar and 2-armed spiral arms as perturbation is used to perform particle test integrations, the pseudo stars get a prominent spiral pattern in the $\langle v_R \rangle$ map in the $x-y$ plane. Additionally, clear breathing mode evidence at a few $\rm{kms}^{-1}$ level can be seen in the $\langle v_z \rangle$ map on the $x-z$ plane, confirming that such breathing modes are non-negligible in the joint presence of a bar and spiral arms. However, no phase-spiral is perceptible in the $(z, v_z)$ plane. When an initial vertical perturbation is added to all pseudo stars and the same orbital integrating scheme is adopted to carry out the simulation, we find that phase spirals can clearly be seen 500~Myr after the perturbation and gradually disappear inside-out; in that case, they are, however, one-armed. Finally, we show as a proof of concept how a toy model of a time-varying non-axisymmetric in-plane perturbation with varying amplitude and pattern speed can, on the other hand, produce a strong two-armed phase-spiral. We conclude that a time-varying strong internal perturbation together with an external vertical perturbation could perhaps explain the transition between the two-armed and one-armed phase-spirals around the Solar radius.

A first solution for the eclipsing binary TYC 4481-358-1 has been found by combining the results of BVIc-photometry with known stellar models and stellar spectral energy distributions (SEDs). The binary shows total and annular eclipses in a circular 90-days orbit. Masses, radii and effective temperatures have been derived: about 3.01 Msun, 14.29 Rsun and 4950 K for the giant, and about 2.39 Msun, 2.39 Rsun and 9600 K for the dwarf. The peculiar dwarf is discolored and shows a Teff - passband dependence from 9000 K in the Ic to 9800 K in the B band. The giant is in its final stage of core helium burning. The interstellar extinction in our line of sight is considerable (Av about 1.44 mag at an elevated Rv of 3.58).

The TIFR Near Infrared Imaging Camera-II (TIRCAM2) is being used at the Devasthal Optical Telescope (DOT) operated by Aryabhatta Research Institute of Observational Sciences (ARIES), Nainital, Uttarakhand, India. In addition to the normal full frame observations, there has been a requirement for high speed sub-array observations for applications such as lunar occultation and star speckle observations. Fast sub-array modes have been implemented in TIRCAM2 with suitable changes in the camera software at the computer and controller DSP code level. Successful observations have been carried out with the fast sub-array mode of observation.

Minkowski's Object and 'Death Star galaxy' are two of the famous cases of rare instances when a radio jet has been observed to directly hit a neighbouring galaxy. RAD12, the RAD@home citizen science discovery with GMRT being presented here, is not only a new system being added to nearly half a dozen rare cases known so far but also the first case where the neighbouring galaxy is not a minor or dwarf companion but a galaxy bigger than the host of the radio jet. Additionally, the jet appears to be one-sided and the jet after interaction completely stops and forms a bubble inflating laterally which is unlike previous cases of minor deviation or loss of collimation. Since the nature of radio jet-ISM coupling is poorly understood so far, more discovery of objects like RAD12 will be important to the understanding of galaxy evolution through merger and AGN feedback.

Radio observations of ultracool dwarfs, objects comprising brown dwarfs and the very lowest mass stars, have mainly focused on analyzing their light-curve and spectral energy distributions providing valuable insights into their magnetic fields. However, spatially-resolved studies of such magnetospheres have been elusive so far. Radio interferometric observations of the brown dwarf LSR J1835+3259 reveal an extended magnetosphere with a morphology compatible with the presence of a radiation belt, similar to that of Jupiter and Earth, consisting of energetic particles confined via magnetic mirroring. Our finding suggests that radio emitting ultracool dwarfs may behave as scaled up versions of Jupiter, validating the connection between dipole-ordered magnetic fields and the presence of belt-like morphologies and aurorae beyond our Solar System.

Connecting the solar wind observed throughout the heliosphere to its origins in the solar corona is one of the central aims of heliophysics. The variability in the magnetic field, bulk plasma, and heavy ion composition properties of the slow wind are thought to result from magnetic reconnection processes in the solar corona. We identify regions of enhanced variability and composition in the solar wind from 2003 April 15 to May 13 (Carrington Rotation 2002), observed by the Wind and Advanced Composition Explorer spacecraft, and demonstrate their relationship to the Separatrix-Web (S-Web) structures describing the corona's large-scale magnetic topology. There are four pseudostreamer (PS) wind intervals and two helmet streamer (HS) heliospheric current sheet/plasma sheet crossings (and an ICME) which all exhibit enhanced alpha-to-proton ratios and/or elevated ionic charge states of carbon, oxygen, and iron. We apply the magnetic helicity-partial variance of increments ($H_m$-PVI) procedure to identify coherent magnetic structures and quantify their properties during each interval. The mean duration of these structures are $\sim$1 hr in both the HS and PS wind. We find a modest enhancement above the power-law fit to the PVI waiting time distribution in the HS-associated wind at the 1.5-2 hr timescales that is absent from the PS intervals. We discuss our results in context of previous observations of the $\sim$90 min periodic density structures in the slow solar wind, further development of the dynamic S-Web model, and future Parker Solar Probe and Solar Orbiter joint observational campaigns.

We present a numerical study on the stability of all fourth- and fifth-order retrograde mean motion resonances (1/3, 3/1, 1/4, 4/1, 2/3, and 3/2) in the 3-body problem composed of a solar mass star, a Jupiter mass planet, and an additional body with zero mass (elliptic restricted problem) or masses corresponding to either Neptune, Saturn, or Jupiter (planetary problem). The fixed point families exist in all cases and are identified through libration of all resonant angles simultaneously. In addition, configurations with libration of a single resonant angle were also observed. Our results for the elliptic restricted 3-body problem are in agreement with previous studies of retrograde periodic orbits, but we also observe new families not previously reported. Our results regarding stable resonant retrograde configurations in the planetary 3-body problem could be applicable to extra-Solar systems.

In recent decades, machine learning has provided valuable models and algorithms for processing and extracting knowledge from time-series surveys. Different classifiers have been proposed and performed to an excellent standard. Nevertheless, few papers have tackled the data shift problem in labeled training sets, which occurs when there is a mismatch between the data distribution in the training set and the testing set. This drawback can damage the prediction performance in unseen data. Consequently, we propose a scalable and easily adaptable approach based on an informative regularization and an ad-hoc training procedure to mitigate the shift problem during the training of a multi-layer perceptron for RR Lyrae classification. We collect ranges for characteristic features to construct a symbolic representation of prior knowledge, which was used to model the informative regularizer component. Simultaneously, we design a two-step back-propagation algorithm to integrate this knowledge into the neural network, whereby one step is applied in each epoch to minimize classification error, while another is applied to ensure regularization. Our algorithm defines a subset of parameters (a mask) for each loss function. This approach handles the forgetting effect, which stems from a trade-off between these loss functions (learning from data versus learning expert knowledge) during training. Experiments were conducted using recently proposed shifted benchmark sets for RR Lyrae stars, outperforming baseline models by up to 3\% through a more reliable classifier. Our method provides a new path to incorporate knowledge from characteristic features into artificial neural networks to manage the underlying data shift problem.

We use the long-slit spectroscopy mode of the VISIR-NEAR experiment to perform diffraction-limited observations of eight nearby Herbig Ae protoplanetary disks. We extract spectra for various locations along the slit with a spectral resolution of R = 300 and perform a compositional fit at each spatial location using spectral templates of silicates and the four PAH bands. This yields the intensity vs. location profiles of each species. Results. We could obtain spatially-resolved intensity profiles of the PAH emission features in the N-band for five objects (AB Aurigae, HD 97048, HD 100546, HD 163296, and HD 169142). We observe two kinds of PAH emission geometry in our sample: centrally-peaked (HD 97048) and ring-like (AB Aurigae, HD 100546, HD 163296, and potentially HD 169142). Comparing the spatial PAH emission profiles with near-infrared scattered light images, we find a strong correlation in the disk sub-structure but a difference in radial intensity decay rate. The PAH emission shows a less steep decline with distance from the star. Finally, we find a correlation between the presence of (sub-) micron-sized silicate grains leading to the depletion of PAH emission within the inner regions of the disks. In this work, we find the following: (1) PAH emission traces the extent of Herbig Ae disks to a considerable radial distance. (2) The correlation between silicate emission within the inner regions of disks and the depletion of PAH emission can result from dust-mixing and PAH coagulation mechanisms and competition over UV photons. (3) For all objects in our sample, PAHs undergo stochastic heating across the entire spatial extent of the disk and are not saturated. (4) The difference in radial intensity decay rates between the PAHs and scattered-light profiles may be attributed to shadowing and dust-settling effects, which affect the scattering grains more than the PAHs.

Starting with models for the compact object merger event rate, the short-duration Gamma-ray Burst (sGRB) luminosity function, and the Swift/BAT detector, we calculate the observed Swift sGRB rate and its uncertainty. Our probabilistic sGRB world model reproduces the observed number distributions in redshift and flux for 123 Swift/BAT detected sGRBs and can be used to predict joint sGRB/LIGO detection rates. We discuss the dependence of the rate predictions on the model parameters and explore how they vary with increasing experimental sensitivity. In particular, the number of bursts in the LIGO volume depends strongly on the parameters that govern sGRB beaming. Our results suggest that nearby sGRBs should be observed to have broader jets on average ($\theta_{\rm jet}\gtrsim 30$ degrees), as compared to the narrowly-beamed appearance of cosmological sGRBs due to detection selection effect driving observed jet angle. Assuming all sGRBs are due to compact object mergers, within a $D < 200$ Mpc aLIGO volume, we predict $0.18^{+0.19}_{-0.08}$ sGRB/GW associations all-sky per year for on-axis events at Swift sensitivities, increasing to $1.2^{+1.9}_{-0.6}$ with the inclusion of off-axis events. We explore the consistency of our model with GW170817/GRB~170817A in the context of structured jets. Predictions for future experiments are made.

The presence of an extra radio background besides the cosmic microwave background has important implications for the observation of the 21-cm signal during the cosmic Dark Ages, Cosmic Dawn, and epoch of Reionization. The strong absorption trough found in the 21-cm global spectrum measured by the EDGES experiment, which has a much greater depth than the standard model prediction, has drawn great interest to this scenario, but more generally it is still of great interest to consider such a cosmic radio background (CRB) in the early Universe. To be effective in affecting the 21-cm signal at early time, such a radio background must be produced by sources which can emit strong radio signals but modest amount of X-rays, so that the gas is not heated up too early. We investigate the scenario that such a radio background is produced by the primordial black holes (PBHs). For PBH with a single mass, we find that if the PBHs' abundance $\log(f_{\rm PBH})$ (ratio of total PBH mass density to total matter density) and mass satisfy the relation $\log(f_{\rm PBH}) \sim -1.8\log(M_\bullet/{\rm M}_{\odot})-3.5$ for $1\,{\rm M}_\odot \lesssim M_\bullet \lesssim 300 {\rm M}_\odot$, and have jet emission, they can generate a CRB required for reproducing the 21-cm absorption signal seen by the EDGES. The accretion rate can be boosted if the PBHs are surrounded by dark matter halos, which permits lower $f_{\rm PBH}$ value to satisfy the EDGES observation. In the latter scenario, since the accretion rate can evolve rapidly during the Cosmic Dawn, the frequency (redshift) and depth of the absorption trough can determine the mass and abundance of the PBHs simultaneously. For absorption trough redshift $\sim$ 17 and depth $\sim -500$ mK, it corresponds to $M_\bullet \sim 1.05\,{\rm M}_{\odot}$ and $f_{\rm PBH}\sim 1.5\times10^{-4}$.

We present HI-based B- and R-band Tully-Fisher relations (TFRs) and the Baryonic TFR (BTFR) at z=0.2 using direct HI detections from the Blind Ultra-Deep HI Environmental Survey (BUDHIES). Deep photometry from the Isaac Newton Telescope was used for 36 out of 166 HI sources, matching the quality criteria required for a robust TFR analysis. Two velocity definitions at 20% and 50% of the peak flux were measured from the global HI profiles and adopted as proxies for the circular velocities. We compare our results with an identically constructed z=0 TFR from the Ursa Major (UMa) association of galaxies. To ensure an unbiased comparison of the TFR, all the samples were treated identically regarding sample selection and applied corrections. We provide catalogues and an atlas showcasing the properties of the galaxies. Our analysis is focused on the zero points of the TFR and BTFR with their slopes fixed to the z=0 relation. Our main results are: (1) The BUDHIES galaxies show more asymmetric HI profiles with shallower wings compared to the UMa galaxies, which is likely due to the environment in which they reside, (2) The luminosity-based z=0.2 TFRs are brighter and bluer than the z=0 TFRs, even when cluster galaxies are excluded from the BUDHIES sample, (3) The BTFR shows no evolution in its zero point over the past 2.5 billion years and does not significantly change on the inclusion of cluster galaxies, and (4) proper sample selection and consistent corrections are crucial for an unbiased analysis of the evolution of the TFR.

We study varying forms of viscous dark matter and try to address the intriguing tensions of the standard model of cosmology with the recent cosmological data, including the Hubble and $S_8$ tensions. We note that assuming the dark matter viscosity depends on the Hubble parameter, dark matter density, or both, one can improve the statistics. Although the models tend to aggravate the Hubble tension a bit, they tend to reduce the $S_8$ tension, even in comparison with the constant viscosity case.

Extremely variable quasars can also show strong changes in broad line emission strength, and are known as Changing Look Quasars (CLQs). To study the CLQ transition mechanism, we present a pilot sample of CLQs with X-ray observations both in the bright and faint states. From a sample of quasars with bright state archival SDSS spectra and (Chandra or XMM-Newton) X-ray data, we identified five new CLQs via optical spectroscopic follow-up, and then obtained new target-of-opportunity X-ray observations with Chandra. No strong absorption is detected in either the bright or the faint state X-ray spectra. The intrinsic X-ray flux generally tracks the optical variability, and the X-ray power-law slope becomes harder in the faint state. Large amplitude mid-infrared variability is detected in all five CLQs. The changing-obscuration model is not consistent with the observed X-ray spectra and spectral energy distribution changes seen in these CLQs. It is highly likely that the observed changes are due to the changing accretion rate of the supermassive black hole, so the multiwavelength emission varies accordingly, with promising analogies to the accretion states of X-ray binaries.

We investigated 634 crater clusters on Mars detected between 2007 and 2021, which represent more than half of all impacts discovered in this period. Crater clusters form when meteoroids in the 10 kg to 10 ton mass range break-up in Mars' atmosphere to produce a few to a few hundred fragments that hit the ground. The properties of the clusters can inform our understanding of meteoroid properties and the processes that govern their fragmentation. We mapped individual craters $>$1 m within each cluster and defined a range of cluster properties based on the spatial and size distributions of the craters. The large data set, with over eight times more cluster observations than previous work, provides a more robust statistical investigation of crater cluster parameters and their correlations. Trends in size, dispersion and large crater fraction with elevation support weak atmospheric filtering of material. The diversity in the number of individual craters within a cluster, and their size-frequency distributions, may reflect either a diversity in fragmentation style, fragility or internal particle sizes.

In this letter, we explore the quiescent lives of central galaxies using the SAGE galaxy model and Uchuu dark matter simulation. We ask three questions: (1) How much of a galaxy's life is spent in quiescence? (2) How often do galaxies transition off the main sequence? (3) What is the typical duration of a quiescent phase? We find low and high-mass galaxies spend the highest fraction of their lives in quiescence: 45 \pm 19% for log10(Mstar) < 9.0 (3.68 \pm 1.80 Gyr) and 26 \pm 25% for log10(Mstar) > 11.5 (3.46 \pm 3.30 Gyr), falling to 7 \pm 13% for galaxies in-between (0.82 \pm 1.57 Gyr). Low mass galaxies move in and out of quiescence frequently, 2.8 \pm 1.3 times on average, though only for short periods, 1.49 \pm 1.04 Gyr. This can be traced to the influence of supernova feedback on their quite stochastic evolution. Galaxies of higher mass have fewer quiescent periods, ~0.7 \pm 0.9, and their length increases with mass, peaking at 1.97 \pm 2.27 Gyr. However, our high-mass population comprises star-forming and quiescent galaxies with diverging evolutionary paths, so the actual time may be even longer. These high-mass trends are driven by radio mode feedback from supermassive black holes, which, once active, tend to remain active for extended periods.

The redshift drift is computed along light rays propagating through a simulated universe based on the Newtonian N-body simulation code GADGET-2 combined with a perturbed Friedmann-Lemaitre-Robertson-Walker metric in the Newtonian gauge. It is found that the mean redshift drift is equal to the drift of the mean redshift to the precision of the numerical computations and that this is due to a high degree of cancellation between two dominant components of the redshift drift. This result is contrary to earlier findings based on inhomogeneous cosmological models exhibiting cosmic backreaction. \newline\indent For simplicity, the results neglect contributions from optical drift. Based on a study of the redshift drift in a Lemaitre-Tolman-Bondi model, the optical drift effects are estimated to be at most of order 10\% of the redshift drift signal. In addition, it is found that the redshift drift contribution from peculiar acceleration of the emitter is negligible in the simulation setup. However, it is expected that the contribution from peculiar acceleration of the emitter is suppressed in the setup due to low resolution of structures and it is hence expected that this contribution will be larger for real observations.

GRB~210704A is a burst of intermediate duration ($T_{90} \sim 1-4$~s) followed by a fading afterglow and an optical excess that peaked about 7 days after the explosion. Its properties, and in particular those of the excess, do not easily fit into the well established classification scheme of GRBs as being long or short, leaving the nature of its progenitor uncertain. We present multi-wavelength observations of the GRB and its counterpart, observed up to 160 days after the burst. In order to decipher the nature of the progenitor system, we present a detailed analysis of the GRB high-energy properties (duration, spectral lag, and Amati correlation), its environment, and late-time optical excess. We discuss three possible scenarios: a neutron star merger, a collapsing massive star, and an atypical explosion possibly hosted in a cluster of galaxies. We find that traditional kilonova and supernova models do not match well the properties of the optical excess, leaving us with the intriguing suggestion that this event was an exotic high-energy merger.

Hubble constant $H_0$ and weighted amplitude of matter fluctuations $S_8$ determinations are biased to higher and lower values, respectively, in the late Universe with respect to early Universe values inferred by the Planck collaboration within flat $\Lambda$CDM cosmology. If these anomalies are physical, i.e. not due to systematics, they naively suggest that $H_0$ decreases and $S_8$ increases with effective redshift. Here, on the assumption that matter density parameter today $\Omega_{m}$ is a constant, we show that $S_8$ determinations from $f \sigma_8(z)$ constraints increase with effective redshift, thereby providing corroborating support for an $S_8$ discrepancy that is physical in origin. Results here and elsewhere suggest that the $\Lambda$CDM cosmological parameters are redshift dependent. Fitting parameters that evolve with redshift is a recognisable hallmark of model breakdown.

By combining high-sensitivity LOFAR 150MHz, uGMRT 400MHz, GMRT 610MHz, and JVLA 5GHz data in the ELAIS-N1 field, we study the radio spectral properties of radio-detected star-forming galaxies (SFGs) at observer-frame frequencies of 150-5000MHz. We select ~3,500 SFGs that have both LOFAR 150MHz and GMRT 610MHz detections by removing AGN from the two radio samples, and obtain a median spectral index of $\alpha_{150}^{610}=-0.51\pm0.01$ with a scatter of $\sigma=0.2$. Due to the relatively lower sensitivity of uGMRT 400MHz data, we apply a flux cut of $S_{610}>300\mu$Jy and obtain the median spectral indices of $\alpha_{150}^{385}=-0.42^{+0.03}_{-0.02}$, $\alpha_{385}^{610}=-0.44^{+0.03}_{-0.04}$, and $\alpha_{150}^{610}=-0.42^{+0.02}_{-0.01}$ for the sample of 258 SFGs that have detections at these three radio frequencies. The JVLA 5GHz observations only cover the central 0.1deg$^{2}$, where ~100 SFGs are selected, for which we obtain median $\alpha_{610}^{5000}=-1.14^{+0.04}_{-0.05}$, $\alpha_{385}^{5000}=-1.08^{+0.01}_{-0.02}$ and $\alpha_{150}^{5000}=-0.87\pm0.01$. Overall, the results show that the radio spectrum is flatter if we include a lower frequency dataset when measuring the radio spectral index at 150-5000MHz. We study the correlations between radio spectral index and physical properties of radio-selected SFGs and find that, on average, the radio spectrum slightly steepens with increasing stellar mass. However, we only find that the radio spectrum flattens with increasing optical depth at V-band at $\nu<\sim1$GHz. We suggest that spectral ageing due to the energy loss of cosmic ray electrons and thermal free-free absorption could be among the most likely physical mechanisms that drive the two correlations respectively. Both of these could be the physical causes of why the radio spectrum is flatter at low frequency than at high frequency.

The role of supersonic turbulence in structuring the interstellar medium (ISM) remains an unsettled question. Here, this problem is investigated using a newexact law of compressible isothermal hydrodynamic turbulence, which involves two-point correlations in physical space. The new law is shown to have a compact expression that contains a single flux term reminiscent of the incompressible case and a source term with a simple expression whose sign is given by the divergence of the velocity. The law is then used to investigate the properties of such a turbulence at integral Mach number $4$ produced by a massive numerical simulation with a grid resolution of $10,048^3$ points. The flux (resp. source) term was found to have positive (resp. negative) contribution to the total energy cascade rate, which is interpreted as a direct cascade amplified by compression, while their sum is constant in the inertial range. Using a local (in space) analysis it is shown that the source is mainly driven by filamentary structures in which the flux is negligible. Taking positive defined correlations reveals the existence of different turbulent regimes separated by the sonic scale, which determines the scale over which the non-negligible source modifies the scaling of the flux. Our study provides new insight into the dynamics and structures of supersonic interstellar turbulence.

Type Ia Supernovae (SNe Ia) are considered the most reliable \textit{standard candles} and they have played an invaluable role in cosmology since the discovery of the Universe's accelerated expansion. During the last decades, the SNe Ia samples have been improved in number, redshift coverage, calibration methodology, and systematics treatment. These efforts led to the most recent \textit{``Pantheon"} (2018) and \textit{``Pantheon +"} (2022) releases, which enable to constrain cosmological parameters more precisely than previous samples. In this era of precision cosmology, the community strives to find new ways to reduce uncertainties on cosmological parameters. To this end, we start our investigation even from the likelihood assumption of Gaussianity, implicitly used in this domain. Indeed, the usual practise involves constraining parameters through a Gaussian distance moduli likelihood. This method relies on the implicit assumption that the difference between the distance moduli measured and the ones expected from the cosmological model is Gaussianly distributed. In this work, we test this hypothesis for both the \textit{Pantheon} and \textit{Pantheon +} releases. We find that in both cases this requirement is not fulfilled and the actual underlying distributions are a logistic and a Student's t distribution for the \textit{Pantheon} and \textit{Pantheon +} data, respectively. When we apply these new likelihoods fitting a flat $\Lambda$CDM model, we significantly reduce the uncertainties on $\Omega_M$ and $H_0$ of $\sim 40 \%$. This boosts the SNe Ia power in constraining cosmological parameters, thus representing a huge step forward to shed light on the current debated tensions in cosmology.

The zodiacal light (ZL) is sunlight scattered by interplanetary dust (IPD) in the optical wavelengths. The spatial distribution of IPD in the Solar system may hold an important key to understanding the evolution of the Solar system and material transportation within it. The IPD number density can be expressed as $n(r) \sim r^{-\alpha}$, and the result of $\alpha \sim 1.3$ was obtained by the previous observations from the interplanetary space by Helios 1/2 and Pioneer 10/11 in the 1970s and 1980s. However, no direct measurements of $\alpha$ based on the ZL observation from the interplanetary space outside the Earth's orbit have been conducted since then. Here we introduce the initial result of the ZL radial profile at optical wavelengths observed at 0.76-1.06 au by ONC-T with Hayabusa2# mission in 2021-2022. The obtained ZL brightness is well reproduced by the model brightness, but there is a small excess of the observed ZL brightness over the model brightness at around 0.9 au. The obtained radial power-law index is $\alpha = 1.30 \pm 0.08$, which is consistent with the previous results based on the ZL observations. The dominant uncertainty source in $\alpha$ arises from the uncertainty in the Diffuse Galactic Light estimation.

Solar jets are well-collimated plasma ejections in the solar atmosphere. They are prevalent in active regions, the quiet Sun, and even coronal holes. They display a range of temperatures, yet the nature of the cool components has not been fully investigated. In this paper, we show the existence of the precursors and quasi-periodic properties for two chromospheric jets, mainly utilizing the He i 10830 {\AA} narrowband filtergrams taken by the Goode Solar Telescope (GST). The extreme ultraviolet (EUV) counterparts present during the eruption correspond to a blowout jet (jet 1) and a standard jet (jet 2), as observed by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO). The high-resolution He i 10830 {\AA} observation captures a long-lasting precursor for jet 1, signified by a series of cool ejections. They are recurrent jet-like features with a quasi-period of about five minutes. On the other hand, the cool components of jet 2, recurrently accompanied by EUV emissions, present a quasi-periodic behavior with a period of about five minutes. Both the EUV brightening and He i 10830 {\AA} absorption show that there was a precursor for jet 2 that occurred about five minutes before its onset. We propose that the precursor of jet 1 may be the consequence of chromospheric shock waves, since the five-minute oscillation from the photosphere can leak into the chromosphere and develop into shocks. Then, we find that the quasi-periodic behavior of the cool components of jet 2 may be related to magnetic reconnections modulated by the oscillation in the photosphere.

A key first step to constrain the impact of energetic particles in exoplanet atmospheres is to detect the chemical signature of ionisation due to stellar energetic particles and Galactic cosmic rays. We focus on GJ$\,$436, a well-studied M dwarf with a warm Neptune-like exoplanet. We demonstrate how the maximum stellar energetic particle momentum can be estimated from the stellar X-ray luminosity. We model energetic particle transport through the atmosphere of a hypothetical exoplanet at orbital distances between $a=0.01-0.2\,$au from GJ$\,$436, including GJ$\,$436$\,$b's orbital distance (0.028$\,$au). For these distances we find that, at top-of-atmosphere, stellar energetic particles ionise molecular hydrogen at a rate of $\zeta_{\rm StEP,H_2} \sim 4\times10^{-10}-2\times10^{-13}\,\mathrm{s^{-1}}$. In comparison, Galactic cosmic rays alone lead to $\zeta_{\rm GCR, H_2}\sim2\times 10^{-20}-10^{-18} \,\mathrm{s^{-1}}$. At 10au we find that ionisation due to Galactic cosmic rays equals that of stellar energetic particles: $\zeta_{\rm GCR,H_2} = \zeta_{\rm StEP,H_2} \sim 7\times10^{-18}\,\rm{s^{-1}}$ for the top-of-atmosphere ionisation rate. At GJ$\,$436$\,$b's orbital distance, the maximum ion-pair production rate due to stellar energetic particles occurs at pressure $P\sim 10^{-3}\,$bar while Galactic cosmic rays dominate for $P>10^2\,$bar. These high pressures are similar to what is expected for a post-impact early Earth atmosphere. The results presented here will be used to quantify the chemical signatures of energetic particles in warm Neptune-like atmospheres.

This paper presents a theoretical investigation on the polarization and magnetic sensitivity of the near-ultraviolet (near-UV) solar spectral lines of Fe II between 250 and 278 nm. In recent years, UV spectropolarimetry has become key to uncover the magnetism of the upper layers of the solar chromosphere. The unprecedented data obtained by the CLASP2 suborbital space experiment across the Mg II h and k lines around 280 nm are a clear example of the capabilities of near-UV spectropolarimetry for the magnetic field diagnostics throughout the whole solar chromosphere. Recent works have pointed out the possible complementary diagnostic potential of the many Fe II lines in the unexplored spectral region between 250 and 278 nm, but no quantitative analysis of the polarization and magnetic sensitivity of those spectral lines has been carried out yet. To study the polarization signals in these spectral lines, we create a comprehensive atomic model including all the atomic transitions resulting in strong spectral lines. We then study the magnetic sensitivity of the linear and circular polarization profiles in a semi-empirical model representative of the quiet sun. We present a selection of Fe II spectral lines with significant linear and circular polarization signals and evaluate their diagnostic capabilities by studying their formation heights and magnetic sensitivity through the action of the Hanle and Zeeman effects. We conclude that when combined with the CLASP2 spectral region these Fe II lines are of interest for the inference of magnetic fields throughout the solar chromosphere, up to near the base of the corona.

The recent $\sim 4 \, \sigma$ Hubble constant, $H_{0}$, tension is observed between the value of $H_{0}$ from the Cosmic Microwave Background (CMB) and Type Ia Supernovae (SNe Ia). It is a decade since this tension is excruciating the modern astrophysical community. To shed light on this problem is key to consider probes at intermediate redshifts between SNe Ia and CMB and reduce the uncertainty on $H_0$. Toward these goals, we fill the redshift gap by employing Gamma-Ray Bursts (GRBs) and Quasars (QSOs), reaching $z=9.4$ and $z=7.6$, respectively, combined with Baryonic Acoustic Oscillations (BAO) and SNe Ia. To this end, we employ the ``Dainotti GRB 3D relation" among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity, and the ``Risaliti-Lusso" QSO relation between ultraviolet and X-ray luminosities. We inquire the commonly adopted Gaussianity assumption on GRBs, QSOs, and BAO. With the joint sample, we fit the flat $\Lambda$ Cold Dark Matter model with both the Gaussian and the newly discovered likelihoods. We also investigate the impact of the calibration assumed for \textit{Pantheon} and \textit{Pantheon +} SNe Ia on this analysis. Remarkably, we show that only GRBs fulfill the Gaussianity assumption. We achieve small uncertainties on the matter density parameter $\Omega_M$ and $H_0$. We find $H_0$ values compatible within 2 $\sigma$ with the one from the Tip of the Red Giant Branch. Finally, we show that the cosmological results are heavily biased against the arbitrary calibration choice for SNe Ia.

We use optical Integral Field Spectroscopy (IFU) to study the gas emission structure and kinematics in the inner 3.4$\times$4.9 kpc$^2$ region of the galaxy UGC 8782 (3C 293), host of a radio loud Active Galactic Nucleus (AGN). The observations were performed with the GMOS-IFU on the Gemini North telescope, resulting in a spatial resolution of $\sim725$ pc at the distance of the galaxy. While the stars present ordered rotation following the orientation of the large scale disc, the gas shows a disturbed kinematics. The emission-line profiles present two kinematic components: a narrow ($\sigma\lesssim200$ km s$^{-1}$) component associated with the gas in the disc of the galaxy and a broad ($\sigma\gtrsim200$ km s$^{-1}$) component produced by gas outflows. Emission-line ratio diagrams indicate that the gas in the disc is excited by the AGN radiation field, while the emission of the outflow includes additional contribution of shock excitation due to the interaction of the radio jet with the environment gas. Deviations from pure rotation, of up to 30 km s$^{-1}$, are observed in the disc component and likely produced by a previous merger event. The broad component is blueshifted by $\sim150-500$ km s$^{-1}$ relative to the systemic velocity of the galaxy in all locations. We construct radial profiles of the mass outflow rate and kinetic power of the ionized gas outflows, which have the maximum values at $\sim1$ kpc from the nucleus with peak values of $\dot{M}_{\rm out,\Delta R}=0.5\pm0.1$ M$_\odot$ yr$^{-1}$ and $\dot{K}_{\rm out,\Delta R} =$(6.8$\pm$1.1)$\times$10$^{41}$ erg s$^{-1}$. The kinetic coupling efficiency of these outflows is in the range of 1$-$3 per cent, indicating that they could be powerful enough to affect the star formation in the host galaxy as predicted by theoretical simulations.

While it is being conjectured that a chromospheric canopy plays a role in penumbra formation, it has been difficult to find observational evidence of the connectivity between the photosphere and the chromosphere. We investigate the existence of a chromospheric canopy as a necessary condition for the formation of a penumbra and aim to find the origin of the inclined magnetic fields. Spectropolarimetric observations of NOAA AR 12776 from the GRIS@GREGOR instrument were analyzed. Atmospheric parameters were obtained from the deep photospheric Ca I 10839 {\AA} line (VFISV inversion code), the mostly photospheric Si I 10827 {\AA} line (SIR inversion code) and the chromospheric He I 10830 {\AA} triplet (HAZEL inversion code). In the deepest atmospheric layers, we find that the magnetic properties (inclination and field strength distribution) measured on the sunspot sector with fully fledged penumbra are similar to those measured on the sector without penumbra. Yet, in higher layers, magnetic properties are different. In the region showing no penumbra, almost vertical chromospheric magnetic fields are observed. Additionally, thin filamentary structures with a maximum width of 0.1 arcsec are seen in photospheric high-resolution TiO-band images in this region. The existence of a penumbra is found to be discriminated by the conditions in the chromosphere. This indicates that a chromospheric canopy is a necessary condition for the formation of a penumbra. However, our results demonstrate that inclined fields in the chromospheric canopy are not needed for the development of inclined fields in the photosphere. We question the `fallen-magnetic-flux-tubes' penumbra formation scenario and favor a scenario, in which inclined fields emerge from below the surface and are blocked by the overlying chromospheric canopy.

One of the main objectives of stage IV galaxy surveys is to constrain gravity on cosmological scales. To this end, it is crucial to make accurate theoretical predictions in the nonlinear regime of structure formation in order to maximise the scientific return. This is possible at a relatively low computational cost thanks to COLA simulations, an approximate and much faster alternative to full $N$-body simulations. In this thesis, we focus on two modified gravity theories, $f(R)$ and nDGP, and present an analysis of how COLA simulations can be combined with empirical models for the galaxy-halo connection to produce realistic mock galaxy catalogues in modified gravity. We then use the resulting galaxy catalogues to study the effects of modified gravity and validate COLA simulations for several probes of the large-scale structure like the galaxy power spectrum, bispectrum and the void-galaxy cross-correlation function. Finally, we propose the COLA method for the extension of cosmological emulators to modified gravity theories and showcase its potential by training an emulator for the nDGP boost factor of the matter power spectrum, i.e., the ratio of the power spectrum in nDGP gravity with that in general relativity.

Observed first by amateur astronomer Stephen J. O'Meara in the 1970s and then subsequently observed by the Voyager Spacecraft flybys in the early 1980s, it was realised that the 'spokes' flare out like spokes on a bicycle wheel. The observed behaviour of the 'spokes' indicates that they are not governed by gravitational interactions with the planets, moons, or ring material. In 2005 the Cassini probe confirmed that the 'spokes' are likely under the influence of the gas giant's global magnetic field. Here we show that the 'spokes' that appear in Saturn's rings consist of grains of silicates coated in pyrolytic carbon through the process of Chemical Vapour Deposition (CVD). Pyrolytic carbon is a highly diamagnetic substance that can levitate above a sufficiently strong magnetic field. The 'spokes' also consist of ice particles that are diamagnetic as well. The photoelectric effect can be used to explain why the silicates coated in pyrolytic carbon return to the main ring structure when exposed to sunlight of a specific frequency. The pyrolytic carbon grains become paramagnetic when some of the unhybridised 2pz orbitals lose their unpaired delocalised electrons, thus collapsing the pi bond molecular orbital structure. The pyrolytic carbon grains are now attracted to the magnetic field emanating above and below the main ring structure. It is suggested that the 'spokes' in Saturn's B ring are always present and no plasma triggering event is required to increase plasma density. The 'spokes', however, are only visible when a favourable viewing angle is allowed, and their visibility is also dependent on the angle of the sunlight hitting Saturn's rings.

HD165052 is a short-period massive eccentric binary system that undergoes apsidal motion. As the rate of apsidal motion is directly related to the internal structure constants of the binary components, its study allows to get insight into the internal structure of the stars. We use medium- and high-resolution spectroscopic observations of HD165052 to provide constraints on the fundamental properties of the binary system and the evolutionary state of its components. We apply a spectral disentangling code to reconstruct artefact-free spectra of the individual stars and derive the radial velocities (RVs) at the times of the observations. We perform the first analysis of the disentangled spectra with the non-LTE model atmosphere code CMFGEN to determine the stellar properties. We derive the first self-consistent orbital solution of all existing RV data, including those reported in the literature, accounting for apsidal motion. We build, for the very first time, dedicated stellar evolution tracks with the Cl\'es code requesting the theoretical effective temperatures and luminosities to match those obtained from our spectroscopic analysis. The binary system HD165052, consisting of an O6.5V((f)) primary and an O7V((f)) secondary, displays apsidal motion at a rate of (11.30+0.64-0.49)$\deg$yr$^{-1}$. Evolutionary masses are compared to minimum dynamical masses to constrain the orbital inclination. Evolutionary masses Mev,P=24.8$\pm$1.0M$_\odot$ and Mev,S=20.9$\pm$1.0M$_\odot$ and radii Rev,P=7.0+0.5-0.4R$_\odot$ and Rev,S=6.2+0.4-0.3R$_\odot$ are derived, and the inclination is constrained to 22.1$\deg\le i\le 23.3\deg$. Theoretical apsidal motion rates, derived assuming an age of 2.0+/-0.5 Myr for the binary, are in agreement with the observational determination. The agreement with theoretical apsidal motion rates enforces the inferred values of the evolutionary stellar masses and radii.

We report comprehensive spectral and temporal properties of the Be/X-ray binary pulsar SMC X-2 using X-ray observations during the 2015 and 2022 outbursts. The pulse profile of the pulsar is unique and strongly luminosity dependent. It evolves from a broad-humped into a double-peaked profile above luminosity 3$\times$10$^{38}$ ergs s$^{-1}$. The pulse fraction of the pulsar is found to be a linear function of luminosity as well as energy. We also studied the spectral evolution of the source during the latest 2022 outburst with NICER. The observed photon index shows a negative and positive correlation below and above the critical luminosity, respectively, suggesting evidence of spectral transition from the sub-critical to super-critical regime. The broadband spectroscopy of four sets of NuSTAR and XRT/NICER data from both outbursts can be described using a cutoff power-law model with a blackbody component. In addition to the 6.4 keV iron fluorescence line, an absorption-like feature is clearly detected in the spectra. The cyclotron line energy observed during the 2015 outburst is below 29.5 keV, however latest estimates in the 2022 outburst suggest a value of 31.5 keV. Moreover, an increase of 3.4 keV is detected in the cyclotron line energy at equal levels of luminosity observed in 2022 with respect to 2015. The observed cyclotron line energy variation is explored in terms of accretion induced screening mechanism or geometrical variation in line forming region.

Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the HI mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 HI maps are less concentrated than their WHISP counterparts (median $\Delta C\approx0.3$), due in part to central HI deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the HI column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total HI abundance in these datasets as well as the postprocessing method. To fully exploit the synergy between cosmological simulations and upcoming deep HI/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps.

Recent progresses of electronics, essentially due to its miniaturization, are opening new fields that were just dreamed of, notably in astronomy. At start in paragraph 3, we introduce the time variation of images expressing the dual nature of the optical signal (ZO) and we expose several useful applications where the optical signal variations are not faster than CCD. However we prefered to initiate the article with a deeper question posed inadvertently in paragraph 2: what causes the rapid, well timed and regular variation of the signals induced in our test setup, which we see in Fig. 1. The answer proposed are two causes: one is a light photon acting indirectly through the induction of a large number of secondary electrons (paragraph 2), the other are the RF photons (subliminal, but acting directly) as detailed in paragraph 4. For both light and RF, using a sum of induced currents instead of a single photon quadri-vector transform the case.

Small Solar System bodies are pristine remnants of Solar System formation, which provide valuable insights for planetary science and astronomy. Their discovery and cataloging also have strong practical implications to life on Earth as the nearest asteroids could pose a serious impact threat. Concurrently with dedicated observational projects, searches for small bodies have been performed on numerous archival data sets from different facilities. Here, we present a framework to increase the scientific return of an exoplanet transit-search survey by recovering serendipitous detections of small bodies in its daily and archival data using a GPU-based synthetic tracking algorithm. As a proof of concept, we analysed $12 \times 12 \mathrm{arcmin^2}$ sky fields observed by the 1-m telescopes of the SPECULOOS survey. We analysed 90 sky fields distributed uniformly across the sky as part of the daily search for small bodies and 21 archival fields located within 5 deg from the ecliptic plane as part of the archival search (4.4 deg$^2$ in total). Overall, we identified 400 known objects of different dynamical classes from Inner Main-belt Asteroids to Jupiter Trojans and 43 potentially new small bodies with no priors on their motion. We were able to reach limiting magnitude for unknown objects of $V$=23.8 mag, and a retrieval rate of $\sim$80% for objects with $V<$22 mag and $V<$23.5 mag for the daily and archival searches, respectively. SPECULOOS and similar exoplanet surveys can thus serve as pencil-beam surveys for small bodies and probe parameter space beyond $V$=22 mag.

Surveys of star-forming regions reveal that the dust mass of protoplanetary discs decreases by several orders of magnitude on a timescale of a few million years. This decrease in the mass budget of solids is likely due to the gas-drag-induced radial drift of mm-sized solids, called pebbles. However, quantifying the evolution of this dust component in young stellar clusters is difficult due to the inherent large spread in stellar masses and formation times. Therefore, we aim to model the collective evolution of a cluster to investigate the effectiveness of radial drift in clearing the discs of mm-sized particles. We use a protoplanetary disc model that numerically solves for disc formation, and the viscous evolution and photoevaporative clearing of the gas component, while also including the drift of particles limited in size by fragmentation. We find that discs are born with dust masses between 50 Earth masses and 1000 Earth masses, for stars with, respectively, masses between 0.1 solar masses and 1 solar masses. The majority of this initial dust reservoir is typically lost through drift before photoevaporation opens a gap in the gas disc for models both with and without strong X-ray-driven mass loss rates. We conclude that the decrease in time of the mass locked in fragmentation-limited pebbles is consistent with the evolution of dust masses and ages inferred from nearby star-forming regions when assuming viscous evolution rates corresponding to mean gas disc lifetimes between 3 Myr and 8 Myr.

We employ a machine learning (ML) algorithm to analyze cosmological background data and the linear red-shift space distortion (RSD) data in a model-independent way, with specific focus on the Hubble expansion rate and the growth of large-scale structure. We find strong evidence that the natural enhancement in the Hubble parameter at low redshifts is due to the underlying phantom nature of dark energy, rather than low matter density. As for the RSD data, we find a higher value of $\sigma^8_{(0)}$ which is consistent with CMB's predictions, but the outcome of low matter density leads to unresolved tension. This might point towards a new physical phenomenon at the perturbative level in the low redshift regime. From a statistical perspective, we have demonstrated that our results hold greater preference compared to those obtained by the $\Lambda$CDM model.

Using resolved optical stellar photometry from the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey, we measured the star formation history (SFH) near the position of 85 supernova remnants (SNRs) in M33. We constrained the progenitor masses for 60 of these SNRs, finding the remaining 25 remnants had no local SF in the last 56 Myr consistent with core-collapse SNe (CCSNe), making them potential Type Ia candidates. We then infer a progenitor mass distribution from the age distribution, assuming single star evolution. We find that the progenitor mass distribution is consistent with being drawn from a power-law with an index of $-2.9^{+1.2}_{-1.0}$. Additionally, we infer a minimum progenitor mass of $7.1^{+0.1}_{-0.2}$ $M_{\odot}$ from this sample, consistent with several previous studies, providing further evidence that stars with ages older than the lifetimes of single 8 $M_{\odot}$ stars are producing supernovae.

We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform simulations for different grid resolutions (768$^3$, 1536$^3$ and 2880$^3$), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial setup increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\,H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers.

An increasing number of Near Earth Asteroids (NEAs) in the range of a few hundred meters to a few kilometres in size have relatively high spin rates, from less than 4 h, down to $\sim$2.2 h, depending on spectral type. For some of these bodies, local acceleration near the equator may be directed outwards so that lift off of near-equatorial material is possible. In particular, this may be the case for asteroid Didymos, the primary of the (65803) Didymos binary system, which is the target of the DART (NASA) and Hera (ESA) space missions. The study of the dynamics of particles in such an environment has been carried out -- in the frame of the Hera mission and the EC-H2020 NEO-MAPP project -- according to the available shape model, known physical parameters and orbital information available before the DART impact. The presence of orbiting particles in the system is likely for most of the estimated range of values for mass and volume. The spatial mass density of ejected material is calculated for different particle sizes and at different heliocentric orbit epochs, revealing that large particles dominate the density distribution and that small particle abundance depends on observation epoch. Estimates of take off and landing areas on Didymos are also reported. Available estimates of the system mass and primary extents, after the DART mission, confirm that the main conclusions of this study are valid in the context of current knowledge.

Cosmic Dawn (CD) and Epoch of Reionization (EoR) are epochs of the universe which host invaluable information about the cosmology and astrophysics of X-ray heating and hydrogen reionization. Radio interferometric observations of the 21-cm line at high redshifts have the potential to revolutionize our understanding of the universe during this time. However, modeling the evolution of these epochs is particularly challenging due to the complex interplay of many physical processes. This makes it difficult to perform the conventional statistical analysis using the likelihood-based Markov-Chain Monte Carlo (MCMC) methods, which scales poorly with the dimensionality of the parameter space. In this paper, we show how the Simulation-Based Inference (SBI) through Marginal Neural Ratio Estimation (MNRE) provides a step towards evading these issues. We use 21cmFAST to model the 21-cm power spectrum during CD-EoR with a six-dimensional parameter space. With the expected thermal noise from the Square Kilometre Array (SKA), we are able to accurately recover the posterior distribution for the parameters of our model with an order of magnitude fewer simulations than the conventional likelihood-based methods. We further show how the same training dataset can be utilized to investigate the sensitivity of the model parameters over different redshifts. Our results support that such efficient and scalable inference techniques enable us to significantly extend the modeling complexity beyond what is currently achievable with conventional MCMC methods.

We have proposed a novel application for cosmic-ray muography, called Magic-{\mu}, which is short for Magnetic field Imaging by Cosmic-ray Muons. The general goal of Magic-{\mu} is to detect the presence of a magnetic field or magnetic flux density whose three-dimensional distribution is unknown. Depending on the application, Magic-{\mu} can have three detection modes. The first is "magnetic field imaging," which is detecting the presence of a magnetic field in specific voxels within a region of space. The other two modes, transmission and deflection, aim not only to detect the presence but also to measure the flux density of the magnetic field. We have performed a feasibility study using the PHITS Monte Carlo simulation code, for strong and weak magnetic fields. In this paper, we first give an overview of the concept and basic principles of magnetic field muography. Then, the results of the feasibility study on magnetic field imaging for a strong magnetic field (more than 500 mT) are presented.

The mechanism for generating neutrino masses remains a puzzle in particle physics. If neutrino masses follow from a Dirac mass term, then neutrino states exist with opposite chirality compared to their weakly-interacting counterparts. These inactive states do not interact with their active counterparts at measurable scales in the standard model. However, the existence of these states can have implications for cosmology as they contribute to the radiation energy density at early times, and the matter energy density at late times. How Dirac neutrinos may populate thermal states via an anomalous magnetic moment operator is the focus of this work. A class of models where all neutrinos have a magnetic moment independent of flavor or chirality is considered. Subsequently, the cross sections for neutrinos scattering on background plasma particles are calculated so that the relic inactive neutrino energy is derived as a function of plasma temperature. To do so, one needs cross sections for scattering on all electrically charged standard-model particles. Therefore, the scattering cross section between a neutrino and $W$-boson via the magnetic moment vertex is derived. Current measurements put a constraint on the size of the neutrino magnetic moment from the cosmological parameter $N_{\rm eff}$ and light-element primordial abundances. Finally, how the extra Dirac states contribute to the matter energy density at late times is investigated by examining neutrino free-streaming.

Gravitational vacuum condensate stars, also known as gravastars, have been proposed as an alternative to black holes. Their interior contains a perfect fluid with an equation of state akin to that of a cosmological constant. For this reason, they have recently been considered as a possible astrophysical source of dark energy. In this work we argue that gravitational vacuum condensate stars cannot be the source of dark energy and highlight that a direct coupling of their mass to the dynamics of the Universe would lead to an additional velocity dependent acceleration, damping their motion with respect to the cosmological frame. We briefly discuss the potential impact of this additional acceleration in the context of a recent proposal that the observed mass growth of compact objects at the core of elliptical galaxies might result from such a cosmological coupling.

This work generalizes the discrete implicit Monte-Carlo (DIMC) method for modeling the radiative transfer equation from a gray treatment to an frequency-dependent one. The classic implicit Monte-Carlo (IMC) algorithm, that has been used for several decades, suffers from a well-known numerical problem, called teleportation, where the photons might propagate faster than the exact solution due to the finite size of the spatial and temporal resolution. The Semi-analog Monte-Carlo algorithm proposed the use of two kinds of particles, photons and material particles that are born when a photon is absorbed. The material particle can `propagate' only by transforming into a photon, due to black-body emission. While this algorithm produces a teleportation-free result, it is noisier results compared to IMC due to the discrete nature of the absorption-emission process. In a previous work [Steinberg and Heizler, ApJS, 258:14 (2022)], proposed a gray version of DIMC, that makes use of two kinds of particles, and therefore has teleportation-free results, but also uses the continuous absorption algorithm of IMC, yielding smoother results. This work is a direct frequency-dependent (energy-dependent) generalization of the DIMC algorithm. We find in several one and two dimensional benchmarks, that the new frequency-dependent DIMC algorithm yields teleportation-free results on one hand, and smooth results with IMC-like noise level.

A non-minimal approximation for effective masses of light and heavy neutrinos in the framework of a type-I seesaw mechanism with three generations of sterile Majorana neutrinos which recover the symmetry between quarks and leptons is considered. The main results are: (a) the next-order corrections to the effective mass matrix of heavy neutrinos due to terms O({\theta}MD) are obtained, which modify the commonly used representation for the effective mass (MD is a Dirac neutrino mass when the electroweak symmetry is spontaneously broken); and (b) the general form of the mixing matrix is found in non-minimal approximation parametrized by a complex 3x3 matrix satisfying a nontrivial constraint. Numerical analysis within the {\nu}MSM framework demonstrates the very small effect of new contributions of direct collider observables as opposed to their possible significance for cosmological models.

We investigate non-radial oscillations of pure and hybrid neutron stars, employing equations of state of nuclear matter from Brueckner-Hartree-Fock theory, and of quark matter from the Dyson-Schwinger quark model, performing a Gibbs construction for the mixed phase in hybrid stars. Characteristic differences between neutron-star and hybrid-star $g_1$-mode oscillation frequencies, damping times, and gravitational wave strains are pointed out. Prospects of observations are also discussed.

We investigate the effect of density perturbations and local anisotropy on the stability of stellar matter structures in general relativity using the concept of cracking. Adopting a core-envelope model of a super-dense star, we examine the properties and stability conditions by introducing anisotropic pressure to the envelope region. Furthermore, we propose self-bound compact stars with an anisotropic envelope as a potential progenitor for starquakes. We show how the difference between sound propagation in radial and tangential directions would be used to identify potentially stable regions within a configuration. Due to an increase in the anisotropic parameter, strain energy accumulates in the envelope region and becomes a potential candidate for building-up quake like situation. This stress-energy stored in the envelope region that would be released during a starquake of a self-bound compact star is computed as a function of the magnitude of anisotropy at the core-envelope boundary. Numerical studies for spherically asymmetric compact stars indicate that the stress-energy can be as high as $10^{50}$ erg if the tangential pressure is slightly more significant than the radial pressure. It is happened to be of the same order as the energy associated with giant $\gamma$-ray bursts. Thus, the present study will be useful for the correlation studies between starquakes and GRBs.

We investigate the evolution of the general relativistic phase of an electromagnetic (EM) wave as it propagates in a weak gravitational field. For that, we introduce coordinate systems used in the Earth's vicinity, the relevant coordinate transformations, and discuss the transformations between proper and coordinate times. Then, by treating the source as an isolated, weakly aspherical gravitating body, we represent its external gravitational field using Cartesian symmetric trace-free (STF) mass multipole moments. The light propagation equation is solvable along the trajectory of a light ray to all STF orders l. We express the STF moments via spherical harmonic coefficients of various degree and order, C_lk, S_lk. Although we focus primarily on the quadrupole (l=2), octupole (l=3), and hexadecapole (l=4) cases, our approach is valid to all orders. The result is the gravitational phase shift expressed in terms of the spherical harmonics. These results are new, useful, and important as increasingly accurate clocks are used for many practical and scientific applications, including space-based time and frequency transfers, relativistic geodesy, and navigation. We also consider tidal contributions and contributions due to the Earth's rotation. We estimate the characteristic numerical magnitudes of each term of the resulting overall gravitational phase shift.