Papers for Thursday, Oct 28 2021

Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$\alpha$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$\alpha$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$) $z\sim 6$ quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within $3$ pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the $\textit{Cosmic Reionization on Computers}$ simulation, and find a good agreement between $1.5$ to $3$pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter $\sigma_8=0.6\pm0.3$. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than $>200$ known quasars at $z\gtrsim 6$, this method will allow us in the future to tighten the constraint on $\sigma_8$ to the percent level. In the region closer to the quasar within $1.5$ pMpc, we find the density is higher than predicted in the simulation by $1.23 \pm 0.17$, suggesting the typical host dark matter halo mass of a bright quasar ($M_{\rm 1450}<-26.5$) at $z\sim 6$ is $\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}$.

Determining the origin of turbulence in galaxy clusters, and quantifying its transport of heat, is an outstanding problem, with implications for our understanding of their thermodynamic history and structure. As the dilute plasma of the intracluster medium (ICM) is magnetized, heat and momentum travel preferentially along magnetic field lines. This anisotropy triggers a class of buoyancy instabilities that destabilize the ICM, and whose turbulent motions can augment or impede heat transport. We focus on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters. We aim to take a fresh look at the problem and construct a general theory that explains the MTI saturation mechanism and provides scalings and estimates for the turbulent kinetic energy, magnetic energy, and heat flux. We simulate MTI turbulence with a Boussinesq code, SNOOPY, which, in contrast to previous work, allows us to perform an extensive sampling of the parameter space. In two dimensions the saturation mechanism involves an inverse cascade carrying kinetic energy from the short MTI injection scales to larger scales, where it is arrested by the stable entropy stratification; at a characteristic "buoyancy scale", the energy is dumped into large-scale g-modes, which subsequently dissipate. Consequently, the entropy stratification sets an upper limit on the size and strength of turbulent eddies. Meanwhile, the MTI conveys a substantial fraction of heat, despite the tangled geometry of the magnetic field. In a companion paper, these results are extended to three-dimensional flows, and compared to real cluster observations.

The formation of the first galaxies during cosmic dawn and reionization (at redshifts $z=5-30$), triggered the last major phase transition of our universe, as hydrogen evolved from cold and neutral to hot and ionized. The 21-cm line of neutral hydrogen will soon allow us to map these cosmic milestones and study the galaxies that drove them. To aid in interpreting these observations, we upgrade the public code 21cmFAST, improving the treatment of feedback in molecular-cooling galaxies. We introduce a new, flexible parametrization of the additive feedback from: (i) an inhomogeneous, $H_2$-dissociating (Lyman-Werner; LW) background; and (ii) dark matter -- baryon relative velocities. We demonstrate that our flexible model can recover results from recent, small-scale hydrodynamical simulations. We perform a large (1.5 comoving Gpc on a side), "best-guess" simulation as the 2021 installment of the Evolution of 21-cm Structure (EOS) project. This improves on the previous EOS release by using an updated galaxy model that reproduces the observed UV luminosity functions (UVLFs), and by including an additional population of molecular-cooling galaxies. The resulting 21-cm global signal and power spectrum are significantly weaker than in the 2016 EOS releases, due to a more rapid evolution of the star-formation rate density required to match the UVLFs. Nevertheless, we forecast high signal-to-noise detections for both HERA and the SKA. We demonstrate how the stellar-to-halo mass relation of the unseen, first galaxies can be inferred from the evolution of 21-cm fluctuations. Finally, we show that the spatial modulation of X-ray heating due to the relative velocities provides a unique acoustic signature that is detectable at $z \approx 10-15$ in our fiducial model. Ours are the first public simulations with joint inhomogeneous LW and relative-velocity feedback across cosmic dawn and reionization.

A simple model for star formation based on supernova (SN) feedback and gravitational heating via the collapse of perturbations in gravitationally unstable disks reproduces the Schmidt-Kennicutt relation between the star formation rate (SFR) per unit area, $\Sigma_{SFR}$, and the gas surface density, $\Sigma_g$, remarkably well. The gas velocity dispersion, $\sigma_g$, is derived self-consistently in conjunction with $\Sigma_{SFR}$ and is found to match the observations. Gravitational instability triggers {``gravito-turbulence"} at the scale of the least stable perturbation mode, boosting $\sigma_g$ at $\Sigma_g> \, \Sigma_g^\textrm{thr}=50\, {\rm M}_\odot\, {\rm pc}^{-2}$, and contributing to the pressure needed to carry the disk weight vertically. $\Sigma_{SFR}$ is reduced to the observed level at $ \Sigma_g > \Sigma_g^\textrm{thr}$, whereas at lower surface densities, SN feedback is the prevailing energy source. Our proposed star formation recipes require efficiencies of order 1\%, and the Toomre parameter, $Q$, for the joint gaseous and stellar disk is predicted to be close to the critical value for marginal stability for $\Sigma_g< \, \Sigma_g^\textrm{thr}$, spreading to lower values and larger gas velocity dispersion at higher $\Sigma_g$.

PLATO (PLAnetary Transits and Oscillations of stars) is an ESA M-class satellite planned for launch by end 2026 and dedicated to the wide-field search of transiting planets around bright and nearby stars, with a strong focus on discovering habitable rocky planets hosted by solar-like stars. The choice of the fields to be pointed at is a crucial task since it has a direct impact on the scientific return of the mission. In this paper we describe and discuss the formal requirements and the key scientific prioritization criteria that have to be taken into account in the Long-duration Observation Phase (LOP) field selection, and apply a quantitative metric to guide us in this complex optimization process. We identify two provisional LOP fields, one for each hemisphere (LOPS1, LOPN1), and discuss their properties and stellar content. While additional fine-tuning shall be applied to LOP selection before the definitive choice (to be made two years before launch), we expect their position will not move by more than a few degrees with respect to what is proposed in this paper.

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high angular resolution (sub-arcsecond) observations of $^13$CH$_{3}$OH $\nu=0$ ($\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$) and CH$_3$CN $\nu_8=1$ ($\sim$0.2$^{\prime\prime}$ and $\sim$0.9$^{\prime\prime}$) line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) to investigate Orion KL's structure on small spatial scales (${\le}350$ au). Gas kinematics, excitation temperatures, and column densities were derived from the molecular emission via a pixel-by-pixel spectral line fitting of the image cubes, enabling us to examine the small-scale variation of these parameters. Sub-regions of the Hot Core have a higher excitation temperature in a 0.2$^{\prime\prime}$ beam than a 0.9$^{\prime\prime}$ beam, indicative of possible internal sources of heating. Furthermore, the velocity field includes a bipolar ${\sim}7{-}8$ km s$^{-1}$ feature with a southeast-northwest orientation against the surrounding ${\sim}4{-}5$ km s$^{-1}$ velocity field, which may be due to an outflow. We also find evidence of a possible source of internal heating toward the Northwest Clump, since the excitation temperature there is higher in a smaller beam versus a larger beam. Finally, the region southwest of the Hot Core (Hot Core-SW) presents itself as a particularly heterogeneous region bridging the Hot Core and Compact Ridge. Additional studies to identify the (hidden) sources of luminosity and heating within Orion KL are necessary to better understand the nebula and its chemistry.

The bispectrum is the leading non-Gaussian statistic in large-scale structure, carrying valuable information on cosmology that is complementary to the power spectrum. To access this information, we need to model the bispectrum in the weakly non-linear regime. In this work we present the first two-loop, i.e., next-to-next-to-leading order perturbative description of the bispectrum within an effective field theory (EFT) framework. Using an analytic expansion of the perturbative kernels up to $F_6$ we derive a renormalized bispectrum that is demonstrated to be independent of the UV cutoff. We show that the EFT parameters associated with the four independent second-order EFT operators known from the one-loop bispectrum are sufficient to absorb the UV sensitivity of the two-loop contributions in the double-hard region. In addition, we employ a simplified treatment of the single-hard region, introducing one extra EFT parameter at two-loop order. We compare our results to N-body simulations using the realization-based grid-PT method and find good agreement within the expected range, as well as consistent values for the EFT parameters. The two-loop terms start to become relevant at $k\approx 0.07h~\mathrm{Mpc}^{-1}$. The range of wavenumbers with percent-level agreement, independently of the shape, extends from $0.08h~\mathrm{Mpc}^{-1}$ to $0.15h~\mathrm{Mpc}^{-1}$ when going from one to two loops at $z=0$. In addition, we quantify the impact of using exact instead of Einstein-de-Sitter kernels for the one-loop bispectrum, and discuss in how far their impact can be absorbed into a shift of the EFT parameters.

We investigate the properties of the dark matter (DM) halo surrounding the nearby galaxy group NGC 1600. Through the use of deep (252 ks) Chandra observations and 64.3 ks of XMM-Newton observations, we construct surface brightness profiles in multiple energy bands in order to perform hydrostatic equilibrium analysis of the hot plasma within NGC 1600. Regardless of the DM model profile assumed, we measure a halo concentration (c$_{200}$) that is an extreme, positive outlier of the $\Lambda$CDM c$_{200}$-M$_{200}$ relation. For a typical NFW DM profile, we measure c$_{200}\!=\!26.7\pm1.4$ and M$_{200}\!=\!(2.0\pm0.2)\times10^{13}$ M$_\odot$; assuming a similar halo mass, the average concentration expected is c$_{200}=6-7$ for the theoretical $\Lambda$CDM c-M relation. Such a high concentration is similar to that of well-known fossil groups MRK 1216 and NGC 6482. While NGC 1600 exhibits some properties of a fossil group, it fails to meet the X-ray luminosity threshold of L$_X>5\times10^{41}$ erg s$^{-1}$. Whether or not it is considered a fossil group, the high concentration value makes it part of a select group of galaxy groups.

Centaurs with high orbital inclinations and perihelia (i > 60 degrees; q >= 15 au) are a small group of poorly understood minor planets that are predicted to enter the giant planet region of the Solar System from the inner Oort Cloud. As such they are one of the few samples of relatively unaltered Oort Cloud material that can currently be directly observed. Here we present two new reflectance spectra of one of the largest of these objects, 2012 DR30, in order to constrain its color and surface composition. Contrary to reports that 2012 DR30 has variable optical color, we find that consistent measurements of its spectral gradient from most new and published datasets at 0.55-0.8 micron agree with a spectral gradient of S ~ 10+/-1 %/0.1 micron within their uncertainties. The spectral variability of 2012 DR30 at Near-UV/blue and Near-Infrared wavelengths, however, is still relatively unconstrained; self-consistent rotationally resolved followup observations are needed to characterise any spectral variation in those regions. We tentatively confirm previous detections of water ice on the surface of 2012 DR30 , and also consistently observe a steady steepening of the gradient of its spectrum from wavelengths of around 0.6 micron towards Near-UV wavelengths. Plausible surface materials responsible for the observed reddening may include ferric oxides contained within phyllosilicates, and aromatic refractory organics.

We report the discovery of a velocity coherent, kpc-scale molecular structure towards the Galactic center region with an angular extent of 30deg and an aspect ratio of 60:1. The kinematic distance of the CO structure ranges between 4.4 to 6.5 kpc. Analysis of the velocity data and comparison with the existing spiral arm models support that a major portion of this structure is either a sub-branch of the Norma arm or an inter-arm giant molecular filament, likely to be a kpc-scale feather (or spur) of the Milky Way, similar to those observed in nearby spiral galaxies. The filamentary cloud is at least 2.0 kpc in extent, considering the uncertainties in the kinematic distances, and it could be as long as 4 kpc. The vertical distribution of this highly elongated structure reveals a pattern similar to that of a sinusoidal wave. The exact mechanisms responsible for the origin of such a kpc-scale filament and its wavy morphology remains unclear. The distinct wave-like shape and its peculiar orientation makes this cloud, named as the Gangotri wave, one of the largest and most intriguing structures identified in the Milky Way.

Context. RR Lyrae stars are useful standard candles allowing one to derive accurate distances for old star clusters. Based on the recent catalogues from OGLE-IV and Gaia Early Data Release 3 (EDR3), the distances can be improved for a few bulge globular clusters. Aims. The aim of this work is to derive an accurate distance for the following six moderately metal-poor, relatively high-reddening bulge globular clusters: NGC 6266, NGC 6441, NGC 6626, NGC 6638, NGC 6642, and NGC 6717. Methods. We combined newly available OGLE-IV catalogues of variable stars containing mean I magnitudes, with Clement's previous catalogues containing mean V magnitudes, and with precise proper motions from Gaia EDR3. Astrometric membership probabilities were computed for each RR Lyrae, in order to select those compatible with the cluster proper motions. Applying luminosity-metallicity relations derived from BaSTI $\alpha$-enhanced models (He-enhanced for NGC 6441 and canonical He for the other clusters), we updated the distances with relatively low uncertainties. Results. Distances were derived with the I and V bands, with a $5-8\%$ precision. We obtained 6.6 kpc, 13.1 kpc, 5.6 kpc, 9.6 kpc, 8.2 kpc, and 7.3 kpc for NGC 6266, NGC 6441, NGC 6626, NGC 6638, NGC 6642, and NGC 6717, respectively. The results are in excellent agreement with the literature for all sample clusters, considering the uncertainties. Conclusions. The present method of distance derivation, based on recent data of member RR Lyrae stars, updated BaSTI models, and robust statistical methods, proved to be consistent. A larger sample of clusters will be investigated in a future work.

Increasing main sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas pressure inversions in 1D models. Additionally, recent photometric surveys (e.g. TESS) have detected excess broadband low frequency variability in power spectra of OB star lightcurves, called stochastic low frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35$\,$M$_\odot$ star envelopes (the outer $\approx$15$\%$ of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron opacity peak convection zone (FeCZ), obliterating any ``quiet" part of the near-surface region and leading to velocities at the photosphere of 10-100$\,$km$\,$s$^{-1}$, directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ ($\approx$3-7$\,$days$^{-1}$). These simulations, which have no free parameters, are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main sequence O star envelopes exhibit SLFV originating from the FeCZ.

We use helioseismic data obtained over two solar cycles to determine whether there are changes in the near-surface shear layer (NSSL). We examine this by determining the radial gradient of the solar rotation rate. The radial gradient itself shows a solar-cycle dependence, and the changes are more pronounced in the active latitudes than at adjoining higher latitudes; results at the highest latitudes (greater than about70 degrees) are unreliable. The pattern changes with depth, even within the NSSL. We find that the near-surface shear layer is deeper at lower latitudes than at high latitudes and that the extent of the layer also shows a small solar-cycle related change.

In this work, we compute the contribution from clusters of galaxies to the diffuse neutrino background. Clusters of galaxies can potentially produce cosmic rays (CRs) up to very-high energies via large-scale shocks and turbulent acceleration. Due to their unique magnetic-field configuration, CRs with energy $\leq 10^{17}$ eV can be trapped within these structures over cosmological time scales, and generate secondary particles, including neutrinos and gamma rays, through interactions with the background gas and photons. We employ three-dimensional cosmological magnetohydrodynamical simulations of structure formation to model the turbulent intergalactic medium. We use the distribution of clusters within this cosmological volume to extract the properties of this population. We propagate CRs in this environment using multi-dimensional Monte Carlo simulations across different redshifts (from $z \sim 5 \; \text{to} \; z =0$), considering all relevant photohadronic, photonuclear, and hadronuclear interactions. We also include the cosmological evolution of the CR sources. We find that, for CRs injected with a spectral index $1.5 - 2.7$ and cutoff energy $E_{max} = 10^{16} - 10^{17}$ eV, clusters contribute to a substantial fraction to the diffuse flux observed by the IceCube Neutrino Observatory, and most of the contribution comes from clusters with $M > 10^{14} \; M_{\odot}$ and redshift $z < 0.3$.

CMB photons redshift and blueshift as they move through gravitational potentials $\Phi$ while propagating across the Universe. If the potential is not constant in time, the photons will pick up a net redshift or blueshift, known as the Integrated Sachs-Wolfe (ISW) effect. In the $z \ll 1000$ universe, $\dot{\Phi}$ is nonzero on large scales when the Universe transitions from matter to dark energy domination. This effect is only detectable in cross-correlation with large-scale structure at $z \sim 1$. In this paper we present a 3.2$\sigma$ detection of the ISW effect using cross-correlations between unWISE infrared galaxies and Planck CMB temperature maps. We use 3 tomographic galaxy samples spanning $0 < z < 2$, allowing us to fully probe the dark energy domination era and the transition into matter domination. This measurement is consistent with $\Lambda$CDM ($A_{\rm ISW} = 0.96 \pm 0.30$). We study constraints on a particular class of dynamical dark energy models (where the dark energy equation of state is different in matter and dark energy domination), finding that unWISE-ISW improves constraints from type Ia supernovae due to improved constraints on the time evolution of dark energy. When combining with BAO measurements, we obtain the tightest constraints on specific dynamical dark energy models. In the context of a phenomenological model for freezing quintessence, the Mocker model, we constrain the dark energy density within 10% at $z < 2$ using ISW, BAO and supernovae. Moreover, the ISW measurement itself provides an important independent check when relaxing assumptions about the theory of gravity, as it is sensitive to the gravitational potential rather than the expansion history.

The Pierre Auger Observatory is the largest facility in the world to study ultra-high-energy cosmic rays. It has a hybrid detection technique that combines the observation of the longitudinal development of extensive air showers and the measurement of their particles at the ground. This capability has opened the possibility to probe hadronic interactions taking place at energies well beyond those accessible by human-made accelerators. In this report, we present a selection of the latest results on hadronic interactions with measurements from the Pierre Auger Observatory. These data span over three decades in energy, showing the tension between data from the muon component of air showers and predictions based on the most updated hadronic interaction models.

Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere.

Any perturbation to a disc galaxy that creates a misalignment between the planes of the inner and outer disc, will excite a slowly evolving bending wave in the outer disc. The torque from the stiff inner disc drives a retrograde, leading-spiral bending wave that grows in amplitude as it propagates outward over a period of several Gyr. This behaviour creates warps that obey the rules established from observations, and operates no matter what the original cause of the misalignment between the inner and outer disc. The part of the disc left behind by the outwardly propagating wave is brought into alignment with the inner disc. Here we confirm that mild warps in simulations of disc galaxies can be excited by shot noise in the halo, as was recently reported. We show that the quadrupole component of the noise creates disc distortions most effectively. Bending waves caused by shot noise in carefully constructed equilibrium simulations of isolated galaxies are far too mild to be observable, but perturbations from halo substructure and galaxy assembly must excite larger amplitude bending waves in real galaxies.

In this review/tutorial we explore planetary nebulae as a stage in the evolution of low-to-intermediate-mass stars, as major contributors to the mass and chemical enrichment of the interstellar medium, and as astrophysical laboratories. We discuss many observed properties of planetary nebulae, placing particular emphasis on element abundance determinations and comparisons with theoretical predictions. Dust and molecules associated with planetary nebulae are considered as well. We then examine distances, binarity, and planetary nebula morphology and evolution. We end with mention of some of the advances that will be enabled by future observing capabilities.

Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk.

Neutron star models with maximum mass close to $2 \ M_{\odot}$ reach high central densities, which may activate nucleonic and hyperon direct Urca neutrino emission. To alleviate the tension between fast theoretical cooling rates and thermal luminosity observations of moderately magnetized, isolated thermally-emitting stars (with $L_{\gamma} \gtrsim 10^{31}$ erg s$^{-1}$ at $t \gtrsim 10^{5.3}$ yr), some internal heating source is required. The power supplied by the internal heater is estimated for both a phenomenological source in the inner crust and Joule heating due to magnetic field decay, assuming different superfluidity models and compositions of the outer stellar envelope. It is found that a thermal power of $W(t) \approx 10^{34}$ erg s$^{-1}$ allows neutron star models to match observations of moderately magnetized, isolated stars with ages $t \gtrsim 10^{5.3}$ yr. The requisite $W(t)$ can be supplied by Joule heating due to crust-confined initial magnetic configurations with (i) mixed poloidal-toroidal fields, with surface strength $B_{\textrm{dip}} = 10^{13}$ G at the pole of the dipolar poloidal component and $\sim 90$ per cent of the magnetic energy stored in the toroidal component; and (ii) poloidal-only configurations with $B_{\textrm{dip}} = 10^{14}$ G.

Nova Her 2021 (V1674 Her), which erupted on 2021 June 12, reached naked-eye brightness and has been detected from radio to $\gamma$-rays. An extremely fast optical decline of 2 magnitudes in 1.2 days and strong Ne lines imply a high-mass white dwarf. The optical pre-outburst detection of a 501.42s oscillation suggests a magnetic white dwarf. This is the first time that an oscillation of this magnitude has been detected in a classical nova prior to outburst. We report X-ray outburst observations from {\it Swift} and {\it Chandra} which uniquely show: (1) a very strong modulation of super-soft X-rays at a different period from reported optical periods; (2) strong pulse profile variations and the possible presence of period variations of the order of 0.1-0.3s; and (3) rich grating spectra that vary with modulation phase and show P Cygni-type emission lines with two dominant blue-shifted absorption components at $\sim 3000$ and 9000 km s$^{-1}$ indicating expansion velocities up to 11000 km s$^{-1}$. X-ray oscillations most likely arise from inhomogeneous photospheric emission related to the magnetic field. Period differences between reported pre- and post-outburst optical observations, if not due to other period drift mechanisms, suggest a large ejected mass for such a fast nova, in the range $2\times 10^{-5}$-$2\times 10^{-4} M_\odot$. A difference between the period found in the {\it Chandra} data and a reported contemporaneous post-outburst optical period, as well as the presence of period drifts, could be due to weakly non-rigid photospheric rotation.

We present the first analysis of internal coronal mass ejection (CME) structure observed very close to the Sun by the Wide-field Imager for Solar PRobe (WISPR) instrument on board Parker Solar Probe (PSP). The transient studied here is a CME observed during PSP's second perihelion passage on 2019 April 2, when PSP was only 40 R_sun from the Sun. The CME was also well observed from 1 au by the STEREO-A spacecraft, which tracks the event all the way from the Sun to 1 au. However, PSP/WISPR observes internal structure not apparent in the images from 1 au. In particular, two linear features are observed, one bright and one dark. We model these features as two loops within the CME flux rope channel. The loops can be interpreted as bundles of field lines, with the brightness of the bright loop indicative of lots of mass being loaded into those field lines, and with the dark loop being devoid of such mass loading. It is possible that these loops are actually representative of two independent flux rope structures within the overall CME outline.

Whilst the underlying assumption of the Friedman-Lema\^itre-Robertson-Walker (FLRW) cosmological model is that matter is homogeneously distributed throughout the universe, gravitational influences over the life of the universe have resulted in mass clustered on a range of scales. Hence we expect that, in our inhomogeneous universe, the view of an observer will be influenced by the location and local environment. Here we analyse the one-point probability distribution functions and angular power spectra of weak-lensing (WL) convergence and magnification numerically to investigate the influence of our local environment on WL statistics in relativistic $N$-body simulations. To achieve this, we numerically solve the null geodesic equations which describe the propagation of light bundles backwards in time from today, and develop a ray-tracing algorithm, and from these calculate various WL properties. Our findings demonstrate how cosmological observations of large-scale structure through WL can be impacted by the locality of the observer. We also calculate the constraints on the cosmological parameters as a function of redshift from the theoretical and numerical study of the angular power spectrum of WL convergence. This study concludes the minimal redshift for the constraint on the parameter $\Omega_m$ ($H_0$) is $z \sim 0.2$ $(z \sim 0.6 )$ beyond which the local environment's effect is negligible and the data from WL surveys are more meaningful above that redshift. The outcomes of this study will have direct consequences for future surveys, where percent-level-precision is necessary.

We present a spectral analysis of four LMC WC-type Wolf-Rayet (WR) stars (BAT99-8, BAT99-9, BAT99-11, and BAT99-52) to shed light on two evolutionary questions surrounding massive stars. The first is: are WO-type WR stars more oxygen enriched than the WC-type stars, indicating further chemical evolution, or are the strong high-excitation oxygen lines in the WO-type stars an indication of higher temperatures. This study will act as a baseline for answering the question of where WO-type stars fall in WR evolution. Each star's spectrum, extending from 1100~\AA\ to 25000~\AA, was modeled using \cmfgen\ to determine the star's physical properties such as luminosity, mass-loss rate, and chemical abundances. The oxygen abundance is a key evolutionary diagnostic, and with higher resolution data and an improved stellar atmosphere code, we found the oxygen abundance to be up to a factor of 5 lower than previous studies. The second evolutionary question revolves around the formation of WR stars: Do they evolve by themselves or is a close companion star necessary for their formation? Using our derived physical parameters, we compared our results to the Geneva single-star evolutionary models and the BPASS binary evolutionary models. We found that both the Geneva solar metallicity models and BPASS LMC metallicity models are in agreement with the four WC-type stars, while the Geneva LMC metallicity models are not. Therefore, these four WC4 stars could have been formed either via binary or single-star evolution.

WD 1856+534 b is a Jupiter-sized, cool giant planet candidate transiting the white dwarf WD 1856+534. Here, we report an optical transmission spectrum of WD 1856+534 b obtained from ten transits using the Gemini Multi-Object Spectrograph. This system is challenging to observe due to the faintness of the host star and the short transit duration. Nevertheless, our phase-folded white light curve reached a precision of 0.12 %. WD 1856+534 b provides a unique transit configuration compared to other known exoplanets: the planet is $8\times$ larger than its star and occults over half of the stellar disc during mid-transit. Consequently, many standard modeling assumptions do not hold. We introduce the concept of a `limb darkening corrected, time-averaged transmission spectrum' and propose that this is more suitable than $(R_{\mathrm{p}, \lambda} / R_{\mathrm{s}})^2$ for comparisons to atmospheric models for planets with grazing transits. We also present a modified radiative transfer prescription. Though the transmission spectrum shows no prominent absorption features, it is sufficiently precise to constrain the mass of WD 1856+534 b to be > 0.84 M$_\mathrm{J}$ (to $2 \, \sigma$ confidence), assuming a clear atmosphere and a Jovian composition. High-altitude cloud decks can allow lower masses. WD 1856+534 b could have formed either as a result of common envelope evolution or migration under the Kozai-Lidov mechanism. Further studies of WD 1856+534 b, alongside new dedicated searches for substellar objects around white dwarfs, will shed further light on the mysteries of post-main sequence planetary systems.

This paper calibrates how metrics derivable from the SAO/NASA Astrophysics Data System can be used to estimate the future impact of astronomy research careers and thereby to inform decisions on resource allocation such as job hires and tenure decisions. Three metrics are used, citations of refereed papers, citations of all publications normalized by the numbers of co-authors, and citations of all first-author papers. Each is individually calibrated as an impact predictor in the book Kormendy (2020), "Metrics of Research Impact in Astronomy" (Publ Astron Soc Pac, San Francisco). How this is done is reviewed in the first half of this paper. Then, I show that averaging results from three metrics produces more accurate predictions. Average prediction machines are constructed for different cohorts of 1990-2007 PhDs and used to postdict 2017 impact from metrics measured 10, 12, and 15 years after the PhD. The time span over which prediction is made ranges from 0 years for 2007 PhDs to 17 years for 1990 PhDs using metrics measured 10 years after the PhD. Calibration is based on perceived 2017 impact as voted by 22 experienced astronomers for 510 faculty members at 17 highly-ranked university astronomy departments world-wide. Prediction machinery reproduces voted impact estimates with an RMS uncertainty of 1/8 of the dynamic range for people in the study sample. The aim of this work is to lend some of the rigor that is normally used in scientific research to the difficult and subjective job of judging people's careers.

Observations of the Sun's off-limb white-light (WL) flares offer rare opportunities to study the energy release and transport mechanisms in flare loops. One of the best such events was SOL2017-09-10, an X8.2 flare that occurred near the Sun's west limb on 2017 September 10 and produced a WL loop system lasting more than 60 minutes and reaching an altitude higher than 30 Mm. The event was well observed by a suite of ground- and space-based instruments, including the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) that captured its off-limb loops in WL continuum near Fe I 6173 A, and the Atmospheric Imager Assembly (SDO/AIA) that observed its ultraviolet (UV) and extreme-ultraviolet (EUV) counterparts. We found quasi-periodic pulsations in the WL and UV emissions at the flare loop-top with a period around 8.0 min. Each pulsation appears to have an EUV counterpart that occurs earlier in time and higher in altitude. Despite many similarities in the WL and UV images and light curves, the WL flux at the loop-top continues to grow for about 16 minutes while the UV fluxes gradually decay. We discuss the implication of these unprecedented observations on the understanding of the enigmatic off-limb WL flare emission mechanisms.

Some studies suggest that the dust temperatures ($T_\mathrm{d}$) in high-redshift ($z\gtrsim 5$) Lyman break galaxies (LBGs) are high. However, possible observational bias in $T_\mathrm{d}$ is yet to be understood. Thus, we perform a simple test using random realizations of LBGs with various stellar masses, dust temperatures, and dust-to-stellar mass ratios, and examine how the sample detected by ALMA is biased in terms of $T_\mathrm{d}$. We show that ALMA tends to miss high-$T_\mathrm{d}$ objects even at total dust luminosity. LBGs are, however, basically selected by the stellar UV luminosity. The dust-temperature bias in a UV-selected sample is complicated because of the competing effects between high $T_\mathrm{d}$ and low dust abundance. For ALMA Band 6, there is no tendency of high-$T_\mathrm{d}$ LBGs being more easily detected in our experiment. Thus, we suggest that the observed trend of high $T_\mathrm{d}$ in $z\gtrsim 5$ LBGs is real. We also propose that the 450 $\mu$m band is useful in further clarifying the dust temperatures. To overcome the current shallowness of 450 $\mu$m observations, we examine a future Antarctic 30-m class telescope with a suitable atmospheric condition for wavelengths $\lesssim 450~\mu$m, where the detection is not confusion-limited. We find that, with this telescope, an $L_\mathrm{IR}$-selected sample with $\log(L_\mathrm{IR}/\mathrm{L}_{\odot})>11$ is constructed for $z\gtrsim 5$, and detection in the intermediate-$M_\star$ (stellar mass) range [$9<\log (M_\star /M_{\odot})<9.5$] is much improved, especially at high $T_\mathrm{d}$.

Cosmic ray muons have emerged as a non-conventional high-energy radiation probe to monitor dense and large objects. Muons are the most abundant cosmic radiation on Earth, however, their flux at sea level is approximately 10,000/min-m2 much less than that of induced radiation, i.e., x-rays or electron beams. Cosmic ray muon flux varies with the particle incident angle and it is frequently approximated using a cosine-squared which can introduce large errors for high zenith angles. However, the cosmic ray muon flux depends on not only the zenith angle but also on the effective solid angle and the geometric characteristics of the detectors. Since the low muon flux typically results in long measurement times, an accurate estimation of the measurable muon counts is important for many muon applications. Here we propose a simple and versatile semi-empirical model to improve the accuracy in muon flux estimation at all zenith angles by incorporating the geometric parameters of detectors. We call this the Effective solid angle model. To demonstrate the functionality of our model, we compare with i) cosmic ray muon measurements, ii) the cosine-squared model, and iii) Monte-Carlo simulations. Our results show that the muon count rate estimation capability is significantly improved resulting in a reduced mean relative error from 30 % (for the cosine-squared model) to less than 15 % for the effective solid angle model. In addition, this model is simple enough and works universally for all detector geometries and configurations. By selecting an appropriate intensity correlation, the model can be easily extended to estimate muon flux at any altitude and underground level. Finally, a simple empirical correlation is derived in order to compute the expected cosmic ray muon counts in a single step.

We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers.

Understanding solar coronal magnetic fields is crucial to address the long-standing mysteries of the solar corona and solar wind. Although routine photospheric magnetic fields (MFs) are available for decades, coronal MFs are rarely reported. Visible Emission Line Coronagraph (VELC) on board Aditya-L1 mission (planned to launch in the near future) can directly measure the MFs in the inner solar corona. This can be achieved with the help of spectropolarimetric observations of the forbidden coronal emission line centered at 1074.7 nm over a field of view 1.05 R$_{\odot}$ - 1.5 R$_{\odot}$. In this article, we summarize various direct and indirect techniques used to estimate the MFs at different wavelength regimes. Further, we summarize the expected accuracies that are required to estimate MFs using VELC's spectropolarimetry channel.

We have conducted mapping observations ($\sim 2'\times2'$) of the Class I protostar L1489 IRS using the 7-m array of the Atacama Compact Array (ACA) and the IRAM-30m telescope in the $\mathrm{C^{18}O}$ 2-1 emission to investigate the gas kinematics on 1000-10,000 au scales. The $\mathrm{C^{18}O}$ emission shows a velocity gradient across the protostar in a direction almost perpendicular to the outflow. The radial profile of the peak velocity was measured from a $\mathrm{C^{18}O}$ position-velocity diagram cut along the disk major axis. The measured peak velocity decreases with radius at a radii of $\sim$1400-2900 au, but increases slightly or is almost constant at radii of $r\gtrsim$2900 au. Disk-and-envelope models were compared with the observations to understand the nature of the radial profile of the peak velocity. The measured peak velocities are best explained by a model where the specific angular momentum is constant within a radius of 2900 au but increases with radius outside 2900 au. We calculated the radial profile of the specific angular momentum from the measured peak velocities, and compared it to analytic models of core collapse. The analytic models reproduce well the observed radial profile of the specific angular momentum and suggest that material within a radius of $\sim$4000-6000 au in the initial dense core has accreted to the central protostar. Because dense cores are typically $\sim$10,000-20,000 au in radius, and as L1489 IRS is close to the end of mass accretion phase, our result suggests that only a fraction of a dense core eventually forms a star.

Comparisons of the alignment of exoplanets with a common host star can be used to distinguish among concurrent evolution scenarios. However, multi-planet systems usually host mini-Neptunes and super-Earths, whose size make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions technique, which can access spin-orbit angles of small planets by exploiting the full information contained in spectral transit time series. We validated the technique on published HARPS-N data of the mini-Neptune HD3167c, refining its high sky-projected spin-orbit angle (-108.9+5.4-5.5 deg), and we applied it to new ESPRESSO observations of the super-Earth HD3167b, revealing an aligned orbit (-6.6+6.6-7.9 deg). Surprisingly different variations in the contrast of the stellar lines occulted by the planets can be reconciled with a latitudinal dependence of the stellar line shape. In this scenario, a joint fit to both datasets constrains the inclination of the star (111.6+3.1-3.3 deg) and the 3D spin-orbit angles of HD3167b (29.5+7.2-9.4 deg) and HD3167c (107.7+5.1-4.9 deg). The projected spin-orbit angles do not depend on the model for the line contrast variations, and so, with a mutual inclination of 102.3+7.4-8.0 deg, we conclude that the two planets are on perpendicular orbits. This could be explained by HD3167b being strongly coupled to the star and retaining its primordial alignment, whereas HD3167c would have been brought to a nearly polar orbit via secular gravitational interactions with an outer companion. Follow-up observations and dynamical evolution simulations are required to search for this companion and explore this scenario. HD3167b is the smallest exoplanet with a confirmed spectroscopic Rossiter-McLaughlin signal. Our new technique opens the way to determining the orbital architectures of the super-Earth and Earth-sized planet populations.

The ionization and thermal state of the intergalactic medium (IGM) during the epoch of reionization has been of interest in recent times because of their close connection to the first stars. We present in this paper a semi-numerical code which computes the large-scale temperature and ionized hydrogen fields in a cosmologically representative volume accounting for the patchiness in these quantities arising from reionization. The code is an extension to a previously developed version for studying the growth of ionized regions, namely, Semi Numerical Code for ReionIzation with PhoTon Conservation (SCRIPT). The main addition in the present version are the inhomogeneous recombinations which are essential for temperature calculations. This extended version of SCRIPT also implements physical consequences of photoheating during reionization, e.g., radiative feedback. These enhancements allow us to predict observables which were not viable with the earlier version. These include the faint-end of the ultra-violet luminosity function of galaxies (which can get affected by the radiative feedback) and the temperature-density relation of the low-density IGM at $z \sim 6$. We study the effect of varying the free parameters and prescriptions of our model on a variety of observables. The conclusion of our analysis is that it should be possible to put constraints on the evolution of thermal and ionization state of the IGM using available observations accounting for all possible variations in the free parameters. A detailed exploration of the parameter space will be taken up in the future.

Radio relics in the outskirts of galaxy clusters imply the diffusive shock acceleration (DSA) of electrons at merger-driven shocks with Mach number $M_{s}\lesssim3-4$ in the intracluster medium (ICM). Recent studies have suggested that electron preacceleration and injection, prerequisite steps for DSA, could occur at supercritical shocks with $M_{s}\gtrsim2.3$ in the ICM, thanks to the generation of multiscale waves by microinstabilities such as the Alfv\'en ion cyclotron (AIC) instability, the electron firehose instability (EFI), and the whistler instability (WI). On the other hand, some relics are observed to have subcritical shocks with $M_{s}\lesssim2.3$, leaving DSA at such weak shocks as an outstanding problem. Reacceleration of preexisting nonthermal electrons has been contemplated as one of possible solutions for that puzzle. To explore this idea, we perform Particle-in-Cell (PIC) simulations for weak quasi-perpendicular shocks in high-$\beta$ ($\beta=P_{\rm gas}/P_{B}$) plasmas with power-law suprathermal electrons in addition to Maxwellian thermal electrons. We find that suprathermal electrons enhance the excitation of electron-scale waves via the EFI and WI. However, they do not affect the ion reflection and the ensuing generation of ion-scale waves via the AIC instability. The presence of ion-scale waves is the key for the preacceleration of electrons up to the injection momentum, thus the shock criticality condition for electron injection to DSA is preserved. Based on the results, we conclude that preexisting nonthermal electrons in the preshock region alone would not resolve the issue of electron preacceleration at subcritical ICM shocks.

Aims. We investigate the synthetic observational signatures from numerical models of coronal heating in an arcade to determine what features are associated with such heating, and what tools can be used to identify them. Methods. We consider two simulations of coronal arcades driven by footpoint motions with different characteristic timescales. Forward modelling is then used, and the synthetic emission data is analysed. Results. The total intensity and Doppler velocities clearly show the magnetic structure of the coronal arcade. Contrasts in the local Doppler shift also highlight the locations of separatrix surfaces. The distinguishing feature of the AC and DC models is that of the frequencies. Through FFT analysis of the Doppler velocities, when short timescale footpoint motions are present, higher frequencies are observed. For longer timescale motions, the dominant signal is that of lower frequencies; however, higher frequencies were also detected, which matched the natural Alfv\'en frequency of the background magnetic field. Alfv\'enic wave signatures were identified in both models with fast wave signatures observable in the AC model. Finally, the estimates of the kinetic energy using the Doppler shifts were found to significantly underestimate within these models. Conclusions. Observables identified within this article were from features such as Alfv\'en waves, fast waves, the arcade structure and separatrix surfaces. The two models were differentiated by examining the frequencies present. The Doppler velocities cannot provide accurate estimates of the total kinetic energy or the component parallel to the LOS. This is due to some plasma outside the formation temperature range of the ion, the multi-directional driver, and cancellation of the velocity along the LOS. The impact each factor has on the estimation is dependent on the set up of the model and the chosen emission line.

We present a study based on the high-resolution spectroscopy and K2 space photometry of five chemically peculiar stars in the region of the open cluster M44. The analysis of the high-precision photometric K2 data reveals that the light variations in HD 73045 and HD 76310 are rotational in nature and caused by spots or cloud-like co-rotating structures, which are non-stationary and short-lived. The time-resolved radial velocity measurements, in combination with the K2 photometry, confirm that HD 73045 does not show any periodic variability on timescales shorter than 1.3 d, contrary to previous reports in the literature. In addition to these new rotational variables, we discovered a new heartbeat system, HD 73619, where no pulsational signatures are seen. The spectroscopic and spectropolarimetric analyses indicate that HD 73619 belongs to the peculiar Am class, with either a weak or no magnetic field considering the 200 G detection limit of our study. The Least-Squares Deconvolution (LSD) profiles for HD 76310 indicate a complex structure in its spectra suggesting that this star is either part of a binary system or surrounded by a cloud shell. When placed in the Hertzsprung-Russell diagram, all studied stars are evolved from main-sequence and situated in the $\delta$ Scuti instability strip. The present work is relevant for further detailed studies of CP stars, such as inhomogeneities (including spots) in the absence of magnetic fields and the origin of the pulsational variability in heartbeat systems.

The orbital motion of the transiting hot Jupiter TrES-5 b was reported to be perturbed by a planetary companion on a nearby orbit. Such compact systems do not frequently occur in nature, and learning their orbital architecture could shed some light on hot Jupiters' formation processes. We acquired fifteen new precise photometric time series for twelve transits of TrES-5 b between June 2019 and October 2020 using 0.9-2.0 m telescopes. The method of precise transit timing was employed to verify the deviation of the planet from the Keplerian motion. Although our results show no detectable short-time variation in the orbital period of TrES-5 b and the existence of the additional nearby planet is not confirmed, the new transits were observed about two minutes earlier than expected. We conclude that the orbital period of the planet could vary in a long timescale. We found that the most likely explanation of the observations is the line-of-sight acceleration of the system's barycentre due to the orbital motion induced by a massive, wide-orbiting companion.

The new generation of CMB polarization experiments will reach limits of sensitivity never achieved before to detect the primordial B-mode signal. However, all these efforts will be futile if we lack a tight control of systematics. Here, we focus on the systematic that arises from the uncertainty on the calibration of polarization angles. Miscalibrated polarization angles induce a mixing of E- and B-modes that obscures the primordial B-mode signal. We introduce an iterative angular power spectra maximum likelihood-based method to calculate polarization angles from the multi-frequency signal. The basis behind this methodology grounds on nulling the EB power spectra. To simplify the likelihood, we assume that the rotation angles are small (<6 deg) and, the maximum likelihood solution for the rotation angles is obtained by applying an iterative process where the covariance matrix does not depend on the angles per iteration, i.e., the rotation angles are fixed to the estimated angles in the previous iteration. With these assumptions, we obtain an analytical linear system which leads to a very fast computational implementation. We show that with this methodology we are able to determine the rotation angle for each frequency with sufficiently good accuracy. To prove the latter point we perform component separation analyses using the parametric component separation method B-SeCRET with two different approaches. In the first approach we apply the B-SeCRET pipeline to the signal de-rotated with the estimation of the angles, while in the second, the rotation angles are treated as model parameters using the estimation of the angles as a prior information. We obtain that the rotation angles estimations improve after applying the second approach, and show that the systematic residuals due to the non-null calibration polarization angles are mitigated to the order of a 1% at the power spectrum level.

There are 23 long-period binary systems discovered to date that contain a B-type hot subdwarf(sdB) whose orbital parameters have been fully solved. They evolve into O-type subdwarfs (sdO) once the helium burning transitions from the core to the He shell. Their study will help constraint parameters on the formation and evolution of these binaries and explain some of their puzzling features. In this study, we aim to solve orbital and atmospheric parameters of two long-period sdO binaries and, for the first time, investigate the chemical composition of their main-sequence (MS) companions. HERMES high-resolution spectra are used to obtain radial velocities and solve their orbits. The Grid Search in Stellar Parameter code (GSSP) is used to derive the atmospheric parameters and photospheric chemical abundances of the MS companions. Stellar evolution models (MIST) are fitted to the companion atmospheric parameters to derive masses. In the bimodal period-eccentricity diagram, the orbital parameters indicate that Feige 80 matches the same correlation as the majority of the systems. The analysis suggests that Feige 80 has a canonical subdwarf mass and followed a standard formation channel. However, BD-11$^{\rm{o}}$162 is an exceptional system with a lower mass. It also shows a carbon overabundance, which could be caused by a higher progenitor mass. The yttrium depletion in both MS companions could indicate the existence of a circumbinary disk in these systems' pasts. Nevertheless, a chemical analysis of a larger sample is necessary to draw strong conclusions.

Meter-sized ground-based telescopes are frequently used today for the follow-up of extrasolar planet candidates. While the transit signal of a Jupiter-sized object can typically be detected to a high level of confidence with small telescope apertures as well, the shallow transit dips of planets with the size of Neptune and smaller are more challenging to reveal. We employ new observational data to illustrate the photometric follow-up capabilities of meter-sized telescopes for shallow exoplanet transits. We describe in detail the capability of distinguishing the photometric signal of an exoplanet transit from an underlying trend in the light curve. The transit depths of the six targets we observed, Kepler-94b, Kepler-63b, K2-100b, K2-138b, K2-138c, and K2-138e, range from 3.9 ppt down to 0.3 ppt. For five targets of this sample, we provide the first ground-based photometric follow-up. We detect or rule out the transit features significantly in single observations for the targets that show transits of 1.3 ppt or deeper. The shallower transit depths of two targets of 0.6 and 0.8 ppt were detected tentatively in single light curves, and were detected significantly by repeated observations. Only for the target of the shallowest transit depth of 0.3 ppt were we unable to draw a significant conclusion despite combining five individual light curves. An injection-recovery test on our real data shows that we detect transits of 1.3 ppt depth significantly in single light curves if the transit is fully covered, including out-of-transit data toward both sides, in some cases down to 0.7 ppt depth. For Kepler-94b, Kepler-63b, and K2-100b, we were able to verify the ephemeris. In the case of K2-138c with a 0.6 ppt deep transit, we were able to refine it, and in the case of K2-138e, we ruled out the transit in the time interval of more than +-1.5 sigma of its current literature ephemeris.

Gamma-ray bursts are the most powerful explosions in the universe and are mainly placed at very large redshifts, up to $z\simeq 9$. In this short review, we first discuss gamma-ray burst classification and morphological properties. We then report the likely relations between gamma-ray bursts and other astronomical objects, such as black holes, supernovae, neutron stars, etc., discussing in detail gamma-ray burst progenitors. We classify long and short gamma-ray bursts, working out their timescales, and introduce the standard fireball model. Afterwards, we focus on direct applications of gamma-ray bursts to cosmology and underline under which conditions such sources would act as perfect standard candles if correlations between photometric and spectroscopic properties were not jeopardized by the \emph{circularity problem}. In this respect, we underline how the shortage of low-$z$ gamma-ray bursts prevents anchor gamma-ray bursts with primary distance indicators. Moreover, we analyze in detail the most adopted gamma-ray burst correlations, highlighting their main differences. We therefore show calibration techniques, comparing such treatments with non-calibration scenarios. For completeness, we discuss the physical properties of the correlation scatters and systematics occurring during experimental computations. Finally, we develop the most recent statistical methods, star formation rate and high-redshift gamma-ray burst excess and show the most recent constraints got from experimental analyses.

Optical photometry and spectra of the luminous type IIn supernova SN~2008iy are analysed in detail with implications for cosmic ray acceleration and the radio emission. The light curve and expansion velocities indicate ejecta with the kinetic energy of $3\times10^{51}$ erg to collide with the $\sim10$ Msun circumstellar envelope. The luminous Ha is explained as originated primarily from circumstellar clouds interacting with the forward shock. For the first time the fluorescent OI 8446A emission is used to demonstrate that the cloud fragmentation cascade spans a scale range > 2.3 dex. The narrow circumstellar Ha permitted us to estimate the acceleration efficiency of cosmic rays. The found value is close to the efficiency inferred in the same way for other two SNe~IIn, SN~1997eg and SN~2002ic. The efficiency of cosmic ray acceleration is utilized to reproduce the radio flux from SN~2008iy for the amplified magnetic field consistent with the saturated turbulent magnetic field in the diffusive shock acceleration mechanism.

GRAND is designed to detect ultra-high-energy cosmic particles -- specially neutrinos, cosmic rays and gamma rays using radio antennas. With $\sim$20 mountainous sites around the world it will cover a total area of 200,000 km$^{2}$. The planned sensitivity of 10$^{-10}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above $5\times10^{17}$ eV will likely ensure the detection of cosmogenic neutrinos predicted by most common scenarios enabling neutrino astronomy. Furthermore, PeV--EeV neutrinos can test particle interactions at energies above those achieved in accelerators. The pathfinder stage GRANDProto300 is planned to start taking data in 2021. We present the current overall status of the project with emphasis on the neutrino physics.

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), that provides twelve photometric optical passbands over 2176 deg2. We analysed 5926 white dwarfs with r <= 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen (H) and helium (He) dominated atmospheres. We estimated the PDF for effective temperature (Teff), surface gravity, parallax, and spectral type; and the probability of having a H-dominated atmosphere (pH) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as reference, and derived a self-consistent prior for the atmospheric composition as a function of Teff. We described the fraction of He-dominated atmosphere white dwarfs (fHe) with a linear function of Teff at 5000 < Teff < 30000 K. We found fHe = 0.24 +- 0.01 at Teff = 10000 K, a change rate along the cooling sequence of 0.14 +- 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 +- 0.02 at the high-temperature end. We tested the obtained pH by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 Msun, and Teff > 6000 K. The result for H-dominated white dwarfs agrees with previous work, with a dominant M = 0.59 Msun peak and the presence of an excess at M ~ 0.8 Msun. This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provides a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 +- 3 % increase in the fraction of He-dominated white dwarfs from Teff = 20000 K to Teff = 5000 K.

Active galactic nucleus (AGN) emission is dominated by stochastic, aperiodic variability which overwhelms any periodic/quasi-periodic signal (QPO) if one is present. The Auto Correlation Function (ACF) and Phase Dispersion Minimization (PDM) techniques have been used previously to claim detections of QPOs in AGN light curves. In this paper we perform Monte Carlo simulations to empirically test QPO detection feasibility in the presence of red noise. Given the community's access to large databases of monitoring light curves via large-area monitoring programmes, our goal is to provide guidance to those searching for QPOs via data trawls. We simulate evenly-sampled pure red noise light curves to estimate false alarm probabilities; false positives in both tools tend to occur towards timescales longer than (very roughly) one-third of the light curve duration. We simulate QPOs mixed with pure red noise and determine the true-positive detection sensitivity; in both tools, it depends strongly on the relative strength of the QPO against the red noise and on the steepness of the red noise PSD slope. We find that extremely large values of peak QPO power relative to red noise (typically $\sim 10^{4-5}$) are needed for a 99.7 per cent true-positive detection rate. Given that the true-positive detections using the ACF or PDM are generally rare to obtain, we conclude that period searches based on the ACF or PDM must be treated with extreme caution when the data quality is not good. We consider the feasibility of QPO detection in the context of highly-inclined, periodically self-lensing supermassive black hole binaries.

We study the leading-order nonlinear gravitational waves (GWs) produced by an electromagnetic (EM) stress in reheating magnetogenesis scenarios. Both nonhelical and helical magnetic fields are considered. By numerically solving the linear and leading-order nonlinear GW equations, we find that the GW energy from the latter is usually larger. We compare their differences in terms of the GW spectrum and parameterize the GW energy difference due to the nonlinear term, $\Delta\mathcal{E}_{\rm GW}$, in terms of EM energy $\mathcal{E}_{\rm EM}$ as $\Delta\mathcal{E}_{\rm GW}=(\tilde p\mathcal{E}_{\rm EM}/k_*)^3$, where $k_*$ is the characteristic wave number, $\tilde p=0.84$ and $0.88$ are found in the nonhelical and helical cases, respectively, with reheating around the QCD energy scale, while $\tilde p=0.45$ is found at the electroweak energy scale. We also compare the polarization spectrum of the linear and nonlinear cases and find that adding the nonlinear term usually yields a decrease in the polarization that is proportional to the EM energy density. We parameterize the fractional polarization suppression as $|\Delta \mathcal{P}_{\rm GW}/\mathcal{P}_{\rm GW}|=\tilde r \mathcal{E}_{\rm EM}/k_*$ and find $\tilde r = 1.2 \times 10^{-1}$, $7.2 \times 10^{-4}$, and $3.2 \times 10^{-2}$ for the helical cases with reheating temperatures $T_{\rm r} = 300 {\rm TeV}$, $8 {\rm GeV}$, and $120 {\rm MeV}$, respectively. Prospects of observation by pulsar timing arrays, space-based interferometers, and other novel detection proposals are also discussed.

We present a comprehensive study of GRS 1915+105 in wide energy band ($0.5-60$ keV) using AstroSat observations during the period of $2016-2019$. The MAXI X-ray lightcurve of the source shows rise and decay profiles similar to canonical outbursting black holes. However, the source does not follow the exemplary 'q'-diagram in the Hardness-Intensity Diagram (HID). Model independent analysis of lightcurves suggests that GRS 1915+105 displays various types of variability classes ($\delta,\chi,\rho,\kappa,\omega$ and $\gamma$). We also report possible transitions from one class to another ($\chi\rightarrow\rho,\rho\rightarrow\kappa$ via an 'unknown' class and $\omega\rightarrow\gamma\rightarrow\omega+\gamma$) within a few hours duration. Broadband energy spectra are well modeled with multi-coloured disc blackbody and Comptonised components. We explore the 'spectro-temporal' features of the source in the different variability classes, transitions between classes, and evolution during $2016-2019$. Detailed analysis indicates a gradual increase in the photon index ($\Gamma$) from $1.83$ to $3.8$, disc temperature ($kT_{in}$) from $1.33$ to $2.67$ keV, and Quasi-periodic Oscillation (QPO) frequency ($\nu$) from $4$ to $5.64$ Hz during the rise, while the parameters decrease to $\Gamma$ ~$1.18$, $kT_{in}$ ~$1.18$ keV, and $\nu$ ~$1.38$ Hz respectively in the decline phase. The source shows maximum bolometric luminosity (L$_{bol}$) during the peak at ~$36$% of Eddington luminosity (L$_{EDD}$), and a minimum of ~$2.4$% L$_{EDD}$ during the decay phase. Further evolution of the source towards an obscured low-luminosity (L$_{bol}$ of ~ 1% L$_{EDD}$) phase, with a decrease in the intrinsic bolometric luminosity of the source due to obscuration, has also been indicated from our analysis. The implication of our results are discussed in the context of accretion disc dynamics around the black hole.

We present a novel approach to the riddle of star cluster multiple populations. Stars form from molecular cores. But not all cores form stars. Following their initial compression, such 'failed' cores re-expand, rather than collapsing. We propose that their formation and subsequent dispersal regulate the gas density of cluster-forming clumps and, therefore, their core and star formation rates. Clumps for which failed cores are the dominant core type experience star formation histories with peaks and troughs. In contrast, too few failed cores results in smoothly decreasing star formation rates. We identify three main parameters shaping the star formation history of a clump: the star and core formation efficiencies per free-fall time, and the time-scale on which failed cores return to the clump gas. The clump mass acts as a scaling factor. We use our model to constrain the density and mass of the Orion Nebula Cluster progenitor clump, and to caution that the star formation histories of starburst clusters may contain close-by peaks concealed by stellar age uncertainties. Our model generates a great variety of star formation histories. Intriguingly, the chromosome maps and O-Na anti-correlations of old globular clusters also present diverse morphologies. This prompts us to discuss our model in the context of globular cluster multiple stellar populations. More massive globular clusters exhibit stronger multiple stellar population patterns, which our model can explain if the formation of the polluting stars requires a given stellar mass threshold.

The Ly$\alpha$ luminosity function (LF) of Ly$\alpha$ emitters (LAEs) has been used to constrain the neutral hydrogen fraction in the intergalactic medium (IGM) and thus the timeline of cosmic reionization. Here we present the results of a new narrow-band imaging survey for $z=7.3$ LAEs in a large area of $\sim 3\ \mathrm{deg}^2$ with Subaru/Hyper Suprime-Cam. No LAEs are detected down to $L_{\mathrm{Ly}\alpha}\simeq 10^{43.2}\ \mathrm{erg\ s^{-1}}$ in an effective cosmic volume of $\sim 2\times 10^6$ Mpc$^3$, placing an upper limit to the bright part of the $z=7.3$ Ly$\alpha$ LF for the first time and confirming a decrease in bright LAEs from $z=7.0$. By comparing this upper limit with the Ly$\alpha$ LF in the case of the fully ionized IGM, which is predicted using an observed $z=5.7$ Ly$\alpha$ LF on the assumption that the intrinsic Ly$\alpha$ LF evolves in the same way as the UV LF, we obtain the relative IGM transmission $T^\mathrm{IGM}_{\mathrm{Ly}\alpha}(7.3)/T^\mathrm{IGM}_{\mathrm{Ly}\alpha}(5.7)<0.77$, and then the volume-averaged neutral fraction $x_\mathrm{HI}(7.3)>0.28$. Cosmic reionization is thus still ongoing at $z=7.3$, being consistent with results from other $x_\mathrm{HI}$ estimation methods. A similar analysis using literature Ly$\alpha$ LFs finds that at $z=6.6$ and 7.0 the observed Ly$\alpha$ LF agrees with the predicted one, consistent with full ionization.

We propose a closed-form normalization method suitable for the study of the secular dynamics of small bodies in heliocentric orbits perturbed by the tidal potential of a planet with orbit external to the orbit of the small body. The method makes no use of relegation, thus, circumventing all convergence issues related to that technique. The method is based on a convenient use of a book-keeping parameter keeping simultaneously track of all the small quantities in the problem. The book-keeping affects both the Lie series and the Poisson structure employed in successive perturbative steps. In particular, it affects the definition of the normal form remainder at every normalization step. We show the results obtained by assuming Jupiter as perturbing planet and we discuss the validity and limits of the method.

The Cherenkov Telescope Array (CTA) will use three telescope sizes to efficiently detect cosmic gamma rays in the energy range from several tens of GeV to hundreds of TeV. The Small-Sized Telescopes (SSTs) will form the largest section of the array, covering an area of many square kilometres on the CTA southern site in Paranal, Chile. Up to 70 SSTs will be implemented by an international consortium of institutes and teams as an in-kind contribution to the CTA Observatory. The SSTs will provide unprecedented sensitivity to gamma rays above 1 TeV and the highest angular resolution of any instrument above the hard X-ray band. CTA has recently finalised the technology that will be used for the SSTs: the telescopes will be a dual-reflector design with a primary reflector of ~4 m diameter, equipped with an SiPM-based camera with full waveform readout from $\sim$2000 channels covering a $\sim$9$^\circ$ field of view. The Schwarzschild-Couder optical configuration leads to a small plate-scale, and consequently a compact, cost-efficient camera. In this contribution, we describe the experience gained operating telescope and camera prototypes during the CTA preparatory phase, and the development of the final SST design.

The elliptical power-law (EPL) model of the mass in a galaxy is widely used in strong gravitational lensing analyses. However, the distribution of mass in real galaxies is more complex. We quantify the biases due to this model mismatch by simulating and then analysing mock {\it Hubble Space Telescope} imaging of lenses with mass distributions inferred from SDSS-MaNGA stellar dynamics data. We find accurate recovery of source galaxy morphology, except for a slight tendency to infer sources to be more compact than their true size. The Einstein radius of the lens is also robustly recovered with 0.1% accuracy, as is the global density slope, with 2.5% relative systematic error, compared to the 3.4% intrinsic dispersion. However, asymmetry in real lenses also leads to a spurious fitted `external shear' with typical strength, $\gamma_{\rm ext}=0.015$. Furthermore, time delays inferred from lens modelling without measurements of stellar dynamics are typically underestimated by $\sim$5%. Using such measurements from a sub-sample of 37 lenses would bias measurements of the Hubble constant $H_0$ by $\sim$9%. Although this work is based on a particular set of MaNGA galaxies, and the specific value of the detected biases may change for another set of strong lenses, our results strongly suggest the next generation cosmography needs to use more complex lens mass models.

Context. ULXs ($L_\textrm{X}>10^{39}$ erg s$^{-1}$) are excellent probes of accretion physics, star formation and IMBHs searches. As the sample size of X-ray data from modern observatories increases, producing extensive catalogues of ULXs and studying their collective properties is a possibility and a priority. Aims. We build a ULX catalogue based on one of the latest XMM-Newton releases, 4XMM-DR9, and the galaxy catalogue HECATE and study whether the properties of the expanded XMM-Newton ULX population are consistent with previous findings. Methods. We perform cross-matching between XMM-Newton sources and HECATE objects to identify host galaxies, and flag known interlopers in external catalogues and databases. Manual inspection of image data from PanSTARRS1 and the NASA/IPAC database is occasionally performed. We use distance and luminosity arguments to identify ULX candidates. Spectral, abundance and variability properties of the candidates are studied from 4XMM-DR9, HECATE and 4XMM-DR9s parameters. Results. We identify 779 ULX candidates. In spiral galaxies the number of ULX candidates per star forming rate is consistent with previous studies, and a significant ULX population in elliptical and lenticular is also found. Candidates hosted by late-type galaxies are more abundant, present harder spectra and undergo more and more extreme inter-observation variability than the early-type hosted ones. Around 30 candidates with $L_\textrm{X}>10^{41}$ erg s$^{-1}$ are also identified. Conclusions. Our results on the spectral and abundance properties of ULXs confirm the findings made by previous studies based on XMM-Newton and Chandra data, while our population-scale study on variability properties is unprecedented. Our study provides limited insight on the properties of the brightest ULX candidates due to the small sample size.

We present a Markov chain Monte Carlo pipeline which can be used for unbiased large-scale tests of gravity using galaxy cluster observations. The pipeline, which currently uses cluster number counts to constrain the present-day background scalar field $f_{R0}$ of Hu-Sawicki $f(R)$ gravity, fully accounts for the effects of the fifth force on cluster properties including the dynamical mass, the halo concentration and the observable-mass scaling relations. This is achieved using general models which have been calibrated over a wide and continuous mass range ($10^{11}M_{\odot}\lesssim M\lesssim10^{15}M_{\odot}$) using a large suite of cosmological simulations, including the first to simultaneously incorporate both screened modified gravity and full baryonic physics. We show, using mock cluster catalogues, that an incomplete treatment of the observable-mass scaling relations in $f(R)$ gravity, which do not necessarily follow the usual power-law behaviour, can lead to unbiased and imprecise constraints. It is therefore essential to fully account for these effects in future cosmological tests of gravity that will make use of vast cluster catalogues from ongoing and upcoming galaxy surveys. Our constraint framework can be easily extended to other gravity models; to demonstrate this, we have carried out a similar modelling of cluster properties in the normal-branch Dvali-Gabadadze-Porrati model (nDGP), which features a very different screening mechanism. Using our full-physics simulations, we also study the angular power spectra of the thermal and kinetic Sunyaev-Zel'dovich effects in $f(R)$ gravity and nDGP, and demonstrate the potential for precise constraints of gravity using data from upcoming CMB experiments. Finally, we present a retuned baryonic physics model, based on the IllustrisTNG model, which can be used for full-physics simulations within large cosmological volumes.

Intrinsic iron abundance spreads in globular clusters, although usually small, are very common, and are signatures of self enrichment: some stars within the cluster have been enriched by supernova ejecta from other stars within the same cluster. We use the Bailin (2018) self enrichment model to predict the relationship between properties of the protocluster -- its mass and the metallicity of the protocluster gas cloud -- and the final observable properties today -- its current metallicity and the internal iron abundance spread. We apply this model to an updated catalog of Milky Way globular clusters where the initial mass and/or the iron abundance spread is known to reconstruct their initial metallicities. We find that with the exception of the known anomalous bulge cluster Terzan 5 and three clusters strongly suspected to be nuclear star clusters from stripped dwarf galaxies, the model provides a good lens for understanding their iron spreads and initial metallicities. We then use these initial metallicities to construct age-metallicity relations for kinematically-identified major accretion events in the Milky Way's history. We find that using the initial metallicity instead of the current metallicity does not alter the overall picture of the Milky Way's history, since the difference is usually small, but does provide information that can help distinguish which accretion event some individual globular clusters with ambiguous kinematics should be associated with, and points to potential complexity within the accretion events themselves.

The Reheating process at the end of inflation is often modeled by an oscillating scalar field which shows a background dust-like behaviour, prompting the analysis of gravitational collapse and black hole formation in this era to be approached by the spherical collapse of standard structure formation. In the scalar field dark matter structure formation process virialized halos halt the direct collapse, resulting in halos with condensed central cores at the de Broglie scale of the dominant scalar field. We show that a similar process can take place during reheating, leading to the formation of primordial black holes (PBHs). We study the formation of PBHs through the gravitational further collapse of structures virialized during reheating, looking at the collapse of either the whole structure, or that of the central core within these configurations. We compute the threshold amplitude for the density contrast to undergo this process, for both free and self-interacting scalar fields. We discuss the relevance of our results for the abundance of PBHs at the lower end of the mass spectrum.

We analyze the uncertainty in grain size estimation of pure methane (CH4) and nitrogen saturated with methane (N2:CH4) ices, the most abundant volatile materials on trans-Neptunian objects (TNOs) and Kuiper belt objects (KBOs). We compare the single scattering albedo, which determines the grain size estimation of outer solar system regolith (Hansen, 2009), of these ices using the Mie scattering model and two other Hapke approximations (Hapke, 1993) in radiative transfer scattering models (RTM) at near-infrared (NIR) wavelengths (1-5 $\mu$m). The equivalent slab (Hapke Slab) approximation model predicts results much closer to Mie scattering over the NIR wavelengths at a wide range of grain sizes. In contrast, even though the internal scattering model (ISM) predicts an approximate particle diameter close to the Mie model for particles with a 10 $\mu$m radii, it exhibits higher discrepancies in the predicted estimation for larger grain sizes (e.g., 100 and 1000 $\mu$m radii). Owing to the Rayleigh effect on single-scattering properties, neither Hapke approximate models could predict an accurate grain size estimation for the small particles (radii $\leq$ 5 $\mu$m). We recommend that future studies should favor the equivalent slab approximation when employing RTMs for estimating grain sizes of the vast number of TNOs and KBOs in the outer solar system.

Juno Microwave Radiometer (MWR) observations of Jupiter's mid-latitudes reveal a strong correlation between brightness temperature contrasts and zonal winds, confirming that the banded structure extends throughout the troposphere. However, the microwave brightness gradient is observed to change sign with depth: the belts are microwave-bright in the $p<5$ bar range and microwave-dark in the $p>10$ bar range. The transition level (which we call the jovicline) is evident in the MWR 11.5 cm channel, which samples the 5-14 bar range when using the limb-darkening at all emission angles. The transition is located between 4 and 10 bars, and implies that belts change with depth from being NH$_3$-depleted to NH$_3$-enriched, or from physically-warm to physically-cool, or more likely a combination of both. The change in character occurs near the statically stable layer associated with water condensation. The implications of the transition are discussed in terms of ammonia redistribution via meridional circulation cells with opposing flows above and below the water condensation layer, and in terms of the `mushball' precipitation model, which predicts steeper vertical ammonia gradients in the belts versus the zones. We show via the moist thermal wind equation that both the temperature and ammonia interpretations can lead to vertical shear on the zonal winds, but the shear is $\sim50\times$ weaker if only NH$_3$ gradients are considered. Conversely, if MWR observations are associated with kinetic temperature gradients then it would produce zonal winds that increase in strength down to the jovicline, consistent with Galileo probe measurements; then decay slowly at higher pressures.

Stimulated decays of axion dark matter, triggered by a source in the sky, could produce a photon flux along the continuation of the line of sight, pointing backward to the source. The strength of this so-called axion "echo" signal depends on the entire history of the source and could still be strong from sources that are dim today but had a large flux density in the past, such as supernova remnants (SNRs). This echo signal turns out to be most observable in the radio band. We study the sensitivity of radio telescopes such as the Square Kilometer Array (SKA) to echo signals generated by SNRs that have already been observed. In addition, we show projections of the detection reach for signals coming from old SNRs and from newly born supernovae that could be detected in the future. Intriguingly, an observable echo signal could come from old "ghost" SNRs which were very bright in the past but are now so dim that they haven't been observed.

Axions are a theoretically promising dark matter (DM) candidate. In the presence of radiation from bright astrophysical sources at radio frequencies, nonrelativistic DM axions can undergo stimulated decay to two nearly back-to-back photons, meaning that bright sources of radio waves will have a counterimage (''gegenschein'') in nearly the exact opposite spatial direction. The counterimage will be spectrally distinct from backgrounds, taking the form of a narrow radio line centered at $\nu = m_a/4\pi$ with a width determined by Doppler broadening in the DM halo, $\Delta \nu/\nu \sim 10^{-3}$. In this work, we show that the axion decay-induced echoes of supernova remnants may be bright enough to be detectable. Their non-detection may be able to set the strongest limits to date on axion DM in the $\sim 1-10 \, \mu$eV mass range where there are gaps in coverage from existing experiments.

We investigate the possibility of keV scale thermal dark matter in a minimal gauged $B - L$ extension of the standard model with three right-handed neutrinos in the context of cosmic inflation. The complex singlet scalar field responsible for the spontaneous breaking of $B-L$ gauge symmetry is non-minimally coupled to gravity and serves the role of inflaton. The lightest right-handed neutrino $N_1$ can be a dark matter candidate if its couplings to leptons are sufficiently suppressed or forbidden. While keV scale $N_1$ gives rise to the possibility of warm dark matter, its thermal production leads to over-abundance. The subsequent entropy dilution due to heavier right-handed neutrino decay can bring the DM abundance within the observed limit. We constrain the model parameters from the requirement of producing sufficient entropy dilution, inflationary observables along with other phenomenological constraints. While these requirements cannot satisfy light neutrino data, suitable extension of the minimal model can accommodate it.

The dynamics of charged particle in the vicinity of event horizon of a rotating charged black hole (BH) with gravitomagnetic charge from a class of vector-tensor theories of modified gravity immersed in an axially symmetric magnetic field (Bfield) has been studied. The presence of Bfield reduces the radii of innermost stable circular orbit (ISCO) of charged particle and the opposite phenomenon occurs due the parameter $\beta$ which describes the deviation of modified gravity considered here from the usual Einstein - Maxwell gravity also acts to increase the radius of ISCO. Angular momentum also plays the role regarding the motion of particle which implies that larger the angular momentum easier for a particle to escape. We study the stability of orbits using Lyapunov characteristic exponent which implies more stability in the presence of Bfield. Generally, in the presence of Bfield it is easier for a particle to escape than its absence. We find twist in the trajectories of charged particles close to event horizon due to spin of BH in the presence of Bfield and there is no twist in the absence of it. Finally, we comment on possible utilization of our findings for the relativistic jets and magnetohydrodynamical out flows.

Simple climate models have been around for more than a century but have recently come back into fashion: they are useful for explaining global warming and the habitability of extrasolar planets. The Climate App (https://www.climateapp.ca) is an interactive web-based application that describes the radiative transfer governing planetary climate. The App is currently available in French and English and is suitable for teaching high-school through college students, or public outreach. The beginner version can be used to explore the greenhouse effect and planetary albedo, sufficient for explaining anthropogenic climate change, the Faint Young Sun Paradox, the habitability of TRAPPIST planets and other simple scenarios. There is also an advanced option with more atmospheric layers and incorporating the absorption and scattering of shortwave radiation for students and educators wishing a deeper dive into atmospheric radiative transfer. A number of pedagogical activities are being beta tested and rolled out.

If entropy production takes place after dark matter particles abundance is frozen out, dark matter density is reduced. We propose two scenarios of such reduction due to entropy production in electrweak phase transition in the early Universe. We study entropy production in the standard model followed by the simplest extension of the standard model, namely two Higgs doublet model (2HDM). We propose the electroweak phase transition (EWPT) is of second order or of smooth crossover in the former scenario while it is of first order in the later. We calculate the entropy release in these scenarios and the corresponding dilution of preexisting dark matter density in the early universe.

In this paper, we consider the possibilities of generating baryon number asymmetry in thermal equilibrium within the frameworks of teleparallel and symmetric teleparallel gravities. Through the derivative couplings of the torsion scalar or the non-metricity scalar to baryons, the baryon number asymmetry is indeed produced in the radiation dominated epoch. For gravitational baryogenesis mechanisms in these two frameworks, the produced baryon-to-entropy ratio is too small to be consistent with observations. But the gravitational leptogenesis models within both frameworks have the possibilities to interpret the observed baryon-antibaryon asymmetry.

Primordial black holes (PBHs) in the mass range $10^{17} - 10^{23}~{\rm gm}$ are considered as possible dark matter candidates as they are not subject to big-bang nucleosynthesis constraints and behave like cold dark matter. If PBHs are indeed dark matter, they cannot be treated as isolated objects in asymptotic flat space-time. Furthermore, when compared to stellar-mass black holes, the rate at which the Hawking particles radiate out from PBHs is significantly faster. In this work, we obtain an exact time-dependent solution that models evaporating black holes in the cosmological background. As a result, the solution considers all three aspects of PBHs -- mass-loss due to Hawking radiation, black hole surrounded by mass distribution, and cosmological background. Furthermore, our model predicts that the decay of PBHs occurs faster for larger masses; however, \emph{the decay rate reduces for lower mass}. Finally, we discuss the implications of theoretical constraints on PBHs as dark matter.

Using a holographic derivation of a quantum effective action for a scalar operator at strong coupling, we compute quasi-equilibrium parameters relevant for the gravitational wave signal from a first order phase transition in a simple dual model. We discuss how the parameters of the phase transition vary with the effective number of degrees of freedom of the dual field theory. Our model can produce an observable signal at LISA if the critical temperature is around a TeV, in a parameter region where the field theory has an approximate conformal symmetry.