Papers for Thursday, Mar 17 2022

The environment of actively repeating fast radio bursts (FRBs) has been shown to be complex and varying. The recently localized FRB 20190520B is extremely active, has the largest confirmed host dispersion measure, and is only the second FRB source associated with a compact, persistent radio source (PRS). The main tracer of the magneto-ionic environments is the rotation measure (RM), a path-integral of the line-of-sight component of magnetic field strength (B) and electron density, which does not allow a direct probe of the B-field configuration. Here we report direct evidence for a B-field reversal based on the observed sign change and extreme variation of FRB 20190520B's RM, which changed from $\sim10000$ rad m$^{-2}$ to $\sim-16000$ rad m$^{-2}$ between June 2021 and January 2022. Such extreme RM reversal has never been observed before in any FRB nor in any astronomical object. The implied short-term change of the B-field configuration in or around the FRB could be due to the vicinity of massive black holes, or a magnetized companion star in binary systems, or a young supernova remnant along the line of sight.

We introduce a new estimator of the CMB lensing power spectrum, together with its likelihood, based on iterative lensing reconstruction. Despite the increased complexity of the lensing maps, this estimator shares similarities with the standard quadratic estimator. Most importantly, it is unbiased towards the assumptions done on the noise and cosmology for the lensing reconstruction. This new spectrum estimator can double the constraints on the lensing amplitude compared to the quadratic estimator, while keeping numerical cost under control and being robust to errors.

The open star cluster NGC 1605 has recently been reported to in fact consist of two clusters (one intermediate-aged and one old) that merged via a flyby capture. Here we show that Gaia data do not support this scenario. We also report the serendipitous discovery of a new open cluster, Can Batll\'o 1, with a similar age and distance.

We present high-cadence optical, ultraviolet (UV), and near-infrared data of the nearby ($D\approx23$ Mpc) Type II supernova (SN) 2021yja. Many Type II SNe show signs of interaction with circumstellar material (CSM) during the first few days after explosion, implying that their red supergiant (RSG) progenitors experience episodic or eruptive mass loss. However, because it is difficult to discover SNe early, the diversity of CSM configurations in RSGs has not been fully mapped. SN 2021yja, first detected within ${\approx}5.4$ hours of explosion, shows some signatures of CSM interaction (high UV luminosity, radio and x-ray emission) but without the narrow emission lines or early light curve peak that can accompany CSM. Here we analyze the densely sampled early light curve and spectral series of this nearby SN to infer the properties of its progenitor and CSM. We find that the most likely progenitor was an RSG with an extended envelope, encompassed by low-density CSM. We also present archival Hubble Space Telescope imaging of the host galaxy of SN 2021yja, which allows us to place a stringent upper limit of ${\lesssim}9\ M_\odot$ on the progenitor mass. However, this is in tension with some aspects of the SN evolution, which point to a more massive progenitor. Our analysis highlights the need to consider progenitor structure when making inferences about CSM properties, and that a comprehensive view of CSM tracers should be made to give a fuller view of the last years of RSG evolution.

We use zoom simulations to show how merger-driven disruption of the gas disc in a galaxy provides its central active galactic nucleus (AGN) with fuel to drive outflows that entrain and expel a significant fraction of the circumgalactic medium (CGM). This in turn suppresses replenishment of the interstellar medium, causing the galaxy to quench up to several Gyr after the merger. We start by performing a zoom simulation of a present-day star-forming disc galaxy with the EAGLE galaxy formation model. Then, we re-simulate the galaxy with controlled changes to its initial conditions, using the genetic modification technique. These modifications either increase or decrease the stellar mass ratio of the galaxy's last significant merger, which occurs at $z\approx 0.74$. The halo reaches the same present-day mass in all cases, but changing the mass ratio of the merger yields markedly different galaxy and CGM properties. We find that a merger can unlock rapid growth of the central supermassive black hole if it disrupts the co-rotational motion of gas in the black hole's vicinity. Conversely, if a less disruptive merger occurs and gas close to the black hole is not disturbed, the AGN does not strongly affect the CGM, and consequently the galaxy continues to form stars. Our result illustrates how a unified view of AGN feedback, the baryon cycle and the interstellar medium is required to understand how mergers and quenching are connected over long timescales.

We model the strong lensing effect in the galaxy cluster PSZ1 G311.65-18.48 (z=0.443) with an improved version of the hybrid method WSLAP+. We extend the number of constraints by including the position of critical points, which are combined with the classic positional constraints of the lensed galaxies. We pay special attention to a transient candidate source (Tr) previously discovered in the giant Sunburst arc (z=2.37). Our lens model predicts Tr to be within a fraction of an arcsecond from the critical curve, having a larger magnification factor than previously found, but still not large enough to explain the observed flux and lack of counterimages. Possible candidate counterimages are discussed that would lower the magnification required to explain Tr, but extreme magnification factors ($\mu>1000$) are still required, even in that case. The presence of a small mass perturber with a mass comparable to a dwarf galaxy ($M\sim 10^8 \,{\rm M}_{\odot}$) near the position of Tr is needed in order to explain the required magnification and morphology of the lensed galaxy. We discuss how the existence of this perturber could potentially be used to constrain models of dark matter. The large apparent brightness and unresolved nature of the magnified object implies a combination of extreme magnification and a very luminous and compact source ($r<0.3$ pc). Possible candidates are discussed, including an hyperluminous star or an accretion disc around an intermediate-mass black hole (IMBH). Based on spectral information, we argue that a luminous blue variable (LBV) star caught during an outburst is the most likely candidate. Owing to the extreme magnification and luminosity of this source we dub it Godzilla.

The abundance of dark matter subhalos orbiting a host galaxy is a generic prediction of the cosmological framework. It is a promising way to constrain the nature of dark matter. Here we describe the challenges of detecting stars whose phase-space distribution may be perturbed by the passage of dark matter subhalos using a machine learning approach. The training data are three Milky Way-like galaxies and nine synthetic Gaia DR2 surveys derived from these. We first quantify the magnitude of the perturbations in the simulated galaxies using an anomaly detection algorithm. We also estimate the feasibility of this approach in the Gaia DR2-like catalogues by comparing the anomaly detection based approach with a supervised classification. We find that a classification algorithm optimized on about half a billion synthetic star observables exhibits mild but nonzero sensitivity. This classification-based approach is not sufficiently sensitive to pinpoint the exact locations of subhalos in the simulation, as would be expected from the very limited number of subhalos in the detectable region. The enormous size of the Gaia dataset motivates the further development of scalable and accurate computational methods that could be used to select potential regions of interest for dark matter searches to ultimately constrain the Milky Way's subhalo mass function.

Periodic nuclear transients have been detected with increasing frequency, with one such system -- ASASSN-14ko -- exhibiting highly regular outbursts on a timescale of $114 \pm 1$ days. It has been postulated that the outbursts from this source are generated by the repeated partial disruption of a star, but how the star was placed onto such a tightly bound orbit about the supermassive black hole remains unclear. Here we use analytic arguments and three-body integrations to demonstrate that the Hills mechanism, where a binary system is destroyed by the tides of the black hole, can lead to the capture of a star on a $\sim 114$ day orbit and with a pericenter distance that is comparable to the tidal radius of one of the stars within the binary. Thus, Hills capture can produce stars on tightly bound orbits that undergo repeated partial disruption, leading to a viable mechanism for generating not only the outbursts detected from ASASSN-14ko, but for periodic nuclear transients in general. We also show that the rate of change of the period of the captured star due to gravitational-wave emission is likely too small to produce the observed value for ASASSN-14ko, indicating that in this system there must be additional effects that contribute to the decay of the orbit. In general, however, gravitational-wave emission can be important for limiting the lifetimes of these systems, and could produce observable period decay rates in future events.

Dynamical interactions in dense star clusters are considered one of the most effective formation channels of binary black holes (BBHs). Here, we present direct $N-$body simulations of two different star cluster families: low-mass ($\sim{500-800}$ M$_\odot$) and high-mass star clusters ($\ge{5000}$ M$_\odot$). We show that the formation channels of BBHs in low- and high-mass star clusters are extremely different and lead to two completely distinct populations of BBH mergers. Low-mass clusters host mainly low-mass BBHs born from binary evolution, while BBHs in high-mass clusters are relatively massive (chirp mass up to $\sim{50}$ M$_\odot$) and driven by dynamical exchanges. Tidal disruption dramatically quenches the formation and dynamical evolution of BBHs in low-mass clusters on a very short timescale ($\lesssim{100}$ Myr), while BBHs in high-mass clusters undergo effective dynamical hardening until the end of our simulations (1.5 Gyr). In high-mass clusters we find 4% of BBHs with primary mass in the pair-instability mass gap, all of them born via stellar collisions, while none of them forms in low-mass clusters. These differences are crucial for the interpretation of the formation channels of gravitational-wave sources.

The prototype Schwarzschild-Couder Telescope (pSCT) is a candidate for a medium-sized telescope in the Cherenkov Telescope Array. The pSCT is based on a novel dual mirror optics design which reduces the plate scale and allows for the use of silicon photomultipliers as photodetectors. The prototype pSCT camera currently has only the central sector instrumented with 25 camera modules (1600 pixels), providing a 2.68$^{\circ}$ field of view (FoV). The camera electronics are based on custom TARGET (TeV array readout with GSa/s sampling and event trigger) application specific integrated circuits. Field programmable gate arrays sample incoming signals at a gigasample per second. A single backplane provides camera-wide triggers. An upgrade of the pSCT camera is in progress, which will fully populate the focal plane. This will increase the number of pixels to 11,328, the number of backplanes to 9, and the FoV to 8.04$^{\circ}$. Here we give a detailed description of the pSCT camera, including the basic concept, mechanical design, detectors, electronics, current status and first light.

Magnetic reconnection in the deep solar atmosphere can give rise to enhanced emission in the Balmer hydrogen lines, a phenomenon referred to as Ellerman bombs. Recent high quality H$\beta$ observations indicate that Ellerman bombs are more common than previously thought and it was estimated that at any time about half a million Ellerman bombs are present in the quiet Sun. We performed an extensive statistical characterization of the quiet Sun Ellerman bombs (QSEBs) in these new H$\beta$ observations. We analyzed a 1 h dataset of quiet Sun observed with the Swedish 1-m Solar Telescope that consists of spectral imaging in the H$\beta$ and H$\alpha$ lines, as well as spectropolarimetric imaging in Fe I 617.3 nm. We detected a total of 2809 QSEBs. The lifetime varies between 9 s and 20.5 min with a median of 1.14 min. The maximum area ranges between 0.0016 and 0.2603 Mm$^2$ with a median of 0.018 Mm$^2$. A subset (14%) of the QSEBs display enhancement of the H$\beta$ line core. On average, the line core brightening appears 0.88 min after the onset of brightening in the wings, and the distance between these brightenings is 243 km. This gives rise to an apparent propagation speed ranging between $-$14.3 and +23.5 km s$^{-1}$, with an average that is upward propagating at +4.4 km $^{-1}$. The average orientation is nearly parallel to the limbward direction. QSEBs are nearly uniformly distributed over the field of view but we find empty areas with the size of mesogranulation. QSEBs are located more frequent near the magnetic network where they are often bigger, longer lived and brighter. We conclude that QSEBs are ubiquitous in quiet Sun and appear everywhere except in areas of mesogranular size with weakest magnetic field ($B_{\rm{LOS}}\lesssim50$~G). Our observations support the interpretation of reconnection along vertically extended current sheets.

We present a spectroscopic analysis of a 1-year intensive monitoring campaign of the 22-Myr old planet-hosting M dwarf AU Mic using the HARPS spectrograph. In a companion paper, we reported detections of the planet radial velocity (RV) signatures of the two close-in transiting planets of the system, with respective semi-amplitudes of 5.8 $\pm$ 2.5 m/s and 8.5 $\pm$ 2.5 m/s for AU Mic b and AU Mic c. Here, we perform an independent measurement of the RV semi-amplitude of AU Mic c using Doppler imaging to simultaneously model the activity-induced distortions and the planet-induced shifts in the line profiles. The resulting semi-amplitude of 13.3 $\pm$ 4.1 m/s for AU Mic c reinforces the idea that the planet features a surprisingly large inner density, in tension with current standard models of core accretion. Our brightness maps feature significantly higher spot coverage and lower level of differential rotation than the brightness maps obtained in late 2019 with the SPIRou spectropolarimeter, suggesting that the stellar magnetic activity has evolved dramatically over a $\sim$1-yr time span. Additionally, we report a 3-$\sigma$ detection of a modulation at 8.33 $\pm$ 0.04 d of the He I D3 (587.562 nm) emission flux, close to the 8.46-d orbital period of AU Mic b. The power of this emission (a few 10$^{17}$ W) is consistent with 3D magnetohydrodynamical simulations of the interaction between stellar wind and the close-in planet if the latter hosts a magnetic field of $\sim$10 G. Spectropolarimetric observations of the star are needed to firmly elucidate the origin of the observed chromospheric variability.

Massive stars lose a large fraction of their mass to radiation-driven winds throughout their entire life. These outflows impact both the life and death of these stars and their surroundings. Theoretical mass-loss rates of hot, massive stars are derived to be used in applications such as stellar evolution. The behaviour of these rates in the OB-star regime is analysed, and their effects on massive-star evolution predictions is studied. Dynamically-consistent models are computed by solving the spherically symmetric, steady-state equation-of-motion for a large grid of hot, massive stars. The radiative acceleration is derived from non-local thermodynamic equilibrium radiative transfer in the co-moving frame. The resulting mass-loss rates are used to derive a simple scaling recipe with stellar parameters and to evaluate some first impacts upon evolution tracks. We provide a new prescription for steady-state, radiation-driven mass-loss from hot, massive stars depending on their fundamental parameters. The rates for O-stars are lower by about a factor ~3 than the rates typically used in previous stellar-evolution calculations, where differences generally decrease with increasing luminosity and temperature. For cooler B-giants/supergiants we find larger discrepancies. This arises because we do not find any systematic increase in mass-loss rates below the so-called bi-stability region; indeed, our results do not show any sign of a significant bi-stability jump within the parameter range covered by the grid. Due to the lower mass-loss rates we find that envelopes are not easily stripped by means of standard steady-state winds, making it difficult to create classical Wolf-Rayet stars via this channel. However, a remaining key uncertainty regarding these predictions concerns unsteady mass loss for very high-luminosity stars close to the Eddington limit as well as the impact of non-line-driven winds.

Small fluctuations around homogeneous and isotropic expanding backgrounds are the main object of study in cosmology. Their origin and evolution is sensitive to the physical processes that happen during inflation and in the late Universe. As such, they hold the key to answering many of the major open questions in cosmology. Given a large separation of relevant scales in many examples of interest, the most natural description of these fluctuations is formulated in terms of effective field theories. This was the main avenue for many of the important modern developments in theoretical cosmology, which provided a unifying framework for a plethora of cosmological models and made a clear connection between the fundamental cosmological parameters and observables. In this review we summarize these results in the context of effective field theories of inflation, large-scale structure, and dark energy.

Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and planet formation. We aim to study the the D/ H ratio in H2CS across the prototypical pre-stellar core L1544. This study allows us to test current gas-dust chemical models involving sulfur in dense cores. We present here single-dish observations of H2CS, HDCS and D2CS with the IRAM 30m telescope. We analyse their column densities and distributions, and compare these observations with gas-grain chemical models. The deuteration maps of H2CS in L1544 are compared with the deuteration maps of methanol, H2CO, N2H+ and HCO+ towards the same source. Furthermore, the single and double deuteration of H2CS towards the dust peak of L1544 is compared with H2CO and c-C3H2. The difference between the deuteration of these molecules in L1544 is discussed and compared with the prediction of chemical models. The maximum deuterium fractionation for the first deuteration of H2CS is N(HDCS)/N(H2CS)$\sim$30$\%$ and is located towards the north-east at a distance of about 10000 AU from the dust peak. While for c-C3H2 the first and second deuteration have a similar efficiency, for H2CS and H2CO the second deuteration is more efficient, leading to D2CX/HDCX$\sim$100$\%$ (with X= O or S). Our results imply that the large deuteration of H2CO and H2CS observed in protostellar cores as well as in comets is likely inherited from the pre-stellar phase. However, the comparison with state-of-the-art chemical models suggests that the reaction network for the formation of the doubly deuterated H2CS and H2CO it is not complete yet.

Dwarf galaxies are missing nearly all of their baryons and metals from the stellar disk, presumed to be in a bound halo or expelled beyond the virial radius. The virial temperature for galaxies with $M_{\rm h} \sim 10^9 - 10^{10}$ $M_{\odot}$ is similar to the collisional ionization equilibrium temperature for the C IV ion. We searched for UV absorption from C IV in six sightlines toward three dwarf galaxies in the anti-M31 direction and at the periphery of the Local Group ($D \approx$ 1.3 Mpc; Sextans A, Sextans B, and NGC 3109). The C IV doublet is detected in only one of six sightlines, toward Sextans A, with $\log N({\rm C IV})$ = $13.06 \pm 0.08$. This is consistent with our gaseous halo models, where the halo gas mass is determined by the cooling rate, feedback, and the star formation rate; the inclusion of photoionization is an essential ingredient. This model can also reproduce the higher detection rate of O VI absorption in other dwarf samples (beyond the Local Group), and with C IV only detectable within $\sim 0.5R_{\rm vir}$.

We study the information content of summary statistics built from the multi-scale topology of large-scale structures on primordial non-Gaussianity of the local and equilateral type. We use halo catalogs generated from numerical N-body simulations of the Universe on large scales as a proxy for observed galaxies. Besides calculating the Fisher matrix for halos in real space, we also check more realistic scenarios in redshift space. Without needing to take a distant observer approximation, we place the observer on a corner of the box. We also add redshift errors mimicking spectroscopic and photometric samples. We perform several tests to assess the reliability of our Fisher matrix, including the Gaussianity of our summary statistics and convergence. We find that the marginalized 1-$\sigma$ uncertainties in redshift space are $\Delta f_{\rm NL}^{\rm loc} \sim 16$ and $\Delta f_{\rm NL}^{\rm equi} \sim 41 $ on a survey volume of $1$ $($Gpc$/h)^3$. These constraints are weakly affected by redshift errors. We close by speculating as to how this approach can be made robust against small-scale uncertainties by exploiting (non)locality.

{"Bare collapse", the collapse of a bare stellar core to a neutron star with a very small mass ejection links two seemingly unrelated phenomena: the formation of binary neutron star (BNS) systems and the observations of fast and luminous optical transients. We carried out calculations of the collapse due to electron-capture of both evolutionary and synthetic isentropic bare stellar cores. We find that the collapse results in {the formation of} a light ~ 1.3 solar mass neutron star and {an} ejection of ~0.1 solar mass at ~0.1c. The outer shell of the ejecta is composed of Ni56 that can power an ultra-stripped supernova. The models we explored can explain most of the observed fast optical flares but not the brightest ones. Collapse of cores surrounded by somewhat more massive envelopes can produce larger amounts of Ni56 and explain brighter flares. Alternatively, those events can arise due to interaction of the very energetic ejecta with winds that were ejected from the progenitor a few days before the collapse.

We present 3D hydrodynamic models of the interaction between the outflows of evolved, pulsating, Asymptotic Giant Branch (AGB) stars and nearby ($< 3$ stellar radii) sub-stellar companions ($M_{\mathrm{comp}} \lesssim 40$ M$_J$). Our models show that due to resonances between the orbital period of the companion and the pulsation period of the AGB star, multiple spiral structures can form; the shocks driven by the pulsations are enhanced periodically in different regions as they encounter the denser material created by the sub-stellar companion's wake. We discuss the properties of these spiral structures and the effect of the companion parameters on them. We also demonstrate that the gravitational potential of the nearby companion enhances the mass loss from the AGB star. For more massive ($M_{\mathrm{comp}} > 40$ M$_J$) and more distant companions ($> 4$ stellar radii), a single spiral arm forms. We discuss the possibility of observing these structures with the new generations of high-resolution, high-sensitivity instruments, and using them to `find' sub-stellar companions around bright, evolved stars. Our results also highlight possible structures that could form in our solar system when the Sun turns into an AGB star.

Discovering radio pulsars, particularly millisecond pulsars (MSPs), is important for a range of astrophysical applications, such as testing theories of gravity or probing the magneto-ionic interstellar medium. We aim to discover pulsars that may have been missed in previous pulsar searches by leveraging known pulsar observables (primarily polarisation) in the sensitive, low-frequency radio images from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS), and have commenced the Targeted search, using LoTSS images, for polarised pulsars (TULIPP) survey. For this survey, we identified linearly and circularly polarised point sources with flux densities brighter than 2 mJy in LoTSS images at a centre frequency of 144 MHz with a 48 MHz bandwidth. Over 40 known pulsars, half of which are MSPs, were detected as polarised sources in the LoTSS images and excluded from the survey. We have obtained beam-formed LOFAR observations of 30 candidates, which were searched for pulsations using coherent de-dispersion. Here, we present the results of the first year of the TULIPP survey. We discovered two pulsars, PSRs J1049+5822 and J1602+3901, with rotational periods of P=0.73 s and 3.7 ms, respectively. We also detected a further five known pulsars (two slowly-rotating pulsars and three MSPs) for which accurate sky positions were not available to allow a unique cross-match with LoTSS sources. This targeted survey presents a relatively efficient method by which pulsars, particularly MSPs, may be discovered using the flexible observing modes of sensitive radio telescopes such as the Square Kilometre Array and its pathfinders/precursors, particularly since wide-area all-sky surveys using coherent de-dispersion are currently computationally infeasible.

Low luminosity active galactic nuclei (LLAGN) probe accretion physics in the low Eddington regime and can provide additional clues about galaxy evolution. AGN variability is ubiquitous and thus provides a reliable tool for finding AGN. We analyze the All-Sky Automated Survey for Supernovae light curves of 1218 galaxies with $g<14$ mag and Sloan Digital Sky Survey spectra in search of AGN. We find 37 objects that are both variable and have AGN-like structure functions, which is about 3% of the sample. The majority of the variability selected AGN are LLAGN with Eddington ratios ranging from $10^{-4}$ to $10^{-2}$. We thus estimate the fraction of LLAGN in the population of galaxies as 2% down to a median Eddington ratio of $2\times10^{-3}$. Combining the BPT line ratio diagnostics and the broad line AGN, up to $\sim$60% of the AGN candidates are confirmed spectroscopically. The BPT diagnostics also classified 10-30% of the candidates as star forming galaxies rather than AGN.

Deep learning methods have been employed in gravitational-wave astronomy to accelerate the construction of surrogate waveforms for the inspiral of spin-aligned black hole binaries, among other applications. We demonstrate, that the residual error of an artificial neural network that models the coefficients of the surrogate waveform expansion (especially those of the phase of the waveform) has sufficient structure to be learnable by a second network. Adding this second network, we were able to reduce the maximum mismatch for waveforms in a validation set by more than an order of magnitude. We also explored several other ideas for improving the accuracy of the surrogate model, such as the exploitation of similarities between waveforms, the augmentation of the training set, the dissection of the input space, using dedicated networks per output coefficient and output augmentation. In several cases, small improvements can be observed, but the most significant improvement still comes from the addition of a second network that models the residual error. Since the residual error for more general surrogate waveform models (when e.g. eccentricity is included) may also have a specific structure, one can expect our method to be applicable to cases where the gain in accuracy could lead to significant gains in computational time.

We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10m Keck telescope to search for Ly$\alpha$ emission from candidate galaxies at z ~9-10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Ly$\alpha$ from one object in ~8.1 hours of integration, at z = 8.665 +/- 0.001 with an integrated signal-to-noise > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field, and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Ly$\alpha$ detected at z=8.683 (Zitrin et al. 2015). The detection of Ly$\alpha$ suggests the existence of large (~1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV Luminosity Function down to M_UV = -13). The Ly$\alpha$ we detect would be able to escape the predominantly-neutral intergalactic medium at this epoch if our galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. 2021. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 which is powering one (or more) ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field.

We perform extensive Monte Carlo simulations to systematically compare the frequentist and Bayesian treatments of the Lomb--Scargle periodogram. The goal is to investigate whether the Bayesian period search is advantageous over the frequentist one in terms of the detection efficiency, how much if yes, and how sensitive it is regarding the choice of the priors, in particular in case of a misspecified prior (whenever the adopted prior does not match the actual distribution of physical objects). We find that the Bayesian and frequentist analyses always offer nearly identical detection efficiency in terms of their tradeoff between type-I and type-II mistakes. Bayesian detection may reveal a formal advantage if the frequency prior is nonuniform, but this results in only $\sim 1$ per cent extra detected signals. In case if the prior was misspecified (adopting nonuniform one over the actual uniform) this may turn into an opposite advantage of the frequentist analysis. Finally, we revealed that Bayes factor of this task appears rather overconservative if used without a calibration against type-I mistakes (false positives), thereby necessitating such a calibration in practice.

The Milky Way is believed to host hundreds of millions of quiescent stellar-mass black holes (BHs). In the last decade, some of these objects have been potentially uncovered via gravitational microlensing events. All these detections resulted in a degeneracy between the velocity and the mass of the lens. This degeneracy has been lifted, for the first time, with the recent astrometric microlensing detection of OB110462. However, two independent studies reported very different lens mass for this event. Sahu et al. (2022) inferred a lens mass of 7.1 $\pm$ 1.3 Msun, consistent with a BH, while Lam et al. (2022) inferred 1.6-4.2 Msun, consistent with either a neutron star or a BH. Here, we study the landscape of isolated BHs formed in the field. In particular, we focus on the mass and center-of-mass speed of four sub-populations: isolated BHs from single-star origin, disrupted BHs of binary-star origin, main-sequence stars with a compact object companion, and double compact object mergers. Our model predicts that most ($\gtrsim$ 50%) isolated BHs in the Milky Way are of binary origin. Moreover, the origin of low-speed (< 50 km/s) isolated BHs depends on their mass: at least $\approx$ 70% of low-mass ($\lesssim$ 10 Msun) BHs are from binary origin. Under the assumption that OB110462 is a free-floating compact object, we conclude that it is more likely to be a BH originally belonging to a binary system. Our results suggest that low-speed BH microlensing events can be useful to understand binary evolution of massive stars in the Milky Way.

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

The missing baryon problem may now be resolved, but the exact location and physical properties of the diffuse component remains unclear. This problem could be tractable, but requires the combination of new galaxy redshift surveys with new X-ray and radio facilities.

A pseudoscalar "axionlike" field, $\phi$, may explain the $3\sigma$ hint of cosmic birefringence observed in the $EB$ power spectrum of the cosmic microwave background (CMB) polarization data. Is $\phi$ dark energy or dark matter? A tomographic approach can answer this question. The effective mass of dark energy field responsible for the accelerated expansion of the Universe today must be smaller than $m_\phi\simeq 10^{-33}$ eV. If $m_\phi \gtrsim 10^{-32}$ eV, $\phi$ starts evolving before the epoch of reionization and we should observe different amounts of birefringence from the $EB$ power spectrum at low ($l\lesssim 10$) and high multipoles. Such an observation, which requires a full-sky satellite mission, would rule out $\phi$ being dark energy. If $m_\phi \gtrsim 10^{-28}$ eV, $\phi$ starts oscillating during the epoch of recombination, leaving a distinct signature in the $EB$ power spectrum at high multipoles, which can be measured precisely by ground-based CMB observations. Our tomographic approach relies on the shape of the $EB$ power spectrum and is less sensitive to miscalibration of polarization angles.

Flat spectrum radio quasar (FSRQ) is the most luminous blazar at the GeV energies. In this paper, we probe the photon-axion-like particle (ALP) oscillation effect on the latest very-high-energy (VHE) $\gamma$-ray observations of the FSRQ 4C+21.35 (PKS 1222+216). The $\gamma$-ray spectra are measured by the collaborations Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC), Very Energetic Radiation Imaging Telescope Array System (VERITAS), and Fermi Large Area Telescope (Fermi-LAT), which cover two activity VHE flares of 4C+21.35 in 2010 and 2014. We show the spectral energy distributions (SEDs) of these two phases under the null and ALP hypotheses, and set the combined limit on the ALP parameter space. The 95% $\rm C.L.$ combined limit set by the FSRQ 4C+21.35 observations measured by MAGIC, VERITAS, and Fermi-LAT in the $m_a-g_{a\gamma}$ plane is roughly at the photon-ALP coupling $g_{a\gamma} \gtrsim 8\times 10^{-12} \rm \, GeV^{-1}$ for the ALP mass $[\,2\times 10^{-10}\, {\rm eV} \lesssim m_a \lesssim 2\times 10^{-8}\, \rm eV\,]$. Compared with the constraint of NGC 1275 set by Fermi-LAT, no stringent limit result is derived with the photon-ALP coupling $g_{a\gamma}$ from the FSRQ 4C+21.35, while this result could slightly broaden the ALP mass $m_a$ limit at the low-mass region.

We study the broadband radiation behavior of a kilonova ejecta-pulsar wind nebula (PWN) system. In this model, we jointly fit the observations of AT 2017gfo in UV-optical-IR bands and the late-time X-ray afterglow of GRB 170817A. Our work shows that a PWN powered by the remnant neutron star (NS) post GW170817 event could affect the optical transient AT 2017gfo and re-brighten the late-time X-ray afterglow of GRB 170817A. The PWN radiation will regulate the trend of future X-ray observations from a flattening to a steep decline until some other sources (e.g., a kilonova afterglow) become dominant. The restricted ranges of the central NS parameters in this work are consistent with the previous works based on the observations of AT 2017gfo only. In addition, the new fitting result indicates that the NS wind is highly magnetized. We point out that the radio and X-ray emission from a kilonova ejecta-PWN system could be an important electromagnetic feature of binary NS mergers when a long-lived remnant NS is formed. Therefore, observations of a kilonova ejecta-PWN system will provide important information to inferring the nature of a merger remnant.

The Onsala twin telescopes (OTT) are two 13 m telescopes located at the Onsala Space Observatory in Sweden. With dual linear polarized broad-band (3-14 GHz) receivers, they are part of the next generation Very Long Baseline Interferometry (VLBI) Global Observing System (VGOS) for geodesy and astrometry. The VGOS network is currently being rolled out globally, and regularly observes hundreds of quasars. In contrast to legacy S/X-band geodetic VLBI, VGOS will eventually - in addition to purely geodetic data products - also produce full-polarisation images and flux densities of all observed radio quasars on a regular basis. This will be a rich archive of broad-band quasar monitoring information, of value both for geodesy and astronomy. We scheduled, correlated and analysed 91 short (<30 min) observing sessions with the OTT spanning 7 months, with 33 sessions placed within a 3 day period to search for short-term variability. We monitored seven quasars (0059+581, 0552+398, 1144+402, 1156+295, 1617+229, 3C418, OJ287) and three absolute flux-density reference calibrators (3C147, 3C286, 3C295). We used the Common Astronomy Software Applications (CASA) package to fringe-fit, bandpass-correct and scale (using auto-correlations, antenna gains and measured system temperatures) the data to obtain simultaneous flux densities in the four standard VGOS bands: A = 3.0-3.5 GHz, B= 5.2-5.7 GHz, C=6.3-6.8 GHz, and D=10.2-10.7 GHz. After correcting for instrumental biases, we find the OTT capable of monitoring flux densities with 5 % uncertainty in all four VGOS bands for Jansky-level sources. Three quasars show significant variability during our observing campaign: 0059+581, OJ287 and 1156+295. We also find a tentative elevation-dependence of about 5 %, suggesting that more detailed characterisation of the antenna gain curves could further improve the accuracy of the monitoring.

In this paper we describe a project we initiated to investigate individual pixels in the downloaded Kepler apertures in order to find objects in the background of the main targets with variable brightness. In the first paper of this series we discovered and investigated 547 short-period eclipsing binaries (Bienias et al. 2021). Here we present the independent discovery of 26 new RR Lyrae stars in the Kepler background pixels obtained during the primary mission, and provide continuous and precise photometry for these objects. Twenty-one of these stars were already noted by Gaia or the Pan-STARRS survey. This new population of dominantly faint and distant RR Lyrae stars increases by 50% and complements nicely the 52 already known main target RR Lyrae stars in the original Kepler field. Despite their faintness, the four-year quasi-uninterrupted light curves of these stars allow an unprecedented view of these faint halo objects. We present an analysis of the light curves of the new RR Lyrae sample, verify their classification using Fourier parameters, and discuss the properties of these newly found pulsating variable stars. Most notably, this is the first time that such faint RR Lyrae stars have been investigated with the help of a photometric data set with outstanding cadence and precision. Interestingly, these objects share the properties of their brighter siblings in terms of sub-class characteristics, additional mode content, and modulation occurrence rates.

We derive the metallicity (traced by the O/H abundance) of the Narrow Line Region ( NLR) of 108 Seyfert galaxies as well as radial metallicity gradients along their galaxy disks and of these of a matched control sample of no active galaxies. In view of that, observational data from the SDSS-IV MaNGA survey and strong emission-line calibrations taken from the literature were considered. The metallicity obtained for the NLRs %each Active Galactic Nucleus (AGN) was compared to the value derived from the extrapolation of the radial oxygen abundance gradient, obtained from \ion{H}{ii} region estimates along the galaxy disk, to the central part of the host galaxies. We find that, for most of the objects ($\sim 80\,\%$), the NLR metallicity is lower than the extrapolated value, with the average difference ($<D>$) between these estimates ranging from 0.16 to 0.30 dex. We suggest that $<D>$ is due to the accretion of metal-poor gas to the AGN that feeds the nuclear supermassive black hole (SMBH), which is drawn from a reservoir molecular and/or neutral hydrogen around the SMBH. Additionally, we look for correlations between $D$ and the electron density ($N_{\rm e}$), [\ion{O}{iii}]$\lambda$5007 and H$\alpha$ luminosities, extinction coefficient ($A_{V})$ of the NLRs, as well as the stellar mass ($M_{*}$) of the host galaxies. Evidences of an inverse correlation between the $D$ and the parameters $N_{\rm e}$, $M_{*}$ and $A_{\rm v}$ were found.

Tumer et al. 2022, ATel #15171, have recently reported the discovery of an X-ray source, NuSTAR J053449+2126.0, during a calibration observation which took place on 25 April 2020. We scan the Zwicky Transient Facility (ZTF) alerts and archival photometry to determine the nature of the source. Palomar Gattini-IR is searched as well. We identify no obvious counterpart candidate. Follow-up X-ray and optical studies are needed to determine the true counterpart.

Meteoritic studies of solar system objects show evidence of nucleosynthetic heterogeneities that are inherited from small presolar grains (< 10 $\mu$m) formed in stellar environments external to our own. The initial distribution and subsequent evolution of these grains are currently unconstrained. Using 3D, gas-dust simulations, we find that isotopic variations on the order of those observed in the solar system can be generated and maintained by drag and viscosity. Small grains are dragged radially outwards without size/density sorting by viscous expansion and backreaction, enriching the outer disc with presolar grains. Meanwhile large aggregates composed primarily of silicates drift radially inwards due to drag, further enriching the relative portion of presolar grains in the outer disc and diluting the inner disc. The late accumulation of enriched aggregates outside Jupiter could explain some of the isotopic variations observed in solar system bodies, such as the enrichment of supernovae derived material in carbonaceous chondrites. We also see evidence for isotopic variations in the inner disc that may hold implications for enstatite and ordinary chondrites that formed closer to the Sun. Initial heterogeneities in the presolar grain distribution that are not continuously reinforced are dispersed by diffusion, radial surface flows, and/or planetary interactions over the entire lifetime of the disc. For younger, more massive discs we expect turbulent diffusion to be even more homogenising, suggesting that dust evolution played a more central role in forming the isotopic anomalies in the solar system than originally thought.

Numerous observational evidence has suggested the presence of active meteorology in the atmospheres of brown dwarfs. A near-infrared brightness variability has been observed. Clouds have a major role in shaping the thermal structure and spectral properties of these atmospheres. The mechanism of such variability is still unclear and both 1D and global circulation model cannot fully study this topics due to resolution. In this study, a convective resolving model is coupled to grey-band radiative transfer in order to study the coupling between the convective atmosphere and the variability of clouds over a large temperature range with a domain of several hundreds of kilometers. Six types of clouds are considered, with microphysics including settling. The clouds are radiatively active using Rosseland mean coefficient. Radiative cloud feedback can drive spontaneous atmospheric variability in both temperature and cloud structure, as modeled for the first time in three dimensions. Silicate clouds have the most effect of the thermal structure with the generation of a secondary convective layer in some cases, depending on the assumed particle size. Iron and Aluminum clouds also have a substantial impact on the atmosphere. Thermal spectra were computed, and we find the strongest effect of clouds is the smoothing of spectral features at optical wavelengths. Compared to observed L and T dwarfs on color-magnitude diagram, the simulated atmospheres are redder for most of the cases. The simulations with the presence of cloud holes are closer to the observations.

The infrared solar spectrum contains a wealth of physical data about the Sun and is being explored using modern detectors and technology with new ground-based solar telescopes. One such instrument will be the ground-based Cryogenic Near-IR Spectro-Polarimeter of the Daniel K. Inouye Solar Telescope that will be capable of sensitive imaging of the faint infrared solar coronal spectra with full Stokes I, Q, U, and V polarization states. Highly ionized magnetic dipole emission lines have been observed in galaxies and the solar corona. Quantifying the accuracy of spectral inversion procedures requires a precise spectroscopic calibration of observations. A careful interpretation of the spectra around prominent magnetic dipole lines is essential for deriving physical parameters, and particularly, for quantifying the off-limb solar coronal observations from DKIST. In this work, we aim to provide an analysis of the spectral regions around the infrared coronal emission lines of Fe XIII 1074.68 nm, Fe XIII 1079.79 nm, Si X 1430.10 nm, and Si IX 3934.34 nm, aligning with the goal of identifying solar photospheric and telluric lines that will help facilitate production of reliable inversions and data products from four sets of solar coronal observations. The outputs will be integrated in the processing pipeline to produce Level-2 science-ready data that will be made available to DKIST observers.

In the framework of compressional magnetohydrodynamics (MHD), we numerically studied the commonly accepted presumption that the Alfv\'enic turbulence is generated by the collisions between counter-propagating Alfv\'en waves (AWs). In the conditions typical for the low-beta solar corona and inner solar wind, we launched two counter-propagating AWs in the three-dimensional simulation box and analyzed polarization and spectral properties of perturbations generated before and after AW collisions. The observed post-collisional perturbations have different polarizations and smaller cross-field scales than the original waves, which supports theoretical scenarios with direct turbulent cascades. However, contrary to theoretical expectations, the spectral transport is strongly suppressed at the scales satisfying the classic critical balance of incompressional MHD. Instead, a modified critical balance can be established by colliding AWs with significantly shorter perpendicular scales. We discuss consequences of these effects for the turbulence dynamics and turbulent heating of compressional plasmas. In particular, solar coronal loops can be heated by the strong turbulent cascade if the characteristic widths of the loop substructures are more than ten times smaller than the loop width. The revealed new properties of AW collisions have to be incorporated in the theoretical models of AW turbulence and related applications.

It has been widely recognized that gamma-ray burst (GRB) afterglows arise from interactions between GRB outflow and circumburst medium, while their evolution follows the behaviors of relativistic shock waves. Assuming the distribution of circumburst medium follows a general power-law form, that is, $n = A_{\ast} R^{-k}$, where $R$ denotes the distance from the burst, it is obvious that the value of density-distribution index $k$ can affect the behaviors of the afterglow. In this paper, we analyze the temporal and spectral behaviors of GRB radio afterglows with arbitrary $k$-values. In the radio band, a standard GRB afterglow produced by forward shock exhibits a late-time flux peak, and the relative peak fluxes as well as peak times at different frequencies show dependencies on $k$. Thus with multi-band radio peak observations, one can determine the density profile of circumburst medium by comparing the relations between peak flux/time and frequency at each observing band. Also, the effects of trans-relativistic shock waves, as well as jets in afterglows are discussed. By analyzing 31 long and 1 short GRBs with multi-band data of radio afterglows, we find that nearly half of them can be explained with uniform interstellar medium ($k=0$), $\sim 1/5$ can be constrained to exhibiting stellar wind environment ($k=2$), while less than $\sim 1/3$ samples show $0< k< 2$.

Galaxy flybys, interactions where two independent halos inter-penetrate but detach at a later time and do not merge, occur frequently at lower redshifts. These interactions can significantly impact the evolution of individual galaxies - from the mass loss and shape transformation to the emergence of tidal features and formation of morphological disc structures. The main focus of this paper is on the dark matter mass loss of the secondary, intruder galaxy, with the goal of determining a functional relationship between the impact parameter and dark matter mass loss. Series of N-body simulations of typical galaxy flybys (10:1 mass ratio) with differing impact parameters show that the dark matter halo leftover mass of the intruder galaxy follows a logarithmic growth law with impact parameter, regardless of the way the total halo mass is estimated. The lost mass then, clearly, follows the exponential decay law. The stellar component stretches faster as the impact parameter decreases, following the exponential decay law with impact parameter. Functional dependence on impact parameter in all cases seems universal, but the fitting parameters are likely sensitive to the interaction parameters and initial conditions (e.g. the mass ratio of interacting galaxies, initial relative velocity of the intruder galaxy, interaction duration). While typical flybys, investigated here, could not be the sole culprit behind the formation of ultra-diffuse or dark matter deficient galaxies, they can still contribute significantly. Rare, atypical and stronger flybys are worth further exploring.

We construct a model for the Milky Way where the interstellar medium (ISM) is equipped with self-consistent dynamics. In simulations a spiral structure emerges from this model that is almost identical with the one in the Milky Way's ISM. Further, the Jeans instability offers an explanation for the observed velocity dispersion of atomic hydrogen in the ISM; this instability vanishes from our model if we choose a velocity dispersion just above the observed one. Surprisingly, our model gets along completely without dark matter. The 'missing mass' distributes uniformly over the baryonic components making it possible to explain the occurring mass gap.

The VISTA Variables in the Via Lactea (VVV) survey obtained near-infrared photometry toward the Galactic bulge and the southern disc plane for a decade (2010 - 2019). We designed a modified Lomb-Scargle method to search for large-amplitude ($\Delta$Ks > 1.5 mag) mid to long-term periodic variables (P > 10 d) in the 2nd version of VVV Infrared Astrometric Catalogue (VIRAC2-$\beta$). In total, 1520 periodic sources were discovered, including 59 candidate periodic outbursting young stellar objects (YSOs), based on the unique morphology of the phase-folded light curves, proximity to Galactic HII regions and mid-infrared colours. Five sources are spectroscopically confirmed as accreting YSOs. Both fast-rise/slow-decay and slow-rise/fast-decay periodic outbursts were found, but fast-rise/slow-decay outbursts predominate at the highest amplitudes. The multi-wavelength colour variations are consistent with a variable mass accretion process, as opposed to variable extinction. The cycles are likely to be caused by dynamical perturbations from stellar or planetary companions within the circumstellar disc. An additional search for periodic variability amongst YSO candidates in published Spitzer-based catalogues yielded a further 71 candidate periodic accretors, mostly with lower amplitudes. These resemble cases of pulsed accretion but with unusually long periods and greater regularity. The majority of other long-period variables are pulsating dusty Miras with smooth and symmetric light curves. We find that some Miras have redder $W3 - W4$ colours than previously thought, most likely due to their surface chemical compositions.

Convective-core overshoot mixing is a significant uncertainty in stellar evolution. Because numerical simulations and turbulent convection models predict exponentially decreasing radial rms turbulent velocity, a popular treatment of the overshoot mixing is to apply a diffusion process with exponentially decreasing diffusion coefficient. It is important to investigate the parameters of the diffusion coefficient because they determine the efficiency of the mixing in the overshoot region. In this paper, we have investigated the effects of the core overshoot mixing on the properties of the core in solar models and have constrained the parameters of the overshoot model by using helioseismic inferences and the observation of the solar 8B neutrino flux. For solar-mass stars, the core overshoot mixing helps to prolong the lifetime of the convective core developed at the ZAMS. If the strength of the mixing is sufficiently high, the convective core in a solar model could survive till the present solar age, leading to large deviations of the sound-speed and density profiles comparing with the helioseismic inferences. The 8B neutrino flux also favours a radiative solar core. Those provide a constraint on the parameters of the exponential diffusion model of the convective overshoot mixing. A limited asteroseismic investigation of 13 Kepler low-mass stars with 1.0 < M < 1.5 shows a mass-dependent range of the overshoot parameter. The overshoot mixing processes for different elements are analyzed in detail. It is found that the exponential diffusion overshoot model leads to different effective overshoot mixing lengths for elements with different nuclear equilibrium timescale.

We present the catalog of Interplanetary Network (IPN) localizations for 199 short-duration gamma-ray bursts (sGRBs) detected by the Konus-Wind (KW) experiment between 2011 January 1 and 2021 August 31, which extends the initial sample of IPN localized KW sGRBs (arXiv:1301.3740) to 495 events. We present the most comprehensive IPN localization data on these events, including probability sky maps in HEALPix format.

An estimate of the expected photon flux above $10^{17}~$eV from the interactions of ultra-high energy cosmic rays with the matter in the Galactic disk is presented. Uncertainties arising from the distribution of the gas in the disk, the absolute level of the cosmic ray flux, and the composition of the cosmic rays are taken into account. Within these uncertainties, the integrated photon flux above $10^{17}~$eV is, averaged out over Galactic latitude less than $5^\circ$, between $\simeq 3.2{\times}10^{-2}~$km$^{-2}~$yr$^{-1}~$sr$^{-1}$ and $\simeq 8.7{\times}10^{-2}~$km$^{-2}~$yr$^{-1}~$sr$^{-1}$. The all-sky average value amounts to $\simeq 1.1{\times}10^{-2}~$km$^{-2}~$yr$^{-1}~$sr$^{-1}$ above $10^{17}~$eV and decreases roughly as $E^{-2}$, making this diffuse flux the dominant one from cosmic-ray interactions for energy thresholds between $10^{17}~$eV and $10^{18}$~eV. Compared to the current sensitivities of detection techniques, a gain between two and three orders of magnitude in exposure is required for a detection below $\simeq 10^{18}$~eV. The implications for searches for photon fluxes from the Galactic center that would be indicative of the decay of super-heavy dark matter particles are discussed, as the photon flux presented in this study can be considered as a floor below which other signals would be overwhelmed.

The large-scale structure of the Universe is a rich source of information to test the consistency of General Relativity on cosmological scales. We briefly describe how the observed distributions of galaxies is affected by redshift-space distortions, but also by gravitational lensing and other relativistic effects. Then, we show how one of this relativistic effects, the gravitational redshift, can be used to build a model independent test that directly measures the anisotropic stress, i.e. the difference between the two gravitational potentials that describe spacetime fluctuations of the geometry.

The current epoch of accelerated cosmic expansion is postulated to be driven by dark energy, which in the standard model takes the form of a cosmological constant with equation of state parameter $w=-1$. We propose an innovative perspective over the nature of dark energy by drawing a parallel with inflation, which we assume to be driven by a single scalar field, the inflaton. The inflaton was not a cosmological constant, as indicated by the fact that inflation ended and by the Planck satellite's constraint of $n_s\neq 1$ at $8\sigma$ confidence. Therefore, it is interesting to verify whether its equation of state parameter was measurably different from $-1$. We analyze this question for a class of single-field slow-roll inflationary models, where the hierarchy of Hubble slow-roll parameters is truncated at different orders. Based on the latest Planck and BICEP2/Keck data, we obtain a $68\%$ upper bound of $1+w < 0.0014$ for the three-parameter model, which gives the best description to the data. This provides a cautionary tale for drawing conclusions about the nature of today's dark energy based upon the non-detection of a deviation from $w=-1$ with current and upcoming cosmological surveys.

Magnetic fields permeate the interstellar medium and are important in the star formation process. Determining the 3D magnetic fields of molecular clouds will allow us to better understand their role in the evolution of these clouds and formation of stars. We fully reconstruct the approximate three-dimensional (3D) magnetic field morphology of the Orion A molecular cloud (on scales of a few to ~100 pc) using Galactic magnetic field models, as well as present line-of-sight and plane-of-sky magnetic field observations. While previous studies identified the 3D magnetic field morphology of the Orion A cloud as an arc shape, in this study we provide the orientation of this arc-shaped field and its plane-of-sky direction, for the first time. We find that this 3D field is a tilted, semi-convex (from our point of view) structure and mostly points in the direction of decreasing latitude and longitude on the plane of the sky, from our vantage point. The previously identified bubbles and events in this region were key in shaping this arc-shaped magnetic field morphology.

We present a detailed study of fossil group candidates identified in "Last Journey", a gravity-only cosmological simulation covering a $(3.4\, h^{-1}\mathrm{Gpc})^3$ volume with a particle mass resolution of $m_p \approx 2.7 \times 10^9\, h^{-1}\mathrm{M}_\odot$. The simulation allows us to simultaneously capture a large number of group-scale halos and to resolve their internal structure. Historically, fossil groups have been characterized by high X-ray brightness and a large luminosity gap between the brightest and second brightest galaxy in the group. In order to identify candidate halos that host fossil groups, we use halo merger tree information to introduce two parameters: a luminous merger mass threshold ($M_\mathrm{LM}$) and a last luminous merger redshift cut-off ($z_\mathrm{LLM}$). The final parameter choices are informed by observational data and allow us to identify a plausible fossil group sample from the simulation. The candidate halos are characterized by reduced substructure and are therefore less likely to host bright galaxies beyond the brightest central galaxy. We carry out detailed studies of this sample, including analysis of halo properties and clustering. We find that our simple assumptions lead to fossil group candidates that form early, have higher concentrations, and are more relaxed compared to other halos in the same mass range.

The Pierre Auger Observatory have reported [1-3] observation of several exotic events (EE) which apparently related to thunderstorms. These events are much larger in size than conventional cosmic ray (CR) events, and they have very distinct timing features. A possible nature of the observed phenomenon is still a matter of active research and debates as many unusual features of these exotic events are hard to explain. In particular, the frequency of appearance of these EE is very low (less than 2 events/year), in huge contrast with a typical rate of a conventional lightning strikes in the area. We propose that the observed EE can be explained within the so-called axion quark nugget (AQN) dark matter model. The idea is that the AQNs may trigger and initiate a special and unique class of lightning strikes during a thunderstorm as a result of ionization of the atmospheric molecules along its path. The corresponding AQN-induced lighting flashes may show some specific features not shared by typical and much more frequent conventional flashes. We support this proposal by demonstrating that the observations[1-3], including the frequency of appearance and time duration are consistent with observations. We also comment on possible relation of AUGER EEs with the Telescope Array bursts and the terrestrial gamma ray flashes (TGF). We list a number of features of the AQN-induced EE (such as specific radio pulses synchronized with these events) which can be directly tested by future experiments. We also suggest to use distributed acoustic sensing (DAS) instruments to detect the acoustic pulses which must be synchronized with AUGER EEs.

The quest to discover the nature of dark matter continues to drive many of the experimental and observational frontiers in particle physics, astronomy, and cosmology. While there are no definitive signatures to date, there exists a rich ecosystem of experiments searching for signals for a broad class of dark matter models, at different epochs of cosmic history, and through a variety of processes with different characteristic energy scales. Given the multitude of candidates and search strategies, effective field theory has been an important tool for parametrizing the possible interactions between dark matter and Standard Model probes, for quantifying and improving model-independent uncertainties, and for robust estimation of detection rates in the presence of large perturbative corrections. This white paper summarizes a wide range of effective field theory applications for connecting dark matter theories to experiments.

An essential goal of the Higgs physics program at the LHC and beyond is to explore the nature of the Higgs potential and shed light on the mechanism of electroweak symmetry breaking. An important class of models defining the strength and order of the electroweak phase transition is driven by the Higgs boson coupling to a light new state. This Snowmass white paper points out the existence of a region of parameter space where a strongly first order electroweak phase transition is compatible with exotic decays of the SM-like Higgs boson. A dedicated search for exotic Higgs decays can actively explore this framework at the Large Hadron Collider (LHC), while future exotic Higgs decay searches at the high-luminosity LHC and future Higgs factories will be vital to conclusively probe the scenario.

The next generation of gravitational-wave observatories can explore a wide range of fundamental physics phenomena throughout the history of the universe. These phenomena include access to the universe's binary black hole population throughout cosmic time, to the universe's expansion history independent of the cosmic distance ladders, to stochastic gravitational-waves from early-universe phase transitions, to warped space-time in the strong-field and high-velocity limit, to the equation of state of nuclear matter at neutron star and post-merger densities, and to dark matter candidates through their interaction in extreme astrophysical environments or their interaction with the detector itself. We present the gravitational-wave detector concepts than can drive the future of gravitational-wave astrophysics. We summarize the status of the necessary technology, and the research needed to be able to build these observatories in the 2030s.

We consider a QCD cold plasma motivated Equation of State (EOS) to examine the impact of an Anomalous Magnetic Moment (AMM) coupling and small shape deformations for static oblate and prolate core shapes of quark stars. Using the Foga\c{c}a QCD motivated EOS which shifts from the high temperature low chemical potential quark gluon plasma environment to the low temperature high chemical potential quark stellar core environment we consider the impact of an AMM coupling with a metric induced shape deformation parameter in the TOV equations. The EOS is developed using a hard gluon and soft gluon decomposition of the gluon field tensor using a mean field effective mass for the gluons. The AMM is considered using the Dirac spin tensor coupled to the EM field tensor with quark flavor based magnetic moments. The shape parameter is introduced in a metric ansatz that represents oblate and prolate static stellar cores for modified TOV equations. These equations are numerically solved for the final mass and radius states representing the core collapse of a massive star with a phase transition leading to an unbound quark-gluon plasma. We find that the combined shape parameter and AMM effects can alter the coupled EOS-TOV equations resulting in an increase in the final mass and a decrease in the final equatorial radius without collapsing the core into a black hole and without violating causality constraints, we find maximum mass values in the range: 2.3 Solar Masses < M < 2.7 Solar Masses. These states are consistent with some astrophysical high mass magnetar/pulsar and gravity wave systems which may provide evidence for a core that has undergone a quark-gluon phase transition such as PSR 0943+10 and the secondary from the GW 190814 event.

The Sanford Underground Research Facility (SURF) has been operating since 2007 supporting underground research in rare-process physics, as well as offering research opportunities in other disciplines. SURF laboratory facilities include a Surface Campus as well as campuses at the 4850-foot level (1500 m, 4300 m.w.e.) that host a range of significant physics experiments, including those studying dark matter, neutrino properties, and nuclear astrophysics topics. SURF is also home to the Long-Baseline Neutrino Facility (LBNF) that will host the international Deep Underground Neutrino Experiment (DUNE). SURF offers an ultra-low background environment, low-background assay capabilities, and electroformed copper is produced at the facility. SURF is proposing additional underground space on the 4850L and 7400L (2300 m, 6500 m.w.e.), and initial engineering designs have been completed. SURF is a dedicated research facility with significant expansion capability, and applications from new experiments are welcome.

The search for particle-like dark matter with meV-to-GeV masses has developed rapidly in the past few years. We summarize the science case for these searches, the recent progress, and the exciting upcoming opportunities. Funding for Research and Development and a portfolio of small dark matter projects will allow the community to capitalize on the substantial recent advances in theory and experiment and probe vast regions of unexplored dark-matter parameter space in the coming decade.

In this paper, with the presence of swampland conjecture, we use the complex form of a scalar field and investigate the anisotropic constant-roll of an inflationary scenario. Recently, the complex quintessence field has been used to describe the accelerating expansion of the universe which has had interesting results connecting with various conditions. Therefore, the Lagrangian density of the quintessence field leads to obtaining the equations of the complex scalar field. Also, we use the anisotropic constant-roll conditions and the field equation and calculate the exact solutions for some cosmology parameters such as Hubble parameter and potential. Using the exact solution of potential and refined swampland conditions, we plot some figures. In that case, the figures lead us to have challenges between three concepts as complex quintessence, anisotropic constant-roll conditions, and swampland conjectures. Then we discuss their compatibility and incompatibility with corresponding results. Finally, we will analyze the complex quintessence in examining the inflationary scenario.

The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-based scenarios are forecasted to probe many square decades of unexplored dark matter parameter space below 5 GeV, covering over 6 decades in mass: 1-100 eV for dark photons and axion-like particles, 1-100 MeV for dark-photon-coupled light dark matter, and 0.05-5 GeV for nucleon-coupled dark matter. They will reach the neutrino fog in the 0.5-5 GeV mass range and test a variety of benchmark models and sharp targets. The novel directions involve greater departures from current SuperCDMS technology but promise even greater reach in the long run, and their development must begin now for them to be available in a timely fashion. The experienced-based upgrade scenarios rely mainly on dramatic improvements in detector performance based on demonstrated scaling laws and reasonable extrapolations of current performance. Importantly, these improvements in detector performance obviate significant reductions in background levels beyond current expectations for the SuperCDMS SNOLAB experiment. Given that the dominant limiting backgrounds for SuperCDMS SNOLAB are cosmogenically created radioisotopes in the detectors, likely amenable only to isotopic purification and an underground detector life-cycle from before crystal growth to detector testing, the potential cost and time savings are enormous and the necessary improvements much easier to prototype.

We address a method of limiting neutron-mirror neutron mixing ($\epsilon_{n n'}$) by analyzing its effect on neutron star (NS) heating. This method employs observational bounds on the surface temperature of NSs to constrain $\epsilon_{n n'}$. The obtained bound is so stringent that it would exclude any discovery of $n-n'$ oscillation by the planned terrestrial experiments. Motivated by this last crucially important fact, we critically analyze this suggestion and note a very interesting new effect present in nearly exact mirror models, which significantly affect this bound. The new element in our discussion is that in mirror models there is the $\beta$ decay $n' \to p'+ e' +\bar{\nu}'_{e}$, which creates a cloud of mirror particles $n'$, $p'$, $e'$ and $D'$ inside the NS. The $e'$ can "rob" the energy generated by the $n \to n'$ transition from the NS. This is achieved via $e-e'$ scattering enabled by the presence of a (minute) milli-charge in mirror particles. This energy is emitted as unobserved mirror photons via fast mirror bremsstrahlung. This will lead to relaxing of the stringent bounds on $\epsilon_{nn'}$.

We study the Vainshtein screening mechanism in Horndeski theories in the presence of a scalar field $\phi$ nonminimally and kinetically coupled to ordinary matter field. A general interacting Lagrangian describing this coupling is characterized by energy transfer $f_1$ and momentum exchange $f_2$. For a spherically symmetric configurations on top of the cosmological background, we investigate the static perturbations in linear and nonlinear regimes with respect to the scalar field perturbation. In the former regime, the PPN parameter generally deviates from unity as long as the matter coupling or $G_{4,\phi}$ exists. On the other hand, in the latter regime, we show that the nonlinear self interaction term of scalar field successfully activates the Vainshtein mechanism even in the presence of the couplings $f_1$ and $f_2$. The gravitational potentials recover the Newtonian behavior deep inside the Vainshtein radius. The bounds on coupling terms not to substantially change the Vainshtein radius are also given.

We apply Bayesian approach to construct a large number of minimally constrained equations of state (EoSs) and study their correlations with a few selected properties of neutron star (NS). Our set of minimal constraints includes a few basic properties of saturated nuclear matter and low density pure neutron matter EoS which is obtained from a precise next-to-next-to-next-to-leading order (N$^{3}$LO) calculation in chiral effective field theory. The tidal deformability and radius of NS with mass $1$-$2M_\odot$ are found to be strongly correlated with the pressure of $\beta$-equilibrated matter as well as with the symmetry energy at densities higher than the saturation density ($\rho_0 = 0.16$ fm$^{-3}$) in a nearly model independent manner. These correlations are employed to parametrize the pressure for $\beta$-equilibrated matter, around 2$\rho_0$, as a function of neutron star mass and the corresponding tidal deformability. The maximum mass of neutron star is also strongly correlated with the pressure of symmetric and $\beta$-equilibrated matter at densities $\sim$ 4.5$\rho_0$.

Axion-photon conversion is a prime mechanism to detect axion-like particles that share a coupling to the photon. We point out that in the vicinity of neutron stars with strong magnetic fields, magnetars, the effective photon mass receives comparable but opposite contributions from free electrons and the radiation field. This leads to an energy-dependent resonance condition for conversion that can be met for arbitrary light axions and leveraged when using systems with detected radio component. Using the magnetar SGR J1745-2900 as an exemplary source, we demonstrate that sensitivity to $|g_{a\gamma}| \sim 10^{-12}\,\rm{GeV^{-1}}$ or better can be gained for $m_a \lesssim 10^{-6}\,\rm eV$, with the potential to improve current constraints on the axion-photon coupling by more than one order of magnitude over a broad mass range. With growing insights into the physical conditions of magnetospheres of magnetars, the method hosts the potential to become a serious competitor to future experiments such as ALPS-II and IAXO in the search for axion-like particles.

Temperature and velocity-dependent $^1$S$_0$ pairing gaps, chemical potentials and entrainment matrix in dense homogeneous neutron-proton superfluid mixtures constituting the outer core of neutron stars, are determined fully self-consistently by solving numerically the time-dependent Hartree-Fock-Bogoliubov equations over the whole range of temperatures and flow velocities for which superfluidity can exist. Calculations have been made for $npe\mu$ in beta-equilibrium using the Brussels-Montreal functional BSk24. The accuracy of various approximations is assessed and the physical meaning of the different velocities and momentum densities appearing in the theory is clarified. Together with the unified equation of state published earlier, the present results provide consistent microscopic inputs for modeling superfluid neutron-star cores.

The concept of a small, single-layer water Cherenkov detector, with three photomultiplier tubes (PMT), placed at its bottom in a $120^{\circ}$ star configuration (\emph{Mercedes} WCD) is presented. The PMTs are placed near the lateral walls of the stations with an adjustable inclination and may be installed inside or outside the water volume. To illustrate the technical viability of this concept and obtain a first-order estimation of its cost, an engineering design was elaborated. The sensitivity of these stations to low energy EAS electrons, photons and muons is discussed, both in compact and sparse array configurations. It is shown that the analysis of the intensity and time patterns of the PMT signals, using Machine Learning techniques, enables the tagging of muons, achieving an excellent gamma/hadron discrimination for TeV showers. This concept minimises the station production and maintenance costs, allowing for a highly flexible and fast installation. Mercedes WCDs are thus well-suited for use in high-altitude large gamma-ray observatories covering an extended energy range from the low energies, closing the gap between satellite and ground-based measurements, to very high energy regions, beyond the PeV scale.