Papers for Friday, May 28 2021

Spin$-$orbit alignment (SOA; i.e., the vector alignment between the halo spin and the orbital angular momentum of neighboring halos) provides an important clue to how galactic angular momenta develop. For this study, we extract virial-radius-wise contact halo pairs with mass ratios between 1/10 and 10 from a set of cosmological $N$-body simulations. In the spin--orbit angle distribution, we find a significant SOA in that 52.7%$\pm$0.2% of neighbors are on the prograde orbit. The SOA of our sample is mainly driven by low-mass target halos ($<10^{11.5}h^{-1}M_{\odot}$) with close merging neighbors, corroborating the notion that the tidal interaction is one of the physical origins of SOA. We also examine the correlation of SOA with the adjacent filament and find that halos closer to the filament show stronger SOA. Most interestingly, we discover for the first time that halos with the spin parallel to the filament experience most frequently the prograde-polar interaction (i.e., fairly perpendicular but still prograde interaction; spin--orbit angle $\sim$ 70$^{\circ}$). This instantly invokes the spin-flip event and the prograde-polar interaction will soon flip the spin of the halo to align it with the neighbor's orbital angular momentum. We propose that the SOA originates from the local cosmic flow along the anisotropic large-scale structure, especially that along the filament, and grows further by interactions with neighbors.

Context. Observations of exoplanets indicate the existence of several correlations in the architecture of planetary systems. Exoplanets within a system tend to be of similar size and mass, evenly spaced and often ordered in size and mass. Small planets are frequently packed in tight configurations, while large planets often have wider orbital spacing. Together, these correlations are called the peas in a pod trends in the architecture of planetary systems. Aims. In this paper, these trends are investigated in theoretically simulated planetary systems and compared with observations. Whether these correlations emerge from astrophysical processes or the detection biases of the transit method is examined. Methods. Using the Generation III Bern Model, synthetic planetary systems are simulated. KOBE, a new computer code, simulates the geometrical limitations of the transit method and applies the detection biases and completeness of the Kepler survey. This allows simulated planetary systems to be confronted with observations. Results. The architecture of synthetic planetary systems, observed via KOBE, show the peas in a pod trends in good agreement with observations. These correlations are also present in the theoretical underlying population, from the Bern Model, indicating that these trends are probably of astrophysical origin. Conclusions. Physical processes, involved in planet formation, are responsible for the emergence of evenly spaced planets with similar sizes and masses. The size/mass similarity trends are primordial and originate from the oligarchic growth of protoplanetary embryos and the uniform growth of planets at early times. Later stages in planet formation allows planets, within a system, to grow at different rates thereby decreasing these correlations. The spacing and packing correlations are absent at early times and arise from dynamical interactions.

The sub-Saturn ($\sim$4--8$R_{\oplus}$) occurrence rate rises with orbital period out to at least $\sim$300 days. In this work we adopt and test the hypothesis that the decrease in their occurrence towards the star is a result of atmospheric mass loss, which can transform sub-Saturns into sub-Neptunes ($\lesssim$4$R_{\oplus}$) more efficiently at shorter periods. We show that under the mass loss hypothesis, the sub-Saturn occurrence rate can be leveraged to infer their underlying core mass function, and by extension that of gas giants. We determine that lognormal core mass functions peaked near $\sim$10--20$M_{\oplus}$ are compatible with the sub-Saturn period distribution, the distribution of observationally-inferred sub-Saturn cores, and gas accretion theories. Our theory predicts that close-in sub-Saturns should be less common and $\sim$30\% more massive around rapidly rotating stars; this should be directly testable for stars younger than $\lesssim$1 Gyr. We also predict that the sub-Jovian desert becomes less pronounced and opens up at smaller orbital periods around M stars compared to solar-type stars ($\sim$0.7 days vs. $\sim$3 days). We demonstrate that exceptionally low-density sub-Saturns, "Super-Puffs", can survive intense hydrodynamic escape to the present day if they are born with even larger atmospheres than they currently harbor; in this picture, Kepler 223 d began with an envelope $\sim$1.5$\times$ the mass of its core and is currently losing its envelope at a rate $\sim$2$\times 10^{-3}M_{\oplus}~\mathrm{Myr}^{-1}$. If the predictions from our theory are confirmed by observations, the core mass function we predict can also serve to constrain core formation theories of gas-rich planets.

Over the past decade, rest-frame color-color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color-color selection tool -- the rest-frame $U-V$ vs. $V-J$ diagram -- by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties in UVJ space as well as the color-color evolution of massive galaxies as they quench at redshifts $z\sim 1-2$. We find that our models broadly reproduce the observed UVJ diagram at $z=1-2$; however, our models do not produce a clear bimodality in UVJ space, largely due to the overpopulation of the green valley in SIMBA. We predict increasing $A_V$ as galaxies move toward redder $U-V$ and $V-J$ colors, with attenuation curves becoming flatter (greyer) with increasing mass. This latter trend results in both relatively muted trends between inferred sSFRs and $U-V$ and $V-J$ colors, as well as a lack of the reddest colors in our most massive systems. When investigating the time evolution of galaxies on the UVJ diagram, we find that the quenching pathway on the UVJ diagram is independent of the quenching timescale, and instead dependent primarily on the average sSFR in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies.

A crucial question in galaxy formation is what role new accretion has in star formation. Theoretical models have predicted a wide range of correlation strengths between halo accretion and galaxy star formation. Previously, we presented a technique to observationally constrain this correlation strength for isolated Milky Way-mass galaxies at $z\sim 0.12$, based on the correlation between halo accretion and the density profile of neighbouring galaxies. By applying this technique to both observational data from the Sloan Digital Sky Survey and simulation data from the UniverseMachine, where we can test different correlation strengths, we ruled out positive correlations between dark matter accretion and recent star formation activity. In this work, we expand our analysis by (1) applying our technique separately to red and blue neighbouring galaxies, which trace different infall populations, (2) correlating dark matter accretion rates with $D_{n}4000$ measurements as a longer-term quiescence indicator than instantaneous star-formation rates, and (3) analyzing higher-mass isolated central galaxies with $10^{11.0} < M_*/M_\odot < 10^{11.5}$ out to $z\sim 0.18$. In all cases, our results are consistent with non-positive correlation strengths with $\gtrsim 85$ per cent confidence, suggesting that processes such as gas recycling dominate star formation in massive $z=0$ galaxies.

We present elemental abundances and stellar population ages for 68 massive quiescent galaxies at $0.59\leq z\leq0.75$ from the LEGA-C survey. The abundance patterns and ages, derived from full-spectrum modeling, are examined as a function of stellar mass ($M_*$) and size (i.e., half-light radius; $R_e$). We find that both [Mg/H] and [Fe/H] do not vary with stellar mass but are correlated with $M_*/R_e$ for quiescent galaxies with $M_*>10^{10.5} M_\odot$. Thus, at fixed mass, compact quiescent galaxies are on average more metal rich. This result reinforces the picture that supernova feedback and gravitational potential regulate chemical enrichment. [Mg/Fe] does not vary with $M_*$ or $M_*/R_e$, but there is a marginal positive relation between age and mass. Our results support low-redshift findings that more massive galaxies form their stars at earlier times. However, in contrast to low-redshift studies, star formation timescale does not appear to depend on mass or size. We also compare the mass-[Fe/H] and mass-[Mg/H] relations to stacks of quiescent galaxies at $z\sim0$ and find that both relations increase by $\sim0.2$ dex over the past 7 Gyr. Furthermore, at $z\sim0.7$ we find a clear trend with age, such that older quiescent galaxies have lower metallicities. Both results can be explained by a chemical evolution model in which galaxies quench via gas removal. Future work, in particular with JWST/NIRSpec, will extend this analysis to higher redshifts, allowing us to fully exploit abundance patterns to study the formation histories of quiescent galaxies.

This paper studied the faint, diffuse extended X-ray emission associated with the radio lobes and the hot gas in the intracluster medium (ICM) environment for a sample of radio galaxies. We used shallow ($\sim 10$ ks) archival Chandra observations for 60 radio galaxies (7 FR I and 53 FR II) with $0.0222 \le z \le 1.785$ selected from the 298 extragalactic radio sources identified in the 3CR catalog. We used Bayesian statistics to look for any asymmetry in the extended X-ray emission between regions that contain the radio lobes and regions that contain the hot gas in the ICM. In the Chandra broadband ($0.5 - 7.0$ keV), which has the highest detected X-ray flux and the highest signal-to-noise ratio, we found that the non-thermal X-ray emission from the radio lobes dominates the thermal X-ray emission from the environment for $\sim 77\%$ of the sources in our sample. We also found that the relative amount of on-jet axis non-thermal emission from the radio lobes tends to increase with redshift compared to the off-jet axis thermal emission from the environment. This suggests that the dominant X-ray mechanism for the non-thermal X-ray emission in the radio lobes is due to the inverse Compton upscattering of cosmic microwave background (CMB) seed photons by relativistic electrons in the radio lobes, a process for which the observed flux is roughly redshift independent due to the increasing CMB energy density with increasing redshift.

We have used the Arecibo Telescope and the Green Bank Telescope to carry out a deep search for H{\sc i}~21\,cm emission from a large sample of Green Pea galaxies, yielding 19 detections, and 21 upper limits on the H{\sc i} mass. We obtain H{\sc i} masses of $\rm M_{HI} \approx (4-300) \times 10^8 \, \rm M_\odot$ for the detections, with a median H{\sc i} mass of $\approx 2.6 \times 10^9 \, \rm M_\odot$; for the non-detections, the median $3\sigma$ upper limit on the H{\sc i} mass is $\approx 5.5 \times 10^8 \, \rm M_\odot$. These are the first estimates of the atomic gas content of Green Pea galaxies. We find that the H{\sc i}-to-stellar mass ratio in Green Peas is consistent with trends identified in star-forming galaxies in the local Universe. However, the median H{\sc i} depletion timescale in Green Peas is $\approx 0.6$~Gyr, an order of magnitude lower than that obtained in local star-forming galaxies. This implies that Green Peas consume their atomic gas on very short timescales. A significant fraction of the Green Peas of our sample lie $\gtrsim 0.6$~dex ($2\sigma$) above the local $\rm M_{HI} - M_B$ relation, suggesting recent gas accretion. Further, $\approx 30$\% of the Green Peas are more than $\pm 2\sigma$ deviant from this relation, suggesting possible bimodality in the Green Pea population. We obtain a low H{\sc i}~21\,cm detection rate in the Green Peas with the highest O32~$\equiv$~[O{\sc iii}]$\lambda$5007/[O{\sc ii}]$\lambda$3727 luminosity ratios, O32~$> 10$, consistent with the high expected Lyman-continuum leakage from these galaxies.

GW190521 is the most massive binary black hole (BBH) merger observed to date, and its primary component lies in the pair-instability (PI) mass gap. Here, we investigate the formation of GW190521-like systems via three-body encounters in young massive star clusters. We performed $10^5$ simulations of binary-single interactions between a BBH and a massive $\geq{60}\,$M$_{\odot}$ black hole (BH), including post-Newtonian terms up to the $2.5$ order and a prescription for relativistic kicks. In our initial conditions, we take into account the possibility of forming BHs in the PI mass gap via stellar collisions. If we assume that first-generation BHs have low spins, $\sim{0.17}\%$ of all the simulated BBH mergers have component masses, effective and precessing spin, and remnant mass and spin inside the $90\%$ credible intervals of GW190521. Half of these systems are first-generation exchanged binaries, while the other half are second-generation BBHs. We estimate a merger rate density $\mathcal{R}_{\rm GW190521}\sim{0.03}\,$Gpc$^{-3}\,$yr$^{-1}$ for GW190521-like binaries formed via binary-single interactions in young star clusters. This rate is extremely sensitive to the spin distribution of first-generation BBHs. Stellar collisions, second-generation mergers and dynamical exchanges are the key ingredients to produce GW190521-like systems in young star clusters.

Observations of hot-Jupiter atmospheres show large variations in the location of the hot spot and the amplitude of spectral features. Atmospheric flow simulations using the commonly-employed forcing and initialization have generally produced a large, monolithic patch of stationary hot area located eastward of the substellar point at $\sim 3\!\times\! 10^{-3}$ MPa pressure level. Here we perform high-resolution (up to T682) pseudospectral simulations that accurately capture small-scale eddies and waves, inherent in hot-Jupiter atmospheres due to ageostrophy. The atmospheres contain a large number of intense storms over a wide range of scales, including the planetary-scale. The latter sized storms dictate the large-scale spatial distribution and temporal variability of hot, as well as cold, regions over the planet. In addition, the large storms exhibit quasi-periodic life cycles within multiple equilibrium states -- all identifiable in the disk-integrated time series of the temperature flux.

Stellar population models produce radiation fields that ionize oxygen up to O+2, defining the limit of standard HII region models (<54.9 eV). Yet, some extreme emission line galaxies, or EELGs, have surprisingly strong emission originating from much higher ionization potentials. We present UV-HST/COS and optical-LBT/MODS spectra of two nearby EELGs that have very-high-ionization emission lines (e.g., HeII 1640,4686 CIV 1548,1550, [FeV] 4227, [ArIV] 4711,4740). We define a 4-zone ionization model that is augmented by a very-high-ionization zone, as characterized by He+2 (>54.4 eV). The 4-zone model has little to no effect on the measured total nebular abundances, but does change the interpretation of other EELG properties: we measure steeper central ionization gradients, higher volume-averaged ionization parameters, and higher central T_e, n_e, and logU values. Traditional 3-zone estimates of the ionization parameter can under-estimate the average log U by up to 0.5 dex. Additionally, we find a model-independent dichotomy in the abundance patterns, where the alpha/H-abundances are consistent but N/H, C/H, and Fe/H are relatively deficient, suggesting these EELGs are alpha/Fe-enriched by >3 times. However, there still is a high-energy ionizing photon production problem (HEIP^3). Even for such alpha/Fe-enrichment and very-high log Us, photoionization models cannot reproduce the very-high-ionization emission lines observed in EELGs.

We present a chronology on the formation and early evolution of the Oort cloud, and test the sequence of events of its formation by simulating the formation process in subsequent amalgamated steps. These simulations start with the Solar system being born with planets and asteroids in a stellar cluster orbiting the Galactic center. Upon ejection from its birth environment, we continue to follow the Solar system's evolution while it sojourns the Galaxy as an isolated planetary system. We conclude that the range in semi-major axis between $\sim 100$\,au and several $\sim 10^3$\,au still bears the signatures of the Sun being born in a $\apgt 1000$\,\Msun/pc$^3$ star cluster, and that most of the outer Oort cloud formed after the Solar system escaped. The escape, we argue, happened between $\sim 20$\,Myr and $50$\,Myr after birth of the Solar system. Today, the bulk of the material in the Oort cloud ($\sim 70$\%) originates from the region in the circumstellar disk that was located between $\sim 15$\,au and $\sim 35$\,au, near the current location of the ice-giants and the Centaur family of asteroids. This population is eradicated if the ice-giant planets were born in orbital resonance. Planet migration or chaotic orbital reorganization, occurring while the Solar system is still a cluster member is, according to our model, inconsistent with the presence of the Oort cloud. About half the inner Oort cloud, between $100$ and $10^4$\,au, and a quarter of the material in the outer Oort cloud $\apgt 10^4$\,au could be non-native to the Solar system but was captured from free-floating derbis in the cluster or from the circumstellar disk of other stars in the birth cluster. Characterizing this population will help us to reconstruct the Solar system's history.

HR 8799 hosts four directly imaged giant planets, but none has a mass measured from first principles. We present the first dynamical mass measurement in this planetary system, finding that the innermost planet HR~8799~e has a mass of $9.6^{+1.9}_{-1.8} \, M_{\rm Jup}$. This mass results from combining the well-characterized orbits of all four planets with a new astrometric acceleration detection (5$\sigma$) from the Gaia EDR3 version of the Hipparcos-Gaia Catalog of Accelerations. We find with 95\% confidence that HR~8799~e is below $13\, M_{\rm Jup}$, the deuterium-fusing mass limit. We derive a hot-start cooling age of $42^{+24}_{-16}$\,Myr for HR~8799~e that agrees well with its hypothesized membership in the Columba association but is also consistent with an alternative suggested membership in the $\beta$~Pictoris moving group. We exclude the presence of any additional $\gtrsim$5-$M_{\rm Jup}$ planets interior to HR~8799~e with semi-major axes between $\approx$3-16\,au. We provide proper motion anomalies and a matrix equation to solve for the mass of any of the planets of HR~8799 using only mass ratios between the planets.

We report spectropolarimetric observations of a supersonic downflow impacting the lower atmosphere within a large sunspot umbra. This work is an extension of Schad et al. 2016 using observations acquired in the He I 10830 Angstrom triplet by the Facility Infrared Spectropolarimeter. Downflowing material accelerating along a cooled coronal loop reaches peak speeds near 200 km s$^{-1}$ and exhibits both high speed emission and absorption within the umbra, which we determine to be a consequence of the strong height dependence of the radiatively-controlled source function above the sunspot umbra. Strong emission profiles close to the rest wavelengths but with long red-shifted tails are also observed at the downflow terminus. From the polarized spectra, we infer longitudinal magnetic field strengths of ${\sim}2.4$ kG in the core portion of the He I strong emission, which we believe is the strongest ever reported in this line. Photospheric field strengths along the same line-of-sight are ${\sim}2.8$ kG as inferred using the Ca I 10839 Angstrom spectral line. The temperatures of the highest speed He I absorption and the near rest emission are similar (${\sim}$10 kK), while a differential emission measure analysis using SDO/AIA data indicates significant increases in radiative cooling for temperatures between $\sim$0.5 and 1 MK plasma associated with the downflow terminus. Combined we interpret these observations in the context of a strong radiative shock induced by the supersonic downflow impacting the low sunspot atmosphere.

We investigate a sample of 6 Herbig Ae/Be stars belonging to the Orion OB1 association, as well as 73 low mass objects, members of the $\sigma$ Orionis cluster, in order to explore the angular momentum evolution at early stages of evolution, and its possible connection with main-sequence Ap/Bp magnetic stars. Using FIES and HECTOCHELLE spectra, we obtain projected rotational velocities through two independent methods. Individual masses, radii, and ages were computed using evolutionary models, distance, and cluster extinction. Under the assumption that similar physical processes operate in both, T Tauri and Herbig Ae/Be stars, we construct snapshots of the protostar's rotation against mass during the first 10 Myr with the aid of a rotational model that includes a variable disc lifetime, changes in the stellar moment of inertia, a dipolar magnetic field with variable strength, and angular momentum loss through stellar winds powered by accretion. We use these snapshots, as well as the rotational data, to infer a plausible scenario for the angular momentum evolution. We find that magnetic field strengths of a few k$G$ at 3 Myr are required to match the rotational velocities of both groups of stars. Models with masses between 2-3 $M_{\odot}$ display larger angular momentum by a factor of $\sim 3$, in comparison to stars of similar spectral types on the main-sequence. Even though some quantitative estimates on this dramatic decrease with age, for Ap/Bp magnetic main-sequence stars are presented, the results obtained for the angular momentum evolution do not explain their low rotation.

Multi band photometry and light curve analysis for two newly recognized contact binary systems, TYC 6995-813-1 and NSVS 13602901 are presented. Both were found to be of extreme low mass ratios 0.11 and 0.17, respectively. The secondary components of both systems show evidence of considerable evolution with elevated densities as well as both luminosity and radii well above their main sequence counterparts. Even in the absence of significant spot activity at least one of the systems, TYC 6995-813-1, shows features of magnetic and chromospheric activity. TYC 6995-813-1 is also determined to be a potential merger candidate with its current separation near the theoretical instability separation.

Redshift-space distortions (RSD) generically affect any spatially-dependent observable that is mapped using redshift information. The effect on the observed clustering of galaxies is the primary example of this. This paper is devoted to another example: the effect of RSD on the apparent peculiar motions of tracers as inferred from their positions in redshift space (i.e. the observed distance). Our theoretical study is motivated by practical considerations, mainly, the direct estimation of the velocity power spectrum, which is preferably carried out using the tracer's redshift-space position (so as to avoid uncertainties in distance measurements). We formulate the redshift-space velocity field and show that RSD enters as a higher-order effect. Physically, this effect may be interpreted as a dissipative correction to the usual perfect-fluid description of dark matter. We show that the effect on the power spectrum is a damping on relatively large, quasilinear scales ($k>0.01\,h\,{\rm Mpc}^{-1}$), as was observed, though unexplained, in $N$-body simulations elsewhere. This paper presents two power spectrum models for the the peculiar velocity field in redshift space, both of which can be considered velocity analogues of existing clustering models. In particular, we show that the "Finger-of-God" effect, while also present in the velocity field, cannot be entirely blamed for the observed damping in simulations. Our work provides some of the missing modelling ingredients required for a density--velocity multi-tracer analysis, which has been proposed for upcoming redshift surveys.

We present here the precise wavelength calibration of a high-resolution spectrum using Uranium lines in the wavelength range of 3809 - 6833 \AA for precision radial velocity measurements for exoplanet detection or related astrophysical sciences. The spectrum is acquired using the PARAS spectrograph (R=67,000) attached to a 1.2 meter telescope at Mount Abu Observatory, India. We identify well-resolved 1540 U lines from a high-resolution spectrum of the UAr hollow cathode lamp (HCL) using PARAS spectrograph in the aforesaid wavelength range. We calculate the Ritz wavelength of U from its known energy levels and compare them with our observed central wavelengths. The comparisons do not show any significant offset in our line list on the absolute wavelength scale. The final U line list has an average measurement uncertainty of 0.28 m\AA. We included these lines to the PARAS data analysis framework to do the wavelength calibration and then calculate the multi-order radial velocity of PARAS spectra. The typical dispersion of residuals around the wavelength solution of a UAr spectrum, using U lines, is found to be 0.8 m\AA. We show our result in the precision radial velocity of an on-sky source (A RV standard star), and an off-sky source (A HCL) observed with PARAS along with UAr HCL. We measure the dispersion in absolute drift difference between two fibers (inter-fiber drift) for a span of 6.5 hours to be 88 cm s$^{-1}$, and the ${\sigma_{RV}}$ for a RV standard star, HD55575 for 450 days to be 3.2 m s$^{-1}$. Comparing these results with the previous ones measured using ThAr HCL shows that the ThAr HCL with about 99\% pure-Th, usually used for precise wavelength calibration, are replaceable with the UAr HCL in the visible-domain for PARAS like spectrograph (R $\leq$ 67,000) to achieve a radial velocity precision of 1-3 m s$^{-1}$.

We reconstruct the Hubble function from cosmic chronometers, supernovae, and baryon acoustic oscillations compiled data sets via the Gaussian process (GP) method and use it to draw out Horndeski theories that are fully anchored on expansion history data. In particular, we consider three well-established formalisms of Horndeski gravity which single out a potential through the expansion data, namely: quintessence potential, designer Horndeski, and tailoring Horndeski. We discuss each method in detail and complement it with the GP reconstructed Hubble function to obtain predictive constraints on the potentials and the dark energy equation of state.

Announcing the release v7.2 of the Milliquas (Million Quasars) catalogue which presents all published quasars to 30 April 2021, including VLASS radio associations for the first time, and concluding the audit of quasars from SDSS-DR16Q and earlier SDSS releases. The totals are 829666 classified type-I QSOs/AGN, 703348 quasar candidates of 60%-100% pQSO, plus type-II objects and blazars which bring the total count to 1573824. Radio and/or X-ray associations, including probable double radio lobes, are shown for 333638 entries. Gaia-DR2 astrometry is given for most objects, as available. The catalogue is available on multiple sites. The inclusion of the SDSS-DR16Q quasars was a complex task with emergent issues which resulted in 13443 DR16Q entries being dropped, 1.79% of their total. There are also 1701 quasars included from earlier visual SDSS releases, as well as 14232 quasars from the SDSS-DR16 pipeline catalogue. All these are explained here, including the validation of 677 additional high-redshift (z>=3.5) SDSS quasars which were not included in DR16Q.

We report the detection, for the first time in space, of cyano thioformaldehyde (HCSCN) and propynethial (HCSCCH) towards the starless core TMC-1. Cyano thioformaldehyde presents a series of prominent a- and b-type lines, which are the strongest previously unassigned features in our Q-band line survey of TMC-1. Remarkably, HCSCN is four times more abundant than cyano formaldehyde (HCOCN). On the other hand, HCSCCH is five times less abundant than propynal (HCOCCH). Surprisingly, we find an abundance ratio HCSCCH/HCSCN of 0.25, in contrast with most other ethynyl-cyanide pairs of molecules for which the CCH-bearing species is more abundant than the CN-bearing one. We discuss the formation of these molecules in terms of neutral-neutral reactions of S atoms with CH2CCH and CH2CN radicals as well as of CCH and CN radicals with H2CS. The calculated abundances for the sulphur-bearing species are, however, significantly below the observed values, which points to an underestimation of the abundance of atomic sulphur in the model or to missing formation reactions, such as ion-neutral reactions.

In a recent MNRAS article, Raposo-Pulido and Pelaez (RPP) designed a scheme for obtaining very close seeds for solving the elliptic Kepler Equation with the classical and the modified Newton-Rapshon methods. This implied an important reduction in the number of iterations needed to reach a given accuracy. However, RPP also made strong claims about the errors of their method that are incorrect. In particular, they claim that their accuracy can always reach the level $\sim5\varepsilon$, where $\varepsilon$ is the machine epsilon (e.g. $\varepsilon=2.2\times10^{-16} $ in double precision), and that this result is attained for all values of the eccentricity $e<1$ and the mean anomaly $M\in[0,\pi]$, including for $e$ and $M$ that are arbitrarily close to $1$ and $0$, respectively. However, we demonstrate both numerically and analytically that any implementation of the classical or modified Newton-Raphson methods for Kepler's equation, including those described by RPP, have a limiting accuracy of the order $\sim\varepsilon/\sqrt{2(1-e)}$. Therefore the errors of these implementations diverge in the limit $e\to1$, and differ dramatically from the incorrect results given by RPP. Despite these shortcomings, the RPP method can provide a very efficient option for reaching such limiting accuracy. We also provide a limit that is valid for the accuracy of any algorithm for solving Kepler equation, including schemes like bisection that do not use derivatives. Moreover, similar results are also demonstrated for the hyperbolic Kepler Equation. The methods described in this work can provide guidelines for designing more accurate solutions of the elliptic and hyperbolic Kepler equations.

We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year old type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29$M_{\odot}$, all evolved with KEPLER at solar metallicity and ignoring rotation. The [OI]$\lambda\lambda$$6300,6364$ line flux generally grows with progenitor mass, and H$\alpha$ exhibits an equally strong and opposite trend. The [CaII]$\lambda\lambda$$7291,\,7323$ strength increases at low $^{56}$Ni mass, low explosion energy, or with clumping. This CaII doublet, which forms primarily in the explosively-produced Si/S zones, depends little on the progenitor mass, but may strengthen if Ca$^+$ dominates in the H-rich emitting zones or if Ca is abundant in the O-rich zones. Indeed, Si-O shell merging prior to core collapse may boost the CaII doublet at the expense of the OI doublet, and may thus mimic the metal line strengths of a lower mass progenitor. We find that the $^{56}$Ni bubble effect has a weak impact, probably because it is too weak to induce much of an ionization shift in the various emitting zones. Our simulations compare favorably to observed SNe II, including SN2008bk (e.g., 9$M_{\odot}$ model), SN2012aw (12$M_{\odot}$ model), SN1987A (15$M_{\odot}$ model), or SN2015bs (25$M_{\odot}$ model with no Si-O shell merging). SNe II with narrow lines and a low $^{56}$Ni mass are well matched by the weak explosion of 9$-$11$M_{\odot}$ progenitors. The nebular-phase spectra of standard SNe II can be explained with progenitors in the mass range 12$-$15$M_{\odot}$, with one notable exception for SN2015bs. In the intermediate mass range, these mass estimates may increase by a few $M_{\odot}$ with allowance for clumping of the O-rich material or CO molecular cooling.

We present a novel approach based on artificial neural networks, so-called geodesyNets, and present compelling evidence of their ability to serve as accurate geodetic models of highly irregular bodies using minimal prior information on the body. The approach does not rely on the body shape information but, if available, can harness it. GeodesyNets learn a three-dimensional, differentiable, function representing the body density, which we call neural density field. The body shape, as well as other geodetic properties, can easily be recovered. We investigate six different shapes including the bodies 101955 Bennu, 67P Churyumov-Gerasimenko, 433 Eros and 25143 Itokawa for which shape models developed during close proximity surveys are available. Both heterogeneous and homogeneous mass distributions are considered. The gravitational acceleration computed from the trained geodesyNets models, as well as the inferred body shape, show great accuracy in all cases with a relative error on the predicted acceleration smaller than 1\% even close to the asteroid surface. When the body shape information is available, geodesyNets can seamlessly exploit it and be trained to represent a high-fidelity neural density field able to give insights into the internal structure of the body. This work introduces a new unexplored approach to geodesy, adding a powerful tool to consolidated ones based on spherical harmonics, mascon models and polyhedral gravity.

Recently it has become apparent that the Galactic center excess (GCE) is spatially correlated with the stellar distribution in the Galactic bulge. This has given extra motivation for the unresolved population of millisecond pulsars (MSPs) explanation for the GCE. However, in the "recycling" channel the neutron star forms from a core collapse supernovae that undergoes a random "kick" due to the asymmetry of the explosion. This would imply a smoothing out of the spatial distribution of the MSPs. We use N-body simulations to model how the MSP spatial distribution changes. We estimate the probability distribution of natal kick velocities using the resolved gamma-ray MSP proper motions, where MSPs have velocities relative to the circular motion of 77 +/- 6 km/s. We find that, due to the natal kicks, there is an approximately 10% increase in each of the bulge MSP spatial distribution dimensions and also the bulge MSP distribution becomes less boxy but is still far from being spherical.

We present an ultra-violet study of two nearby dwarf irregular galaxies WLM and IC~2574, using the Far-UV and Near-UV data from the Ultra-Violet Imaging Telescope (UVIT). We used the F148W band Far-UV images and identified 180 and 782 young star-forming clumps in WLM and IC~2574, respectively. The identified clumps have sizes between 7 - 30 pc in WLM and 26 - 150 pc in IC~2574. We noticed more prominent hierarchical splitting in the structure of star-forming regions at different flux levels in IC~2574 than WLM. We found that the majority of the clumps have elongated shapes in the sky plane with ellipticity ($\epsilon$) greater than 0.6 in both the galaxies. The major axis of the identified clumps is found to show no specific trend of orientation in IC~2574, whereas in WLM the majority are aligned along south-west to north-east direction. We estimated (F148W$-$N242W) colour for the clumps identified in WLM and noticed that the younger ones (with (F148W$-$N242W) $<-0.5$) are smaller in size ($<10$ pc) and are located mostly in the southern half of the galaxy between galactocentric radii 0.4 - 0.8 kpc.

We present a search for dark photon dark matter that could couple to gravitational-wave interferometers using data from Advanced LIGO and Virgo's third observing run. To perform this analysis, we use two methods, one based on cross-correlation of the strain channels in the two nearly aligned LIGO detectors, and one that looks for excess power in the strain channels of the LIGO and Virgo detectors. The excess power method optimizes the Fourier Transform coherence time as a function of frequency, to account for the expected signal width due to Doppler modulations. We do not find any evidence of dark photon dark matter with a mass between $m_{\rm A} \sim 10^{-14}-10^{-11}$ eV/$c^2$, which corresponds to frequencies between 10-2000 Hz, and therefore provide upper limits on the square of the minimum coupling of dark photons to baryons, i.e. $U(1)_{\rm B}$ dark matter. For the cross-correlation method, the best median constraint on the squared coupling is $\sim1.31\times10^{-47}$ at $m_{\rm A}\sim4.2\times10^{-13}$ eV/$c^2$; for the other analysis, the best constraint is $\sim 1.2\times 10^{-47}$ at $m_{\rm A}\sim 5.7\times 10^{-13}$ eV/$c^2$. These limits improve upon those obtained in direct dark matter detection experiments by a factor of $\sim100$ for $m_{\rm A}\sim [2-4]\times 10^{-13}$ eV/$c^2$.

The low water content of the terrestrial planets in the solar system suggests that the protoplanets formed within the water snow line. Accurate prediction of the snow line location moving with time provides a clue to constrain the formation process of the planets. In this paper, we investigate the migration of the snow line in protoplanetary disks whose accretion is controlled by laminar magnetic fields, which have been proposed by various nonideal magnetohydrodynamic (MHD) simulations. We propose an empirical model of the disk temperature based on our nonideal MHD simulations, which show that the accretion heating is significantly less efficient than in turbulent disks, and calculate the snow line location over time. We find that the snow line in the magnetically accreting laminar disks moves inside the current Earth's orbit within 1 Myr after star formation, whereas the time for the conventional turbulent disk is much longer than 1 Myr. This result suggests that either the rocky protoplanets formed in such an early phase of the disk evolution, or the protoplanets moved outward to the current orbits after they formed close to the protosun.

Fast radio bursts (FRBs) are mysterious radio transients with millisecond durations. Recently, a periodic activity of $\sim$16 day and a possible periodicity of $\sim$159 day were detected to arise from FRB 180916.J0158+65 and FRB 121102, respectively, and the spin period of a slow-rotation magnetar was further considered to be one of possibilities to explain the periodic activities of repeating FRBs. For isolated neutron stars, the spin evolution suggests that it's difficult to reach several hours. In this work, we mainly focus on the possible maximum spin period of isolated NSs / magnetars dominated by an interaction between star's magnetic field and the disk. We find that the disk wind plays an important role in spin evolution, whose influence varies the power law index in the evolution equation of mass flow rate. For a magnetar without disk wind, the longest spin period is tens of hours. When the disk wind with a classical parameter is involved, the maximum spin period can reach hundreds of hours. But for a much extremely large index of mass flow rate due to disk wind or other angular momentum extraction processes, a spin period of $\sim$(16-160) days is still possible.

We present and perform a detailed analysis of multi-wavelength observations of \thisgrb, an optical bright GRB with an observed reverse shock (RS) signature. Observations of this GRB were acquired with the BOOTES-4 robotic telescope, the \fermi, and the \swift missions. Time-resolved spectroscopy of the prompt emission shows that changes to the peak energy (\Ep) tracks intensity and the low-energy spectral index seems to follow the intensity for the first episode, whereas this tracking behavior is less clear during the second episode. The fit to the afterglow light curves shows that the early optical afterglow can be described with RS emission and is consistent with the thin shell scenario of the constant ambient medium. The late time afterglow decay is also consistent with the prediction of the external forward shock (FS) model. We determine the properties of the shocks, Lorentz factor, magnetization parameters, and ambient density of \thisgrb, and compare these parameters with another 12 GRBs, consistent with having RS produced by thin shells in an ISM-like medium. The value of the magnetization parameter ($R_{\rm B} \approx 18$) indicates a moderately magnetized baryonic dominant jet composition for \thisgrb. We also report the host galaxy photometric observations of \thisgrb obtained with 10.4m GTC, 3.5m CAHA, and 3.6m DOT telescopes and find the host (photo $z$ = $2.8^{+0.7}_{-0.9}$) to be a high mass, star-forming galaxy with a star formation rate of $20 \pm 10 \msun$ $\rm yr^{-1}$.

The advent of multi-messenger astronomy has allowed for new types of source searches by neutrino detectors. We present the results of the first search for 0.5 - 5 GeV astrophysical neutrinos emitted from all compact binary mergers, i.e., binary black hole, neutron star black, mass gap and binary neutron star mergers, detected by the LIGO and Virgo interferometers during their three first runs of observation. We use an innovative approach that lowers the energy threshold from ~10 GeV to ~0.5 GeV and searches for an excess of GeV-scale events during astrophysical transient events. No significant excess was found from the studied mergers, and there is currently no hint of a population of GeV neutrino emitters found in the IceCube data.

The study of non-thermal processes such as synchrotron emission, inverse Compton scattering, bremsstrahlung and pion production is crucial to understand the properties of the Galactic cosmic-ray population, to shed light on their origin and confinement mechanisms, and to assess the significance of exotic signals possibly associated to new physics. We present a public code called HERMES aimed at generating sky maps associated to a variety of multi-messenger and multi-wavelength radiative processes, spanning from the radio domain all the way up to high-energy gamma-ray and neutrino production. We describe the physical processes under consideration, the code concept and structure, and the user interface, with particular focus on the python-based interactive mode. We especially present the modular and flexible design that allows to easily further extend the numerical package according to the user's needs. In order to demonstrate the capabilities of the code, we describe in detail a comprehensive set of sky maps and spectra associated to all physical processes included in the code. We comment in particular on the radio, gamma-ray, and neutrino maps, and mention the possibility to study signals stemming from dark matter annihilation. HERMES can be successfully applied to constrain the properties of the Galactic cosmic-ray population, improve our understanding of the diffuse Galactic radio, gamma--ray, and neutrino emission, and search for signals associated to particle dark matter annihilation or decay.

The present work is an effort to investigate possible radial variations in the solar coronal rotation by analyzing the solar radio emission data at 15 different frequencies (275-1755 MHz) for the period starting from July 1994 to May 1999. We used a time series of disk-integrated radio flux recorded daily at these frequencies through radio telescopes situated at Astronomical Observatory of the Jagellonian University in Cracow. The different frequency radiation originates from different heights in the solar corona. Existing models, indicate its origin at the height range from nearly $\sim12,000$ km (for emission at 275 MHz), below up to $\sim2,400$ km (for emission at 1755 MHz). There are some data gaps in the time series used for the study, so we used statistical analysis using the Lomb-Scargle Periodogram (LSP). This method has successfully estimated the periodicity present in time series even with such data gaps. The rotation period estimated through LSP shows variation in rotation period, which is compared with the earlier reported estimate using auto correlation technique. The present study indicates some similarity as well as contradiction with studies reported earlier. The radial and temporal variation in solar rotation period are presented and discussed for the whole period analyzed.

We present a possible correlation between the properties of scattered and thermal radiation from dust and the principal dust characteristics responsible for this relationship. To this end, we use the NASA/PDS archival polarimetric data on cometary dust in the Red (0.62--0.73 $\mu$m) and K (2.00--2.39 $\mu$m) domains to leverage the relative excess of the polarisation degree of a comet to the average trend at the given phase angle ($P_{\rm excess}$) as a metric of the dust's scattered light characteristics. The flux excess of silicate emissions to the continuum around 10 $\mu$m ($F_{\rm Si}/F_{\rm cont}$) is adopted from previous studies as a metric of the dust's MIR feature. The two metrics show a positive correlation when $P_{\rm excess}$ is measured in the K domain. No significant correlation was identified in the Red domain. The gas-rich comets have systematically weaker $F_{\rm Si}/F_{\rm cont}$ than the dust-rich ones, yet both groups retain the same overall tendency with different slope values. The observed positive correlation between the two metrics indicates that composition is a peripheral factor in characterising the dust's polarimetric and silicate emission properties. The systematic difference in $F_{\rm Si}/F_{\rm cont}$ for gas-rich versus dust-rich comets would rather correspond with the difference in their dust size distribution. Hence, our results suggest that the current MIR spectral models of cometary dust should prioritise the dust size and porosity over the composition. With light scattering being sensitive to different size scales in two wavebands, we expect the K-domain polarimetry to be sensitive to the properties of dust aggregates, such as size and porosity, which might have been influenced by evolutionary processes. On the other hand, the Red-domain polarimetry reflects the characteristics of sub-$\mu$m constituents in the aggregate.

When modified theories of gravity are considered, at most six gravitational wave polarization modes are allowed and classified in tensor modes, the only ones predicted by General Relativity (GR), along with additional vector and scalar modes. Therefore, gravitational waves represent a powerful tool to test alternative theories of gravitation. In this paper, we forecast the sensitivity of third-generation ground-based interferometers, Einstein Telescope and Cosmic Explorer, to non-GR polarization modes focusing on the stochastic gravitational wave background. We consider the latest technical specifications of the two independent detectors and the full network in order to estimate both the optimal signal-to-noise ratio and the detectable energy density limits relative to all polarization modes in the stochastic background for several locations on Earth and orientations of the two observatories. By considering optimal detector configurations, we find that in 5 years of observation the detection limit for tensor and extra polarization modes could reach $h_0^2\Omega^{T,V,S}_{GW} \approx 10^{-12}-10^{-11}$, depending on the network configuration and the stochastic background (i.e., if only one among vector and scalar modes exists or both are present). This means that the network sensitivity to different polarization modes can be approximately improved by a factor $10^3$ with respect to second-generation interferometers. We finally discuss the possibility of breaking the scalar modes degeneracy by considering both detectors angular responses to sufficiently high gravitational wave frequencies.

Adopting various unified equations of state (EOSs), we examine the quasinormal modes of gravitational waves from cold neutron stars. We focus on the fundamental ($f$-), 1st pressure ($p_1$-), and 1st spacetime ($w_1$-) modes, and derive the empirical formulae for the frequencies and damping rate of those modes. With the resultant empirical formulae, we find that the value of $\eta$, which is a specific combination of the nuclear saturation parameters, can be estimated within $\sim 30 \%$ accuracy, if the $f$-mode frequency from the neutron star whose mass is known would be observed or if the $f$- and $p_1$-mode frequencies would be simultaneously observed, even though this estimation is applicable only for the low-mass neutron stars. Additionally, we find that the mass and radius of canonical neutron stars can be estimated within a few per cent accuracy via the simultaneous observations of the $f$- and $w_1$-mode frequencies. We also find that, if the $f$-, $p_1$-, and $w_1$-mode frequencies would be simultaneously observed, the mass of canonical neutron stars can be estimated within $2\%$ accuracy, while the radius can be estimated within $1\%$ for the neutron star with $M\ge 1.6M_\odot$ or within $0.6\%$ for the neutron star with $M\ge 1.4M_\odot$ constructed with the EOS constrained via the GW170817 event. Furthermore, we find the strong correlation between the maximum $f$-mode frequency and the neutron star radius with the maximum mass, between the minimum $w_1$-mode frequency and the maximum mass, and between the minimum damping rate of the $w_1$-mode and the stellar compactness for the neutron star with the maximum mass.

We present a photometric redshift (photo-$z$) estimation technique for galaxies in the P\lowercase{an}-STARRS1 (PS1) $3\pi $ survey. Specifically, we train and test a regression and a classification Random-Forest (RF) models using photometric features (magnitudes, colors and moments of the radiation intensity) from the optical PS1 data release 2 (PS1-DR2) and from the AllWISE/unWISE infrared source catalogs. The classification RF model ($RF_{clas}$) has better performance in the local universe ($z\lesssim 0.1$), while the second one ($RF_{reg}$) is on average better for $0.1 \lesssim z\lesssim1$. We adopt as labels the spectroscopic redshift of the galaxies from the Sloan Digital Sky Survey (SDSS) data release 16 (SDSS-DR16). We find that the combination of AllWISE/unWISE and PS1-DR2 features leads to an average bias of $\overline{\Delta z_{norm}}=1\times 10^{-3}$, a standard deviation $\sigma(\Delta z_{norm})=0.0225$, (where $\Delta z_{norm} \equiv (z_{phot}-z_{spec})/(1+z_{spec})$), and an outlier rate of $P_0=1.48 \%$ in the test set for the $RF_{clas}$ model. In the low-redshift Universe ($z<0.1$) that is of primary interest to many astronomical transient studies, our model produces an error estimate on the inferred magnitude of an object of $\le$1 mag in 87\% of the test sample.

We present here results from a survey of intervening C IV absorbers at $z < 0.16$ conducted using 223 sightlines from the Hubble Spectroscopic Legacy Archive. Most systems (83%) out of the total sample of 69 have simple kinematics with 1 or 2 C IV components. In the 22 C IV systems with well constrained H I column densities, the temperatures from the $b$-values imply predominantly photoionized plasma ($T\leq 10^5$ K) and non-thermal dynamics. These systems also have solar or higher metallicities. We obtain a C IV line density of $d\mathcal{N}/dX = 5.1\pm 1.0$ for $\log [N(C~IV)~(cm^{-2})]\geq12.9$, and $\Omega_{C~IV}=(8.01\pm 1.62) \times 10^{-8}$ for $12.9 \leq \log [N(C~IV)~(cm^{-2})] \leq 15.0$. The C IV bearing diffuse gas in the $z < 0.16$ Universe has a metallicity of $(2.07~{\pm}~0.43)~\times~10^{-3}$ Z$_{\odot}$, an order of magnitude more than the metal abundances in the IGM at high redshifts ($z \gtrsim 5$), and consistent with the slow build-up of metals in the diffuse circum/intergalactic space with cosmic time. For $z<0.015$ (complete above $L>0.01L^\star$), the Sloan Digital Sky Survey provides a tentative evidence of declining covering fraction for strong C IV ($N>10^{13.5}~cm^{-2}$) with $\rho$ (impact parameter) and $\rho/R_\mathrm{vir}$. However, the increase at high separations suggests that strong systems are not necessarily coincident with such galaxies. We also find that strong C IV absorption at $z<0.051$ is not coincident with galaxy over-dense regions complete for $L>0.13L^\star$

Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles, enhancing the volatile content of the disk. We investigate in this work how the inward drift and evaporation of volatile rich pebbles influences the atmospheric C/O ratio and heavy element content of giant planets growing by pebble and gas accretion. To achieve this goal, we perform semi analytical 1D models of protoplanetary disks including the treatment of viscous evolution and heating, pebble drift and simple chemistry to simulate the growth of planets from planetary embryos to Jupiter mass objects by accretion of pebbles and gas while they migrate through the disk. Our simulations show that the composition of the planetary gas atmosphere is dominated by the accretion of vapor, originating from inward drifting evaporating pebbles. This process allows the giant planets to harbour large heavy element contents, in contrast to models that do not take pebble evaporation into account. In addition, our model reveals that giant planets originating further away from the central star have a higher C/O ratio on average due to the evaporation of methane rich pebbles in the outer disk. However, planets formed in the outer disk harbor a smaller heavy element content, due to a smaller vapor enrichment of the outer disk compared to the inner disk. Our model predicts that giant planets with low/large atmospheric C/O should harbour a large/low total heavy element content. We further conclude that the inclusion of pebble evaporation is a key ingredient to determine the heavy element content and composition of giant planets.

Gravitational waves are a radically new way to peer into the darkest depths of the cosmos. Pulsars can be used to make direct detections of gravitational waves through precision timing. When a gravitational wave passes between a pulsar and the Earth, it stretches and squeezes the intermediate space-time, leading to deviations of the measured pulse arrival times away from model expectations. Combining the data from many Galactic pulsars can corroborate such a signal, and enhance its detection significance. This technique is known as a Pulsar Timing Array (PTA). Here I provide an overview of PTAs as a precision gravitational-wave detection instrument, then review the types of signal and noise processes that we encounter in typical pulsar data analysis. I take a pragmatic approach, illustrating how searches are performed in real life, and where possible directing the reader to codes or techniques that they can explore for themselves. The goal is to provide theoretical background and practical recipes for data exploration that allow the reader to join in the exciting hunt for very low frequency gravitational waves.

We present the results of a systematic search for new rapidly oscillating Ap (roAp) stars using the 2-min cadence data collected by the Transiting Exoplanet Survey Satellite (TESS) during its Cycle 1 observations. We identify 12 new roAp stars. Amongst these stars we discover the roAp star with the longest pulsation period, another with the shortest rotation period, and six with multiperiodic variability. In addition to these new roAp stars, we present an analysis of 44 known roAp stars observed by TESS during Cycle 1, providing the first high-precision and homogeneous sample of a significant fraction of the known roAp stars. The TESS observations have shown that almost 60 per cent (33) of our sample of stars are multiperiodic, providing excellent cases to test models of roAp pulsations, and from which the most rewarding asteroseismic results can be gleaned. We report four cases of the occurrence of rotationally split frequency multiplets that imply different mode geometries for the same degree modes in the same star. This provides a conundrum in applying the oblique pulsator model to the roAp stars. Finally, we report the discovery of non-linear mode interactions in $\alpha$ Cir (TIC 402546736, HD 128898) around the harmonic of the principal mode -- this is only the second case of such a phenomenon.

The origin of the diverse light-curve shapes of Type II supernovae (SNe), and whether they come from similar or distinct progenitors, has been actively discussed for decades. Here we report spectropolarimetry of two fast declining Type II (Type IIL) SNe: SN 2013ej and SN 2017ahn. SN 2013ej exhibited high continuum polarization from very soon after the explosion to the radioactive tail phase with time-variable polarization angles. The origin of this polarimetric behavior can be interpreted as the combination of two different aspherical structures, namely an aspherical interaction of the SN ejecta with circumstellar matter (CSM) and an inherently aspherical explosion. Aspherical explosions are a common feature of slowly declining Type II (Type IIP) SNe. By contrast, SN 2017ahn showed low polarization not only in the photospheric phase but also in the radioactive tail phase. This low polarization in the tail phase, which has never before been observed in other Type IIP/L SNe, suggests that the explosion of SN 2017ahn was nearly spherical. These observations imply that Type IIL SNe have, at least, two different origins: they result from stars that have different explosion properties and/or different mass-loss processes. This fact might indicate that 13ej-like Type IIL SNe originate from a similar progenitor to those of Type IIP SNe accompanied by an aspherical CSM interaction, while 17ahn-like Type IIL SNe come from a more massive progenitor with less hydrogen in its envelope.

Through a new spectral-fitting technique on a Crab giant-pulse dataset at 400-800 MHz with the 46m dish at the Algonquin Radio Observatory, a sub-population of 6 narrow-banded out of 1578 giant pulses was discovered. The narrow-banded giant pulses are detected in both the main-pulse and inter-pulse, thereby being more likely to be caused by an intrinsic emission mechanism as opposed to a propagation effect. Fast Radio Bursts (FRBs) have demonstrated similar narrow-banded behaviour while only little has been observed in the giant pulses of pulsars. The narrow-banded giant pulses reported here yield $\Delta \nu / \nu$ on the order of 0.1, which is close to the value of 0.05 reported for the repeater FRB20190711A. Hence, the connection between FRBs and giant pulses of pulsars is further established.

We revisit moduli stabilisation for type IIB flux compactifications that include a warped throat region corresponding to a warped deformed conifold, with an anti-D3-brane sitting at its tip. The warping induces a coupling between the conifold's deformation modulus and the bulk volume modulus in the K\"ahler potential. Previous works have studied the scalar potential assuming a strong warping such that this coupling term dominates, and found that the anti-D3-brane uplift may destabilise the conifold modulus and/or volume modulus, unless flux numbers within the throat are large, which makes tadpole cancellation a challenge. We explore the regime of parameter space corresponding to a weakly-but-still warped throat, such that the coupling between the conifold and volume moduli is subdominant. We thus discover a new metastable de Sitter solution within the four-dimensional effective field theory. We discuss the position of this de Sitter vacuum in the string theory landscape and swampland.

Ever since the discovery of neutrinos, one question has daunted us, are neutrinos their own antiparticles? One remarkable possibility is that neutrinos have a pseudo-Dirac nature, truly Majorana neutrinos which behave, for all practical purposes, as Dirac fermions, only distinguishable by tiny mass-squared differences. Such mass differences would induce oscillations that could only be conspicuous over astrophysical baselines. We analyze the neutrino data from SN1987A in the light of these active-sterile oscillations and find a mild preference ($\Delta\chi^2\approx 3$) for a non-zero quadratic mass difference $\delta m^2=6.31\times 10^{-20}~{\rm eV}^2$. Notably, the same data is able to exclude $\delta m^2\sim [2.55,3.01]\times 10^{-20}~{\rm eV}^2$ with $\Delta\chi^2> 9$, the tiniest mass differences constrained so far. We further consider the future sensitivity of next-generation experiments like the Deep Underground Neutrino Experiment (DUNE) and Hyper-Kamiokande (HK) and demonstrate that, for a future galactic SN occurring at $10~{\rm kpc}$, mass-squared differences as small as $\sim 10^{-20}~{\rm eV}^2$ could be explored.

We show that monopoles can be pair produced by cosmological magnetic fields in the early universe. The pair production gives rise to relic monopoles, and at the same time induces a self-screening of the magnetic fields. By studying these effects we derive limits on the monopole mass, and also on the initial amplitude of primordial magnetic fields. Monopoles of GUT scale mass can even be produced if primordial magnetic fields exist at sufficiently high redshifts.

At any epoch, particle physics must be open to completely unexpected discoveries, and that is reason enough to extend the reach of searches for ultra-high energy (UHE) photons. The observation of a population of photons with energies $E \gtrsim 100$ EeV would for example imply the existence of either a completely new physical phenomena, or particle acceleration mechanisms heretofore never seen or imagined. But as we outline in this Letter of Interest, there are also good arguments for super-heavy dark matter (SHDM) in a parameter range such that it could be discovered via its decays to, in particular, UHE photons. Only ultra-high energy cosmic ray observatories have capabilities to detect UHE photons. We first investigate how current and future observations can probe and constrain SHDM models in important directions, and then outline some of the scenarios that motivate such searches. We also discuss connections between constraints on SHDM and on the parameter values of cosmological models.

Anomalies in the cosmic microwave background (CMB) refer to features that have been observed, mostly at large angular scales, and which show some tension with the statistical predictions of the standard $\Lambda$CDM model. In this work, we focus our attention on power suppression, dipolar modulation, a preference for odd parity, and the tension in the lensing parameter $A_L$. Though the statistical significance of each individual anomaly is inconclusive, collectively they are significant, and could indicate new physics beyond the $\Lambda$CDM model. In this article, we present a brief, but pedagogical introduction to CMB anomalies and propose a common origin in the context of loop quantum cosmology.

The linearized dynamical equation for metric perturbations in a fully general, non-vacuum, background geometry is obtained from the Hamilton variational principle applied to the action up to second order. We specialize our results to the case of traceless and transverse metric fluctuations, and we discuss how the intrinsic properties of the matter stress tensor can affect (and modify) the process of gravity wave propagation even in most conventional geometric scenarios, like (for instance) those described by a FLRW metric background. We provide explicit examples for fluid, scalar field and electromagnetic field sources.

When the inflaton couples to photons and amplifies electric fields, charged particles produced via the Schwinger effect can dominate the universe after inflation, which is dubbed as the Schwinger preheating. Using the hydrodynamic approach for the Boltzmann equation, we numerically study two cases, the Starobinsky inflation model with the kinetic coupling and the anisotropic inflation model. The Schwinger preheating is not observed in the latter model but occurs for a sufficiently large inflaton-photon coupling in the first model. We analytically address its condition and derive a general attractor solution of the electric fields. The occurrence of the Schwinger preheating in the first model is determined by whether the electric fields enter the attractor solution during inflation or not.

MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC).

Gravitational waves (GWs) have opened a new window on fundamental physics in a number of important ways. The next generation of GW detectors may reveal more information about the polarization structure of GWs. Additionally, there is growing interest in theories of gravity beyond GR. One such theory which remains viable within the context of recent measurements of the speed of propagation of GWs is the teleparallel analogue of Horndeski gravity. In this work, we explore the polarization structure of this newly proposed formulation of Horndeski theory. In curvature-based gravity, Horndeski theory is almost synonymous with extensions to GR since it spans a large portion of these possible extensions. We perform this calculation by taking perturbations about a Minkowski background and consider which mode propagates. The result is that the polarization structure depends on the choice of model parameters in the teleparallel Horndeski Lagrangian with a maximum of seven propagating degrees of freedom. While the curvature-based Horndeski results follows as a particular limit within this setup, we find a much richer structure of both massive and massless cases which produce scalar--vector--tensor propagating degrees of freedom. We also find that the GW polarization that emerges from the teleparallel analogue of Horndeski gravity results in analogous massive and massless modes which take on at most four polarizations in the massless sector and two scalar ones in the massive sector. In none of the cases do we find vector polarizations.

During a solar eclipse the solar irradiance reaching the top-of-atmosphere (TOA) is reduced in the Moon shadow. The solar irradiance is commonly measured by Earth observation satellites before the start of the solar eclipse and is not corrected for this reduction, which results in a decrease of the computed TOA reflectances. Consequently, air quality products that are derived from TOA reflectance spectra, such as the ultraviolet (UV) Absorbing Aerosol Index (AAI), are distorted or undefined in the shadow of the Moon. The availability of air quality satellite data in the penumbral and antumbral shadow during solar eclipses, however, is of particular interest to users studying the atmospheric response to solar eclipses. Given the time and location of a point on the Earth's surface, we explain how to compute the obscuration during a solar eclipse taking into account wavelength-dependent solar limb darkening. With the calculated obscuration fractions, we restore the TOA reflectances and the AAI in the penumbral shadow during the annular solar eclipses on 26 December 2019 and 21 June 2020 measured by the TROPOMI/S5P instrument. In the corrected products, the signature of the Moon shadow disappeared, but only if wavelength-dependent solar limb darkening is taken into account. We conclude that the correction method of this paper can be used to detect real AAI rising phenomena during a solar eclipse and has the potential to restore any other product that is derived from TOA reflectance spectra. This would resolve the solar eclipse anomalies in satellite air quality measurements and would allow for studying the effect of the eclipse obscuration on the composition of the Earth's atmosphere from space.

We consider matter density effects in theories with a false ground state. Large and dense systems, such as stars, can destabilize a metastable minimum and allow for the formation of bubbles of the true minimum. We derive the conditions under which these bubbles form, as well as the conditions under which they either remain confined to the dense region or escape to infinity. The latter case leads to a phase transition in the universe at star formation. We explore the phenomenological consequences of such seeded phase transitions.