Papers for Friday, Sep 10 2021

The observed flat rotation curves of galaxies are among a number of astrophysical phenomena which require a larger acceleration than can be provided by the Newtonian gravity of the detected baryons. The main proposed solutions are additional undetected mass in the form of dark matter, or a low-acceleration modification to Newtonian gravity known as Milgromian dynamics (MOND). It is possible to directly test MOND using wide binary stars in the Solar neighbourhood, as these systems should contain a dynamically insignificant amount of dark matter even if it comprises most of the Galaxy. However, local wide binaries in MOND should orbit each other $\approx 20\%$ faster than in Newtonian dynamics. We describe the detailed plan for how this wide binary test will be conducted, focusing especially on stages with a high numerical cost. The computational costs and memory requirements are estimated for the main stages in the plan. Our overall assessment is that the critically important cost function can be evaluated deterministically at a marginal cost of a few seconds, giving the absolute binomial likelihood of a model. This will allow the cost function to be embedded within a Markov Chain Monte Carlo sampler, or a less expensive gradient descent stage designed to reveal the best-fitting model parameters. Therefore, the wide binary test should be feasible using currently available technology.

We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11}$ $M_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D-1D approach to fit 3D kinematic models to long-slit observations of the H$\alpha$-[NII] emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our sample seem to reach a flat part of the rotation curve within the outermost optical radius. We use the derived kinematics to study the high-mass end of the two most important scaling relations for spiral galaxies: the stellar/baryonic mass Tully-Fisher relation and the Fall (mass-angular momentum) relation. All galaxies in our sample, with the possible exception of the two fastest rotators, lie comfortably on both these scaling relations determined at lower masses, without any evident break or bend at the high-mass regime. When we combine our high-mass sample with lower-mass data from the Spitzer Photometry & Accurate Rotation Curves catalog, we find a slope of $\alpha=4.25\pm0.19$ for the stellar Tully-Fisher relation and a slope of $\gamma=0.64\pm0.11$ for the Fall relation. Our results indicate that most, if not all, of these rare, giant spiral galaxies are scaled up versions of less massive discs and that spiral galaxies are a self-similar population of objects up to the very high-mass end.

We present a detailed spectroscopic analysis of a galaxy at $z \simeq 4.88$ that is, by chance, magnified $\sim 30 \times$ by gravitational lensing. Only three sources at $z \gtrsim 5$ are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width CIV $\lambda$ 1549 {\AA} emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionisation. Using UV nebular line ratio diagnostics, VLT/X-shooter observations rule out strong AGN activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of [NeIII] $\lambda$ 3870 {\AA} and use the [NeIII]/[OII] line ratio to constrain the ionisation parameter and gas-phase metallicity. Closely related to the commonly used [OIII]/[OII] ratio, our [NeIII]/[OII] measurement shows this source is similar to local "Green Pea" galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the Fundamental Metallicity Relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of MgII $\lambda$ 2796 {\AA}, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting MgII absorption. Strong MgII emission has been observed in most $z \sim 0$ Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionising radiation. In conclusion, this strongly lensed galaxy, observed just 300 Myr after reionisation ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the $\mathit{JWST}$/ELT era.

The bipolar, nonthermal, high-latitude lobes known as the Fermi bubbles (FBs) are thought to originate from a massive energy release near the Galactic centre (GC). We constrain the FB engine and the circumgalactic medium (CGM) by analytically and numerically modeling the FB edges as strong forward shocks, as inferred from recent observations. A non-directed energy release produces shocks too spherical to account for observations even for a maximally massive Galactic disc, critical CGM rotation, or injection effectively offset from the GC. In contrast, collimated injection nearly perpendicular to the disc can account for observations in both ballistic (free expansion) and slowdown regimes, as we show using a simple stratified evolution model verified by hydrodynamic simulations. FBs still in their ballistic regime require injection (at $z\simeq100$ pc heights in our model) with a half-opening angle $\theta\simeq4^\circ$, a normalized velocity $\beta_{-2}\equiv v/(0.01c)\gtrsim 0.4$, and an energy $E\gtrsim2\beta_{-2}^2\times 10^{55}$ erg, launched $\mathbb{T}\simeq 3.3\beta_{-2}^{-1}$ Myr ago, showing a distinctive low-pressure region behind the bubble head. Slowing-down (mass accumulated) FBs require a faster injection, a thinner jet, a smaller $E/(\beta_{-2}\theta)^{2}$, and a comparable $\mathbb{T}$, and follow a ballistic stage that must reach a height $z_{s}\gtrsim 5$ kpc.

Chemical tagging is a powerful tool to reveal the origin of stars and globular clusters (GCs), especially when dynamics alone cannot provide robust answers. So far, mostly $\alpha$- and neutron capture elements have been used to distinguish stars born in the Milky Way (MW) from those born in external environments such as that of dwarf galaxies. Here, instead, we use iron-peak elements abundances to investigate the origin of a sample of metal-rich globular clusters. By homogeneously analyzing high-resolution UVES spectra of giant stars belonging to four metal-rich GCs (namely NGC 5927, NGC 6388, NGC 6441, NGC 6496), we find that while the $\alpha$-elements Si and Ca have similar abundance ratios for all the four GCs, and Ti and neutron capture elements (La, Ba and Eu) only show a marginal discrepancy, a stark difference is found when considering the abundances of some iron-peak elements (Sc, V and Zn). In particular, NGC 6388 and NGC 6441 have abundance ratios for these iron-peak elements significantly lower (by ~ 0.5 dex) than those measured in NGC 5927 and NGC 6496, which are clearly identified as born in-situ MW clusters through an analysis of their orbital properties. These measurements indicate that the environment in which these clusters formed is different, and provide robust evidence supporting an accreted origin from the same progenitor for NGC 6388 and NGC 6441.

In the framework of a systematic ALMA study of IR-selected main-sequence and starburst galaxies at z~1-1.7 at typical ~1" resolution, we report on the effects of mid-IR- and X-ray-detected active galactic nuclei (AGN) on the reservoirs and excitation of molecular gas in a sample of 55 objects. We find detectable nuclear activity in ~30% of the sample. The presence of dusty tori influences the IR SED of galaxies, as highlighted by the strong correlation among the AGN contribution to the total IR luminosity budget (fAGN = LIR,AGN/LIR), its hard X-ray emission, and the Rayleigh-Jeans to mid-IR (S1.2mm/S24um) observed color, with consequences on the empirical SFR estimates. Nevertheless, we find only marginal effects of AGN on the CO (J=2,4,5,7) or neutral carbon ([CI](1-0), [CI](2-1)) line luminosities and on the derived molecular gas excitation as gauged by line ratios and the full SLEDs. The [CI] and CO emission up to J=5,7 thus primarily traces the properties of the host in typical IR luminous galaxies. However, we highlight the existence of a large variety of line luminosities and ratios despite the homogeneous selection. In particular, we find a sparse group of AGN-dominated sources with the highest LIR,AGN/LIR,SFR ratios, >3, that are more luminous in CO(5-4) than what is predicted by the L'CO(5-4)-LIR,SFR relation, which might be the result of the nuclear activity. For the general population, our findings translate into AGN having minimal effects on quantities such as gas and dust fractions and SFEs. If anything, we find hints of a marginal tendency of AGN hosts to be compact at far-IR wavelengths and to display 1.8x larger dust optical depths. In general, this is consistent with a marginal impact of the nuclear activity on the gas reservoirs and star formation in average star-forming AGN hosts with LIR>5e11 Lsun, typically underrepresented in surveys of quasars and SMGs.

We perform an aperture-matched analysis of dust-corrected H$\alpha$ and UV SFRs using 303 star-forming galaxies with spectroscopic redshifts $1.36<z_\text{spec}<2.66$ from the MOSFIRE Deep Evolution Field (MOSDEF) survey. By combining H$\alpha$ and H$\beta$ emission line measurements with multi-waveband resolved CANDELS/3D-HST imaging, we directly compare dust-corrected H$\alpha$ and UV SFRs, inferred assuming a fixed attenuation curve shape and constant SFHs, within the spectroscopic aperture. Previous studies have found that H$\alpha$ and UV SFRs inferred with these assumptions generally agree for typical star-forming galaxies, but become increasingly discrepant for galaxies with higher SFRs ($\gtrsim$100 M$_\odot$ yr$^{-1}$), with H$\alpha$-to-UV SFR ratios being larger for these galaxies. Our analysis shows that this trend persists even after carefully accounting for the apertures over which H$\alpha$ and UV-based SFRs (and the nebular and stellar continuum reddening) are derived. Furthermore, our results imply that H$\alpha$ SFRs may be higher in the centers of large galaxies (i.e., where there is coverage by the spectroscopic aperture) compared to their outskirts, which could be indicative of inside-out galaxy growth. Overall, we suggest that the persistent difference between nebular and stellar continuum reddening and high H$\alpha$-to-UV SFR ratios at the centers of large galaxies may be indicative of a patchier distribution of dust in galaxies with high SFRs.

The damped random walk (DRW) process is one of the most commonly used and simplest stochastic models to describe variability of active galactic nuclei (AGN). An AGN light curve can be converted to just two DRW model parameters - the signal decorrelation timescale $\tau$ and the asymptotic amplitude ${SF}_{\infty}$. By simulation means, we have recently shown that in order to measure the decorrelation timescale accurately, the experiment or the light curve length must be at least 10 times the underlying decorrelation timescale. In this paper, we investigate the origin of this requirement and find that typical AGN light curves do not sufficiently represent the intrinsic stationary process. We simulated extremely long (10,000$\tau$) AGN light curves using DRW, and then measured the variance and the mean of short light curves spanning 1-1000$\tau$. We modeled these light curves with DRW to obtain both the signal decorrelation timescale $\tau$ and the asymptotic amplitude ${SF}_{\infty}$. The variance in light curves shorter than $\approx30\tau$ is smaller than that of the input process, as estimated by both a simple calculation from the light curve and by DRW modeling. This means that while the simulated stochastic process is intrinsically stationary, short light curves do not adequately represent the stationary process. Since the variance and timescale are correlated, underestimated variances in short light curves lead to underestimated timescales as compared to the input process. It seems, that a simulated AGN light curve does not fully represent the underlying DRW process until its length reaches even $\approx30$ decorrelation timescales. Modeling short AGN light curves with DRW leads to biases in measured parameters of the model - the amplitude being too small and the timescale being too short.

We present a co-scaling grid formalism and its implementation in the magnetohydrodynamics code Athena++. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The grid is evolved at the same time order as the fluid variables. The user specifies grid evolution laws, which can be independent of the fluid motion. Applying our implementation to standard hydrodynamic test cases leads to improved results and higher efficiency, compared to the fixed-grid solutions.

We report long-termed UBVRIRcIc photometry of the highly eccentric 12.04 day detached eclipsing binary V490 Sct (V =13.1, B9.5+A0, e = 0.40), which we use to determine its relative and absolute parameters. The absolute masses, radii, and temperatures are Ma = 2.33+/-0.1 Msun, Ra = 1.91+/-0.04 Rsun, and Ta = 9960+/-60 K for the primary and Mb = 2.24+/-0.1 Msun, Rb = 1.86+/-0.04 Rsun, and Tb = 9700+/-80 K for the secondary. The system displays a slow periastron advance that is dominated by general relativity (GR). Our measurement, dw/dt = 0.86 deg/century, is 32% less then the expected rate, dw/dt = 1.24 deg/century, which has an 83% contribution from GR. A comparison with current stellar evolution models shows a good match to the measured properties at an age of about 130 mln. years and Solar abundance. The photometrical parallax of the system pi = 0.77 +/- 0.02 mas, matches quite well the GAIA DR2 value, pi = 0.76 +/- 0.04 mas.

A rotating star with a monopole (or split monopole) magnetic field gives the simplest, prototype model of a rotationally driven stellar wind. Winds from compact objects, in particular neutron stars, carry strong magnetic fields with modest plasma loading, and develop ultra-relativistic speeds. We investigate the relativistic wind launched from a dense, gravitationally bound, atmosphere on the stellar surface. We first examine the problem analytically and then perform global kinetic plasma simulations. Our results show how the wind acceleration mechanism changes from centrifugal (magnetohydrodynamic) to electrostatic (charge-separated) depending on the parameters of the problem. The two regimes give winds with different angular distributions and different scalings with the magnetization parameter.

We present a global kinetic plasma simulation of an axisymmetric pulsar magnetosphere with self-consistent $e^\pm$ pair production. We use the particle-in-cell method and log-spherical coordinates with a grid size $4096\times 4096$. This allows us to achieve a high voltage induced by the pulsar rotation and investigate pair creation in a young pulsar far from the death line. We find the following. (1) The energy release and $e^\pm$ creation are strongly concentrated in the thin, Y-shaped current sheet, with a peak localized in a small volume at the Y-point. (2) The Y-point is shifted inward from the light cylinder by $\sim 15\%$, and "breathes" with a small amplitude. (3) The dense $e^\pm$ cloud at the Y-point is in ultra-relativistic rotation, which we call super-rotation, because it exceeds co-rotation with the star. The cloud receives angular momentum flowing from the star along the poloidal magnetic lines. (4) Gamma-ray emission peaks at the Y-point and is collimated in the azimuthal direction, tangent to the Y-point circle. (5) The separatrix current sheet between the closed magnetosphere and the open magnetic field lines is sustained by the electron backflow from the Y-point cloud. Its thickness is self-regulated to marginal charge starvation. (6) Only a small fraction of dissipation occurs in the separatrix inward of the Y-point. A much higher power is released in the equatorial plane, especially at the Y-point where the created dense $e^\pm$ plasma is spun up and intermittently ejected through the nozzle between the two open magnetic fluxes.

The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($\sigma(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning $\sim$1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($\sigma(\mathrm{[Ba~II/H]}) = {0.03}^{+0.10}_{-0.02}$). This abundance signature is best interpreted as slow neutron capture element production from a massive fast-rotating metal-poor star ($15-20 \mathrm{M}_\odot$, $v_\mathrm{ini}/v_\mathrm{crit} = 0.4$, $[\mathrm{Fe/H}] = -3.8$). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A(\mathrm{Li}) = 3.1 \pm 0.1$). At $[\mathrm{Fe/H}] = -2.93$ it is among the most metal-poor lithium-rich giants known.

We employ ionization-parameter mapping (IPM) to infer the optical depth of HII regions in the northern half of M33. We construct [OIII]$\lambda 5007$/[OII]$\lambda 3727$ and [OIII]$\lambda 5007$/[SII]$\lambda 6724$ ratio maps from narrow-band images continuum-subtracted in this way, from which we classify the HII regions by optical depth to ionizing radiation, based on their ionization structure. This method works relatively well in the low metallicity regime, $12 + \log(\rm O/H) \leq 8.4$, where [OIII]$\lambda\lambda4949,5007$ is strong. However, at higher metallicities, the method breaks down due to the strong dependence of the [OIII]$\lambda\lambda4959,5007$ emission lines on the nebular temperature. Thus, although O$^{++}$ may be present in metal-rich HII regions, these commonly used emission lines do not serve as a useful indicator of its presence, and hence, the O ionization state. In addition, IPM as a diagnostic of optical depth is limited by spatial resolution. We also report a region of highly excited [OIII] extending over an area $\sim$ 1 kpc across and [OIII]$\lambda5007$ luminosity of $4.9\pm 1.5\times10^{38}$ erg/s, which is several times higher than the ionizing budget of any potential sources in this portion of the galaxy. Finally, this work introduces a new method for continuum subtraction of narrow-band images based on the dispersion of pixels around the mode of the diffuse-light flux distribution. In addition to M33, we demonstrate the method on C III]$\lambda$1909 imaging of Haro~11, ESO 338-IG004, and Mrk~71.

ASASSN-V J205543.90+240033.5 has been suggested to be a white dwarf pulsar by Kato (2021, arXiv:2108.09060). We obtained time-resolved photometry and identified the orbital and spin periods to be 0.523490(1) d and 0.00678591(1) d = 9.77 min, respectively. These values strengthen the similarity of this object with AR Sco. We estimated that the strength of the spin pulse is 3.6 times smaller than in AR Sco.

Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respectively. In this work, we use high-resolution (sub-pc) ALMA $^{12}$CO ($J$ = $2-1$), $^{13}$CO ($J$ = $2-1$), and C$^{18}$O ($J$ = $2-1$) observations of filamentary molecular clouds in the star-forming region N159E (the Papillon Nebula) in the Large Magellanic Cloud (LMC) to provide the first measurement of turbulence driving parameter in an extragalactic region. We use a non-local thermodynamic equilibrium (NLTE) analysis of the CO isotopologues to construct a gas density PDF, which we find to be largely log-normal in shape with some intermittent features indicating deviations from lognormality. We find that the width of the log-normal part of the density PDF is comparable to the supersonic turbulent Mach number, resulting in $b \approx 0.9$. This implies that the driving mode of turbulence in N159E is primarily compressive. We speculate that the compressive turbulence could have been powered by gravo-turbulent fragmentation of the molecular gas, or due to compression powered by \ion{H}{i} flows that led to the development of the molecular filaments observed by ALMA in the region. Our analysis can be easily applied to study the nature of turbulence driving in resolved star-forming regions in the local as well as the high-redshift Universe.

A major discovery of Parker Solar Probe (PSP) was the presence of large numbers of localized increases in the radial solar wind speed and associated sharp deflections of the magnetic field - switchbacks (SB). A possible generation mechanism of SBs is through magnetic reconnection between open and closed magnetic flux near the solar surface, termed interchange reconnection that leads to the ejection of flux ropes (FR) into the solar wind. Observations also suggest that SBs undergo merging, consistent with a FR picture of these structures. The role of FRs merging in controlling the structure of SB in the solar wind is explored through direct observations, through analytic analysis, and numerical simulations. Analytic analysis reveals key features of the structure of FR and their scaling with the heliocentric distance R that are consistent with observations and that reveal the critical role of merging in controlling the SB structure. FR merging is shown to be energetically favorable to reduce the strength of the wrapping magnetic field and drive the observed elongation of SBs. A further consequence is the resulting dominance of the axial magnetic field within SBs that leads to the characteristic sharp rotation of the magnetic field into the axial direction at the SB boundary that is revealed in observations. Finally, the radial scaling of the SB area in the FR model of SBs suggests that the observational probability of SB identification should be insensitive to R, which is consistent with the most recent statistical analysis of SB observations from PSP.

Basing on the large volume \textit{Gaia} Early Data Release 3 and LAMOST Data Release 5 data, we estimate the bias-corrected binary fractions of the field late G and early K dwarfs. A stellar locus outlier method is used in this work, which works well for binaries of various periods and inclination angles with single epoch data. With a well-selected, distance-limited sample of about 90 thousand GK dwarfs covering wide stellar chemical abundances, it enables us to explore the binary fraction variations with different stellar populations. The average binary fraction is 0.42$\pm$0.01 for the whole sample. Thin disk stars are found to have a binary fraction of 0.39$\pm$0.02, thick disk stars own a higher one of 0.49$\pm$0.02, while inner halo stars possibly own the highest binary fraction. For both the thin and thick disk stars, the binary fractions decrease toward higher [Fe/H], [$\alpha$/H], and [M/H] abundances. However, the suppressing impacts of the [Fe/H], [$\alpha$/H], and [M/H] are more significant for the thin disk stars than those for the thick disk stars. For a given [Fe/H], a positive correlation between [$\alpha$/Fe] and the binary fraction is found for the thin disk stars. However, this tendency disappears for the thick disk stars. We suspect that it is likely related to the different formation histories of the thin and thick disks. Our results provide new clues for theoretical works on binary formation.

Supernova remnants (SNRs) are important objects in terms of their connections with supernova (SN) explosion mechanism(s), progenitor stars, and cosmic-ray acceleration. Non-thermal emission from SNRs is an effective probe of the structure of their surrounding circumstellar media (CSM), which can in turn shed lights on mechanism and history of the elusive mass-loss of massive stars. In this work, we calculate the time evolution of broadband non-thermal emission from SNRs originating from Type II SNe embedded in a CSM environment linked to the mass loss history of the progenitor. Our results predict that Type II SNRs experience a prolonged period of weak radio and $\gamma$-ray emission if they run into a spatially extended bubble of low density and high temperature created by the stellar wind during main sequence. For a typical red supergiant progenitor evolved within an average interstellar medium (ISM), this "dark age" corresponds to a range of SNR ages spanning from ~1000 to 5000 yrs old. This result suggests that a majority of Type II SNRs are too faint to be detected, which may help explain why the number of known Galactic SNRs is significantly less than what we expect from the SN rate in our Galaxy.

'Main-sequence radio pulse-emitters' (MRPs) are magnetic early-type stars from which periodic radio pulses, produced via electron cyclotron maser emission (ECME), are observed. Despite the fact that these stars can naturally offer suitable conditions for triggering ECME, only seven such stars have been reported so far within a span of more than two decades. In this paper, we report the discovery of eight more MRPs, thus more than doubling the sample size of such objects. These discoveries are the result of our sub-GHz observation program using the Giant Metrewave Radio Telescope over the years 2015-2021. Adding these stars to the previously known MRPs, we infer that at least 32 percent of the magnetic hot stars exhibit this phenomenon, thus suggesting that observation of ECME is not a rare phenomenon. The significantly larger sample of MRPs allows us for the first time to perform a statistical analysis comparing their physical properties. We present an empirical relation that can be used to predict whether a magnetic hot star is likely to produce ECME. Our preliminary analysis suggests that the physical parameters that play the primary role in the efficiency of the phenomenon are the maximum surface magnetic field strength and the surface temperature. In addition, we present strong evidence of the influence of the plasma density distribution on ECME pulse profiles. Results of this kind further motivate the search for MRPs as a robust characterization of the relation between observed ECME properties and stellar physical parameters can only be achieved with a large sample.

Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$\alpha$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. We thus discuss chemodynamic selections (such as eccentricity and upper limits on [Na/Fe]).

The third observing run of the LIGO-Virgo collaboration has resulted in about hundred gravitational-wave triggers including the binary neutron star merger GW190425. However, none of these events have been accompanied with an electromagnetic transient found during extensive follow-up searches. In this article, we perform new numerical-relativity simulations of binary neutron star and black hole - neutron star systems that have a chirp mass consistent with GW190425. Assuming that the GW190425's sky location was covered with sufficient accuracy during the electromagnetic follow-up searches, we investigate whether the non-detection of the kilonova is compatible with the source parameters estimated through the gravitational-wave analysis and how one can use this information to place constraints on the properties of the system. Our simulations suggest that GW190425 is incompatible with an unequal mass binary neutron star merger with a mass ratio $q<0.8$ when considering stiff or moderately stiff equations of state if the binary was face-on and covered by the observation. Our analysis shows that a detailed observational result for kilonovae will be useful to constrain the mass ratio of binary neutron stars in future events.

Thermally emitting neutron stars represent a promising environment for probing the properties of axion-like particles. Due to the strong magnetic fields of these sources, surface photons may partially convert into such particles in the large magnetospheric region surrounding the stars, which will result in distinctive signatures in their spectra. However, the interaction depends on the polarization state of the radiation and is rather weak due to the low experimentally allowed values of the coupling constant $g_{\gamma a}$. In this work, we compute the degree of photon-axion transition in the case of 100% O-mode polarization and spectral energy distribution of an isotropic blackbody with uniform surface temperature. The stellar magnetic field is assumed to be dipolar. We show that with the maximum effect reached for the magnetic fields $\sim10^{13}$ - $10^{14}$ G (typical for X-ray dim isolated neutron stars) and $g_{\gamma a} = 2 \times 10^{-11}$ GeV$^{-1}$, the optical flux is reduced by 30 - 40%, while the high-energy part of the spectrum is not affected. The low-energy decrease exceeds 5% at $g_{\gamma a} \geq 2 \times 10^{-12}$ GeV$^{-1}$ and $m_a \leq 2\times 10^{-6}$ eV, which is below the present experimental and astrophysical limits on axion parameters. To obtain the actual observational constraints, rigorous treatment of the radiative surface layers is required.

The detection of the high-energy neutrino IC-170922A in spatial (within the error region) and temporal flare activity correlation with the blazar TXS 0506+056 allowed these objects to be considered as progenitor sources of neutrinos. Besides this, no more detection of this kind was reported. Some other neutrinos detected by IceCube show a spatial correlation (within the error region) from other Fermi-LAT detected sources. However, these objects did not show a flare activity like TXS 0506+056. Assuming a lepto-hadronic scenario through p$\gamma$ interactions, this work describes the SED in some objects from the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi Gamma-ray Space Telescope Large Area Telescope (4LAC) sources, which are in spatial correlation with neutrinos detected by IceCube. Additionally, we estimate the corresponding neutrino flux counterpart from these sources.

We report a $5\sigma$ detection of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect, combining galaxy clusters in DESI imaging surveys and the Planck temperature maps. The detection is facilitated by both improvements in the data and in the analysis method. For the data, we adopt the recently released galaxy group catalog (Y20: \cite{yang2020extended}) with $\sim 10^6$ robustly-identified groups, and construct various galaxy cluster samples for the kSZ measurement. The Y20 catalogue also provides estimation of halo mass, which further improves the kSZ measurement by $\sim 10\%$. For the analysis method, we derive an optimal estimator of pairwise kSZ through the maximum likelihood analysis. It also handles potential systematic errors self-consistently. The baseline cluster sample, containing the $1.2\times 10^5$ richest galaxy clusters of typical mass ~$ 10^{14} M_{\odot}/h$ at typical redshift $0.2$-$0.5$, rules out the null hypothesis at $5\sigma$. When fitting with a pairwise kSZ template from simulations, the signal is detected at $4.7\sigma$ and the average optical depth is constrained as $\bar{\tau}_e=(1.66\pm 0.35)\times 10^{-4}$. We perform various internal checks, with different cluster selection criteria, different sky coverage and redshift range, different CMB maps, different filter sizes, different treatments of potential systematics and the covariance matrix. The kSZ effect is consistently detected with $2.5\leq $S/N$\leq 5.6$ and acceptable $\chi^2_{\rm min}$, across a variety of cluster samples. The S/N is limited by both the Planck resolution and the photo-z accuracy, and therefore can be significant improved with DESI spectroscopic redshift information and with other CMB experiments.

(Shortened) Luminous blue variable stars (LBVs) form dust as a result of episodic, violent mass loss. To investigate their contribution as dust producers in the Magellanic Clouds, we analyse 31 LBVs from a recent census. We built a maximally complete multi-wavelength dataset of these sources from archival data from near-IR to millimetre wavelengths. We review the LBV classification on the basis of the IR SED. To derive characteristic dust parameters, we fitted the photometry resulting from a stacking analysis. For comparison we also stacked the images of low- and intermediate-mass evolved stars in the LMC. We find four classes of sources: 1) LBVs showing mid-IR dust emission plus near-IR free-free emission from an ionised stellar wind (Class 1a) or only mid-IR dust emission (Class 1b); 2) LBVs with a near-IR excess due to free-free emission only (Class 2); 3) objects with an sgB[e] classification; and 4) objects with no detected stellar winds and no circumstellar matter in their SEDs. From the stacking analysis of the 18 Class 1 and 2 objects in the LMC, we derived an integrated dust mass of $0.11^{+0.06}_{-0.03} M_\odot$. This is two orders of magnitude larger than the value inferred from stacking 1342 extreme-AGB stars. The dust mass of individual LBVs does not correlate with the stellar parameters, possibly suggesting that the dust production mechanism is independent of the initial stellar mass or that the stars have different evolutionary histories. The total dust yield from LBVs over the age of the LMC is $\sim 10^4-10^5 M_\odot$. LBVs are potentially the second most important source of dust in normal galaxies. The role of dust destruction in LBV nebulae by a possible subsequent SN blast wave has yet to be determined. Recent theoretical developments in the field of dust processing by SN shocks highlight the potential survival of dust from the pre-existing circumstellar nebula.

A wide range of data formats and proprietary software have traditionally been used in gamma-ray astronomy, usually developed for a single specific mission or experiment. However, in recent years there has been an increasing effort towards making astronomical data open and easily accessible. Within the gamma-ray community this has translated to the creation of a common data format across different gamma-ray observatories: the "gamma-astro-data-format" (GADF). Based on a similar premise, open-source analysis packages, such as Gammapy, are being developed and aim to provide a single, robust tool which suits the needs of many experiments at once. In this contribution we show that data from the High-Altitude Water Cherenkov (HAWC) observatory can be made compatible with the GADF and present the first GADF-based production of event lists and instrument response functions for a ground-based wide-field instrument. We use these data products to reproduce with excellent agreement the published HAWC Crab spectrum using Gammapy. Having a common data format and analysis tools facilitates joint analysis between different experiments and effective data sharing. This will be especially important for next-generation instruments, such as the proposed Southern Wide-field Gamma-ray Observatory (SWGO) and the planned Cherenkov Telescope Array (CTA).

The presence of muons in air-showers initiated by cosmic ray protons and nuclei is well established as a powerful tool to separate such showers from those initiated by gamma-rays. However, so far this approach has been fully exploited only for ground level particle detecting arrays. In this contribution, we explore the feasibility of using Cherenkov light from muons as a background rejection tool for imaging atmospheric Cherenkov telescope arrays at the highest energies. We adopt an analytical model of the Cherenkov light from individual muons to allow rapid simulation of a large number of showers in a hybrid mode. This allows exploration of the very high background rejection power regime at acceptable cost in terms of computing time. We find that for very large telescopes ($\gtrsim$20 m diameter), efficient identification of muons would provide a major improvement with respect to standard background rejection techniques at energies above several tens of TeVs.

The energy spectra of primary and secondary cosmic rays (CR) generally harden at several hundreds of GeV, which can be naturally interpreted by propagation effects. We adopt a spatially dependent CR propagation model to fit the spectral hardening, where a slow-diffusion disk (SDD) is assumed near the Galactic plane. We aim to constrain the propagation parameters with the Bayesian parameter estimation based on a Markov chain Monte Carlo sampling algorithm. The latest precise measurements of carbon spectrum and B/C ratio are adopted in the Bayesian analysis. The $\rm{^{10}Be/^{9}Be}$ and Be/B ratios are also included to break parameter degeneracies. The fitting result shows that all the parameters are well constrained. Especially, the thickness of the SDD is limited to 0.4-0.5 kpc above and below the Galactic plane, which could be the best constraint for the slow-diffusion region among similar works. The $\bar{p}/p$ ratio and amplitude of CR anisotropy predicted by the SDD model are consistent with the observations, while the predicted high-energy electron and positron fluxes are slightly and significantly lower than the observations, respectively, indicating the necessity of extra sources.

Young solar-type stars are known to be strong X-ray emitters and their X-ray spectra have been widely studied. X-rays from the central star may play a crucial role in the thermodynamics and chemistry of the circumstellar material as well as in the atmospheric evolution of young planets. In this paper we present model spectra based on spectral parameters derived from the observations of young stars in the Orion Nebula Cluster from the Chandra Orion Ultradeep Project (COUP). The spectra are then used to calculate new photoevaporation prescriptions that can be used in disc and planet population synthesis models. Our models clearly show that disc wind mass loss rates are controlled by the stellar luminosity in the soft (100 eV - 1 keV) X- ray band. New analytical relations are provided for the mass loss rates and profiles of photoevaporative winds as a function of the luminosity in the soft X-ray band. The agreement between observed and predicted transition disc statistics moderately improved using the new spectra, but the observed population of strongly accreting large cavity discs can still not be reproduced by these models. Furthermore, our models predict a population of non-accreting transition discs that are not observed. This highlights the importance of considering the depletion of millimeter-sized dust grains from the outer disc, which is a likely reason why such discs have not been detected yet.

We present a polarimetric study of the pulsar wind nebula (PWN) in supernova remnant G21.5$-$0.9 using archival Very Large Array (VLA) data. The rotation measure (RM) map of the PWN shows a symmetric pattern that aligns with the presumed pulsar spin axis direction, implying a significant contribution of RM from the nebula. We suggest that the spatial variation of the internal RM is mostly caused by non-uniform distribution of electrons originated from the supernova ejecta. Our high-resolution radio polarization map reveals an overall radial $B$-field. We construct a simple model with an overall radial $B$-field and turbulence in small scale. The model can reproduce many of the observed features of the PWN, including the polarization pattern and polarized fraction. The results also reject a large-scale toroidal $B$-field which implies that the toroidal field observed in the inner PWN cannot propagate to the entire nebula.

Understanding and modelling astrophysical plasmas on atomic levels while taking into account various assumptions (for example, collisional ionisation equilibrium or photoionisation equilibrium) became essential with the progress of high-resolution X-ray spectroscopy. In order to prepare for the upcoming X-ray spectroscopy missions such as XRISM or Athena, the plasma codes with their models and the atomic databases need to be up to date and accurate. One such update for the plasma code SPEX is presented in this paper where we focus on the radiative loss due to collisional excitation in the low-density, optically thin regime. We also update the atomic data for neutral hydrogen and include the contribution of the dielectronic recombination. With all these updates being implemented in SPEX we finally present the new cooling curve. We include the comparison to other plasma codes (MEKAL, APEC, Cloudy) and other atomic databases (CHIANTI, ADAS). We show how the updated cooling impacts the stability curve for photoionised plasmas and find a new stable branch.

We studied the low-frequency $\lesssim 0.5\;$h$^{-1}$ (long-period $\gtrsim 2\;$h) oscillations of active regions (ARs). The investigation is based on an analysis of a time series built from Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) photospheric magnetograms and comprises case studies of several types of AR structures. The main goals are to investigate whether ARs can be engaged in long-period oscillations as unified oscillatory entities and, if so, to determine the spectral pattern of such oscillations. Time series of characteristic parameters of the ARs, such as, the total area, total unsigned radial magnetic flux, and tilt angle, were measured and recorded using the image moment method. The power spectra were built out of Gaussian-apodised and zero-padded datasets. There are long-period oscillations ranging from 2 to 20 h, similarly to the characteristic lifetimes of super-granulation, determined from the datasets of the AR total area and radial magnetic flux, respectively. However, no periodicity in tilt angle data was found. Whatever nature these oscillations have, they must be energetically supported by convective motions beneath the solar surface. The possible interpretations can be related to different types of magnetohydrodynamic (MHD) oscillations of the multi-scale structure of the AR magnetic field, which is probably linked with the characteristic turnover timescales of the super-granulation cells. The presence of oscillations in the radial magnetic flux data may be connected to periodic flux emergence or cancellation processes.

We demonstrate that future radius measurement of the NICER mission have the potential to reveal the existence of a strong phase transition in dense neutron star matter by confirming the existence of so called twin stars, compact star configurations with the same mass but different radii. The latest radius constraints from NICER for the pulsars J0740+6620 as well as J0030+0451 are discussed using relativistic mean field equations of state with varying stiffness, connected with a first order phase transition to quark matter. We show, that twin star solutions are compatible with the new radius constraint but are located at radii \emph{below} the present constraints from NICER serving as a smoking gun for a strong phase transition in neutron star matter. This scenario is realized if a strong phase transition takes place in neutron stars of the first branch with masses above $2\,M_\odot$.

The ungrouped iron meteorite Nedagolla is the first meteorite with bulk Mo, Ru, and Ni isotopic compositions that are intermediate between those of the non-carbonaceous (NC) and carbonaceous (CC) meteorite reservoirs. The Hf-W chronology of Nedagolla indicates that this mixed NC-CC isotopic composition was established relatively late, more than 7 million years after Solar System formation. The mixed NC-CC isotopic composition is consistent with the chemical composition of Nedagolla, which combines signatures of metal segregation under more oxidizing conditions (relative depletions in Mo and W), characteristic for CC bodies, and more reducing conditions (high Si and Cr contents), characteristic for some NC bodies, in a single sample. These data combined suggest that Nedagolla formed as the result of collisional mixing of NC and CC core material, which partially re-equilibrated with silicate mantle material that predominantly derives from the NC body. These mixing processes might have occurred during a hit-and-run collision between two differentiated bodies, which also provides a possible pathway for Nedagolla's extreme volatile element depletion. As such, Nedagolla provides the first isotopic evidence for early collisional mixing of NC and CC bodies that is expected as a result of Jupiter's growth.

The Mid-Infrared Instrument (MIRI) on-board JWST will provide imaging, coronagraphy, low-resolution spectroscopy and medium-resolution spectroscopy at unprecedented sensitivity levels in the mid-infrared wavelength range. The Medium-Resolution Spectrometer (MRS) of MIRI is an integral field spectrograph that provides diffraction-limited spectroscopy between 4.9 and 28.3 um, within a FOV varying from 13 to 56" square. From ground testing, we calculate the physical parameters essential to general observers and calibrating the wavelength solution and resolving power of the MRS is critical for maximising the scientific performance of the instrument. We have used ground-based observations of discrete spectral features in combination with Fabry-Perot etalon spectra to characterize the wavelength solution and spectral resolving power of the MRS. We present the methodology used to derive the MRS spectral characterisation, which includes the precise wavelength coverage of each MRS sub-band, computation of the resolving power as a function of wavelength, and measuring slice-dependent spectral distortions. The resolving power varies from R3500 in channel 1 to R1500 in channel 4. Based on the ground test data, the wavelength calibration accuracy is estimated to be below one tenth of a pixel, with small systematic shifts due to the target position within a slice for unresolved sources, that have a maximum amplitude of about 0.25 spectral resolution elements. Based on ground test data, the MRS complies with the spectral requirements for both the R and wavelength accuracy for which it was designed. We also present the commissioning strategies and targets that will be followed to update the spectral characterisation of the MRS.

We present a novel 2D flux density model for observed HI emission lines combined with a Bayesian stacking technique to measure the baryonic Tully-Fisher relation below the nominal detection threshold. We simulate a galaxy catalogue, which includes HI lines described either with Gaussian or busy function profiles, and HI data cubes with a range of noise and survey areas similar to the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. With prior knowledge of redshifts, stellar masses and inclinations of spiral galaxies, we find that our model can reconstruct the input baryonic Tully-Fisher parameters (slope and zero point) most accurately in a relatively broad redshift range from the local Universe to $z = 0.3$ for all the considered levels of noise and survey areas, and up to $z = 0.55$ for a nominal noise of $90\,\mu$Jy/channel over 5 deg$^{2}$. Our model can also determine the $M_{\rm HI} - M_{\star}$ relation for spiral galaxies beyond the local Universe, and account for the detailed shape of the HI emission line, which is crucial for understanding the dynamics of spiral galaxies. Thus, we have developed a Bayesian stacking technique for measuring the baryonic Tully-Fisher relation for galaxies at low stellar and/or HI masses and/or those at high redshift, where the direct detection of HI requires prohibitive exposure times.

We compute the desorption rate of icy mantles on dust grains as a function of the size and composition of both the grain and the mantle. We combine existing models of cosmic ray (CR) related desorption phenomena with a model of CR transport to accurately calculate the desorption rates in dark regions of molecular clouds. We show that different desorption mechanisms dominate for grains of different sizes, and in different regions of the cloud. We then use these calculations to investigate a simple model of the growth of mantles, given a distribution of grain sizes. We find that modest variations of the desorption rate with grain size lead to a strong dependence of mantle thickness on grain size. Furthermore, we show that freeze-out is almost complete in the absence of an external UV field, even when photodesorption from CR produced UV is taken into consideration. Even at gas densities of $10^4$ ${\rm cm^{-3}}$, less than 30\% of the CO remains in the gas phase after $3\times 10^5$ years for standard values of the CR ionization rate.

Coronal Holes (CHs) have subdued intensity and net blueshifts when compared to Quiet Sun (QS) at coronal temperatures. At transition region temperatures, such differences are obtained for regions with identical photospheric absolute magnetic flux density ($\vert$B$\vert$). In this work, we use spectroscopic measurements of the \car 1334~{\AA} line from Interface Region Imaging Spectrograph (IRIS), formed at chromospheric temperatures, to investigate the intensity, Doppler shift, line width, skew, and excess kurtosis variations with $\vert$B$\vert$. We find the intensity, Doppler shift, and line widths to increase with $\vert$B$\vert$ for CHs and QS. The CHs show deficit in intensity and excess total widths over QS for regions with identical $\vert$B$\vert$. For pixels with only upflows, CHs show excess upflows over QS, while for pixels with only downflows, CHs show excess downflows over QS that cease to exist at $\vert$B$\vert$ $\le$ 40. Finally, the spectral profiles are found to be more skewed and flatter than a Gaussian, with no difference between CH and QS. These results are important in understanding the heating of the atmosphere in CH and QS, including solar wind formation, and provide further constraints on the modeling of the solar atmosphere.

Semi-regular variables (SRVs) are similar to Miras in brightness, and they also follow PLRs, though not necessarily the same as Miras. As potential standard candles they are more challenging than Miras due to their smaller variability amplitudes and less regular light curves, but they are substantially more numerous and especially promising to probe old stellar populations. We aim to characterize the variability of SRVs, with focus on their connection with Miras, in order to prepare the ground for investigating their potential as distance indicators. We examine SRVs and Miras in the Magellanic Clouds from OGLE-III observations, with data from Gaia and 2MASS. After cleaning the sample from variability periods unrelated to pulsation, we classify each source by chemical type and combination of pulsation modes. We examine the results in terms of global photometric and pulsation properties. We identify four SRVs groups that fit the general evolutionary scenario predicted by theory. SRVs dominated by fundamental-mode pulsation are very similar to Miras, especially if mono-periodic. They further split into two sub-groups, one of which follows the same sequence as Miras in the period-luminosity and period-amplitude diagram, without discontinuity. The similarities between Miras and SRVs suggest that the latter can be adopted as complementary distance indicators to the former, thereby at least doubling the available number of LPVs suitable for use as distance indicators. The traditional amplitude-based separation between Miras and SRVs is not necessarily appropriate, and a more physically sound criterion should also involve pulsation periods. While this would require comparatively longer time series, they are expected to become accessible in the coming years even for weak sources thanks to current and future large-scale surveys. The table of reclassified LPVs is made public.

Most of the observed and simulated galactic bars are symmetric in the face-on view. However, there are indeed cases of bars that are off-center with respect to the disk or have an asymmetric shape. The only well-known example showing both these features is the Large Magellanic Cloud. We report on the identification of several lopsided galactic bars in the Illustris TNG100 simulation found among a sample of elongated bar-like galaxies studied in the past. The bars show a clear asymmetry in the face-on view, which is in the shape of a footprint. We measured the evolution of the different parameters of the bars' shape and asymmetry as a function of time and find that the asymmetry is preserved for a few Gyr. It can grow together with a bar or appear later, after bar formation. We considered two scenarios leading to the formation of lopsided bars using controlled simulations. In the first, a Milky Way-like galaxy interacts with a massive companion placed on a radial orbit in the plane of the disk and perpendicular to the orientation of the bar at the time of the first passage. In the second, the galaxy initially has an off-center disk and the growth of the bar and its asymmetry is more similar to the one found in IllustrisTNG galaxies, where it is also preceded by the presence of an asymmetric disk. It is thus possible that lopsided bars are born in lopsided disks, although in some cases, the time difference between the occurrence of the asymmetry in the two components is quite large.

Megaconstellations of thousands to tens of thousands of artificial satellites (satcons) are rapidly being developed and launched. These satcons will have negative consequences for observational astronomy research, and are poised to drastically interfere with naked-eye stargazing worldwide should mitigation efforts be unsuccessful. Here we provide predictions for the optical brightnesses and on-sky distributions of several satcons, including Starlink, OneWeb, Kuiper, and StarNet/GW, for a total of 65,000 satellites on their filed or predicted orbits. We develop a simple model of satellite reflectivity, which is calibrated using published Starlink observations. We use this model to estimate the visible magnitudes and on-sky distributions for these satellites as seen from different places on Earth, in different seasons, and different times of night. For latitudes near 50 degrees North and South, satcon satellites make up a few percent of all visible point sources all night long near the summer solstice, as well as near sunrise and sunset on the equinoxes. Altering the satellites' altitudes only changes the specific impacts of the problem. Without drastic reduction of the reflectivities, or significantly fewer total satellites in orbit, satcons will significantly change the night sky worldwide.

The realtime follow-up of neutrino events is a promising approach to search for astrophysical neutrino sources. It has so far provided compelling evidence for a neutrino point source: the flaring gamma-ray blazar TXS 0506+056 observed in coincidence with the high-energy neutrino IceCube-170922A detected by IceCube. The detection of very-high-energy gamma rays (VHE, $\mathrm{E} > 100\,\mathrm{GeV}$) from this source helped establish the coincidence and constrained the modeling of the blazar emission at the time of the IceCube event. The four major imaging atmospheric Cherenkov telescope arrays (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-opportunity observations of neutrino alerts sent by IceCube. This program has two main components. One are the observations of known gamma-ray sources around which a cluster of candidate neutrino events has been identified by IceCube (Gamma-ray Follow-Up, GFU). Second one is the follow-up of single high-energy neutrino candidate events of potential astrophysical origin such as IceCube-170922A. GFU has been recently upgraded by IceCube in collaboration with the IACT groups. We present here recent results from the IACT follow-up programs of IceCube neutrino alerts and a description of the upgraded IceCube GFU system.

Recent observations of the ultra-hot Jupiter WASP-76b have revealed a diversity of atmospheric species. Here we present new high-resolution transit spectroscopy of WASP-76b with GRACES at the Gemini North Observatory, serving as a baseline for the Large and Long Program "Exploring the Diversity of Exoplanet Atmospheres at High Spectral Resolution" (Exoplanets with Gemini Spectroscopy, or ExoGemS for short). With a broad spectral range of $400 - 1050$ nm, these observations allow us to search for a suite of atomic species. We recover absorption features due to neutral sodium (Na I), and report a new detection of the ionized calcium (Ca II) triplet at ~ $850$ nm in the atmosphere of WASP-76b, complementing a previous detection of the Ca II H & K lines. The triplet has line depths of $0.295 \pm 0.034$% at ~ $849.2$ nm, $0.574 \pm 0.041$% at ~ $854.2$ nm, and $0.454 \pm 0.024$% at ~ $866.2$ nm, corresponding to effective radii close to (but within) the planet's Roche radius. These measured line depths are significantly larger than those predicted by model LTE and NLTE spectra obtained on the basis of a pressure-temperature profile computed assuming radiative equilibrium. The discrepancy suggests that the layers probed by our observations are either significantly hotter than predicted by radiative equilibrium and/or in a hydrodynamic state. Our results shed light on the exotic atmosphere of this ultra-hot world, and will inform future analyses from the ExoGemS survey.

K2-25b is a Neptune-sized exoplanet (3.45 Earth radii) that orbits its M4.5 host with a period of 3.48 days. Due to its membership in the Hyades Cluster, the system has a known age (727 +/- 75 Myr). K2-25b's youth and its similarities with Gl 436b suggested that K2-25b could be undergoing strong atmospheric escape. We observed two transits of K2-25b at Lyman-alpha using HST/STIS in order to search for escaping neutral hydrogen. We were unable to detect an exospheric signature, but placed an upper limit of (R_p/R_s) < 0.56 at 95% confidence by fitting the light curve of the Lyman-alpha red-wing, or < 1.20 in the blue-wing. We reconstructed the intrinsic Lyman-alpha profile of K2-25 to determine its Lyman-alpha flux, and analyzed XMM-Newton observations to determined its X-ray flux. Based on the total X-ray and extreme ultraviolet irradiation of the planet (8763 +/- 1049 erg/s/cm^2), we estimated the maximum energy-limited mass loss rate of K2-25b to be 10.6 x 10^10 g/s (0.56 Earth masses per 1 Gyr), five times larger than the similarly estimated mass loss rate of Gl 436b (2.2 x 10^10 g/s). The photoionization time is about 3 hours, significantly shorter than Gl 436b's 14 hours. A non-detection of a Lyman-alpha transit could suggest K2-25b is not significantly losing its atmosphere, or factors of the system are resulting in the mass loss being unobservable (e.g., atmosphere composition or the system's large high energy flux). Further observations could provide more stringent constraints.

The mechanism accelerating cosmic rays in the Milky Way Galaxy and galaxy clusters is identified and described. The acceleration of cosmic rays is a purely electrostatic process which operates up to the maximum energies of $10^{23}$ eV in galaxy clusters. Galactic cosmic rays are accelerated in a pervasive electrostatic field active in the whole Galaxy except in restricted regions shielded by interstellar and stellar plasmas as, for instance, the region occupied by the solar system. It is proved that the energy spectrum of the cosmic radiation in the Milky Way Galaxy, in the region where the solar system resides, has a constant spectral index comprised between 2.64-2.68 and the maximum energies of Galactic protons are $3.0 \times 10^{19}$ eV. The agreement of these results with the experimental data is discussed in detail and underlined. The various physical processes that maintain the stability of the electrostatic structure in the Milky Way Galaxy are the same that generate the Galactic magnetic field. Accordingly, the intensity, orientation and direction of the Galactic magnetic field are evaluated. The results of the calculation are compared with the observational data, optical and mostly radio astronomy data. The accord of the intensity, orientation and direction of the observed magnetic field with calculation is excellent.

Particles of the cosmic radiation, electrons and nuclei, transport a dominant positive electric charge. A tiny fraction of these particles of extremely high energies in favorable conditions overflow from galaxies. The overflowing of positively charged cosmic nuclei into the intergalactic space uncovers an equal amount of negative charge in the parent galaxy. Negative charge is mainly stored by quiescent electrons. After adequate particle propagation neither the negative electric charge located in the galaxies nor the positive electric charge of the overflowed cosmic nuclei can be neutralized due to the enormous distances. In several ways it is proved that the total electric charge retained by clusters of galaxies after an appropriate time interval generate a repulsive force between clusters which overwhelms gravity. After a few billions years of electrostatic repulsion, peripheral clusters attain relativistic velocities and their mutual distances increase accordingly. Several facts suggest that the expansion of the universe, as determined by optical observations since a century, has been caused by the electrostatic repulsion of the positively charged cosmic nuclei overflowed from galaxy clusters.

The Planck satellite has identified more than 2000 protocluster candidates with extreme star formation rates (SFRs). Here, we present the spectroscopic identification of a Planck-selected protocluster located in the Cosmos field, PHz G237.01+42.50 (G237). G237 contains a galaxy overdensity of 31 spectroscopically identified galaxies at z~2.16 (significant at 5.4 sigma) in a 10'x11' region. The overdensity contains two substructures or protoclusters at <z>~2.16 and 2.195 with estimated halo masses at z=0 of ~(5-6)x10^14 Msun. The overdensity total SFR, ~4000 Msun/yr, is higher than predicted by simulations but much smaller than the SFR derived from the Planck data. The analysis of the Herschel data, in combination with the available ancillary data, shows that such a difference is due to an effect of source alignment along the line of sight that produces a 5 sigma overdensity of red Herschel sources in the field. We analyze the members' UV spectra and UV-far-infrared spectral energy distributions to derive their SFR, stellar mass, and metallicity. Galaxy members include blue star-forming galaxies and AGN with SFRs and stellar masses consistent with the main sequence. AGN, identified through optical spectroscopy or X-ray data, represent a significant fraction (20+/-10%) of all members of the protocluster at z=2.16, and they are powerful enough to produce radiative feedback. The core of this protocluster, besides being denser, includes members that are, on average, more massive and star-forming and contains a larger fraction of AGN and Herschel-detected galaxies than the full sample, suggesting an environmental effect on galaxy growth. A comparison between G237 and other protoclusters in the literature at similar redshifts reveals some common traits and differences that reflect both observational biases and a diversity in intrinsic properties that is not yet fully understood.

We exploit a suite of large N-body simulations (up to N=$4096^3$) performed with Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum in dark matter only cosmological N-body simulations. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At the $1\%$ level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, $\sim 2 k_\Lambda$ at the very latest times we evolve to. At the $5\%$ level accuracy extends up to slightly larger wavenumbers, of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple $n$-dependence, allowing a straightforward extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale free models like LCDM. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we note that studies of the matter power spectrum in the literature have often used data at larger wavenumbers, where convergence to the physical result is poor. Even qualitative conclusions about clustering at small scales, e.g concerning the validity of the stable clustering approximation, may need revision in light of our results.

We re-analyze the Cepheid data used to infer the value of $H_0$ by calibrating SnIa. We do not enforce a universal value of the empirical Cepheid calibration parameters $R_W$ (Cepheid Wesenheit color-luminosity parameter) and $M_H^{W}$ (Cepheid Wesenheit H-band absolute magnitude). Instead, we allow for variation of either of these parameters for each individual galaxy. We also consider the case where these parameters have two universal values: one for low galactic distances $D<D_c$ and one for high galactic distances $D>D_c$ where $D_c$ is a critical transition distance. We find hints for a $3\sigma$ level mismatch between the low and high galactic distance parameter values. We then use AIC and BIC criteria to compare and rank the following types of models: Base models: Universal values for $R_W$ and $M_H^{W}$ (no parameter variation), I Individual fitted galactic $R_W$ with a universal fitted $M_H^{W}$, II Universal fixed $R_W$ with individual fitted galactic $M_H^{W}$, III Universal fitted $R_W$ with individual fitted galactic $M_H^{W}$, IV Two universal fitted $R_W$ (near and far) with one universal fitted $M_H^{W}$, V Universal fitted $R_W$ with two universal fitted $M_H^{W}$ (near and far), VI Two universal fitted $R_W$ with two universal fitted $M_H^{W}$ (near and far). We find that the AIC and BIC criteria consistently favor model IV instead of the commonly used Base model where no variation is allowed for the Cepheid empirical parameters. The best fit value of the SnIa absolute magnitude $M_B$ and of $H_0$ implied by the favored model IV is consistent with the inverse distance ladder calibration based on the CMB sound horizon $H_0=67.4\pm 0.5\,km\,s^{-1}\,Mpc^{-1}$. Thus in the context of the favored model IV the Hubble crisis is not present. This model may imply the presence of a fundamental physics transition taking place at a time more recent than $100\,Myrs$ ago.

We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive constraints on both cosmology and baryon feedback for the first time from this data, marginalising over redshift uncertainties, intrinsic alignment of galaxies, and contamination by the cosmic infrared background (CIB). We find our results to be insensitive to the CIB, while intrinsic alignment provides a small but significant contribution to the lensing--tSZ cross-correlation. The cosmological constraints are consistent with those of other low-redshift probes and prefer strong baryon feedback. The inferred amplitude of the lensing--tSZ cross-correlation signal, which scales as $\sigma_8(\Omega_\mathrm{m}/0.3)^{0.2}$, is low by $\sim 2\,\sigma$ compared to the primary cosmic microwave background constraints by Planck. The lensing--tSZ measurements are then combined with pseudo-$C_\ell$ measurements of KiDS-1000 cosmic shear into a novel joint analysis, accounting for the full cross-covariance between the probes, providing tight cosmological constraints by breaking parameter degeneracies inherent to both probes. The joint analysis gives an improvement of 40% on the constraint of $S_8=\sigma_8\sqrt{\Omega_\mathrm{m}/0.3}$ over cosmic shear alone, while providing constraints on baryon feedback consistent with hydrodynamical simulations, demonstrating the potential of such joint analyses with baryonic tracers like the tSZ effect. We discuss remaining modelling challenges that need to be addressed if these baryonic probes are to be included in future precision-cosmology analyses.

Variational principles for magnetohydrodynamics (MHD) were in\-troduced by previous authors both in Lagrangian and Eulerian form. In this paper we introduce simpler Eulerian variational principles from which all the relevant equations of non-barotropic stationary magnetohydrodynamics can be derived for certain field topologies. The variational principle is given in terms of three independent functions for stationary non-barotropic flows. This is a smaller number of variables than the eight variables which appear in the standard equations of non-barotropic magnetohydrodynamics which are the magnetic field $\vec B$ the velocity field $\vec v$, the entropy $s$ and the density $\rho$. We further investigate the case in the flow along magnetic lines is not ideal.

The most sensitive haloscopes that search for axion dark matter through the two photon electromagnetic anomaly, convert axions into photons through the mixing of axions with a large DC magnetic field. In this work we apply Poynting theorem to the resulting axion modified electrodynamics and identify two possible Poynting vectors, one similar to the Abraham Poynting vector and the other to the Minkowski Poynting vector in electrodynamics. The latter picks up the extra non-conservative terms while the former does not. To understand the source of energy conversion and power flow in the detection systems, we apply the two Poynting theorems to axion modified electrodynamics, for both the resonant cavity and broadband low-mass axion detectors. We show that both Poynting theorems give the same sensitivity for a resonant cavity axion haloscope, but predict markedly different sensitivity for a low-mass broadband capacitive haloscope. Hence we ask the question, can understanding which one is the correct one for axion dark matter detection, be considered under the framework of the Abraham-Minkowski controversy? In reality, this should be confirmed by experiment when the axion is detected. However, many electrodynamic experiments have ruled in favour of the Minkowski Poynting vector when considering the canonical momentum in dielectric media. In light of this, we show that the axion modified Minkowski Poynting vector should indeed be taken seriously for sensitivity calculation for low-mass axion haloscope detectors in the quasi static limit, and predict orders of magnitude better sensitivity than the Abraham Poynting vector equivalent.

We discuss the most general field equations of cosmological spacetimes in gravity theories based on non-linear extensions of the non-metricity scalar and the torsion scalar, which do not only involve the well-known general metric describing the gravitational field, but also an affine connection subject to flatness and metric-compatability or symmetry conditions, respectively. While in the simplest cases the connection disappears from the field equations and one obtains the Friedmann equations of General Relativity, for more general functions the connection modifies the metric equations of motion and come along with its own equations of motion. Here, we derive the most general Ans\"atze for both metric and connections and derive the field equations, and while recovering previously known ones we also obtain new connections. In particular, in non-linear extensions of the non-metricity scalar the connection becomes truly dynamical, leading to cosmological connection hairs. We also give some example spacetimes showing both General Relativity solutions and beyond using the dynamical connection. Such a dynamical connection does not appear in the cosmology of non-linear extensions of the the torsion scalar, showing that non-metricity cosmologies are richer. In some specific limits both cosmologies agree at the background level in the spatially flat case.

A fraction of light scalar dark matter, especially axions, may organize into Bose-Einstein condensates, gravitationally bound clumps, "boson stars", and be present in large number in galactic halos today. We compute the expected number of gravitational microlensing events of clumps composed of the ordinary QCD axion and axion-like-particles and derive microlensing constraints from the EROS-2 survey and the Subaru Hyper Suprime-Cam observation. We perform a detailed lensing calculation, including the finite lens and source size effects in our analysis. We constrain the axion mass in terms of the fraction of dark matter collapsed into clumps, the individual clump densities, and the axion self-coupling. We also consider and constrain clumps composed of a generic scalar dark matter candidate with repulsive self-interactions. Our analysis opens up a new window for the potential discovery of dark matter.