Papers for Friday, Apr 09 2021

The total extragalactic $\gamma$-ray flux provides a powerful probe into the origin and evolution of the highest energy processes in our universe. An important component of this emission is the isotropic $\gamma$-ray background (IGRB), composed of sources that cannot be individually resolved by current experiments. Previous studies have determined that the IGRB can be dominated by either misaligned active galactic nuclei (mAGN) or star-forming galaxies (SFGs). However, these analyses are limited, because they have utilized binary source classifications and examined only one source class at a time. We perform the first combined joint-likelihood analysis that simultaneously correlates the $\gamma$-ray luminosity of extragalactic objects with both star-formation and mAGN activity. We find that SFGs produce 48$^{+33}_{-20}$\% of the total IGRB at~1 GeV and 56$^{+40}_{-23}$\% of the total IGRB at 10 GeV. The contribution of mAGN is more uncertain, but can also be significant. Future work to quantify the radio and infrared properties of nearby galaxies could significantly improve these constraints.

We present updated results of the Large Binocular Telescope Search for Failed Supernovae. This search monitors luminous stars in 27 nearby galaxies with a current baseline of 11~yr of data. We re-discover the failed supernova (SN) candidate N6946-BH1 as well as a new candidate, M101-OC1. M101-OC1 is a blue supergiant that rapidly disappears in optical wavelengths with no evidence for significant obscuration by warm dust. While we consider other options, a good explanation for the fading of M101-OC1 is a failed SN, but follow-up observations are needed to confirm this. Assuming only one clearly detected failed SN, we find a failed SN fraction $f = 0.16^{+0.23}_{-0.12}$ at 90 per~cent confidence. We also report on a collection of stars that show slow ($\sim$decade), large amplitude ($\Delta L/L > 3$) luminosity changes.

The LOFAR Two-metre Sky Survey (LoTSS) is ongoing and plans to map the complete Northern sky in the future. The source catalogue from the public LoTSS first data release covers 1% of the sky and is known to show some correlated noise or fluctuations of the flux density calibration over a few degree scale. Due to its unique and excellent design, observations from LOFAR are expected to be an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe in the future. We explore the LoTSS DR1 to understand the survey systematics and data quality of its very first data release. We produce catalog mocks to determine error estimates and with our detailed and careful analysis, we successfully recover the angular clustering statistics of LoTSS galaxies, which fits the $\Lambda$CDM cosmology reasonably well. We employ a Markov chain Monte Carlo (MCMC) based Bayesian analysis and recover the best galaxy biasing scheme for LoTSS galaxies and also constrain the radial distribution of LoTSS DR1. After masking some noisy and uneven patches and with reasonable flux cuts, the LOFAR survey appears qualified for large-scale cosmological studies. The upcoming data releases from LOFAR are expected to be deeper and wider, thus will be more suitable for drawing cosmological implications.

The detection of redshifted hyperfine line of neutral hydrogen (HI) is the most promising probe of the Epoch of Reionization (EoR). We report an analysis of 55 hours of Murchison Widefield Array (MWA) Phase II drift scan EoR data. The data correspond to a central frequency $\nu_0 = 154.24 \, \rm MHz$ ($z\simeq 8.2$ for the redshifted HI hyperfine line) and bandwidth $B = 10.24 \, \rm MHz$. As one expects greater system stability in a drift scan, we test the system stability by comparing the extracted power spectra from data with noise simulations and show that the power spectra for the cleanest data behave as thermal noise. We compute the HI power spectrum as a function of time in one and two dimensions. The best upper limit on the one-dimensional power spectrum are: $\Delta^2(k) \simeq (1000~\rm mK)^2$ at $k \simeq 0.2$$h~{\rm Mpc}^{-1}$ and at $k \simeq 1$$h~{\rm Mpc}^{-1}$. The cleanest modes, which might be the most suited for obtaining the optimal signal-to-noise, correspond to $k \gtrsim 1$$h~{\rm Mpc}^{-1}$. We also study the time-dependence of the foreground-dominated modes in a drift scan and compare with the expected behaviour.

The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing amongst dark matter models. We use the COCO $N$-body simulations of cold dark matter (CDM) and 3.3keV thermal relic warm dark matter (WDM) to predict the satellite spatial distribution. We demonstrate that the radial distributions of CDM and 3.3keV-WDM luminous satellites are identical if the minimum pre-infall halo mass to form a galaxy is $>10^{8.5}$$\mathrm{M}_{\odot}$ The distribution of dark subhaloes is significantly more concentrated in WDM due to the absence of low mass, recently accreted substructures that typically inhabit the outer parts of a MW halo in CDM. We show that subhaloes of mass $[10^{7},10^{8}]$$\mathrm{M}_{\odot}$ and within 30kpc of the centre are the stripped remnants of larger haloes in both models. Therefore their abundance in WDM is $3\times$ higher than one would anticipate from the overall WDM subhalo population. We estimate that differences between CDM and WDM concentration--mass relations can be probed for subhalo--stream impact parameters $<2$kpc. Finally, we find that the impact of WDM on planes of satellites is likely negligible. Precise predictions will require further work with high resolution, self-consistent hydrodynamical simulations.

R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy.

Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects of such fossil groups to be particularly strong. This paper explores the exceptional system $\mathrm{J}143454.4+522850$. We combine gravitational lensing with stellar population-synthesis to separate the total mass of the lens into stars and dark matter. The enclosed mass profiles are contrasted with state-of-the-art galaxy formation simulations, to conclude that SW05 is likely a fossil group with a high stellar to dark matter mass fraction $0.027\pm0.003$ with respect to expectations from abundance matching $0.012\pm0.004$, indicative of a more efficient conversion of gas into stars in fossil groups.

The chemical inventory of planets is determined by the physical and chemical processes that govern the early phases of star formation. The aim is to investigate N-bearing complex organic molecules towards two Class 0 protostars (B1-c and S68N) at millimetre wavelengths with ALMA. Next, the results of the detected N-bearing species are compared with those of O-bearing species for the same and other sources. ALMA observations in Band 6 ($\sim$ 1 mm) and Band 5 ($\sim$ 2 mm) are studied at $\sim$ 0.5" resolution, complemented by Band 3 ($\sim$ 3 mm) data in a $\sim$ 2.5" beam. NH2CHO, C2H5CN, HNCO, HN13CO, DNCO, CH3CN, CH2DCN, and CHD2CN are identified towards the investigated sources. Their abundances relative to CH3OH and HNCO are similar for the two sources, with column densities that are typically an order of magnitude lower than those of O-bearing species. The largest variations, of an order of magnitude, are seen for NH2CHO abundance ratios with respect to HNCO and CH3OH and do not correlate with the protostellar luminosity. In addition, within uncertainties, the N-bearing species have similar excitation temperatures to those of O-bearing species ($\sim$ 100 $\sim$ 300 K). The similarity of most abundances with respect to HNCO, including those of CH2DCN and CHD2CN, hints at a shared chemical history, especially the high D/H ratio in cold regions prior to star formation. However, some of the variations in abundances may reflect the sensitivity of the chemistry to local conditions such as temperature (e.g. NH2CHO), while others may arise from differences in the emitting areas of the molecules linked to their different binding energies in the ice. The two sources discussed here add to the small number of sources with such a detailed chemical analysis on Solar System scales. Future JWST data will allow a direct comparison between the ice and gas abundances of N-bearing species.

Galaxy mergers are traditionally one of the favoured mechanisms for quenching star formation. To test this paradigm in the context of modern cosmological simulations, we use the IllustrisTNG simulation to investigate the impact of individual merger events on quenching (i.e. star formation rate at least 3sig below the star-forming main sequence) within 500Myr after the coalescence phase.The rate of quenching amongst recently merged galaxies is compared with a control sample that is matched in redshift, stellar mass, star formation rate (SFR), black hole mass and environment.We find quenching to be uncommon among the descendants of post-merger galaxies, with only 5% of galaxies quenching within 500 Myr after the merger.Despite this low absolute rate, we find that quenching occurs in post-mergers at twice the rate of the control galaxies.The fraction of quenched post-merger descendants 1.5 Gyr after the merger becomes statistically indistinguishable from that of non-post-mergers, suggesting that mergers could speed up the quenching process in those post-mergers whose progenitors had physical conditions able to sustain effective active galactic nuclei (AGN) kinetic feedback, thus capable of removing gas from galaxies.Our results indicate that although quenching does not commonly occur promptly after coalescence, mergers nonetheless do promote the cessation of star formation in some post-mergers. We find that, in IllustrisTNG, it is the implementation of the AGN kinetic feedback that is responsible for quenching post-mergers, as well as non-post-merger controls.As a result of the released kinetic energy, galaxies experience gas loss and eventually, they will quench.Galaxies with an initially low gas fraction show a preferable pre-disposition towards quenching.The primary distinguishing factor between quenched and star-forming galaxies is gas fraction, with a sharp boundary at fgas=0.1 in TNG.

The $\Lambda$CDM prediction of $S_8\equiv\sigma_8(\Omega_m/0.3)^{0.5}$ -- where $\sigma_8$ is the root mean square of matter fluctuations on a 8 $h^{-1}$Mpc scale -- once calibrated on Planck CMB data is $2-3\sigma$ lower than its direct estimate by a number of weak lensing surveys. In this paper, we explore the possibility that the '$S_8$-tension' is due to a non-thermal hot dark matter (HDM) fractional contribution to the universe energy density leading to a power suppression at small-scales in the matter power spectrum. Any HDM models can be characterized by its effective mass $ m_{sp}^{\rm eff}$ and its contribution to the relativistic degrees of freedom at CMB decoupling $\Delta N_{\rm eff}$. Taking the specific example of a sterile particle produced from the decay of the inflaton during a matter dominated era, we find that from Planck only the tension can be reduced below $2\sigma$, but Planck does not favor a non-zero ${m_{sp}^{\rm eff},\Delta N_{\rm eff}}$. In combination with a measurement of $S_8$ from KIDS1000+BOSS+2dfLenS, the $S_8$-tension would hint at the existence of a particle of mass $ m_{sp}^{\rm eff} \simeq 0.67_{-0.48}^{+0.26}$ ${\rm eV}$ with a contribution to $\Delta N_{\rm eff} \simeq0.06\pm0.05$. However, Pantheon and BOSS BAO/$f\sigma_8$ data restricts the particle mass to $m_{sp}^{\rm eff} \simeq 0.48_{-0.36}^{+0.17}$ and contribution to $\Delta N_{\rm eff} \simeq 0.046_{-0.031}^{+0.004}$. We discuss implications of our results for other canonical non-thermal HDM models -- the Dodelson-Widrow model and a thermal sterile particle with a different temperature in the hidden sector. We report competitive results on such hidden sector temperature which might have interesting implications for particle physics model building, in particular connecting the $S_8$-tension to the longstanding short baseline oscillation anomaly.

At second-order, scalar perturbations can source traceless and transverse perturbations to the metric, called induced gravitational waves (IGW). The apparent gauge-dependence of the IGW obscures the interpretation of the stochastic gravitational-wave signal. To elucidate the gauge dependence, we study the IGW from manifestly gauge-invariant scalar perturbations, namely, the relative velocity between the baryon and cold dark matter. From this relative velocity perturbation, we compute the dimensionless gravitational wave power spectrum and the corresponding expected angular power spectrum of the B-mode polarization of the cosmic microwave background. Although the effect turns out to be unobservably small, the calculation demonstrates both the importance of using observable quantities to remove the gauge ambiguity and the observable consequences of tensor perturbations which are not propagating gravitational waves.

We report 1.2 mm polarized continuum emission observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the high-mass star formation region G5.89-0.39. The observations show a prominent 0.2 pc north-south filamentary structure. The UCHII in G5.89-0.39 breaks the filament in two pieces. Its millimeter emission shows a dusty belt with a mass of 55-115 M$_{\odot}$ and 4,500 au in radius, surrounding an inner part comprising mostly ionized gas with a dust emission only accounting about 30% of the total millimeter emission. We also found a lattice of convex arches which may be produced by dragged dust and gas from the explosive dispersal event involving the O5 Feldt's star. The north-south filament has a mass between 300-600 M$_{\odot}$ and harbours a cluster of about 20 millimeter envelopes with a median size and mass of 1700 au and 1.5 M$_{\odot}$, respectively, some of which are already forming protostars. We interpret the polarized emission in the filament as mainly coming from magnetically aligned dust grains. The polarization fraction is ~4.4% in the filaments and 2.1% at the shell. The magnetic fields are along the North Filament and perpendicular to the South Filament. In the Central Shell, the magnetic fields are roughly radial in a ring surrounding the dusty belt between 4,500 and 7,500 au, similar to the pattern recently found in the surroundings of Orion BN/KL. This may be an independent observational signpost of explosive dispersal outflows and should be further investigated in other regions.

The Ly$\alpha$ emission line is one of the most promising probes of cosmic reionisation but isolating the signature of a change in the ionisation state of the IGM is challenging because of intrinsic evolution and internal radiation transfer effects. We present the first study of the evolution of Ly$\alpha$ emitters (LAE) during the epoch of reionisation based on a full radiation-hydrodynamics cosmological simulation that is able to capture both the large-scale process of reionisation and the small-scale properties of galaxies. We predict the Ly$\alpha$ emission of galaxies in the $10^3$ cMpc$^3$ SPHINX simulation at $6\leq z\leq9$ by computing the full Ly$\alpha$ radiation transfer from ISM to IGM scales. SPHINX is able to reproduce many observational constraints such as the UV/Ly$\alpha$ luminosity functions and stellar mass functions at z $\geq$ 6 for the dynamical range probed by our simulation ($M_{\rm 1500}\gtrsim-18$, $L_{\rm Ly\alpha}\lesssim10^{42}$ erg/s, $M_{\star}\lesssim10^9$ M$_{\odot}$). As intrinsic Ly$\alpha$ emission and internal Ly$\alpha$ escape fractions barely evolve from $z=6$ to 9, the observed suppression of Ly$\alpha$ luminosities with increasing redshift is fully attributed to IGM absorption. For most observable galaxies ($M_{\rm 1500}\lesssim-16$), the Ly$\alpha$ line profiles are slightly shifted to the red due to internal radiative transfer effects which mitigates the effect of IGM absorption. Overall, the enhanced Ly$\alpha$ suppression during reionisation traces the IGM neutral fraction $x_{\rm HI}$ well but the predicted amplitude of this reduction is a strong function of the Ly$\alpha$ peak shift, which is set at ISM/CGM scales. We find that a large number of LAEs could be detectable in very deep surveys during reionisation when $x_{\rm HI}$ is still $\approx 50\%$.

The atmospheres of gaseous giant exoplanets orbiting close to their parent stars (hot Jupiters) have been probed for nearly two decades. They allow us to investigate the chemical and physical properties of planetary atmospheres under extreme irradiation conditions. Previous observations of hot Jupiters as they transit in front of their host stars have revealed the frequent presence of water vapour and carbon monoxide in their atmospheres; this has been studied in terms of scaled solar composition under the usual assumption of chemical equilibrium. Both molecules as well as hydrogen cyanide were found in the atmosphere of HD 209458b, a well studied hot Jupiter (with equilibrium temperature around 1,500 kelvin), whereas ammonia was tentatively detected there and subsequently refuted. Here we report observations of HD 209458b that indicate the presence of water (H2O), carbon monoxide (CO), hydrogen cyanide (HCN), methane (CH4), ammonia (NH3) and acetylene (C2H2), with statistical significance of 5.3 to 9.9 standard deviations per molecule. Atmospheric models in radiative and chemical equilibrium that account for the detected species indicate a carbon-rich chemistry with a carbon-to-oxygen ratio close to or greater than 1, higher than the solar value (0.55). According to existing models relating the atmospheric chemistry to planet formation and migration scenarios, this would suggest that HD 209458b formed far from its present location and subsequently migrated inwards. Other hot Jupiters may also show a richer chemistry than has been previously found, which would bring into question the frequently made assumption that they have solar-like and oxygen-rich compositions.

Some hydrogen poor superluminous supernovae (SLSNe) exhibit bumps in the tails of their light-curves associated with hydrogen features in their late time spectra. Here we use the explosion parameters of one such SLSN -- SN 2017gci -- to search the stellar models of the Binary Population And Spectral Synthesis (BPASS) code for potential progenitors. We find good matches for a 30 $M_{\odot}$ progenitor star in a binary system and no matches from single star models. Common envelope and mass transfer after the giant branch, combined with increased mass loss from strong stellar winds immediately before death, allow the progenitor to lose its hydrogen envelope decades before the explosion. This results in a hydrogen poor SLSN and allows for delayed interaction of the ejecta with the lost stellar material

The goal of this paper is to study the smallest brightening events observed in the EUV quiet Sun. We use commissioning data taken by the EUI instrument onboard the recently launched Solar Orbiter mission. On 2020 May 30, EUI was situated at 0.556AU from the Sun. Its HRIEUV telescope 17.4nm passband reached an exceptionally high two-pixel spatial resolution of 400km. The size and duration of small-scale structures is determined in the HRIEUV data, while their height is estimated from triangulation with the simultaneous SDO/AIA data. This is the first stereoscopy of small scale brightenings at high resolution. We observed small localised brightenings ("campfires") in a quiet Sun region with lengthscales between 400km and 4000km and durations between 10 and 200s. The smallest and weakest of these HRIEUV brightenings have not been observed before. Simultaneous HRILYA observations do not show localised brightening events, but the locations of the HRIEUV events correspond clearly to the chromospheric network. Comparison with simultaneous AIA images shows that most events can also be identified in the 17.1nm, 19.3nm, 21.1nm, and 30.4nm passbands of AIA, although they appear weaker and blurred. DEM analysis indicates coronal temperatures peaking at log(T)~6.1-6.15. We determined the height of a few campfires, which is between 1000 and 5000km above the photosphere. We conclude that "campfires" are mostly coronal in nature and are rooted in the magnetic flux concentrations of the chromospheric network. We interpret these events as a new extension to the flare/microflare/nanoflare family. Given their low height, the EUI "campfires" could be a new element of the fine structure of the transition region/low corona: apexes of small-scale loops that are internally heated to coronal temperatures.

We present theoretical predictions for the free-free emission at cm wavelengths obtained from photoevaporation and MHD wind disk models adjusted to the case of the TW Hydrae young stellar object. For this system, disk photoevaporation with heating due to the high-energy photons from the star has been proposed as a possible mechanism to open the gap observed in the dust emission with ALMA. We show that the photoevaporation disk model predicts a radial profile for the free-free emission that is made of two main spatial components, one originated from the bound disk atmosphere at 0.5-1 au from the star, and another more extended component from the photoevaporative wind at larger disk radii. We also show that the stellar X-ray luminosity has a significant impact on both these components. The predicted radio emission from the MHD wind model has a smoother radial distribution which extends to closer distances to the star than the photoevaporation case. We also show that a future radio telescope such as the \textit{Next Generation Very Large Array} (ngVLA) would have enough sensitivity and angular resolution to spatially resolve the main structures predicted by these models.

Systematic errors in the galaxy redshift distribution $n(z)$ can propagate to systematic errors in the derived cosmology. We characterize how the degenerate effects in tomographic bin widths and galaxy bias impart systematic errors on cosmology inference using observational data from the Deep Lens Survey. For this we use a combination of galaxy clustering and galaxy-galaxy lensing. We present two end-to-end analyses from the catalogue level to parameter estimation. We produce an initial cosmological inference using fiducial tomographic redshift bins derived from photometric redshifts, then compare this with a result where the redshift bins are empirically corrected using a set of spectroscopic redshifts. We find that the derived parameter $S_8 \equiv \sigma_8 (\Omega_m/.3)^{1/2}$ goes from $.841^{+0.062}_{-.061}$ to $.739^{+.054}_{-.050}$ upon correcting the n(z) errors in the second method.

As astronomical photometric surveys continue to tile the sky repeatedly, the potential to pushdetection thresholds to fainter limits increases; however, traditional digital-tracking methods cannotachieve this efficiently beyond time scales where motion is approximately linear. In this paper weprototype an optimal detection scheme that samples under a user defined prior on a parameterizationof the motion space, maps these sampled trajectories to the data space, and computes an optimalsignal-matched filter for computing the signal to noise ratio of trial trajectories. We demonstrate thecapability of this method on a small test data set from the Dark Energy Camera. We recover themajority of asteroids expected to appear and also discover hundreds of new asteroids with only a fewhours of observations. We conclude by exploring the potential for extending this scheme to larger datasets that cover larger areas of the sky over longer time baselines.

We conduct the first 3D hydrodynamic simulations of oxygen-neon white-dwarf-neutron star/black hole mergers (ONe WD-NS/BH mergers). Such mergers constitute a significant fraction, and may even dominate, the inspiral rates of all WD-NS binaries. We post-process our simulations to obtain the nuclear evolution of these systems and couple the results to a supernova spectral synthesis code to obtain the first light curves and spectra for these transients. We find that the amount of $^{56}$Ni synthesised in these mergers grows as a strong function of the WD mass, reaching typically $0.05$ and up to $0.1\,{\rm M}_\odot$ per merger. Photodisintegration leads to similar amounts of $^4$He and about a ten times smaller amount of $^1$H. The nuclear yields from these mergers, in particular those of $^{55}$Mn, may contribute significantly to Galactic chemical evolution. The transients expected from ONe WD-NS mergers are dominantly red/infrared, evolve on month-long timescales and reach bolometric magnitudes of up to -16.5. The current surveys must have already detected these transients or are, alternatively, putting strong constraints on merger scenarios. The properties of the expected transients from WD-NS mergers best agree with faint type Iax supernovae. The Vera Rubin Observatory (LSST) will be detecting up to hundreds of merging ONe WD-NS systems per year. We simulate a subset of our models with 2D axisymmetric FLASH code to investigate why they have been challenging for previous studies. We find that the likely main challenge has been effectively modelling the nuclear statistical equilibrium regime in such mergers.

The relationship between galaxy characteristics and the reionization of the universe remains elusive, mainly due to the observational difficulty in accessing the Lyman continuum (LyC) at these redshifts. It is thus important to identify low-redshift LyC-leaking galaxies that can be used as laboratories to investigate the physical processes that allow LyC photons to escape. The weakness of the [S II] nebular emission lines relative to typical star-forming galaxies has been proposed as a LyC predictor. In this paper, we show that the [S II]-deficiency is an effective method to select LyC-leaking candidates using data from the Low-redshift LyC Survey, which has detected flux below the Lyman edge in 35 out of 66 star-forming galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We show that LyC leakers tend to be more [S II]-deficient and that the fraction of their detections increases as [S II]-deficiency becomes more prominent. Correlational studies suggest that [S II]-deficiency complements other LyC diagnostics (such as strong Lyman-$\alpha$ emission and high [O III]/[O II]). Our results verify an additional technique by which reionization-era galaxies could be studied.

This paper reviews improved calibration methods for the Lunar Reconnaissance Orbiter Lunar Exploration Neutron Detector. We cross calibrated the set of LEND observations and models of its detectors physical geometry and composition against the McKinney Apollo 17 era measured neutron flux, Lunar Prospector Neutron Spectrometer epithermal neutron observations, Earth based Galactic Cosmic Ray observations and altitude dependent models of the Moon neutron emission flux. Our neutron transport modeling of the LEND system with the Geant4 software package allows us to fully decompose the varying contributions of lunar, spacecraft and instrument dependent sources of neutrons and charged particles during the LEND mission. With this improved calibration, we can now fully predict every observation from the eight helium 3 detectors and the expected total and partial count rates of neutrons and charged particles for the entirety of LEND now ten plus year observation campaign at the Moon. The study has resulted in an improved calibration for all detectors. The high spatial resolution of LEND collimated and uncollimated sensors are illustrated using the neutron suppression region associated with the south polar Cabeus permanent shadowed region.

Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio wavelengths, finding enhancement of the X-ray emission by $3.8\pm0.7\%$ (a 5.4$\sigma$ detection) coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts.

The critical component of radio astronomy radiometers built to detect redshifted 21-cm signals from Cosmic Dawn is the antenna element. We describe the design and performance of an octave bandwidth cone disc antenna built to detect this signal in the band 40 to 90 MHz. The Cosmic Dawn signal is predicted to be a wideband spectral feature orders of magnitude weaker than sky and ground radio brightness. Thus, the engineering challenge is to design an antenna at low frequencies that is able to provide with high fidelity the faint cosmological signal, along with foreground sky, to the receiver. The antenna characteristics must not compromise detection by imprinting any confusing spectral features on the celestial radiation, ground emission or receiver noise. An innovation in the present design is making the antenna electrically smaller than half wavelength and operating it on the surface of a sufficiently large water body. The homogeneous and high permittivity medium beneath the small cone-disc antenna results in an achromatic beam pattern, high radiation efficiency and minimum unwanted confusing spectral features. The antenna design was optimized in WIPL-D and FEKO. A prototype was constructed and deployed on a lake to validate its performance with field measurements. Index Terms: Antenna measurements, radio astronomy, reflector antennas.

Size distribution of sunspots provides key information about the generation and emergence processes of the solar magnetic field. Previous studies on the size distribution have primarily focused on either the whole group or individual spot areas. In this paper, we investigate the organization of spot areas within sunspot groups. In particular, we analyze the ratio, $\rm{R}$, of the area of the biggest spot ($\rm{A_{big\_spot}}$) inside a group, to the total area of that group ($\rm{A_{group}}$). We use sunspot observations from Kislovodsk, Pulkovo and Debrecen observatories, together covering solar cycles 17 to 24. We find that at the time when the group area reaches its maximum, the single biggest spot in a group typically occupies about 60% of the group area. For half of all groups, $\rm R$ lies in the range between roughly 50% and 70%. We also find R to change with the group area, $\rm{A_{group}}$, such that $\rm{R}$ reaches a maximum of about 0.65 for groups with $\rm{A_{group}}\approx 200\mu$Hem and then remains at about 0.6 for lager groups. Our findings imply a scale invariant emergence pattern, providing an observational constraint on the emergence process. Furthermore, extrapolation of our results to larger sunspot groups may have a bearing on the giant unresolved starspot features found in Doppler images of highly active sun-like stars. Our results suggest that such giant features are composed of multiple spots, with the largest spot occupying roughly 55--75% of the total group area (i.e. of the area of the giant starspots seen in Doppler images).

High velocity clumps joined to the outflow source by emission with a "Hubble law" ramp of linearly increasing radial velocity vs. distance are observed in some planetary nebulae and in some outflows in star formation regions. We propose a simple model in which a "clump" is ejected from a source over a period $\tau_0$, with a strong axis to edge velocity stratification. This non-top hat cross section results in the production of a highly curved working surface (initially being pushed by the ejected material, and later coasting along due to its inertia). From both analytic models and numerical simulations we find that this working surface has a linear velocity vs. position ramp, and therefore reproduces in a qualitative way the "Hubble law clumps" in planetary nebulae and outflows from young stars.

We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like structure and a two-arm/bi-symmetric spiral pattern within it, surrounded by a starbursting ring. The CND has a relatively low $\sigma/V\approx0.35$ ($r\sim0.5$") and $\sim0.19$ ($r>0.5"$), suggesting that the gas is dynamically settled and suitable for dynamically deriving the mass of its central source. As is expected from X-ray dominated region (XDR) effects that dramatically increase an atomic carbon abundance by dissociating CO molecules, we suggest that the atomic [CI](1-0) emission is a better probe of SMBH masses than CO emission in AGNs. Our dynamical model using the ${\rm [CI]}$(1-0) kinematics yields a $M_{\rm BH}=1.78^{+2.69}_{-1.10}\times10^7$M$_\odot$ and $M/L_{\rm F547M}=2.25^{+0.40}_{-0.43}$ (M$_\odot$/L$_\odot$). The model using the CO(1-0) kinematics also gives a consistent $M_{\rm BH}$ with a larger uncertainty, up to an order of magnitude, i.e.\ $M_{\rm BH}=1.60^{+11.52}_{-1.45}\times10^7$M$_\odot$. This newly dynamical $M_{\rm BH}$ is $\approx$ 2 times higher than the mass determined from the reverberation mapped (RM) method using emissions arising in the unresolved broad-line region (BLR). Given this new $M_{\rm BH}$, we are able to constrain the specific RM dimensionless scaling factor of $f=7.2^{+4.2}_{-3.4}$ for the AGN BLR in NGC7469. The gas within the unresolved BLR thus has a Keplerian virial velocity component and the inclination of $i\approx11.0^\circ$$_{-2.5}^{+2.2}$, confirming its face-on orientation in a Seyfert 1 AGN by assuming a geometrically thin BLR model.

The probability distribution of temperature of a blackbody can be determined from its power spectrum. This technique is called blackbody radiation inversion. In the present paper blackbody radiation inversion is applied on the spectrum of the Sun. The probability distribution of temperature and the mean temperature of the Sun are calculated without assuming a homogenous temperature and without using Stefan-Boltzmann law. Different properties of this distribution are characterized. This paper presents the very first mention and investigation of the distortions present within the Sun's spectrum.

We investigate how the choice of equation of state (EOS) and resolution conspire to affect the outcomes of giant impact (GI) simulations. We focus on the simple case of equal mass collisions of two Earth-like $0.5\,M_\oplus$ proto-planets showing that the choice of EOS has a profound impact on the outcome of such collisions as well as on the numerical convergence with resolution. In simulations where the Tillotson EOS is used, impacts generate an excess amount of vapour due to the lack of a thermodynamically consistent treatment of phase transitions and mixtures. In oblique collisions this enhances the artificial angular momentum (AM) transport from the planet to the circum-planetary disc reducing the planet's rotation period over time. Even at a resolution of $1.3 \times 10^6$ particles the result is not converged. In head-on collisions the lack of a proper treatment of the solid/liquid-vapour phase transition allows the bound material to expand to very low densities which in turn results in very slow numerical convergence of the critical specific impact energy for catastrophic disruption $Q_{RD}^*$ with increasing resolution as reported in prior work. The simulations where ANEOS is used for oblique impacts are already converged at a modest resolution of $10^5$ particles, while head-on collisions converge when they evidence the post-shock formation of a dense iron-rich ring, which promotes gravitational re-accumulation of material. Once sufficient resolution is reached to resolve the liquid-vapour phase transition of iron in the ANEOS case, and this ring is resolved, the value of $Q_{RD}^*$ has then converged.

Magnetic clouds (MCs) are transient structures containing large-scale magnetic flux ropes from solar eruptions. The twist of magnetic field lines around the rope axis reveals information about flux rope formation processes and geoeffectivity. During propagation, MC flux ropes may erode via reconnection with the ambient solar wind. Any erosion reduces the magnetic flux and helicity of the ropes, and changes their cross-sectional twist profiles. This study relates twist profiles in MC flux ropes observed at 1 AU to the amount of erosion undergone by the MCs in interplanetary space. The twist profiles of two well-identified MC flux ropes associated with the clear appearance of post eruption arcades in the solar corona are analysed. To infer the amount of erosion, the magnetic flux content of the ropes in the solar atmosphere is estimated, and compared to estimates at 1 AU. The first MC shows a monotonically decreasing twist from the axis to periphery, while the second displays high twist at the axis, rising twist near the edges, and lower twist in between. The first MC displays a larger reduction in magnetic flux between the Sun and 1 AU, suggesting more erosion, than that seen in the second MC. In the second cloud, rising twist at the rope edges may have been due to an envelope of overlying coronal field lines with relatively high twist, formed by reconnection beneath the erupting flux rope in the low corona. This high-twist envelope remained almost intact from the Sun to 1 AU due to the low erosion levels. In contrast, the high-twist envelope of the first cloud may have been entirely peeled away via erosion by the time it reaches 1 AU.

The discovery of gravitational wave radiation from merging black holes (BHs) also uncovered BHs with masses in the range of ~20-90 Msun, which upon their merger became even more massive ones. In contrast, the most massive Galactic stellar-mass BH currently known has a mass ~21 Msun. While low-mass X-ray binaries (LMXBs) will never independently evolve into a binary BH system, and binary evolution effects can play an important role explaining the different BH masses found through studies of LMXBs, high-mass X-ray binaries, and gravitational wave events, (electromagnetic) selection effects may also play a role in this discrepancy. Assuming BH LMXBs originate in the Galactic Plane where massive stars are formed, we show that both the spatial distribution of the current sample of 20 Galactic LMXBs with dynamically confirmed BH masses, and that of candidate BH LMXBs, are both strongly biased to sources that lie at a large distance from the Galactic Plane. Specifically, most of the confirmed and candidate BH LMXBs are found at a Galactic height larger than 3 times the scale height for massive star formation. In addition, the confirmed BHs in LMXBs are found at larger distances to the Galactic Center than the candidate BH LMXBs. Interstellar absorption makes candidate BH X-ray binaries in the Plane and those in the Bulge close to the Galactic Center too faint for a dynamical mass measurement using current instrumentation. Given the observed and theoretical evidence for BH natal and/or Blaauw kicks, their relation with BH mass and binary orbital period, and the relation between outburst recurrence time and BH mass, the observational selection effects imply that the current sample of confirmed BH LMXBs is biased against the most massive BHs.

The EAS Cherenkov light array Tunka-133, with $\sim$ 3 km$^2$ geometric area, is taking data since 2009.The array permits a detailed study of energy spectrum and mass composition of cosmic rays in the energy range from $6\cdot 10^{15}$ to $10^{18}$ eV. We describe the methods of time and amplitude calibration of the array and the methods of EAS parameters reconstruction. We present the all-particle energy spectrum, based on 7 seasons of operation.

The Rasnik system is a 3-point optical displacement monitor with sub-nanometer precision. The CCD-Rasnik alignment system was developed in 1993 for the monitoring of the alignment of the muon chambers of the ATLAS Muon Spectrometer at CERN. Since then, the development has continued as new CMOS imaging pixel chips became available. The system's processes and parameters that limit the precision have been studied in detail. We conclude that only the quantum fluctuations to which the light level content of sensor pixels are subject to, is limiting the spatial resolution. The results of two Rasnik systems are compared to results from simulations, which are in good agreement: the best reached precision of $\SI{7}{pm/\sqrt{Hz}}$ is reported. Finally, some applications of high-precision Rasnik systems are set out.

We present multi-scale and multi-wavelength data of the Galactic HII region G25.4-0.14 (hereafter G25.4NW, distance ~5.7 kpc). The SHARC-II 350 micron continuum map displays a hub-filament configuration containing five parsec scale filaments and a central compact hub. Through the 5 GHz radio continuum map, four ionized clumps (i.e., Ia-Id) are identified toward the central hub, and are powered by massive OB-stars. The Herschel temperature map depicts the warm dust emission (i.e., Td ~23-39 K) toward the hub. High resolution Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum map (resolution ~0".82 X 0".58) reveals three cores (c1-c3; mass ~80-130 Msun) toward the ionized clumps Ia, and another one (c4; mass ~70 Msun) toward the ionized clump Ib. A compact near-infrared (NIR) emission feature (extent ~0.2 pc) is investigated toward the ionized clump Ia excited by an O8V-type star, and contains at least three embedded K-band stars. In the direction of the ionized clump Ia, the ALMA map also shows an elongated feature (extent ~0.2 pc) hosting the cores c1-c3. All these findings together illustrate the existence of a small cluster of massive stars in the central hub. Considering the detection of the hub-filament morphology and the spatial locations of the mm cores, a global non-isotropic collapse (GNIC) scenario appears to be applicable in G25.4NW, which includes the basic ingredients of the global hierarchical collapse and clump-fed accretion models. Overall, the GNIC scenario explains the birth of massive stars in G25.4NW.

Milli-second pulsars with highly stable periods can be considered as very precise clocks and can be used for pulsar timing array (PTA) which attempts to detect nanoheltz gravitational waves (GWs) directly. Main sources of nanoheltz GWs are supermassive black hole (SMBH) binaries which have sub-pc-scale orbits. On the other hand, a SMBH binary which is in an earlier phase and has pc-scale orbit emits ultra-low-frequency ($\lesssim 10^{-9}\,\mathrm{Hz}$) GWs cannot be detected with the conventional methodology of PTA. Such binaries tend to obtain high eccentricity, possibly $\sim 0.9$. In this paper, we develop a formalism for extending constraints on GW amplitudes from single sources obtained by PTA toward ultra-low frequencies considering the waveform expected from an eccentric SMBH binary. GWs from an eccentric binaries are contributed from higher harmonics and, therefore, have a different waveform those from a circular binary. Furthermore, we apply our formalism to several hypothetical SMBH binaries at the center of nearby galaxies, including M87, using the constraints from NANOGrav's 11-year data set. For a hypothetical SMBH binary at the center of M87, the typical upper limit on the mass ratio is $0.16$ for eccentricity of $0.9$ and semi-major axis of $a=1~\mathrm{pc}$, assuming the binary phase to be the pericenter.

In this work, we use $N-$body simulations performed on GPU to trace the past 10 Gyr dynamical history of a globular cluster (GC) similar to NGC 6397 in the tidal field of a Milky Way-like galaxy and we compare our simulated GCs with data from the third Gaia early data release. Our simulations predict, in contrast to what is deduced from the data, that such a cluster should present strong and extended tidal tails by more than 6 Gyr ago (right after the first third of its life), exceeding 1 kpc of length, and should be roughly disrupted by current time. We analyzed each of our initial conditions, such as initial mass and density parameters, as well as the dark matter shape, and we argue that the most likely reason for such discrepancy between the data and our simulations is related to the fact that we consider a purely baryonic cluster in the beginning of each model we test. We discuss that if our globular cluster was initially embedded in a dark matter minihalo, the latter could act as a protecting envelope, which prevents the tidal stripping of the luminous matter, while being itself gradually disrupted and removed in the course of the cluster evolution. This could explain why an insignificant amount of dark matter is required to describe the velocity dispersion in NGC 6397, up to at least a few half-mass radii.

Swing amplification is a model of spiral arm formation in disk galaxies. Previous $N$-body simulations show that the epicycle phases of stars in spiral arms are synchronized. However, the elementary process of the phase synchronization is not well understood. In order to investigate phase synchronization, we investigate the orbital evolution of stars due to gravitational scattering by a perturber under the epicycle approximation and its dependence on orbital elements and a disk parameter. We find that gravitational scattering by the perturber can cause phase synchronization of stellar orbits. The epicycle phases are better synchronized for smaller initial epicycle amplitudes of stars and larger shear rates of galactic disks. The vertical motion of stars does not affect the phase synchronization. The phase synchronization forms trailing dense regions, which may correspond to spiral arms.

In this paper, we reanalyze the M1.2 confined flare with a large extreme-ultraviolet (EUV) late phase on 2011 September 9, focusing on its energy partition. The radiation ($\sim$5.4$\times$10$^{30}$ erg) in 1$-$70 {\AA} is nearly eleven times larger than the radiation in 70$-$370 {\AA}, and is nearly 180 times larger than the radiation in 1$-$8 {\AA}. The peak thermal energy of the post-flare loops is estimated to be (1.7$-$1.8)$\times$10$^{30}$ erg based on a simplified schematic cartoon. Based on previous results of Enthalpy-Based Thermal Evolution of Loops (EBTEL) simulation, the energy inputs in the main flaring loops and late-phase loops are (1.5$-$3.8)$\times$10$^{29}$ erg and 7.7$\times$10$^{29}$ erg, respectively. The nonthermal energy ((1.7$-$2.2)$\times$10$^{30}$ erg) of the flare-accelerated electrons is comparable to the peak thermal energy and is sufficient to provide the energy input of the main flaring loops and late-phase loops. The magnetic free energy (9.1$\times$10$^{31}$ erg) before flare is large enough to provide the heating requirement and radiation, indicating that the magnetic free energy is adequate to power the flare.

Observations of Type II supernovae imply that a large fraction of its progenitors experience enhanced mass loss years to decades before core collapse, creating a dense circumstellar medium (CSM). Assuming that the CSM is produced by a single mass eruption event, we analytically model the density profile of the resulting CSM. We find that a double power-law profile gives a good fit to the CSM profile obtained using radiation hydrodynamical simulations. With our profile the CSM is well described by just two parameters, the transition radius $r_*$ and density at $r=r_*$ (alternatively $r_*$ and the total CSM mass). We encourage future studies to include this profile, if possible, when modelling emission from interaction-powered transients.

(abridged) Within the Orion A molecular cloud, the integral-shaped filament (ISF) is a prominent, degree-long structure of dense gas and dust, with clear signs of recent and on-going high-mass star formation. We used the ArTeMiS bolometer camera at APEX to map a 0.6x0.2 deg^2 region covering OMC-1, OMC-2, OMC-3 at 350 and 450 micron. We combined these data with Herschel-SPIRE maps to recover extended emission. The combined Herschel-ArTeMiS maps provide details on the distribution of dense, cold material, with a high spatial dynamic range, from our 8'' resolution (0.016 pc) up to the size of the map ~10-15 deg. By combining Herschel and ArTeMiS data at 160, 250, 350 and 450 micron, we constructed high-resolution temperature and H2 column density maps. We extracted radial profiles from the column density map in several, representative portions of the ISF, that we fitted with Gaussian and Plummer models to derive their intrinsic widths. We also compared the distribution of material traced by ArTeMiS with that seen in the higher density tracer N2H+(1-0) recently observed with the ALMA interferometer. All the radial profiles that we extracted show clear deviation from a Gaussian, with evidence for an inner plateau, previously not seen using Herschel-only data. We measure intrinsic half-power widths in the range 0.06 to 0.11 pc. This is significantly larger than the Gaussian widths measured for fibers seen in N2H+, which probably traces only the dense innermost regions of the large-scale filament. These half-power widths are within a factor of two of the value of 0.1 pc found for a large sample of nearby filaments in various low-mass star-forming regions, which tends to indicate that the physical conditions governing the fragmentation of prestellar cores within transcritical or supercritical filaments are the same over a large range of masses per unit length.

We present the photometric and spectroscopic evolution of supernova (SN) 2019cad during the first $\sim100$ days from explosion. Based on the light curve morphology, we find that SN 2019cad resembles the double-peaked type Ib/c SN 2005bf and the type Ic PTF11mnb. Unlike those two objects, SN 2019cad also shows the initial peak in the redder bands. Inspection of the g-band light curve indicates the initial peak is reached in $\sim8$ days, while the r band peak occurred $\sim15$ days post-explosion. A second and more prominent peak is reached in all bands at $\sim45$ days past explosion, followed by and fast decline from $\sim60$ days. During the first 30 days, the spectra of SN 2019cad show the typical features of a type Ic SN, however, after 40 days, a blue continuum with prominent lines of Si II ${\lambda}6355$ and C II ${\lambda}6580$ is observed again. Comparing the bolometric light curve to hydrodynamical models, we find that SN 2019cad is consistent with a pre-SN mass of 11 M$_{\odot}$, and an explosion energy of $3.5\times 10^{51}$ erg. The light curve morphology can be reproduced either by a double-peaked $^{56}$Ni distribution with an external component of 0.041 M$_{\odot}$ and an internal component of 0.3 M$_{\odot}$ or a double-peaked $^{56}$Ni distribution plus magnetar model (P $\sim11$ ms and B $\sim26\times 10^{14}$ G). If SN 2019cad were to suffer from significant host reddening (which cannot be ruled out), the $^{56}$Ni model would require extreme values, while the magnetar model would still be feasible.

The Tibet ASgamma experiment just reported their measurement of sub-PeV diffuse gamma ray emission from the Galactic disk, with the highest energy up to 957 TeV. These gamma-rays are most likely the hadronic origin by cosmic ray interaction with interstellar gas in the Galaxy. This measurement provides direct evidence to the hypothesis that the Galactic cosmic rays can be accelerated beyond PeV energies. In this work, we try to explain the sub-PeV diffuse gamma-ray spectrum within cosmic rays diffusive propagation model. We find there is a tension between the sub-PeV diffuse gamma rays and the local cosmic ray spectrum. To describe the sub-PeV diffuse gamma-ray flux, it generally requires larger local cosmic-ray flux than measurement in the knee region. We further calculate the PeV neutrino flux from the cosmic ray propagation model. Even all of these sub-PeV diffuse gamma rays originate from the propagation, the Galactic neutrinos only account for less than ~15% of observed flux, most of which are still from extragalactic sources.

Euler springs are used for vertical suspension and vibration isolation as they provide a large static supporting force with low spring-rate and use minimal spring material. To date multiple single-width, vertically-stacked, rectangular blades of uniform thickness have been used in the post buckled state, with half of the number buckling in each of opposing directions. This structure requires clamping the ends of the blades giving stick-slip problems. In this study we investigate the benefits of forming multiple, oppositely buckling blades, side- by-side from a single monolithic sheet of spring material. We investigate how to distribute the stress evenly along the length of the blade by contouring its width, as well as finding the optimal joining contour to distribute the stress evenly around the tearing joints between oppositely bending blades.

We present here the results of broadband spectral analysis of low-mass X-ray binary and a black hole candidate 4U 1957+11. The source was observed nine times with the Nuclear Spectroscopic Telescope Array (NuSTAR) between 2018 September and 2019 November. During these observations, the spectral state of 4U 1957+11 evolved marginally. The disc dominant spectra are well described with a hot, multicolour disc blackbody with disc temperature varying in the range $kT_{\rm in} \sim$ 1.35--1.86 keV and a non-thermal component having a steep slope ($\Gamma =$ 2--3). A broad Fe emission line feature (5--8 keV) was observed in the spectra of all the observations. The relativistic disc model was used to study the effect of distance, inclination, and the black hole mass on its spin. Simulations indicate a higher spin for smaller distances and lower black hole mass. At smaller distances and higher mass, spin is maximum and almost independent of the distance. An inverse correlation exists between the spin and the spectral hardening factor for all the cases. The system prefers a moderate spin of about 0.85 for black hole masses between 4--6 M_sun for a 7 kpc distance.

Context. While the shapes of many observed bow shocks can be reproduced by simple astrosphere models, more elaborate approaches have recently been used to explain differing observable structures. Aims. By placing perturbations of an otherwise homogeneous interstellar medium in front of the astrospheric bow shock of the runaway blue supergiant $\lambda$ Cephei, the observable structure of the model astrosphere is significantly altered, providing insight into the origin of perturbed bow shock images. Methods. Three-dimensional single-fluid magnetohydrodynamic (MHD) models of stationary astrospheres were subjected to various types of perturbations and simulated until stationarity was reached again. As examples, simple perturbations of the available MHD parameters (number density, bulk velocity, temperature, and magnetic field) as well as a more complex perturbation were chosen. Synthetic observations were generated by line-of-sight integration of the model data, producing H$\alpha$, $70\,\mu$m dust emission, and bremsstrahlung maps of the perturbed astrosphere's evolution. Results. The resulting shock structures and observational images differ strongly depending on the type of the injected perturbation and the viewing angles, forming arc-like protrusions or bifurcations of the bow shock structure, as well as rings, arcs, and irregular structures detached from the bow shock.

We report results of proper motions of 25 known Small Magellanic Cloud (SMC) clusters (ages ~ 1 - 10 Gyr old) derived from Gaia EDR3 data sets. When these mean proper motions are gathered with existent radial velocity measurements to compose the clusters' velocity vectors, we found the parameter values of a rotation disk that best resemble their observed motions, namely: central coordinates and distance, inclination and position angle of the line-of-node, proper motion in right ascension and declination and systemic velocity, rotation velocity and velocity dispersion. The SMC cluster rotation disk seems to be at some level kinematically synchronized with the rotation of field red giants recently modeled using DR2 data sets. Such a rotation disk is seen in the sky as a tilted edge-on disk, with a velocity dispersion perpendicular to it twice as big as that in the plane of the disk. Because the direction perpendicular to the disk is nearly aligned with the Magellanic Bridge, we interpret the larger velocity dispersion as a consequence of the SMC velocity stretching caused by the tidal interaction with the Large Magellanic Cloud. Rotation alone would not seem sufficient to explain the observed kinematic behaviors in the SMC.

The Astrometric Gravitation Probe mission is a modern version of the 1919 Dyson-Eddington-Davidson experiment, based on a space-borne telescope with a permanent built-in eclipse, provided by a coronagraphic system. The expected improvement on experimental bounds to General Relativity and competing gravitation theories is by at least two orders of magnitude. The measurement principle is reviewed, in particular the principle of Fizeau-like combination of a set of individual inverted coronagraphs simultaneously feeding a common high resolution telescope. Also, the payload has a dual field of view property, in order to support simultaneous observations of stellar fields either very close, or far away, from the Sun, i.e. fields affected by either high or low light bending. We discuss a set of solutions introduced in the optical design to improve on technical feasibility and robustness of the optical performance against perturbations, in particular induced by manufacturing and alignment tolerances, and launch stresses.

Predicting how a young planet shapes the gas and dust emission of its parent disc is key to constraining the presence of unseen planets in protoplanetary disc observations. We investigate the case of a 2 Jupiter mass planet that becomes eccentric after migrating into a low-density gas cavity in its parent disc. Two-dimensional hydrodynamical simulations are performed and post-processed by three-dimensional radiative transfer calculations. In our disc model, the planet eccentricity reaches ~0.25, which induces strong asymmetries in the gas density inside the cavity. These asymmetries are enhanced by photodissociation and form large-scale asymmetries in 12CO J=3-2 integrated intensity maps. They are shown to be detectable for an angular resolution and a noise level similar to those achieved in ALMA observations. Furthermore, the planet eccentricity renders the gas inside the cavity eccentric, which manifests as a narrowing, stretching and twisting of iso-velocity contours in velocity maps of 12CO J=3-2. The planet eccentricity does not, however, give rise to detectable signatures in 13CO and C18O J=3-2 inside the cavity because of low column densities. Outside the cavity, the gas maintains near-circular orbits, and the vertically extended optically thick CO emission displays a four-lobed pattern in integrated intensity maps for disc inclinations > 30 degrees. The lack of large and small dust inside the cavity in our model further implies that synthetic images of the continuum emission in the sub-millimetre, and of polarized scattered light in the near-infrared, do not show significant differences when the planet is eccentric or still circular inside the cavity.

In this research note, we outline the extension of the grid of upper atmosphere models first presented in Kubyshkina et al. (2018a). The original grid is based on a 1D hydrodynamic model and consists of about 7000 models covering planets of the size from Earth to twice Neptune at orbits corresponding to equilibrium temperatures between 300 and 2000 K around solar-like (0.4 to 1.3 solar mass) stars. The extended and revised grid of models consists of 10235 points and covers a planetary mass range of up to 109 Earth masses, which allows one to outline the transition between low- and high-gravity hot planets in short orbital separations. We prepared the interpolation tool allowing one to use the grid to define the mass-loss of a planet that falls into the parameter range of the grid. We provide a comparison of our results to common analytical models.

We study a scenario by which a giant wide tertiary star engulfs and forces a tight binary system of a white dwarf (WD) and a main sequence (MS) star to enter a common envelope evolution (CEE) with each other, and then unbinds the WD-MS common envelope. The WD-MS binary system, now with the WD inside the MS envelope, does not have sufficient orbital energy to unbind their common envelope. However, as they approach the center of the giant star Roche lobe overflow to the core of the giant star and/or merger of the WD with the core remove a large fraction of the WD-MS common envelope or all of it. Namely, the energy source for unbinding the WD-MS tight common envelope is the triple-star CEE. For that we term this scenario a parasite CEE. Overall, the destruction of the MS star absorbs energy from the triple-star system, a process that might lead to WD-core merger during the triple-star CEE. The parasite CEE leaves behind either one massive WD that in some cases might explode as a peculiar type Ia supernova or two close WDs that at later time might explode as a type Ia supernova. We very crudely estimate the rate of the parasite CEE to be a fraction of about 0.001 out of all evolved triple stars.

We report deep imaging observations with DOLoRes@TNG of an ultra-faint dwarf satellite candidate of the Triangulum galaxy (M33) found by visual inspection of the public imaging data release of the Dark Energy Camera Legacy Survey. Pisces VII/Triangulum (Tri) III is found at a projected distance of 72 kpc from M33, and using the tip of the red giant branch method we find a distance to this faint system of D=820^(+195, -190) kpc. We estimate an absolute magnitude of M_V=-4.4^(+0.8,-0.7) and a half-light radius of r_half =100 +/-14 pc for the galaxy, consistent with similarly faint galaxies around the Milky Way. As the tip of the red giant branch is sparsely populated, constraining a precision distance is difficult, but if Pisces VII/Tri III can be confirmed as a true satellite of M33, it is a significant finding. Firstly, it would be the faintest dwarf galaxy detected to-date outside of the Milky Way. With only one potential satellite detected around M33 previously (Andromeda XXII/Tri I), it lacks a significant satellite population in stark contrast to the similarly massive Large Magellanic Cloud. The detection of more satellites in the outskirts of M33 could help to better illuminate if this discrepancy between expectation and observations is due to a poor understanding of the galaxy formation process, or if it is due to the low luminosity and surface brightness of the M33 satellite population which has thus far fallen below the detection limits of previous surveys.

I consider a sample of eight pressure-supported low-surface brightness galaxies (seven nearby dwarfs and one ultra-diffuse object) in terms of Milgrom's modified Newtonian dynamics (MOND). These objects are modelled as Milgromian isotropic isothermal spheres characterised by two parameters that are constrained by observations: the constant line-of-sight velocity dispersion and the central surface density. The velocity dispersion determines the total mass, and, with the implied mass-to-light ratio, the central surface brightness. This then specifies the radial run of surface brightness over the entire isothermal sphere. For the objects in this sample the predicted radial distribution of surface brightness is shown to be entirely consistent with observations which constitutes a success for MOND that is independent of the reduced dynamical mass.

HD49798 / RXJ0648.0-4418 is the only confirmed X-ray binary in which the mass donor is a hot subdwarf star of O spectral type and, most likely, it contains a massive white dwarf (1.28$\pm$0.05 M$_{\rm SUN}$) with a very fast spin period of 13.2 s. Here we report the results of new XMM-Newton pointings of this peculiar binary, carried out in 2018 and in 2020, together with a reanalysis of all the previous observations. The new data indicate that the compact object is still spinning-up at a steady rate of $(-2.17\pm0.01)\times10^{-15}$ s s$^{-1}$, consistent with its interpretation in terms of a young contracting white dwarf. Comparison of observations obtained at similar orbital phases, far from the ecplise, shows evidence for long term variability of the hard ($>$0.5 keV) spectral component at a level of $\sim$(70$\pm$20)\%, suggesting the presence of time-dependent inhomogeneities in the weak stellar wind of the HD49798 subdwarf. To investigate better the soft spectral component that dominates the X-ray flux from this system, we computed a theoretical model for the thermal emission expected from an atmosphere with element abundances and surface gravity appropriate for this massive white dwarf. This model gives a best fit with effective temperature of T$_{\rm eff}$=2.25$\times$10$^5$ K and an emitting area with radius of $\sim$1600 km, larger than that found with blackbody fits. This model also predicts a contribution of the pulsed emission from the white dwarf in the optical band significantly larger than previously thought and possibly relevant for optical variability studies of this system.

The core mass of galaxy clusters is both an important anchor of the radial mass distribution profile and probe of structure formation. With thousands of strong lensing galaxy clusters being discovered by current and upcoming surveys, timely, efficient, and accurate core mass estimates are needed. We assess the results of two efficient methods to estimate the core mass of strong lensing clusters: the mass enclosed by the Einstein radius ($M_{corr}(<\theta_E)$; Remolina Gonz\'{a}lez et al. 2020), and single-halo lens model ($M_{\rm{SHM}}(<\rm{e}\theta_{\rm{E}})$; Remolina Gonz\'{a}lez et al. 2021), against measurements from publicly available detailed lens models ($M_{\rm{DLM}}$) of the same clusters. We use publicly available lens models from the Sloan Giant Arc Survey, the Reionization Lensing Cluster Survey, the \Hubble\ Frontier Fields, and the Cluster Lensing and Supernova Survey with \Hubble. We find a scatter of $18.3\%$ ($8.4\%$) with a bias of $-7.5\%$ ($0.4\%$) between $M_{corr}(<\theta_E)$ ($M_{\rm{SHM}}(<\rm{e}\theta_{\rm{E}})$) and $M_{\rm{DLM}}$. Last, we compare the statistical uncertainties measured in this work to those from simulations. This work demonstrates the successful application of these methods to observational data. As the effort to efficiently model the mass distribution of strong lensing galaxy clusters continues, we are in need of fast and reliable methods to advance the field.

With a metallicity of 12 + Log(O/H) $\approx$ 7.1-7.2, I Zw 18 is a canonical low-metallicity blue compact dwarf (BCD) galaxy. A growing number of BCDs, including I Zw 18, have been found to host strong, narrow-lined, nebular He II ($\lambda$4686) emission with enhanced intensities compared to H$\beta$ (e.g., He II($\lambda$4686)/H$\beta$ > 1%). We present new observations of I Zw 18 using the Keck Cosmic Web Imager. These observations reveal two nebular He II emission regions (or He III regions) northwest and southeast of the He III region in the galaxy's main body investigated in previous studies. All regions exhibit He II($\lambda4686$)/Hbeta greater than 2%. The two newly resolved He III regions lie along an axis that intercepts the position of I Zw 18's ultraluminous X-ray (ULX) source. We explore whether the ULX could power the two He III regions via shock activity and/or beamed X-ray emission. We find no evidence of shocks from the gas kinematics. If the ULX powers the two regions, the X-ray emission would need to be beamed. Another potential explanation is that a class of early-type nitrogen-rich Wolf-Rayet stars with low winds could power the two He III regions, in which case the alignment with the ULX would be coincidental.

Gravitational wave science is a pioneering field with rapidly evolving data analysis methodology currently assimilating and inventing deep learning techniques. The bulk of the sophisticated flagship searches of the field rely on the time-tested matched filtering principle within their core. In this paper, we make a key observation on the relationship between the emerging deep learning and the traditional techniques: matched filtering is formally equivalent to a particular neural network. This means that a neural network can be constructed analytically to exactly implement matched filtering, and can be further trained on data or boosted with additional complexity for improved performance. This fundamental equivalence allows us to define a "complexity standard candle" allowing us to characterize the relative complexity of the different approaches to gravitational wave signals in a common framework. Additionally it also provides a glimpse of an intriguing symmetry that could provide clues on how neural networks approach the problem of finding signals in overwhelming noise. Moreover, we show that the proposed neural network architecture can outperform matched filtering, both with or without knowledge of a prior on the parameter distribution. When a prior is given, the proposed neural network can approach the statistically optimal performance. We also propose and investigate two different neural network architectures MNet-Shallow and MNet-Deep, both of which implement matched filtering at initialization and can be trained on data. MNet-Shallow has simpler structure, while MNet-Deep is more flexible and can deal with a wider range of distributions. Our theoretical findings are corroborated by experiments using real LIGO data and synthetic injections. Finally, our results suggest new perspectives on the role of deep learning in gravitational wave detection.

The discovery of the Higgs boson has opened a new avenue for exploring physics beyond the Standard Model. In this review, we discuss cosmological aspects of the simplest Higgs couplings to the hidden sector, known as the Higgs portal. We focus on implications of such couplings for inflation, vacuum stability and dark matter, with the latter including both the traditional weakly interacting massive particles as well as feebly interacting scalars. The cosmological impact of the Higgs portal can be important even for tiny values of the couplings.

Are some cosmologists trying to return human beings to the centre of the cosmos? In the view of some critics, the so-called "anthropic principle" is a desperate attempt to salvage a scrap of dignity for our species after a few centuries of demotion at the hands of science. It is all things archaic and backwards - teleology, theology, religion, anthropocentrism - trying to sneak back in scientific camouflage. We argue that this is a mistake. The anthropic principle is not mere human arrogance, nor is it religion in disguise. It is a necessary part of the science of the universe.

We study the possibility of probing new physics accounting for $(g-2)_\mu$ anomaly and gravitational waves with pulsar timing array measurements. The model we consider is either a light gauge boson or neutral scalar interacting with muons. We show that the parameter spaces of dark $U(1)$ model with kinetic mixing explaining $(g-2)_\mu$ anomaly can realize a first-order phase transition, and the yield-produced gravitational wave may address the common red noise observed in the NANOGrav 12.5-yr dataset.

It is a long-standing debate as to whether or not the annual modulation in the event rate observed by the DAMA sodium iodide experiment is caused by the interaction of dark matter particles. To resolve this issue, several groups have been working to develop new experiments with the aim of reproducing or refuting DAMA's results using the same sodium iodide target medium. The COSINE-100 experiment is one of these that is currently operating with 106 kg of low-background sodium iodide crystals at the Yangyang underground laboratory. Analysis of the initial 59.5 days of COSINE-100 data showed that the annual modulation signal reported by DAMA is inconsistent with explanation using spin-independent interaction of weakly interacting massive particles (WIMPs), a favored candidate of dark matter particles, with sodium or iodine nuclei in the context of the standard halo model. However, this first result left open interpretations using certain alternative dark matter models, dark matter halo distributions, and detector responses that could allow room for consistency between DAMA and COSINE-100. Here we present new results from over 1.7 years of COSINE-100 operation with improved event selection and energy threshold reduced from 2 keV to 1 keV. We find an order of magnitude improvement in sensitivity, sufficient for the first time to strongly constrain these alternative scenarios, as well as to further strengthen the previously observed inconsistency with the WIMP-nucleon spin-independent interaction hypothesis.

Solar Orbiter was launched on February 10, 2020 with the purpose of investigating solar and heliospheric physics using a payload of instruments designed for both remote and in-situ sensing. Similar to the recently launched Parker Solar Probe, and unlike earlier missions, Solar Orbiter carries instruments designed to measure the low frequency DC electric fields. In this paper we assess the quality of the low-frequency DC electric field measured by the Radio and Plasma Waves instrument (RPW) on Solar Orbiter. In particular we investigate the possibility of using Solar Orbiter's DC electric and magnetic field data to estimate the solar wind speed. We use deHoffmann-Teller (HT) analysis based on measurements of the electric and magnetic fields to find the velocity of solar wind current sheets which minimizes a single component of the electric field. By comparing the HT velocity to proton velocity measured by the Proton and Alpha particle Sensor (PAS) we develop a simple model for the effective antenna length, $L_\text{eff}$ of the E-field probes. We then use the HT method to estimate the speed of the solar wind. Using the HT method, we find that the observed variations in $E_y$ are often in excellent agreement with the variations in the magnetic field. The magnitude of $E_y$, however, is uncertain due to the fact that the $L_\text{eff}$ depends on the plasma environment. We derive an empirical model relating $L_\text{eff}$ to the Debye length, which we can use to improve the estimate of $E_y$ and consequently the estimated solar wind speed. The low frequency electric field provided by RPW is of high quality. Using deHoffmann-Teller analysis, Solar Orbiter's magnetic and electric field measurements can be used to estimate the solar wind speed when plasma data is unavailable.

It has been argued that supergravity models of inflation with vanishing sound speeds, $c_s$, lead to an unbounded growth in the production rate of gravitinos. We consider several models of inflation to delineate the conditions for which $c_s = 0$. In models with unconstrained superfields, we argue that the mixing of the goldstino and inflatino in a time-varying background prevents the uncontrolled production of the longitudinal modes. This conclusion is unchanged if there is a nilpotent field associated with supersymmetry breaking with constraint ${\bf S^2} =0$, i.e. sgoldstino-less models. Models with a second orthogonal constraint, ${\bf S(\Phi-\bar{\Phi})} =0$, where $\bf{\Phi}$ is the inflaton superfield, which eliminates the inflatino, may suffer from the over-production of gravitinos. However, we point out that these models may be problematic if this constraint originates from a UV Lagrangian, as this may require using higher derivative operators. These models may also exhibit other pathologies such as $c_s > 1$, which are absent in theories with the single constraint or unconstrained fields.

The effect of energy absorption during the binary evolution of Exotic-compact-objects (ECOs) is extensively studied. We review the underlying mechanism that provides the energy dissipation in material objects - tidal friction. We show that unlike typical astrophysical objects, where absorption due to viscosity is negligible, in ECOs, absorption could potentially mimic the analogous effect of black-holes (BHs) - tidal heating. We stand for their differences and similarities in the context of energy dissipation during the inspiral. Inspired by the membrane paradigm and recent studies, we demonstrate how viscosity is a defining feature that quantifies how close is the ECO absorption to that of a classical BH absorption. We show that for ECOs, viscosity can induce significant modifications to the GW waveform, which in some favorable scenarios of super-massive binaries of equal mass and spin, enables the measurement of the ECO absorption in the future precision gravitational-wave (GW) observations. Finally, we discuss the implications on the ECO reflection coefficient and the relation to the universal viscosity to volume entropy bound.

The effective approach is applied to the analysis of inflationary magnetogenesis. Rather than assuming a particular underlying description, all the generally covariant terms potentially appearing with four space-time derivatives in the effective action have been included and weighted by inflaton-dependent couplings. The higher derivatives are suppressed by the negative powers of a typical mass scale whose specific values ultimately depend on the tensor to scalar ratio. During a quasi-de Sitter stage the corresponding corrections always lead to an asymmetry between the hypermagnetic and the hyperelectric susceptibilities. After presenting a general method for the estimate of the gauge power spectra, the obtained results are illustrated for generic models and also in the case of some non-generic scenarios where either the inflaton has some extra symmetry or the higher-order terms are potentially dominant.