Papers for Friday, Nov 19 2021

A new instrumented baffle was installed in Spring 2021 at Virgo surrounding the suspended mirror in the input mode cleaner triangular cavity. It serves as a demonstrator of the technology designed to instrument the baffles in the main arms in the near future. We present, for the first time, results on the measured scattered light distribution inside the cavity as determined by the new device using data collected between May and July 2021, with Virgo in commissioning phase and operating with an input laser power in the cavity of 28.5~W. The sensitivity of the baffle is discussed and the data is compared to scattered light simulations.

Mass and spin are two fundamental properties of astrophysical black holes. While some established, indirect methods are adopted to measure both these properties of active galactic nuclei (AGN) when viewed relatively face-on, very few suggested methods exist to measure these properties when AGN are viewed highly inclined and potentially obscured by large amounts of gas. In this context, we explore the accuracy and performance of a recently proposed method to estimate the spin of AGN through fitting their accretion disk spectral energy distribution, when adapted for highly inclined and obscured systems, and in particular to a sample of six, local water megamasers. For these sources, both the accretion rate and inclination angle are known, allowing us to rely only on the AGN bolometric luminosity to infer their spin. Using the bolometric luminosity as a proxy for the accretion disk peak luminosity, we derive the expected bolometric luminosity as a function of spin. Then, we measure the bolometric luminosity of each source through X-ray spectroscopy, and compare it with the expected value to constrain the spin of the AGN. The quality of the constraints depend critically on the accuracy of the measured bolometric luminosity, which is difficult to estimate in heavily obscured systems. Three out of six sources do not show consistency between the expected and measured bolometric luminosities, while other three (four, when considering the [OIII] line as tracer of the bolometric luminosity) are formally consistent with high spin values. Our results suggest that this method, although promising (and possibly considered as a future calibrator for other methods) needs better observational data and further theoretical modeling to be successfully applied to obscured AGN and to infer robust results.

I describe and demonstrate a new approach to using spectroscopic data to exploit Poisson sampling fluctuations in unresolved stellar populations. The method is introduced using spectra predicted for independent samples of stars from a 10 Gyr population a using a simple stochastic spectral synthesis model. A principal components analysis shows that >99 per cent of the spectral variation in the red-optical can be attributed to just three "fluctuation eigenspectra", which can be related to the number of giant stars present in each sample, and their distribution along the isochrone. The first eigenspectrum effectively encodes the spectrum of the coolest giant branch stars, and is equivalent to the ratio between high- and low-flux pixels discussed in previous literature. The second and third eigenspectra carry higher-order information from which the giant-star spectral sequence can in principle be reconstructed. I demonstrate the method in practice using observations of part of NGC 5128, obtained with the MUSE narrow-field adaptive optics mode. The expected first eigenspectrum is easily recovered from the data, and closely matches the model results except for small differences around the Ca II triplet. The second eigenspectrum is below the noise level of the present observations. A future application of the method would be to the cores of giant ellipticals to probe the spectra of cool giant stars at high metallicity and with element abundance patterns not accessible in the Milky Way.

Using stellar population synthesis models to infer star formation histories (SFHs), we analyse photometry and spectroscopy of a large sample of quiescent galaxies which are members of Sunyaev-Zel'dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < z < 0.9) and 3 clusters in the SPT Hi-z cluster sample (1.25 < z < 1.4). We analyse spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our dataset sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, t$_{50}$) of $z=2.8\pm0.5$, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age-stellar mass relations for the full sample (age of the Universe at $t_{50}$ (Gyr) = $2.52 (\pm0.04) - 1.66 (\pm0.11)$ log$_{10}(M/10^{11} M\odot))$ and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate $\sim0.75$ Gyr earlier than lower mass galaxies. We also find marginally steeper age-mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass.

Blazar flares have been suggested as ideal candidates for enhanced neutrino production. While the neutrino signal of $\gamma$-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Here, we compute the predicted neutrino signal from X-ray flares detected in 66 blazars observed more than 50 times with the X-ray Telescope (XRT) on board the Neil Gehrels Swift Observatory. We consider a scenario where X-ray flares are powered by synchrotron radiation of relativistic protons, and neutrinos are produced through photomeson interactions between protons with their own synchrotron X-ray photons. Using the 1 keV X-ray light curves for flare identification, the 0.5-10 keV fluence of each flare as a proxy for the all-flavour neutrino fluence, and the IceCube point-source effective area for different detector configurations, we calculate the number of muon and antimuon neutrinos above 100 TeV expected for IceCube from each flaring source. The bulk of the neutrino events from the sample originates from flares with durations $\sim 1-10$ d. Accounting for the X-ray flare duty cycle of the sources in the sample, which ranges between $\sim2$ and 24 per cent, we compute an average yearly neutrino rate for each source. The median of the distribution (in logarithm) is $\sim0.03$ yr$^{-1}$, with Mkn 421 having the highest predicted rate $1.2\pm 0.3$ yr$^{-1}$, followed by 3C 273 $(0.33\pm0.03$ yr$^{-1})$ and PG 1553+113 ($0.25\pm0.02$ yr$^{-1}$). Next-generation neutrino detectors together with regular X-ray monitoring of blazars could constrain the duty cycle of hadronic X-ray flares.

Dark matter (DM) from freeze-in or superWIMP production is well known to imprint non-cold DM signatures on cosmological observables. We derive constraints from Lyman-$\alpha$ forest observations for both cases, basing ourselves on a reinterpretation of the existing Lyman-$\alpha$ limits on thermal warm DM. We exclude DM masses below 15 keV for freeze-in, in good agreement with previous literature, and provide a generic lower mass bound for superWIMPs that depends on the mother particle decay width. Special emphasis is placed on the mixed scenario, where contributions from both freeze-in and superWIMP are similarly important. In this case, the imprint on cosmological observables can deviate significantly from thermal warm DM. Furthermore, we provide a modified version of the Boltzmann code CLASS, analytic expressions for the DM distributions, and fits to the DM transfer functions that account for both mechanisms of production. Moreover, we also derive generic constraints from $\Delta N_\mathrm{eff}$ measurements and show that they cannot compete with those arising from Lyman-$\alpha$ observations. For illustration, we apply the above generic limits to a coloured $t$-channel mediator DM model, in which case contributions from both freeze-in through scatterings and decays, as well as superWIMP production can be important. We map out the entire cosmologically viable parameter space, cornered by bounds from Lyman-$\alpha$ observations, the LHC, and Big Bang Nucleosynthesis.

Supersonic gas turbulence is a ubiquitous property of the interstellar medium. The level of turbulence, quantified by the gas velocity dispersion ($\sigma_{\rm g}$), is observed to increase with the star formation rate (SFR) rate of a galaxy, but it is yet not established whether this trend is driven by stellar feedback or gravitational instabilities. In this work we carry out hydrodynamical simulations of entire disc galaxies, with different gas fractions, to understand the origins of the SFR-$\sigma_{\rm g}$ relation. We show that disc galaxies reach the same levels of turbulence regardless of the presence of stellar feedback processes, and argue that this is an outcome of the way disc galaxies regulate their gravitational stability. The simulations match the SFR-$\sigma_{\rm g}$ relation up to SFRs of the order of tens of M$_\odot$ yr$^{-1}$ and $\sigma_{\rm g}\sim 50$ km s$^{-1}$ in neutral hydrogen and molecular gas, but fail to reach the very large values ($>100$ km s$^{-1}$) reported in the literature for rapidly star forming galaxies. We demonstrate that such high values of $\sigma_{\rm g}$ can be explained by 1) insufficient beam smearing corrections in observations, and 2) stellar feedback being coupled to the ionised gas phase traced by recombination lines. Given that the observed SFR-$\sigma_{\rm g}$ relation is composed of highly heterogeneous data, with $\sigma_{\rm g}$ at high SFRs almost exclusively being derived from H$\alpha$ observations of high redshift galaxies with complex morphologies, we caution against analytical models that attempt explain the SFR-$\sigma_{\rm g}$ relation without accounting for these effects.

Using action-based distribution function for the dynamical model of the Milky Way we have estimated its total mass and its density profile. Constraints are coming from the globular cluster proper motions from Gaia EDR3, from the rotation curve based on Gaia DR2 data, and from the vertical force data. We use Bayesian MCMC method to explore the parameters, for which the globular cluster distribution function and the Galactic potential are fully constrained. Numerical simulations are used to study the uncertainties on the potential constraint if considering a possible massive Large Magellanic Could (LMC). We found that a massive LMC (1.5$\times10^{11}$ M$_{\odot}$) will affect the MW mass measurement at large radius, which includes both the Milky Way and the LMC. We also use the FIRE2 Latte cosmological hydrodynamic simulations to make mock data set from a Milky-Way like galaxy that includes many unrelaxed substructures. We test the effect of these unrelaxed substructures on the final results, and found that the measured rotation curve fluctuated around input value within 5 percent. By keeping a large freedom in choosing a priori mass profile for both baryonic and dark matter leads a total mass of the MW that ranges from $5.36_{-0.68}^{+0.81}\times10^{11}$ M$_{\odot}$ to $7.84_{-1.97}^{+3.08} \times 10^{11}$ M$\odot$. This includes the contribution of a putative massive LMC and significantly narrows the MW total mass range published earlier. Such total mass leads to dark matter density at solar position of $0.34_{-0.02}^{+0.02}$ GeV cm$^{-3}$.

The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Juri\'c (2019) to study the mass distribution of the Milky Way halo; the new model allows for the use of all available 6D phase-space measurements. We use kinematic data of halo stars out to $142~{\rm kpc}$, obtained from the H3 Survey and $\textit{Gaia}$ EDR3, to infer the mass of the Galaxy. Inference is carried out with the No-U-Turn sampler, a fast and scalable extension of Hamiltonian Monte Carlo. We report a median mass enclosed within $100~{\rm kpc}$ of $\rm M(<100 \; kpc) = 0.69_{-0.04}^{+0.05} \times 10^{12} \; M_\odot$ (68% Bayesian credible interval), or a virial mass of $\rm M_{200} = M(<216.2_{-7.5}^{+7.5} \; kpc) = 1.08_{-0.11}^{+0.12} \times 10^{12} \; M_\odot$, in good agreement with other recent estimates. We analyze our results using posterior predictive checks and find limitations in the model's ability to describe the data. In particular, we find sensitivity with respect to substructure in the halo, which limits the precision of our mass estimates to $\sim 15\%$.

In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations ($\sim 10-20\%$ $1\sigma$-level) in the extinction, corresponding to $\sim 0.1-1\,$mag level in the stellar types above, on timescales of order months to a year. The fluctuations are surprisingly robust to finite source-size effects as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.

We present the results of ALMA $\sim$2 mm, $\lesssim$1$''$-resolution observations of ten (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity $\gtrsim$ 10$^{11.7}$L$_{\odot}$) at $z <$ 0.15, targeting dense ($>$10$^{4}$ cm$^{-3}$) molecular (HCN, HCO$^{+}$, and HNC J=2-1) and 183 GHz H$_{2}$O 3$_{1,3}$-2$_{2,0}$ emission lines. Higher HCN to HCO$^{+}$ J=2-1 flux ratios are observed in some, but not all, AGN-important ULIRGs than in starburst-classified sources. We detect 183 GHz H$_{2}$O emission in almost all AGN-important ULIRGs, and elevated H$_{2}$O emission is found in two sources with elevated HCN J=2-1 emission, relative to HCO$^{+}$ J=2-1. Except one ULIRG (the Superantennae), the H$_{2}$O emission largely comes from the entire nuclear regions ($\sim$1 kpc), rather than AGN-origin megamaser at the very center ($<<$1 kpc). Nuclear ($\sim$1 kpc) dense molecular gas mass derived from HCO$^{+}$ J=2-1 luminosity is $\gtrsim$a few $\times$ 10$^{8}$M$_{\odot}$, and its depletion time is estimated to be $\gtrsim$10$^{6}$ yr in all sources. Vibrationally excited J=2-1 emission lines of HCN and HNC are detected in a few (U)LIRGs, but those of HCO$^{+}$ are not. It is suggested that in mid-infrared-radiation-exposed innermost regions around energy sources, HCO$^{+}$ and HNC are substantially less abundant than HCN. In our ALMA $\sim$2 mm data of ten (U)LIRGs, two continuum sources are serendipitously detected within $\sim$10$''$, which are likely to be an infrared luminous dusty galaxy at $z >$ 1 and a blazar.

The Enthalpy Based Thermal Evolution of Loops (EBTEL) approximate model for static and dynamic coronal loops is developed to include the effect of a loop cross-sectional area which increases from the base of the transition region (TR) to the corona. The TR is defined as the part of a loop between the top of the chromosphere and the location where thermal conduction changes from an energy loss to an energy gain. There are significant differences from constant area loops due to the manner in which the reduced volume of the TR responds to conductive and enthalpy fluxes from the corona. For static loops with modest area variation the standard picture of loop energy balance is retained, with the corona and TR being primarily a balance between heating and conductive losses in the corona, and downward conduction and radiation to space in the TR. As the area at the loop apex increases, the TR becomes thicker and the density in TR and corona larger. For large apex areas, the coronal energy balance changes to one primarily between heating and radiation, with conduction playing an increasingly unimportant role, and the TR thickness becoming a significant fraction of the loop length. Approximate scaling laws are derived that give agreement with full numerical solutions for the density, but not the temperature. For non-uniform areas, dynamic loops have a higher peak temperature and are denser in the radiative cooling phase by of order 50% than the constant area case for the examples considered. They also show a final rapid cooling and draining once the temperature approaches 1 MK. Although the magnitude of the emission measure will be enhanced in the radiative phase, there is little change in the important observational diagnostic of its temperature dependence.

We present updated cosmological constraints from measurements of the gas mass fractions ($f_{gas}$) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus Cluster at low redshifts, two new clusters at redshifts $z>0.97$, and significantly longer observations of four clusters at $0.6<z<0.9$. Our low-redshift ($z<0.16$) $f_{gas}$ data, combined with the cosmic baryon fraction measured from the cosmic microwave background (CMB), imply a Hubble constant of $h = 0.722 \pm 0.067$. Combining the full $f_{gas}$ data set with priors on the cosmic baryon density and the Hubble constant, we constrain the dark energy density to be $\Omega_\Lambda = 0.865 \pm 0.119$ in non-flat $\Lambda$CDM (cosmological constant) models, and its equation of state to be $w = -1.13_{-0.20}^{+0.17}$ in flat, constant-w models, respectively 41 and 29 per cent tighter than our previous work, and comparable to the best constraints available from other probes. Combining $f_{gas}$, CMB, supernova, and baryon acoustic oscillation data, we also constrain models with global curvature and evolving dark energy. For the massive, relaxed clusters employed here, we find the scaling of $f_{gas}$ with mass to be consistent with a constant, with an intrinsic scatter that corresponds to just 3 per cent in distance.

The Photo-z InfraRed Telescope (PIRT) is an instrument on the Gamow Explorer, currently proposed for a NASA Astrophysics Medium Explorer. PIRT works in tandem with a companion wide-field instrument, the Lobster Eye X-ray Telescope (LEXT), that will identify x-ray transients likely to be associated with high redshift gamma-ray bursts (GRBs) or electromagnetic counterparts to gravitational wave (GW) events. After receiving an alert trigger from LEXT, the spacecraft will slew to center the PIRT field of view on the transient source. PIRT will then begin accumulating data simultaneously in five bands spanning 0.5 - 2.5 microns over a 10 arc-minute field of view. Each PIRT field will contain many hundreds of sources, only one of which is associated with the LEXT transient. PIRT will gather the necessary data in order to identify GRB sources with redshift $z > 6$, with an expected source localization better than 1 arcsec. A near real-time link to the ground will allow timely follow-up as a target of opportunity for large ground-based telescopes or the James Webb Space Telescope (JWST). PIRT will also allow localization and characterization of GW event counterparts. We discuss the instrument design, the on-board data processing approach, and the expected performance of the system.

We present two new rapidly oscillating Ap (roAp) stars, TIC 198781841 and TIC 229960986, discovered in TESS photometric data. The periodogram of TIC 198781841 has a large peak at 166.506 d$^{-1}$ (1.93 mHz), with two nearby peaks at 163.412 d$^{-1}$ (1.89 mHz) and 169.600 d$^{-1}$ (1.96 mHz). These correspond to three independent high-overtone pressure modes, with alternating even and odd $\ell$ values. TIC 229960986 has a high-frequency triplet centered at 191.641 d$^{-1}$ (2.218 mHz), with sidebands at 191.164 d$^{-1}$ (2.213 mHz) and 192.119 d$^{-1}$ (2.224 mHz). This pulsation appears to be a rotationally split dipole mode, with sideband amplitudes significantly larger than that of the central peak; hence, both pulsation poles are seen over the rotation cycle. Our photometric identification of two new roAp stars underscores the remarkable ability of TESS to identify high-frequency pulsators without spectroscopic observations.

The first mass-estimate of an exoplanet around a Sun-like star, 51 Peg b and the first radius measurement of an exoplanet, HD209458b pointed to the challenges of understanding the atmosphere, interior, and evolution of exoplanets including the possibility of mass loss of planets on close-orbits that are exposed to strong irradiation. These discoveries raised the question of heating and inflation mechanisms, and of the nature of these objects in terms of composition compared to the known planets in the Solar system. The field of exoplanet interior modeling was born. Here, we outline and discuss current big science questions: (i) What is the amount of heavy elements in a planet and do all planets possess an iron-rock core? We suggest that a promising and novel approach for exoplanets can be measuring their tidal response in form of the Love numbers h2 and k2. (ii) How much and through what mechanisms are the interiors of planets heated or delayed from cooling? Many strongly irradiated gaseous planets require an additional heat source to explain their large radii. (iii) What is the origin of the observed populations in the radius-period diagram? Objects in and along the radius valley are excellent targets to study planetary formation and evaporation. (iv) What does the composition of rocky planets tell us about their formation? Planets more iron-rich than Mercury are found, as well as planets that if rocky, are depleted in iron with respect to Earth. We do not have yet a reliable formation theory that would explain their existence.

Using Bayesian analyses we study the solar electron density with the NANOGrav 11-year pulsar timing array (PTA) dataset. Our model of the solar wind is incorporated into a global fit starting from pulse times-of-arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into existing PTA software. We perform an ab initio recovery of various solar wind model parameters. We then demonstrate the richness of information about the solar electron density, $n_E$, that can be gleaned from PTA data, including higher order corrections to the simple $1/r^2$ model associated with a free-streaming wind (which are informative probes of coronal acceleration physics), quarterly binned measurements of $n_E$ and a continuous time-varying model for $n_E$ spanning approximately one solar cycle period. Finally, we discuss the importance of our model for chromatic noise mitigation in gravitational-wave analyses of pulsar timing data and the potential of developing synergies between sophisticated PTA solar electron density models and those developed by the solar physics community.

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $\alpha$ in our sample is =1.73 with a $\sigma$ of 0.15. We find that $\alpha$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $\alpha$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate modes that are included. Five different approaches - namely the HLL, HLLC, HLLD, HLLEM and GFORCE schemes - are compared in terms of accuracy and robustness against one- and multi-dimensional standard numerical benchmarks. Our results demonstrate that - for weak or moderate magnetizations - the HLLD Riemann solver yields the most accurate results, followed by HLLC solver(s). The GFORCE approach provides a valid alternative to the HLL solver being less dissipative and equally robust for strongly magnetized environments. Finally, our tests show that the HLLEM Riemann solver is not cost-effective in improving the accuracy of the solution and reducing the numerical dissipation.

The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources.

High-energy neutrinos have thus far been observed in coincidence with time-variable emission from three different accreting black holes: a gamma-ray flare from a blazar (TXS 0506+056), an optical transient following a stellar tidal disruption (AT2019dsg), and an optical outburst from an active galactic nucleus (AT2019fdr). Here we present a unified explanation for the latter two of these sources: accretion flares that reach the Eddington limit. A signature of these events is a luminous infrared reverberation signal from circumnuclear dust that is heated by the flare. Using this property we construct a sample of similar sources, revealing a third event coincident with a PeV-scale neutrino. This sample of three accretion flares is correlated with high-energy neutrinos at a significance of 3.7 sigma. Super-Eddington accretion could explain the high particle acceleration efficiency of this new population.

WR22 = HD 92740 is a bright (V = 6.4 $mag$), intrinsically luminous, double-line WN7h + O9III-V binary exhibiting one sharp 8\% deep eclipse near periastron in its elliptical (e = 0.6) 80-day orbit, when the WR-star passes in front of the O star, with no secondary eclipse. We apply two models (L96, A13) to probe the optical space-based light curves from {\em BRITE-Constellation}, including three separate, complete eclipses, that show increased (o-c) scatter compared to the rest of the observations outside the eclipses, likely due to O-star light encountering WR wind-clumps. L96 is a simple atmospheric-eclipse model, often applied to close WR+O binaries, where the O-star is considered a point-source. A13 considers a finite-disk O-star and allows for atmospheric, photospheric and reflection components to the eclipse, permitting a better characterization of its shape through a more physically realistic description of the structures for both stars in WR22. Nevertheless, A13 is still susceptible to uncertainties in the luminosity of the O-star before unique values for the orbital inclination and WR mass-loss rate can be estimated. We present solutions for the two extremes of the O-star, O9V and O9III. As photometry alone cannot allow us to discriminate between these, we compared our results to the spectral models found in the literature and determined the correct solution to be O9V. Our best-fit A13 Model 1 gives $i = 83.5 \pm 0.4^{\circ}$, with $\dot M_{\rm WR} = (1.86 \pm 0.2) \times 10^{-5} \dot M_{\odot}/yr$. The flux ratio in the red {\em BRITE} band in this model is $F_{\rm O}/F_{\rm WR} = 0.064\pm 0.002$.

Black hole components in LIGO/Virgo mergers are more massive than black holes observed in X-ray binaries. Black--hole spins deduced from model--dependent fits to X-ray binary observations give larger values than those deduced for LIGO/Virgo black holes. These two statements have been used to conclude that LIGO/Virgo black holes and X-ray binary black holes form two distinct populations (apples and oranges). Here we show why this conclusion is not secure.

4FGL J0822.8-4207 is a point source found in the 4FGL-DR2 catalog by the gamma-ray observatory Fermi-LAT and has no known association. We carry out X-ray observations of 4FGL J0822.8-4207 to help understand its nature. We explore two scenarios for the origin of 4FGL J0822.8-4207. In the first case we study the possibility that cosmic rays from the supernova remnant (SNR) Puppis A, seen nearby in the sky, reach the dense gas at the location of the source and produce the gamma-rays through inelastic proton-proton collisions. We apply a standard model for particle diffusion in the interstellar medium and derive the required physical parameters. We find that this scenario for the gamma rays is possible if the gas is located at a distance that is not higher than ~40 pc from Puppis A, unless the SNR is older than 7 kyr or the diffusion coefficient is higher than typical Galactic values, and relatively low-energy cosmic rays are currently escaping from the SNR. In the second scenario, we consider the protostellar jet HH219 as the origin of the GeV source and find the very interesting possibility that particles could be accelerated up to energies of at least several TeV in HH219. This would make this system the first known of its kind to produce gamma-ray emission extending up to hundreds of GeV without any apparent cutoff and an excellent laboratory to study the process of particle acceleration.

Correlations in and with the flux transmission of the Lyman$-\alpha$ (Ly$\alpha$) forest in the spectra of high-redshift quasars are powerful cosmological tools, yet these measurements can be compromised if the intrinsic quasar continuum is significantly uncertain. One particularly problematic case is broad absorption line (BAL) quasars, which exhibit blueshifted absorption associated with many spectral features that are consistent with outflows of up to $\sim0.1c$. As these absorption features can both fall in the forest region and be difficult to distinguish from Ly$\alpha$ absorption, cosmological analyses eliminate the 12 - 16% of quasars that exhibit BALs. In this paper we explore an alternate approach that includes BALs in the Ly$\alpha$ auto correlation function, with the exception of the expected locations of the BAL absorption troughs. This procedure returns over 95% of the pathlength that is lost by the exclusion of BALs, as well as increases the density of sightlines. We show that including BAL quasars reduces the fractional uncertainty in the covariance matrix and correlation function and does not significantly change the shape of the correlation function relative to analyses that exclude BAL quasars. We also evaluate different definitions of BALs, masking strategies, and potential differences in the quasar continuum in the forest region for BALs with different amounts of absorption.

In this Letter we report a discovery of a prominent flash of a peculiar overluminous Type Ia supernova, SN 2020hvf, in about 5 hours of the supernova explosion by the first wide-field mosaic CMOS sensor imager, the Tomo-e Gozen Camera. The fast evolution of the early flash was captured by intensive intranight observations via the Tomo-e Gozen high-cadence survey. Numerical simulations show that such a prominent and fast early emission is most likely generated from an interaction between $0.01~M_{\odot}$ circumstellar material (CSM) extending to a distance of $\sim$$10^{13}~\text{cm}$ and supernova ejecta soon after the explosion, indicating a confined dense CSM formation at the final evolution stage of the progenitor of SN 2020hvf. Based on the CSM-ejecta interaction-induced early flash, the overluminous light curve, and the high ejecta velocity of SN 2020hvf, we suggest that the SN 2020hvf may originate from a thermonuclear explosion of a super-Chandrasekhar-mass white dwarf ("super-$M\rm_{Ch}$ WD"). Systematical investigations on explosion mechanisms and hydrodynamic simulations of the super-$M\rm_{Ch}$ WD explosion are required to further test the suggested scenario and understand the progenitor of this peculiar supernova.

We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg (ATLAS20nuc) discovered in NGC 5012 at $\sim 40$ Mpc. With peak absolute magnitude $M_g \approx -14.9 \pm 0.2$, the bolometric light curve requires only $\approx 7 \times 10^{-3}$ \msol\ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ \msol. The photospheric velocity around maximum is $\sim 4500$ \kms\, implying a low kinetic energy of $E_{51} \approx 0.05 \pm 0.02$ erg. A self-consistent model constructed using the 1D radiative transfer code TARDIS for the early optical spectra is dominated by carbon, oxygen, neon and intermediate mass elements like silicon, sulphur and magnesium. In order to constrain the rates of SNe Iax in the local Universe, we construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) SNe within 60 Mpc at $12^{+13}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+15}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for less than $1\%$ of the Ia rate. We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates.

Low and intermediate mass stars with super solar metallicities comprise a known portion of the universe. Yet yields for asymptotic giant branch (AGB) stars with metallicities greater than $Z=0.04$ do not exist in the literature. This contributes a significant uncertainty to galactic chemical evolution simulations. We present stellar yields of AGB stars for $M=1-8$$M_\odot$ and $Z=0.04-0.10$. We also weight these yields to represent the chemical contribution of a metal-rich stellar population. We find that as metallicity increases, the efficiency of the mixing episodes (known as the third dredge up) on the thermally pulsing AGB (TP-AGB) decrease significantly. Consequently, much of the nucleosynthesis that occurs on the TP-AGB is not represented on the surface of very metal-rich stars. It instead remains locked inside the white dwarf remnant. The temperatures at the base of the convective envelope also decrease with increasing metallicity. For the intermediate mass models, this results in the occurrence of only partial hydrogen burning at this location, if any burning at all. We also investigate heavy element production via the slow neutron capture process (s-process) for three 6$M_\odot$ models: $Z=0.04, 0.05$ and $0.06$. There is minor production at the first s-process peak at strontium, which decreases sharply with increasing metallicity. We find the chemical contributions of our models are dominated by proton capture nucleosynthesis, mixed to the surface during first and second dredge up events. This conclusion is mirrored in our stellar population yields, weighted towards the lower mass regime to reflect the mass distribution within a respective galaxy.

Supernova remnants (SNRs) are important objects in investigating the links among supernova (SN) explosion mechanism(s), progenitor stars, and cosmic-ray acceleration. Non-thermal emission from SNRs is an effective and promising tool for probing their surrounding circumstellar media (CSM) and, in turn, the stellar evolution and mass-loss mechanism(s) of massive stars. In this work, we calculate the time evolution of broadband non-thermal emissions from Type Ib/c SNRs whose CSM structures are derived from the mass-loss history of their progenitors. Our results predict that Type Ib/c SNRs make a transition of brightness in radio and $\gamma$-ray bands from an undetectable dark for a certain period to a re-brightening phase. This transition originates from their inhomogeneous CSM structures in which the SNRs are embedded within a low-density wind cavity surrounded by a high-density wind shell and the ambient interstellar medium (ISM). The "resurrection" in non-thermal luminosity happens at an age of ~1,000 yrs old for a Wolf-Rayet star progenitor evolved within a typical ISM density. Combining with the results of Type II SNR evolution recently reported by Yasuda et al. (2021), this result sheds light on a comprehensive understanding of non-thermal emissions from SNRs with different SN progenitor types and ages, which is made possible for the first time by the incorporation of realistic mass-loss histories of the progenitors.

We revisit the big bang nucleosynthesis (BBN) limits on primordial magnetic fields and/or turbulent motions accounting for the decaying nature of turbulent sources between the time of generation and BBN. This leads to larger estimates for the gravitational wave (GW) signal than previously expected. We address the detection prospects through space-based interferometers (for GWs generated around the electroweak energy scale) as well as pulsar timing arrays and astrometric missions (for GWs generated around the quantum chromodynamics energy scale).

Galactic bars are important drivers of galactic evolution, and yet how they impact the interstellar medium and correspondingly star formation, remains unclear. We present simulation results for two barred galaxies with different formation mechanisms, bars formed in isolation or via a tidal interaction, to consider the spatially and temporally varying trends of star formation. We focus on the early (< 1Gyr) epoch of bar formation so that the interaction is clearly identifiable. The nearby NGC 4303 (isolated) and NGC 3627 (interaction history) are selected as observational analogues to tailor these simulations. Regardless of formation mechanism, both models show similar internal dynamical features, although the interaction appears to promote bar-arm disconnection in the outer disc velocity structure. Both bars trigger similar boosts in star formation (79%; 66%), while the interaction also triggers an earlier 31% burst. Significant morphological dependence is observed in the relation between surface gas and star formation rate. In both cases, the bar component is notably steepest; the arm is similar to the overall disc average; and the inter-arm clearly the shallowest. A distinguishable feature of the tidal disc is the presence of moderately dense, inefficiently star forming gas mostly confined to tidal debris outside the optical disc. The tidal disc also exhibits a unique trend of radially increasing star formation efficiency and a clear dearth of star formation which persists along the bar between the centre and bar ends. These are potential signatures for identifying a barred system post-interaction.

The lateral density distributions (LDD) of inclined cosmic ray air shower are asymmetric and the corresponding iso-density contours are of increasing eccentric ellipses with zenith angles of different showers. The polar asymmetry of the iso-density contours introduces a significant shift of the EAS core, which is quantitatively expressed as a gap length (GL) parameter between the EAS core and the center of the modified density pattern consisting of several equi-density ellipses. The LDD of EAS particles is usually approximated by a particular type of lateral density function (LDF) which is generally assumed to be polar symmetric about the EAS axis, and cannot describe the asymmetric LDDs accurately. A polar angle-dependent modified lateral density function of EASs has been derived analytically by considering the effect of attenuation of EAS particles in the atmosphere. From the simulation studies, it has been found that the GL manifests sensitivity to the cosmic ray mass composition. The cosmic ray mass sensitivity of the lateral shower age is also re-examined by applying the modified LDF to the simulated data.

We further explored the radius-to-frequency mapping in cases of FRBs. An analytical treatment of Lyutikov (2020) is presented. The frequency dependence of the drifting rate and the drifting timescale are obtained. The aberration effect and the twist of the magnetic field lines may result in drifting in both directions. For one FRB, the burst width is larger at lower frequency according to the radius-to-frequency mapping. For the FRB population, the magnetic fields of the repeaters may be larger than that of the non-repeaters. Then, according to the radius-to-frequency mapping, the burst widths of the repeaters will be wider than that of the apparent non-repeaters. If similar window function (or emission cones) like that of pulsars and magnetars is also at work in the case of FRBs, then the window function may explain the single or multiple components of FRB profiles. The radius-to-frequency mapping modeling is to some degree independent of the underlying radio emission mechanism.

The telescopes and the infrastructures may alter the local wind environment around the observatory and further affect the observing environment. After the completion of site testing, it is necessary to analyze the wind environment of the entire site and plan the telescope layout to make use of the excellent conditions scientifically and rationally. Taking a typical observatory as an example, the effect of topographical features on wind environment and the mutual interference between telescope enclosures are analyzed by using Computational Fluid Dynamics (CFD) method. The CFD simulations are compared with the seeing data from Differential Image Motion Monitor (DIMM), the results are in good agreement, which verifies the effectiveness of the CFD method. The results of wind environment analysis can provide reasonable suggestions for site layout and construction, improving the observing environment and the image quality.

The binary status of gamma-Cas stars has been discussed while theoretically examining the origin of their peculiar X-ray emission. However, except in two cases, no systematic radial velocity monitoring of these stars had been undertaken yet to clarify their status. We now fill this gap using TIGRE, CARMENES, and UVES high-resolution spectroscopy. Velocities were determined for 16 stars, revealing shifts and/or changes in line profiles. The orbit of six new binaries could be determined: the long periods (80-120d) and small velocity amplitudes (5-7km/s) suggest low mass companions (0.6-1 M$_{\odot}$). The properties of the known gamma-Cas binaries appear similar to those of other Be systems, with no clear-cut separation between them. One of the new systems is a candidate for a rare case of quadruple system involving a Be star. Five additional gamma-Cas stars display velocity variations compatible with the presence of companions, but no orbital solution could yet be formally established for them hence they only receive the status of "binary candidate".

Testing the DAMA/LIBRA annual modulation result independently of dark matter particle and halo models has been a challenge for twenty years. Using the same target material, NaI(Tl), is required and presently two experiments, ANAIS-112 and COSINE-100, are running for such a goal. A precise knowledge of the detector response to nuclear recoils is mandatory because this is the most likely channel to find the dark matter signal. The light produced by nuclear recoils is quenched with respect to that produced by electrons by a factor that has to be measured experimentally. However, current quenching factor measurements in NaI(Tl) crystals disagree within the energy region of interest for dark matter searches. To disentangle whether this discrepancy is due to intrinsic differences in the light response among different NaI(Tl) crystals, or has its origin in unaccounted for systematic effects will be key in the comparison among the different experiments. We present measurements of the quenching factors for five small NaI(Tl) crystals performed in the same experimental setup to control systematics. Quenching factor results are compatible between crystals and no clear dependence with energy is observed from 10 to 80 keVnr.

The Low Frequency Array (LOFAR) is an international radio telescope array, consisting of 38 stations in the Netherlands and 14 international stations spread over Europe. Here we present an observation method to study the jovian decametric radio emissions from several LOFAR stations (here DE604, FR606 and IE613), at high temporal and spectral resolution. This method is based on prediction tools, such as radio emission simulations and probability maps, and data processing. We report an observation of Io-induced decametric emission from June 2021, and a first case study of the substructures that compose the macroscopic emissions (called millisecond bursts). The study of these bursts make it possible to determine the electron populations at the origin of these emissions. We then present several possible future avenues for study based on these observations. The methodology and study perspectives described in this paper can be applied to new observations of jovian radio emissions induced by Io, but also by Ganymede or Europa, or jovian auroral radio emissions.

We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5\sigma$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM.

Focusing on the electron and positron spectrum measured with CALET, which shows characteristic structures, we calculate flux contributions of cosmic rays escaped from supernova remnants, which were randomly born. We adopt a Monte Carlo method to take into account the stochastic property of births of nearby sources. We find that without a complicated energy dependence of the diffusion coefficient, simple power-law diffusion coefficients can produce spectra similar to the CALET spectrum even with a dispersion in the injection index. The positron component measured with AMS-02 is consistent with a bump-like structure around 300 GeV in the CALET spectrum. One to three nearby supernovae can contribute up to a few tens of percent of the CALET flux at 2--4 TeV, ten or more unknown and distant ($\gtrsim 500$ pc) supernovae account for the remaining several tens of percent of the flux. The CALET spectrum, showing a sharp drop at $\sim 1$ TeV, allows for a contribution of cosmic rays from an extraordinary event which occured $\sim 400$ kyr ago. This type of event releases electrons/positrons with a total energy more than 10 times the average energy for usual supernovae, and its occurrence rate is lower than $1/300$ of the usual supernova rate.

The population of known low- to intermediate-mass exoplanets shows a large spread in densities, which is believed to be due to the diversity of planetary atmospheres and thus controlled by planetary atmospheric mass loss. One of the main drivers of long-term atmospheric escape is the absorption of high-energy XUV radiation from the host star. For main sequence solar-like stars, rotation and XUV radiation are closely connected, with faster rotating stars being XUV brighter and with both rotation and XUV decreasing with time. This evolution, however, does not follow a unique path, as stars born with the same mass and metallicity can have widely different initial rotation rates. This non-uniqueness holds up to about 1 Gyr, while atmospheric escape from exoplanets is strongest. The atmospheric mass loss through this period is often deciding the future of the planet and its position in the observed population. Therefore, the diversity of possible stellar histories can be an uncertain factor affecting the predictions of population studies. Here, I explore its relevance for different planets and different host stars.

We propose a new method to determine the shape of the gravitational potential of the dark matter (DM) halo of the Milky Way (MW) with the galactocentric tangential velocities of a sample of hypervelocity stars (HVSs). We compute the trajectories of different samples of HVSs in a MW where the baryon distribution is axisymmetric and the DM potential either is spherical or is spheroidal or triaxial with radial-dependent axis ratios. We determine the shape of the DM potential with the distribution of the latitudinal velocity $|v_{\vartheta}|$ in axisymmetric Galactic potentials, or with the distribution of $|v_{\vartheta}|$ and of a function $\bar v_{\varphi}$ of the azimuthal velocity in non-axisymmetric Galactic potentials. We recover the correct shape of the DM potential by comparing the distribution of $|v_{\vartheta}|$ and $\bar v_{\varphi}$ against the corresponding distributions of mock samples of HVSs that traveled in DM halos of different shapes. We use the largest possible sample of $\sim 800$ HVSs of $4~M_\odot$ ejected with the Hills mechanism at a rate $\sim 10^{-4}$ yr$^{-1}$, currently outgoing, and located at more than 10 kpc from the Galactic center. In the ideal case of galactocentric velocities with null uncertainties, our method recovers the correct shape of the DM potential with a success rate $S\gtrsim 89\%$ in axisymmetric Galactic potentials, and $S > 96\%$ in the explored non-axisymmetric cases. The unsuccessful cases yield axis ratios of the DM potential that are off by $\pm 0.1$. The success rate decreases with decreasing sample size: for $\sim 80$ HVSs, close to the current number of HVS candidates, $S$ is in the range $\sim 40\%-60\%$, depending on the actual shape of the DM halo. A robust determination of the shape of the DM potential thus requires increasing by a factor $\sim 10$ the size of the sample of genuine HVSs with measured galactocentric velocity.

The discovery of a population of close-orbiting giant planets ($\le$ 1 au) has raised a number of questions about their origins and dynamical histories. These issues have still not yet been fully resolved, despite over 20 years of exoplanet detections and a large number of discovered exoplanets. In particular, it is unclear whether warm Jupiters (WJs) form in situ, or whether they migrate from further outside and are even currently migrating to form hot Jupiters (HJs). Here, we report the possible discovery and characterization of the planets in a highly mutually-inclined ($I_{\rm mut}\simeq 45^\circ$), compact two-planet system (KOI-984), in which the newly discovered warm Jupiter KOI-984$c$ is on a 21.5-day, moderately eccentric ($e\simeq 0.4$) orbit, in addition to a previously known 4.3-day planet candidate KOI-984$b$. Meanwhile, the orbital configuration of a moderately inclined ($I_{\rm mut}\simeq 15^\circ$), low-mass ($m_{c}\simeq 24 M_{\oplus}$;$P_b\simeq 8.6$ days) perturbing planet near 1:2 mean motion resonace with KOI-984$b$ could also well reproduce observed transit timing variations and transit duration variations of KOI-984$b$. Such an eccentric WJ with a close-in sibling would pose a challenge to the proposed formation and migration mechanisms of WJs, if the first scenario is supported with more evidences in near future; this system with several other well-measured inclined WJ systems (e.g., Kepler-419 and Kepler-108) may provide additional clues for the origin and dynamical histories of WJs.

We present modelling of ~0.1arcsec resolution Atacama Large Millimetre/sub-millimeter Array imaging of seven strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Four of these systems are galaxy-galaxy scale strong lenses, with the remaining three being group-scale lenses. Through careful modelling of visibilities, we infer the mass profiles of the lensing galaxies and by determining the magnification factors, we investigate the intrinsic properties and morphologies of the lensed sub-millimetre sources. We find that these sub-millimetre sources all have ratios of star formation rate to dust mass that is consistent with or in excess of the mean ratio for high-redshift sub-millimetre galaxies and low redshift ultra-luminous infrared galaxies. The contribution to the infrared luminosity from possible AGN is not quantified and so could be biasing our star formation rates to higher values. The majority of our lens models have mass density slopes close to isothermal, but some systems show significant differences.

It is widely anticipated that the James Webb Space Telescope (JWST) will be transformative for exoplanet studies. It has even been suggested that JWST could provide the first opportunity to search for biosignatures in an alien atmosphere using transmission spectroscopy. This claim is investigated, specifically for the proposed anoxic biosignature pair CH4-CO2. The most favourable known target is adopted (TRAPPIST-1e), with an assumed atmospheric composition similar to the Archean Earth. Compared to previous studies, a more systematic investigation of the effect that cloud/haze-layers have on the detectability of CH4 and CO2 is performed. In addition to a clear atmosphere scenario, cloud/haze-layers are considered at eight pressure levels between 600mbar and 1mbar. These pressures cover a plausible range for H2O cloud and photochemical haze, based on observations of solar system atmospheres and physical models of tidally-locked planets such as TRAPPIST-1e, although no assumptions regarding the cloud/haze-layer composition are made in this study. For the clear atmosphere and cloud/haze-layer pressures of 600-100mbar, strong (5-sigma) detections of both CH4 and CO2 are found to be possible with approximately 5-10 co-added transits measured using the Near Infrared Spectrograph (NIRSpec) prism, assuming a dry stratosphere. However, approximately 30 co-added transits would be required to achieve the same result if a cloud/haze-layer is present at 10mbar. A cloud/haze-layer at 1mbar would prevent the detection of either molecule with the NIRSpec prism for observing programs up to 50 transits (>200 hours of JWST time), the maximum considered.

Coronal mass ejections (CMEs) are more energetic than any other class of solar phenomena. They arise from the rapid release of up to $10^{33}$ erg of magnetic energy mainly in the form of particle acceleration and bulk plasma motion. Their stellar counterparts, presumably involving much larger energies, are expected to play a fundamental role in shaping the environmental conditions around low-mass stars, in some cases perhaps with catastrophic consequences for planetary systems due to processes such as atmospheric erosion and depletion. Despite their importance, the direct observational evidence for stellar CMEs is almost non-existent. In this way, numerical simulations constitute extremely valuable tools to shed some light on eruptive behavior in the stellar regime. Here we review recent results obtained from realistic modeling of CMEs in active stars, highlighting their key role in the interpretation of currently available observational constraints. We include studies performed on M-dwarf stars, focusing on how emerging signatures in different wavelengths related to these events vary as a function of the magnetic properties of the star. Finally, the implications and relevance of these numerical results are discussed in the context of future characterization of host star-exoplanet systems.

The cumulative effect of the magnetized stellar winds on exoplanets dominates over other forms of star-planet interactions. When combined with photoevaporation, these winds will lead to atmospheric erosion. This is directly connected with the concept of Habitable Zone (HZ) planets around late-type stars. Our knowledge of these magnetized winds is limited, making numerical models useful tools to explore them. In this preliminary study, we focus on solar-like stars exploring how different stellar wind properties scale with one another. We used one of the most detailed physics-based models, the 3D Alfv\'en Wave Solar Model part of the Space Weather ModelingFramework, and applied it to the stellar winds domain. Our simulations showed that the magnetic field topology on the star surface plays a fundamental role in shaping the different stellar wind properties (wind speed, mass loss rate, angular momentum loss rate). We conclude that a characterization of the Alfv\'en surface is crucial when studying star-planet interaction as it can serve as an inner-boundary of the HZ

Suprathermal ions in the corona are thought to serve as seed particles for large gradual solar energetic particle (SEP) events associated with fast and wide coronal mass ejections (CMEs). A better understanding of the role of suprathermal particles as seed populations for SEP events can be made by using observations close to the Sun. We study a series of SEP events observed by the Integrated Science Investigation of the Sun (IS$\odot$IS) suite on board the Parker Solar Probe (PSP) from 2020 May 27 to June 2, during which PSP was at heliocentric distances between $\sim$0.4 and $\sim$0.2 au. These events were also observed by the Ahead Solar TErrestrial RElations Observatory (STEREO-A) near 1 au, but the particle intensity magnitude was much lower than that at PSP. We find that the SEPs should have spread in space as their source regions were distant from the nominal magnetic footpoints of both spacecraft, and the parent CMEs were slow and narrow. We study the decay phase of the SEP events in the $\sim$1-2 MeV proton energy range at PSP and STEREO-A, and discuss their properties in terms of both continuous injections by successive solar eruptions and the distances where the measurements were made. This study indicates that seed particles can be continuously generated by eruptions associated with slow and narrow CMEs, spread over a large part of the inner heliosphere, and remain there for tens of hours, even if minimal particle intensity enhancements were measured near 1 au.

We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restricted to individual clouds. We use this method to simulate star and star cluster formation in spiral arms, investigating massive GMCs and spiral arm regions with lower mass clouds, from two galaxy models with different spiral potentials. After selecting these regions from pre-run galaxy simulations, we re-sample the particles to obtain a higher resolution. We then re-simulate these regions for 3 Myr to study where and how star clusters form. We analyse the early evolution of the embedded star clusters in these regions. We find that the massive GMC regions, which are more common with stronger spiral arms, form more massive clusters than the sections of spiral arms containing lower mass clouds. Clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.

Magnitudes for the VisorSat and Original-design types were analyzed separately and by time. Mean values are compared with those from other large-scale photometric studies, and some signficant differences are noted. The illumination phase functions for Starlink satellites indicate strong forward scattering of sunlight. They are also time-dependent on a scale of months and years. These phase functions improve the predictability of satellite magnitudes. A Starlink Brightness Function tailored to the satellite shape also improves magnitude predictions. Brightness flares lasting a few seconds are characterized and the mean rate of magnitude variation during a pass is determined. Observation planning tools, including graphs and statistics of predicted magnitudes, are discussed and illustrated.

We report on simultaneous 30 - 60 MHz LOFAR / AARTFAAC12 radio observations and CAMS low-light video observations of +4 to -10 magnitude meteors at the peak of the Perseid meteor shower on August 12/13, 2020. 204 meteor trains were imaged in both the radio and optical domain. Aside from scattered artificial radio sources, we identify broadband radio emission from many persistent trains, one of which lingered for up to 6 minutes. Unexpectedly, fewer broadband radio meteor trains were recorded when the experiment was repeated during the 2020 Geminids and 2021 Quadrantids. Intrinsic broadband radio emission was reported earlier by the Long Wavelength Array, but for much brighter meteors and observed with lower spatial resolution. The new results offer insight into the unknown radio emission mechanism.

The consequences of complex disturbed environments in the vicinity of a supermassive black hole are not well represented by standard statistical models of optical variability in active galactic nuclei (AGN). Thus, developing new methodologies for investigating and modeling AGN light curves is crucial. Conditional Neural Processes (CNPs) are nonlinear function models that forecast stochastic time-series based on a finite amount of known data without the use of any additional parameters or prior knowledge (kernels). We provide a CNP algorithm that is specifically designed for simulating AGN light curves. It was trained using data from the All-Sky Automated Survey for Supernovae, which included 153 AGN. We present CNP modeling performance for a subsample of five AGNs with distinctive difficult-to-model properties. The performance of CNP in predicting temporal flux fluctuation was assessed using a minimizing loss function, and the results demonstrated the algorithm's usefulness. Our preliminary parallelization experiments show that CNP can efficiently handle large amounts of data. These results imply that CNP can be more effective than standard tools in modeling large volumes of AGN data (as anticipated from time-domain surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time).

Short-lived radioactive nulcei (half-life $\tau_{1/2}\sim1$ Myr) influence the formation of stars and planetary systems by providing sources of heating and ionization. Whereas many previous studies have focused on the possible nuclear enrichment of our own Solar System, the goal of this paper is to estimate the distributions of short-lived radionuclides (SLRs) for the entire population of stars forming within a molecular cloud. Here we focus on the nuclear species $^{60}$Fe and $^{26}$Al, which have the largest impact due to their relatively high abundances. We construct molecular cloud models and include nuclear contributions from both supernovae and stellar winds. The resulting distributions of SLRs are time dependent with widths of $\sim3$ orders of magnitude and mass fractions $\rho_{\scriptstyle SLR}/\rho_\ast\sim10^{-11}-10^{-8}$. Over the range of scenarios explored herein, the SLR distributions show only modest variations with the choice of cloud structure (fractal dimension), star formation history, and cluster distribution. The most important variation arises from the diffusion length scale for the transport of SLRs within the cloud. The expected SLR distributions are wide enough to include values inferred for the abundances in our Solar System, although most of the stars are predicted to have smaller enrichment levels. In addition, the ratio of $^{60}$Fe/$^{26}$Al is predicted to be greater than unity, on average, in contrast to Solar System results. One explanation for this finding is the presence of an additional source for the $^{26}$Al isotope.

The characteristic age of PSR J1734$-$3333 estimated from its current spin down rate implies that it is a young pulsar ($\tau_c<10$ kyr). But the time derivative of its spin down rate differs markedly from that assumed for normal radio pulsars, meaning its actual age is uncertain. G354.8$-$0.8 is a supernova remnant (SNR) whose centre is located 21 arcmin away of the pulsar, and with a morphology that suggests an association with the pulsar. We want to assess the likelihood of the association between PSR J1734$-$3333 and G354.8$-$0.8 or other nearby supernova remnants quantitatively, with the objective of shedding light on the real age of this pulsar. Observations with the Karl G. Jansky Very Large Array were carried out in 2015 and 2019 that allow precise astrometric measurements and consequently a proper motion estimate for the pulsar. The proper motion was found to be $\mu_\alpha=10\pm10$ mas yr$^{-1}$ and $\mu_\delta=-29\pm11$ mas yr$^{-1}$ (error bars are $1$-$\sigma$). Though marginal, this detection rules out the association with G354.8$-$0.8 because it means the pulsar is not moving away from the centre of the SNR. No SNR consistent with the measured proper motion and an age $\sim\tau_c$ could be found. We also present the first measurement of the spectral index for this pulsar, $\alpha=-1.1\pm0.3$, measured between $1.5$ and $3.0$ GHz. The SNR produced by the birth supernova of PSR J1734$-$3333 could have already faded to undetectable brightness, for which estimates suggest timescales of $10$-$100$ kyr. This and other considerations lead us to conclude that the pulsar is possibly older than $45$-$100$ kyr. PSR J1734$-$3333 is a pulsar with rotational properties that place it between standard radio pulsars and magnetars, and we interpret our result in the context of a possible future life as a magnetar for this pulsar.

Context. There is currently a niche for providing high-cadence, high resolution, time-series optical spectroscopy from space, which can be filled by using a low-cost cubesat mission. The Belgian-led ESA CubeSpec mission is specifically designed to provide space-based, low-cost spectroscopy with specific capabilities that can be optimised for a particular science need. Approved as an ESA in-orbit demonstrator, the CubeSpec satellite's primary science objective will focus on obtaining high-cadence, high resolution optical spectroscopic data to facilitate asteroseismology of pulsating massive stars. Aims. In this first paper, we aim to search for pulsating massive stars suitable for the CubeSpec mission, specifically $\beta$ Cep stars, which typically require time series spectroscopy to identify the geometry of their pulsation modes. Methods. Based on the science requirements needed to enable asteroseismology of massive stars with the capabilities of CubeSpec's spectrograph, we combine a literature study for pulsation with the analysis of recent high-cadence time series TESS photometry to classify the variability for stars brighter than V < 4 mag and between O9 and B3 in spectral type. Results. Among the 90 stars that meet our magnitude and spectral type requirements, we identify 23 promising $\beta$ Cep stars with high-amplitude (non-)radial pulsation modes with frequencies below 7 d$^{-1}$. Using further constraints on projected rotational velocities, pulsation amplitudes and number of pulsation modes, we devise a prioritised target list for the CubeSpec mission according to its science requirements and the potential of the targets for asteroseismology. The full target catalogue further provide a modern TESS-based review of line profile and photometric variability properties among bright O9-B3 stars.

With the increased sensitivity and resolution of radio interferometry the study of the collimation and acceleration region of jets in Active Galactic Nuclei (AGN) has come into focus within the last years. Whereas a large fraction of AGN jets reveal a change from parabolic to conical collimation around the Bondi radius, a small number of sources deviate from this standard picture, including the radio galaxy NGC1052. We study the jet width profile, which provides valuable information about the interplay between the central engine and accretion disk system and the collimation and acceleration zone of the jets. We observed the double-sided active galaxy NGC1052 at six frequencies with the VLBA in 2017 and at 22GHz with RadioAstron in 2016. These data are combined with archival 15, 22, and 43 GHz multi-epoch VLBA observations. From ridge-line fitting we obtained width measurements along the jet and counter-jet which were fitted with broken power-laws. We find a break in the jet collimation profile at ~10^4 R_s (Schwarzschild radii). Downstream of the break the collimation is conical with a power-law index of 1.0 - 1.2 (cylindrical 0; parabolic 0.5; conical 1). The upstream power-law index of 0.36 for the approaching jet is neither cylindrical nor parabolic and for the receding jet with 0.16 close-to cylindrical. Both jets have an opening angle of ~30 degree at a distance of ~10^3 R_S and are well collimated with an opening angle of <10 degrees downstream of the break. There are significant differences in the upstream collimation profile between approaching (Eastern) and receding (Western) jet. Absorption or scattering in the surrounding torus as well as an accretion wind may mimic a cylindrical profile. We need to increase the observing frequencies, which do not suffer from absorption to find the true jet collimation profile upstream of 10^4 R_s.

Most baryonic matter in the universe exists in gaseous form and can be found in structures such as galactic halos and the low-density intergalactic medium. proposed-ray spectroscopy missions such as Athena, Arcus, and Lynx will have the capability to identify absorption lines in spectra toward bright active galactic nuclei (AGNs), which can be used as a tool to probe this missing matter. In this study, we examine the optical fields surrounding 15 primary observing targets and identify the foreground galaxies and galaxy groups that are potential hosts of absorption. We record the basic properties of the potential host and their angular and physical separation from the AGN line of sight. This process is done by marking the location of various galaxies and groups in optical images of the field surrounding the target and plotting their angular separation vs. redshift to gauge physical proximity to the background source. We identify the surrounding objects according to those which have measured redshifts and those that require them.

The magnification of galaxies in modern galaxy surveys induces additional correlations in the cosmic shear, galaxy-galaxy lensing and clustering observables used in modern lensing "3x2pt" analyses, due to sample selection. In this paper, we emulate the magnification contribution to all three observables utilising the SLICS simulations suite, and test the sensitivity of the cosmological model, galaxy bias and redshift distribution calibration to un-modelled magnification in a Stage-IV-like survey using Monte-Carlo sampling. We find that magnification cannot be ignored in any single or combined observable, with magnification inducing $>1\sigma$ biases in the $w_0-\sigma_8$ plane, including for cosmic shear and 3x2pt analyses. Significant cosmological biases exist in the 3x2pt and cosmic shear from magnification of the shear sample alone. We show that magnification induces significant biases in the mean of the redshift distribution where a position sample is analysed, which may potentially be used to identify contamination by magnification.

We report results of our analysis of the Kepler superaperture LC data of the open cluster NGC6791 to search for pulsating sdB stars. We checked all pixels and we found only three sdB stars to be pulsating, KIC2569576 (B3), KIC2438324 (B4) and KIC2437937 (B5). These stars were known to be pulsators before, though we extended data coverage detecting more frequencies and features in their amplitude spectra, i.e. new multiplets and more complete period spacing sequences that we used for identifying geometry of the pulsation modes. The multiplet splittings were also used to derive rotation periods. The remaining known sdBs do not show any pulsation-related light variation down to our detection thresholds. We analyzed already existing spectroscopic observations taken with the HECTOSPEC at the MMT telescope in Smithsonian Arizona and with the GMOS at the Gemini North telescope, and fitted atmospheric parameters using the Balmer lines. Four stars, B3-B6, show atmospheric parameters that are consistent with g-mode dominated sdBs. We detected hints of radial velocity variability in B3, B5, and B6, indicating these three stars may be in binaries.

We propose a novel class of composite models that feature both a technicolor and a composite Higgs vacuum limit, resulting in an asymmetric dark matter candidate. These Techni-Composite Higgs models are based on an extended left-right electroweak symmetry with a pseudo-Nambu Goldstone boson Higgs and stable dark matter candidates charged under a global $\mathrm{U}(1)_X$ symmetry, connected to the baryon asymmetry at high temperatures via the $SU(2)_{\rm R}$ sphaleron. We consider, as explicit examples, four-dimensional gauge theories with fermions charged under a new confining gauge group $G_{\rm HC} $.

Non-trivial inflaton self-interactions can yield calculable signatures of primordial non-Gaussianity that are measurable in cosmic surveys. Surprisingly, we find that the phase transition to slow-roll eternal inflation is often incalculable in the same models. Instead, this transition is sensitive to the non-Gaussian tail of the distribution of scalar fluctuations, which probes physics inside the horizon, potentially beyond the cutoff scale of the Effective Field Theory of Inflation. We demonstrate this fact directly by calculating non-Gaussian corrections to Stochastic Inflation within the framework of Soft de Sitter Effective Theory, from which we derive the associated probability distribution for the scalar fluctuations. We find parameter space consistent with current observations and weak coupling at horizon crossing in which the large fluctuations relevant for eternal inflation can only be determined by appealing to a UV completion. We also show this breakdown of the perturbative description is required for the de Sitter entropy to reflect the number of de Sitter microstates.

We explore the two-dimensional motion of relativistic electrons when they are trapped in magnetic fields having spatial power-law variation. Its impacts include lifting of degeneracy that emerged in the case of the constant magnetic field, special alignment of Landau levels of spin-up and spin-down electrons depending on whether the magnetic field is increasing or decreasing from the centre, splitting of Landau levels of electrons with zero angular momentum from that of positive one and the change in the equation of state of matter. Landau quantization (LQ) in variable magnetic fields has interdisciplinary applications in a variety of disciplines ranging from condensed matter to quantum information. As examples, we discuss the increase in quantum speed of the electron in presence of spatially increasing magnetic field; and the attainment of super Chandrasekhar mass of white dwarfs by taking into account LQ and Lorentz force simultaneously.

The data from the event horizon telescope (EHT) have provided a novel view of the vicinity of the horizon of a black hole (BH), by imaging the region around the light-ring. They have also raised hopes for measuring in the near future, features of the image (or the shadow) related to higher order effects of photons traveling in these regions, such as the appearance of higher order bright rings produced by more than one windings of photons around the light-ring. While the prospect of measuring these fine features of Kerr BHs is very interesting in itself, there are some even more intriguing prospects for observing novel features of possible non-Kerr objects, in the case that the subjects of our images are not the BH solutions of general relativity. In the hope of sufficient resolution being available in the future, we explore in this work the structure and properties of null geodesics around a Hartle-Thorne spacetime that includes a deformation from the Kerr spacetime characterised by the quadrupole deformation $\delta q$. These spacetimes have been found to exhibit a bifurcation of the equatorial light-ring to two off-equatorial light-rings in a range of $\delta q$s and spin parameters. In addition to this, there is a range of parameters where both the equatorial and the off-equatorial light-rings are present. This results in the formation of a pocket that can trap photon orbits. We investigate the properties of these trapped orbits and find that chaotic behaviour emerges. Some of these chaotic orbits are additionally found to be "sticky" and get trapped close to periodic orbits for long times. We also explore how these novel features affect the shadow and find that the off-equatorial light-rings produce distinctive features that deform its circular shape, while the chaotic behaviour associated to the pocket creates features with fractal structure.

Data from current gravitational wave detectors contains a high rate of transient noise (glitches) that can trigger false detections and obscure true astrophysical events. Existing noise-detection algorithms largely rely on model-based methods that may miss noise transients unwitnessed by auxiliary sensors or with exotic morphologies. We propose the Unicorn Multi-window Anomaly-detection Pipeline (UniMAP): a model-free algorithm to identify and characterize transient noise leveraging the Temporal Outlier Factor (TOF) via a multi-window data-resampling scheme. We show this windowing scheme extends the anomaly detection capabilities of the TOF algorithm to resolve noise transients of arbitrary morphology and duration. We demonstrate the efficacy of this pipeline in detecting glitches during LIGO and Virgo's third observing run, and discuss potential applications.

We investigate how inflationary predictions are affected by the difference in the number of $e$-folds between the Jordan and Einstein frames. We study several test models in relation to a Jordan frame defined by the common Higgs-inflation-like non-minimal coupling to gravity and consider two different formulations of gravity: metric and Palatini. We find that the difference is quite contained in case of a metric Jordan frame while can be quite remarkable in case of a Palatini Jordan frame. We also discuss a way to overcome the discrepancy by introducing a frame invariant physical distance and a consequently frame invariant number of $e$-folds.

This article summarizes the talks in the session GW2 of the Sixteenth Marcel Grossmann Meeting on Recent Developments in Theoretical and Experimental General Relativity, Gravitation, and Relativistic Field Theories, 5-10 July, 2021, on Mid-frequency (0.1-10 Hz) gravitational waves: Sources and detection methods with a review on strain power spectral density amplitude of various mid-frequency gravitational wave projects/concepts and with extended summaries on the progress of ZAIGA project and on the conceptual study of AMIGO.