Papers for Monday, May 31 2021

In previous work we showed that a semi-empirical model in which galaxies in host dark matter haloes are assigned stellar masses via a stellar mass-halo mass (SMHM) relation and sizes ($R_e$) via a linear and tight $R_e-R_h$ relation, can faithfully reproduce the size function of local SDSS central galaxies and the strong size evolution of massive galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$). In this third paper of the series, we focus on the population of satellite MGs. We find that without any additional calibration and irrespective of the exact SMHM relation, fraction of quenched galaxies or level of stellar stripping, the same model is able to reproduce the local size function of quiescent satellite MGs in SDSS. In addition, the same model can reproduce the puzzling weak dependence of mean size on host halo mass for both central and satellite galaxies. The model also matches the size function of starforming satellite MGs, after assuming that some of them transform into massive lenticulars in a few Gyr after infalling in the group/cluster environment. However, the vast majority of satellite lenticulars is predicted to form before infall. The $R_e-R_h$ appears to be fundamental to connect galaxies and their host haloes.

We present results from the Large Adaptive optics Survey for Substellar Objects (LASSO), where the goal is to directly image new substellar companions (<70 M$_{Jup}$) at wide orbital separations ($\gtrsim$50 AU) around young ($\lesssim$300 Myrs), nearby (<100 pc), low-mass ($\approx$0.1-0.8 M$_{\odot}$) stars. We report on 427 young stars imaged in the visible (i') and near-infrared (J or H) simultaneously with Robo-AO on the Kitt Peak 2.1-m telescope and later the Maunakea University of Hawaii 2.2-m telescope. To undertake the observations, we commissioned a new infrared camera for Robo-AO that uses a low-noise high-speed SAPHIRA avalanche photodiode detector. We detected 121 companion candidates around 111 stars, of which 62 companions are physically associated based on Gaia DR2 parallaxes and proper motions, another 45 require follow-up observations to confirm physical association, and 14 are background objects. The companion separations range from 2-1101 AU and reach contrast ratios of 7.7 magnitudes in the near infrared compared to the primary. The majority of confirmed and pending candidates are stellar companions, with ~5 being potentially substellar and requiring follow-up observations for confirmation. We also detected a 43$\pm$9 M$_{Jup}$ and an 81$\pm$5 M$_{Jup}$ companion that were previously reported. We found 34 of our targets have acceleration measurements detected using Hipparcos-Gaia proper motions. Of those, 58$^{+12}_{-14}$% of the 12 stars with imaged companion candidates have significant accelerations ($\chi^2 >11.8$), while only 23$^{+11}_{-6}$% of the remaining 22 stars with no detected companion have significant accelerations. The significance of the acceleration decreases with increasing companion separation. These young accelerating low-mass stars with companions will eventually yield dynamical masses with future orbit monitoring.

We investigate how cosmic web environment impacts the average properties of central galaxies in the Sloan Digital Sky Survey (SDSS). We analyse how the average specific star-formation rate, stellar age, metallicity and element abundance ratio [$\alpha$/Fe] of SDSS central galaxies depend on distance from the cosmic web nodes, walls and filaments identified by DisPerSE. In our approach we control for galaxy stellar mass and local density differentiated between field and group environment. Our results confirm the known trend whereby galaxies exhibit lower specific star-formation rates with decreasing distance to the cosmic web features. Furthermore, we show that centrals closer to either nodes, walls or filaments are on average older, metal richer and $\alpha$-enhanced compared to their equal mass counterparts at larger distances. The identified property gradients appear to have the same amplitude for central galaxies in the field as for those in groups. Our findings support a cosmic web quenching that stems from nurture effects, such as ram pressure stripping and strangulation, and/or nature effects linked to the intrinsic properties of the cosmic web.

The Taurus-Auriga complex is the prototypical low-mass star forming region, and provides a unique testbed of the star formation process, which left observable imprints on the spatial, kinematic, and temporal structure of its stellar population. Taurus's rich observational history has uncovered peculiarities that suggest a complicated star forming event, such as members at large distances from the molecular clouds and evidence of an age spread. With Gaia, an in-depth study of the Taurus census is possible to confirm membership, identify substructure, and reconstruct its star formation history. We have compiled an expansive census of the greater Taurus region, identifying spatial subgroups and confirming that Taurus is substructured across stellar density. There are two populations of subgroups: clustered groups near the clouds and sparse groups spread throughout the region. The sparse groups comprise Taurus's distributed population, which is on average older than the population near the clouds, and hosts sub-populations up to 15 Myr old. The ages of the clustered groups increase with distance, suggesting that the current star formation was triggered from behind. Still, the region is kinematically coherent, and its velocity structure reflects an initial turbulent spectrum similar to Larson's Law that has been modified by dynamical relaxation. Overall, Taurus has a complicated star formation history, with at least two epochs of star formation featuring both clustered and distributed modes. Given the correlations between age and spatial distribution, Taurus might be part of a galaxy-scale star forming event that can only begin to be understood in the Gaia era.

The formation of the most massive quasars observed at high redshifts requires extreme inflows of gas down to the length scales of the central compact object. Here, we estimate the maximum inflow rate allowed by gravity down to the surface of supermassive stars, the possible progenitors of these supermassive black holes. We use the continuity equation and the assumption of free-fall to derive maximum allowed inflow rates for various density profiles. We apply our approach to the mass-radius relation of rapidly accreting supermassive stars to estimate an upper limit to the accretion rates allowed during the formation of these objects. We find that the maximum allowed rate $\dot M_{\rm max}$ is given uniquely by the compactness of the accretor. For the compactness of rapidly accreting supermassive stars, $\dot M_{\rm max}$ is related to the stellar mass $M$ by a power-law $\dot M_{\rm max}\propto M^{3/4}$. The rates of atomically cooled halos (0.1 -- 10 M$_\odot$ yr$^{-1}$) are allowed as soon as $M\gtrsim1$ M$_\odot$. The largest rates expected in galaxy mergers ($10^4-10^5$ M$_\odot$ yr$^{-1}$) become accessible once the accretor is supermassive ($M\gtrsim10^4$ M$_\odot$). These results suggest that supermassive stars can accrete up to masses $>10^6$ M$_\odot$ before they collapse via the general-relativistic instability. At such masses, the collapse is expected to lead to the direct formation of a supermassive black hole even within metal-rich gas, resulting in a black hole seed that is significantly heavier than in conventional direct collapse models for atomic cooling halos.

We use the red tunable filter of the instrument OSIRIS at GTC to obtain the pseudo spectra of emission line sources in the OTELO field. From these pseudo spectra we identify the objects with H$\beta$ emission. As the resolution of the pseudo spectra allows us to separate H$\beta$ from [OIII], we are able to derive the H$\beta$ flux without contamination from its adjacent line. Using data from the extended OTELO catalogue we discriminate AGNs and study the star formation rate, the stellar mass and the morphology of the star-forming galaxies. We find that our sample is located on the main sequence of star-forming galaxies. The sources are morphologically classified mostly as disc-like galaxies (76%) and 90% of the sample as low-mass galaxies ($M_*<10^{10}\;\mathrm{M}_\odot$). The low-mass star-forming galaxies at $z \sim 0.9$ detected by OTELO present similar properties to low-mass star-forming galaxies in the local universe, suggesting that this kind of object have not a favorite epoch of formation and star formation enhancement from $z \sim 1$ to now. Our sample of 40 H$\beta$ star-forming galaxies includes the faintest H$\beta$ emitters detected so far. This allows us to constrain the faint end of the luminosity function for the H$\beta$ line alone with a minimum luminosity of $\log L = 39 \;\mathrm{erg\,s}^{-1}$, a hundred times fainter than previous surveys. The dust-corrected OTELO H$\beta$ luminosity function established the faint-end slope as $\alpha =-1.36\pm 0.15$. We increased the scope of the analysis to the bright end by adding ancillary data from the literature, non dust-corrected in this case. The obtained slope for this extended luminosity function is $\alpha = -1.43\pm 0.12$.

An analysis of 7.5 years of data in the high-energy starting event sample has been recently published by the IceCube collaboration. The hottest spot in a search for neutrino sources was found far above the Galactic plane and is thus, at first sight, difficult to reconcile with a Galactic origin. In this work, we calculate the cosmic ray (CR) density around nearby, young supernova remnants assuming anisotropic diffusion. Combining the obtained CR densities with the matter distribution deduced from extinction maps, we find two prominent hot spots: One of them is close to the most significant point in the IceCube search for point sources, with an intensity corresponding to two neutrino events.

The age-metallicity relation is a fundamental tool for constraining the chemical evolution of the Galactic disc. In this work we analyse the observational properties of this relation using binary stars that have not interacted consisting of a white dwarf - from which we can derive the total age of the system - and a main sequence star - from which we can derive the metallicity as traced by the [Fe/H] abundances. Our sample consists of 46 widely separated, but unresolved spectroscopic binaries identified within the Sloan Digital Sky Survey, and 189 white dwarf plus main sequence common proper motion pairs identified within the second data release of Gaia. This is currently the largest white dwarf sample for which the metallicity of their progenitors have been determined. We find a flat age-metallicity relation displaying a scatter of [Fe/H] abundances of approximately 0.5 dex around the solar metallicity at all ages. This independently confirms the lack of correlation between age and metallicity in the solar neighbourhood that is found in previous studies focused on analysing single main sequence stars and open clusters.

X-ray bursts have recently been discovered in the Cepheids $\delta$ Cep and $\beta$ Dor modulated by the pulsation cycle. We have obtained an observation of the Cepheid $\eta$ Aql with the XMM-Newton satellite at the phase of maximum radius, the phase at which there is a burst of X-rays in $\delta$ Cep. No X-rays were seen from the Cepheid $\eta$ Aql at this phase, and the implications for Cepheid upper atmospheres are discussed. We have also used the combination of X-ray sources and Gaia and 2MASS data to search for a possible grouping around the young intermediate mass Cepheid. No indication of such a group was found.

Solar Radio Bursts (SRBs) are generally observed in dynamic spectra and have five major spectral classes, labelled Type I to Type V depending on their shape and extent in frequency and time. Due to their complex characterisation, a challenge in solar radio physics is the automatic detection and classification of such radio bursts. Classification of SRBs has become fundamental in recent years due to large data rates generated by advanced radio telescopes such as the LOw-Frequency ARray, (LOFAR). Current state-of-the-art algorithms implement the Hough or Radon transform as a means of detecting predefined parametric shapes in images. These algorithms achieve up to 84% accuracy, depending on the Type of radio burst being classified. Other techniques include procedures that rely on Constant-FalseAlarm-Rate detection, which is essentially detection of radio bursts using a de-noising and adaptive threshold in dynamic spectra. It works well for a variety of different Types of radio bursts and achieves an accuracy of up to 70%. In this research, we are introducing a methodology named You Only Look Once v2 (YOLOv2) for solar radio burst classification. By using Type III simulation methods we can train the algorithm to classify real Type III solar radio bursts in real-time at an accu

We have fit the far-ultraviolet (FUV) to mid-infrared (MIR) spectral energy distributions (SEDs) for several nearby galaxies ($<$ 20 Mpc). Global, radial, and local photometric measurements are explored to better understand how SED-derived star formation histories (SFHs) and classic star formation rate (SFR) tracers manifest at different scales. Surface brightness profiles and radial SED fitting provide insight into stellar population gradients in stellar discs and haloes. A double exponential SFH model is used in the SED fitting to better understand the distributions of young vs. old populations throughout these galaxies. Different regions of a galaxy often have undergone very different SFHs, either in strength, rate, timing, or some combination of all these factors. An analysis of individual stellar complexes within these galaxies shows a relationship between the ages of stellar clusters and how these clusters are distributed throughout the galaxy. These star formation properties are presented alongside previously published HI observations to provide a holistic picture of a small sample of nearby star-forming galaxies. The results presented here show that there is a wide variety of star formation gradients and average stellar age distributions that can manifest in a $\Lambda$CDM universe.

We explore the presence of active galactic nuclei (AGN)/black hole (BH) in Green Pea galaxies (GPs), motivated by the presence of high ionization emission lines such as HeII and [NeIII] in their optical spectra. In order to identify AGN candidates, we used mid-infrared (MIR) photometric observations from the all-sky Wide-field Infrared Survey Explorer (WISE) mission for a sample of 516 GPs. We select 58 GPs as candidate AGN based on a stringent 3-band WISE color diagnostic. Using multi-epoch photometry of W1 and W2 bands from the WISE/NEOWISE-R observations, we find 38 GPs showing significant variability in both the WISE bands. Four of these were selected as AGN by the WISE 3-band color diagnostic as well. Interestingly, we find a high fraction of MIR variable sources among GPs which demonstrates the uniqueness and importance of studying these extreme objects. Through this work, we demonstrate that photometric variability is a promising tool to select AGN that may be missed by other selection techniques (including optical emission-line ratios and X-ray emission) in star-formation dominated, low-mass, low-metallicity galaxies.

Time series of spectroscopic, speckle-interferometric, and optical long-baseline-interferometric observations confirm that $\nu$ Gem is a hierarchical triple system. It consists of an inner binary composed of two B-type stars and an outer classical Be star. Several photospheric spectral lines of the inner components were disentangled, revealing two stars with very different rotational broadening ($\sim$260 and $\sim$140 kms$^{-1}$, respectively), while the photospheric lines of the Be star remain undetected. From the combined spectroscopic and astrometric orbital solution it is not possible to unambiguously cross-identify the inner astrometric components with the spectroscopic components. In the preferred solution based on modeling of the disentangled line profiles, the inner binary is composed of two stars with nearly identical masses of 3.3 M$_\odot$ and the more rapidly rotating star is the fainter one. These two stars are in a marginally elliptical orbit ($e$ = 0.06) about each other with a period of 53.8 d. The third star also has a mass of 3.3 M$_\odot$ and follows a more eccentric ($e$ = 0.24) orbit with a period of 19.1 yr. The two orbits are co-directional and, at inclinations of 79$^{\circ}$ and 76$^{\circ}$ of the inner and the outer orbit, respectively, about coplanar. No astrometric or spectroscopic evidence could be found that the Be star itself is double. The system appears dynamically stable and not subject to eccentric Lidov-Kozai oscillations. After disentangling, the spectra of the components of the inner binary do not exhibit peculiarities that would be indicative of past interactions. Motivations for a wide range of follow-up studies are suggested.

Axion-like particles (ALPs) are promising dark matter candidates. A classical field description is typically employed, motivated by large phase space occupation numbers. Here we show that such a description is accompanied by a quantum effect: squeezing due to gravitational self-interactions. For a typical QCD axion today, the onset of squeezing is reached on microsecond-scales and grows over millennia. Thus within the usual models based on the classical Schr\"odinger-Poisson equation, a type of Gross-Pitaevskii equation, any viable ALP is nonclassical. We also show that squeezing may be relevant on scales of axion haloscopes, or within galactic solitonic cores. Conversely, our results highlight the incompleteness and limitations of the typically employed classical single field description of ALPs.

Feedback from active galactic nuclei (AGN) has proven to be a critical ingredient in the current picture of galaxy assembly and growth. However, observational constraints on AGN-driven outflows face technical challenges and as a result, the cold molecular gas outflow properties of type-1 AGN are not well known. We present new IRAM NOrthern Extended Milimeter Array (NOEMA) observations of CO (1-0) in F07599+6508, Z11598$-$0112, F13342+3932, and PG1440+356, all nearby type-1 AGN and ultraluminous infrared galaxies (ULIRGs). We achieve spatial resolution of 1-3 arcsec corresponding to physical scales of 2--8~kpc and spectral resolution of 15-60 km s$^{-1}$ which enables updated CO (1-0) redshifts and a detailed morphological view of the cold molecular gas in these sources. The CO (1-0) luminosities, $L_{CO}^{\prime}$, are in the range 2-12$\times 10^9$ K km s$^{-1}$ pc$^{2}$ and inferred molecular gas masses, $M$(H$_2$), are in the range 2-9$\times 10^9$ M$_{\odot}$. The velocity fields and gas distributions do not unambiguously identify any of these sources as having outflows. However, Z11598$-$0112 has signs of infalling material and after the subtraction of a rotating disk model PG 1440+356 shows complex kinematics in the residuals that may indicate an outflow or warped disk.

Stellar activity poses one of the main obstacles for the detection and characterisation of small exoplanets around cool stars, as it can induce radial velocity (RV) signals that can hide or mimic the presence of planetary companions. Several indicators of stellar activity are routinely used to identify activity-related signals in RVs, but not all indicators trace exactly the same activity effects, nor are any of them always effective in all stars. We evaluate the performance of a set of spectroscopic activity indicators for M dwarf stars with different masses and activity levels with the aim of finding a relation between the indicators and stellar properties. In a sample of 98 M dwarfs observed with CARMENES, we analyse the temporal behaviour of RVs and nine spectroscopic activity indicators: cross-correlation function (CCF) full-width-at-half-maximum (FWHM), contrast, and bisector inverse slope (BIS), chromatic index (CRX), differential line width (dLW), and indices of the chromospheric lines H$\alpha$ and calcium infrared triplet. A total of 56 stars of the initial sample show periodic signals related to activity in at least one of these ten parameters. RV is the parameter for which most of the targets show an activity-related signal. CRX and BIS are effective activity tracers for the most active stars in the sample, especially stars with a relatively high mass, while for less active stars, chromospheric lines perform best. FWHM and dLW show a similar behaviour in all mass and activity regimes, with the highest number of activity detections in the low-mass, high-activity regime. Most of the targets for which we cannot identify any activity-related signals are stars at the low-mass end of the sample. These low-mass stars also show the lowest RV scatter, which indicates that ultracool M dwarfs could be better candidates for planet searches than earlier types, which show larger RV jitter.

Outflows from super-massive black holes (SMBHs) play an important role in the co-evolution of themselves, their host galaxies, and the larger scale environments. Such outflows are often characterized by emission and absorption lines in various bands and in a wide velocity range blueshifted from the systematic redshift of the host quasar. In this paper, we report a strong broad line region (BLR) outflow from the z~4.7 quasar BR 1202-0725 based on the high-resolution optical spectrum taken with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph installed on the 6.5m Magellan/Clay telescope, obtained from the `Probing the He II re-Ionization ERa via Absorbing C IV Historical Yield' (HIERACHY) project. This rest-frame ultraviolet (UV) spectrum is characterized by a few significantly blueshifted broad emission lines from high ions; the most significant one is the C IV line at a velocity of -6500 km/s relative to the H{\alpha} emission line, which is among the highest velocity BLR outflows in observed quasars at z > 4. The measured properties of UV emission lines from different ions, except for O I and Ly{\alpha}, also follow a clear trend that higher ions tend to be broader and outflow at higher average velocities. There are multiple C IV and Si IV absorbing components identified on the blue wings of the corresponding emission lines, which may be produced by either the outflow or the intervening absorbers.

We quantitatively investigate the dependence of central galaxy HI mass ($M_{\rm HI}$) on the stellar mass ($M_\ast$), halo mass ($M_{\rm h}$), star formation rate (SFR), and central stellar surface density within 1 kpc ($\Sigma_1$), taking advantage of the HI spectra stacking technique using both the Arecibo Fast Legacy ALFA Survey and the Sloan Digital Sky Survey. We find that the shapes of $M_{\rm HI}$-$M_{\rm h}$ and $M_{\rm HI}$-$M_\ast$ relations are remarkably similar for both star-forming and quenched galaxies, with massive quenched galaxies having constantly lower HI masses of around 0.6 dex. This similarity strongly suggests that neither halo mass nor stellar mass is the direct cause of quenching, but rather the depletion of HI reservoir. While the HI reservoir for low-mass galaxies of $M_\ast<10^{10.5}M_\odot$ strongly increases with $M_{\rm h}$, more massive galaxies show no significant dependence of $M_{\rm HI}$ on $M_{\rm h}$, indicating the effect of halo to determine the smooth cold gas accretion. We find that the star formation and quenching of central galaxies are directly regulated by the available HI reservoir, with an average relation of ${\rm SFR}\propto M_{\rm HI}^{2.75}/M_\ast^{0.40}$, implying a quasi-steady state of star formation. We further confirm that galaxies are depleted of their HI reservoir once they drop off the star-formation main sequence and there is a very tight and consistent correlation between $M_{\rm HI}$ and $\Sigma_1$ in this phase, with $M_{\rm HI}\propto\Sigma_1^{-2}$. This result is in consistent with the compaction-triggered quenching scenario, with galaxies going through three evolutionary phases of cold gas accretion, compaction and post-compaction, and quenching.

Short-orbit gas giant planet formation/evolution mechanisms are still not well understood. One promising pathway to discriminate between mechanisms is to constrain the occurrence rate of these peculiar exoplanets at the earliest stage of the system's life. However, a major limitation when studying newly born stars is stellar activity. This cocktail of phenomena triggered by fast rotation, strong magnetic fields and complex internal dynamics, especially present in very young stars, compromises our ability to detect exoplanets. In this paper, we investigated the limitations of such detections in the context of already acquired data solely using radial velocity data acquired with a non-stabilised spectrograph. We employed two strategies: Doppler Imaging and Gaussian Processes and could confidently detect Hot Jupiters with semi-amplitude of 100 $m.s^{-1}$ buried in the stellar activity. We also showed the advantages of the Gaussian Process approach in this case. This study serves as a proof of concept to identify potential candidates for follow-up observations or even discover such planets in legacy datasets available to the community.

We present reconstructed convergence maps, \textit{mass maps}, from the Dark Energy Survey (DES) third year (Y3) weak gravitational lensing data set. The mass maps are weighted projections of the density field (primarily dark matter) in the foreground of the observed galaxies. We use four reconstruction methods, each is a \textit{maximum a posteriori} estimate with a different model for the prior probability of the map: Kaiser-Squires, null B-mode prior, Gaussian prior, and a sparsity prior. All methods are implemented on the celestial sphere to accommodate the large sky coverage of the DES Y3 data. We compare the methods using realistic $\Lambda$CDM simulations with mock data that are closely matched to the DES Y3 data. We quantify the performance of the methods at the map level and then apply the reconstruction methods to the DES Y3 data, performing tests for systematic error effects. The maps are compared with optical foreground cosmic-web structures and are used to evaluate the lensing signal from cosmic-void profiles. The recovered dark matter map covers the largest sky fraction of any galaxy weak lensing map to date.

In this work we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, MagLim, and a sample of luminous red galaxies selected with the redMaGiC algorithm. MagLim / redMaGiC sample contains over 10 million / 2.5 million galaxies and is divided into six / five photometric redshift bins spanning the range $z\in[0.20,1.05]$ / $z\in[0.15,0.90]$. Both samples cover 4143 deg$^2$ over which we perform our analysis blind, measuring the angular correlation function with a S/N $\sim 63$ for both samples. In a companion paper (DES Collaboration et al. 2021)), these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy-galaxy lensing of each sample to place cosmological constraints with a 3$\times$2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using log-normal mocks, showing that our decontamination procedure induces biases no greater than $0.5\sigma$ in the $(\Omega_m, b)$ plane, where $b$ is galaxy bias. We demonstrate that failure to remove the artificial clustering would introduce strong biases up to $\sim 7 \sigma$ in $\Omega_m$ and of more than $4 \sigma$ in galaxy bias.

We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg$^2$. These galaxy-galaxy measurements are used in the DES Y3 3$\times$2pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redMaGic sample, which span the redshift range $\sim 0.2-1$ with 10.7 M and 2.6 M galaxies respectively. For the source catalog, we use the Metacalibration shape sample, consisting of $\simeq$100 M galaxies separated into 4 tomographic bins. Our galaxy-galaxy lensing estimator is the mean tangential shear, for which we obtain a total S/N of $\sim$148 for MagLim ($\sim$120 for redMaGic), and $\sim$67 ($\sim$55) after applying the scale cuts of 6 Mpc/$h$. Thus we reach percent-level statistical precision, which requires that our modeling and systematic-error control be of comparable accuracy. The tangential shear model used in the 3$\times$2pt cosmological analysis includes lens magnification, a five-parameter intrinsic alignment model (TATT), marginalization over a point-mass to remove information from small scales and a linear galaxy bias model validated with higher-order terms. We explore the impact of these choices on the tangential shear observable and study the significance of effects not included in our model, such as reduced shear, source magnification and source clustering. We also test the robustness of our measurements to various observational and systematics effects, such as the impact of observing conditions, lens-source clustering, random-point subtraction, scale-dependent Metacalibration responses, PSF residuals, and B-modes.

Using the first three years of data from the Dark Energy Survey, we use ratios of small-scale galaxy-galaxy lensing measurements around the same lens sample to constrain source redshift uncertainties, intrinsic alignments and other nuisance parameters of our model. Instead of using a simple geometric approach for the ratios, we use the full modeling of the galaxy-galaxy lensing measurements, including the corresponding integration over the power spectrum and the contributions from intrinsic alignments and lens magnification. We perform extensive testing of the small-scale shear ratio (SR) modeling by studying the impact of different effects such as the inclusion of baryonic physics, non-linear biasing, halo occupation distribution descriptions and lens magnification, among others, and using realistic $N$-body simulations. We validate the robustness of our constraints in the data by using two independent lens samples, and by deriving constraints using the corresponding large-scale ratios for which the modeling is simpler. The DES Y3 results demonstrate how the ratios provide significant improvements in constraining power for several nuisance parameters in our model, especially on source redshift calibration and intrinsic alignments (IA). For source redshifts, SR improves the constraints from the prior by up to 38\% in some redshift bins. Such improvements, and especially the constraints it provides on IA, translate to tighter cosmological constraints when SR is combined with cosmic shear and other 2pt functions. In particular, for the DES Y3 data, SR improves $S_8$ constraints from cosmic shear by up to 31\%, and for the full combination of probes (3$\times$2pt) by up to 10\%. The shear ratios presented in this work are used as an additional likelihood for cosmic shear, 2$\times$2pt and the full 3$\times$2pt in the fiducial DES Y3 cosmological analysis.

This work, together with its companion paper, Secco and Samuroff et al. (2021), presents the Dark Energy Survey Year 3 cosmic shear measurements and cosmological constraints based on an analysis of over 100 million source galaxies. With the data spanning 4143 deg$^2$ on the sky, divided into four redshift bins, we produce the highest significance measurement of cosmic shear to date, with a signal-to-noise of 40. We conduct a blind analysis in the context of the $\Lambda$CDM model and find a 3% constraint of the clustering amplitude, $S_8\equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.759^{+0.025}_{-0.023}$. A $\Lambda$CDM-Optimized analysis, which safely includes smaller scale information, yields a 2% precision measurement of $S_8= 0.772^{+0.018}_{-0.017}$ that is consistent with the fiducial case. The two low-redshift measurements are statistically consistent with the Planck Cosmic Microwave Background result, however, both recovered $S_8$ values are lower than the high-redshift prediction by $2.3\sigma$ and $2.1\sigma$ ($p$-values of 0.02 and 0.05), respectively. The measurements are shown to be internally consistent across redshift bins, angular scales and correlation functions. The analysis is demonstrated to be robust to calibration systematics, with the $S_8$ posterior consistent when varying the choice of redshift calibration sample, the modeling of redshift uncertainty and methodology. Similarly, we find that the corrections included to account for the blending of galaxies shifts our best-fit $S_8$ by $0.5\sigma$ without incurring a substantial increase in uncertainty. We examine the limiting factors for the precision of the cosmological constraints and find observational systematics to be subdominant to the modeling of astrophysics. Specifically, we identify the uncertainties in modeling baryonic effects and intrinsic alignments as the limiting systematics.

This work and its companion paper, Amon et al. (2021), present cosmic shear measurements and cosmological constraints from over 100 million source galaxies in the Dark Energy Survey (DES) Year 3 data. We constrain the lensing amplitude parameter $S_8\equiv\sigma_8\sqrt{\Omega_\textrm{m}/0.3}$ at the 3% level in $\Lambda$CDM: $S_8=0.759^{+0.025}_{-0.023}$ (68% CL). Our constraint is at the 2% level when using angular scale cuts that are optimized for the $\Lambda$CDM analysis: $S_8=0.772^{+0.018}_{-0.017}$ (68% CL). With cosmic shear alone, we find no statistically significant constraint on the dark energy equation-of-state parameter at our present statistical power. We carry out our analysis blind, and compare our measurement with constraints from two other contemporary weak-lensing experiments: the Kilo-Degree Survey (KiDS) and Hyper-Suprime Camera Subaru Strategic Program (HSC). We additionally quantify the agreement between our data and external constraints from the Cosmic Microwave Background (CMB). Our DES Y3 result under the assumption of $\Lambda$CDM is found to be in statistical agreement with Planck 2018, although favors a lower $S_8$ than the CMB-inferred value by $2.3\sigma$ (a $p$-value of 0.02). This paper explores the robustness of these cosmic shear results to modeling of intrinsic alignments, the matter power spectrum and baryonic physics. We additionally explore the statistical preference of our data for intrinsic alignment models of different complexity. The fiducial cosmic shear model is tested using synthetic data, and we report no biases greater than 0.3$\sigma$ in the plane of $S_8\times\Omega_\textrm{m}$ caused by uncertainties in the theoretical models.

We constrain cosmological parameters and galaxy-bias parameters using the combination of galaxy clustering and galaxy-galaxy lensing measurements from the Dark Energy Survey Year-3 data. We describe our modeling framework and choice of scales analyzed, validating their robustness to theoretical uncertainties in small-scale clustering by analyzing simulated data. Using a linear galaxy bias model and redMaGiC galaxy sample, we obtain constraints on the matter content of the universe to be $\Omega_{\rm m} = 0.325^{+0.033}_{-0.034}$. We also implement a non-linear galaxy bias model to probe smaller scales that includes parameterizations based on hybrid perturbation theory, and find that it leads to a 17% gain in cosmological constraining power. Using the redMaGiC galaxy sample as foreground lens galaxies, we find the galaxy clustering and galaxy-galaxy lensing measurements to exhibit significant signals akin to decorrelation between galaxies and mass on large scales, which is not expected in any current models. This likely systematic measurement error biases our constraints on galaxy bias and the $S_8$ parameter. We find that a scale-, redshift- and sky-area-independent phenomenological decorrelation parameter can effectively capture this inconsistency between the galaxy clustering and galaxy-galaxy lensing. We perform robustness tests of our methodology pipeline and demonstrate stability of the constraints to changes in the theory model. After accounting for this decorrelation, we infer the constraints on the mean host halo mass of the redMaGiC galaxies from the large-scale bias constraints, finding the galaxies occupy halos of mass approximately $1.5 \times 10^{13} M_{\odot}/h$.

Two of the most sensitive probes of the large scale structure of the universe are the clustering of galaxies and the tangential shear of background galaxy shapes produced by those foreground galaxies, so-called galaxy-galaxy lensing. Combining the measurements of these two two-point functions leads to cosmological constraints that are independent of the galaxy bias factor. The optimal choice of foreground, or lens, galaxies is governed by the joint, but conflicting requirements to obtain accurate redshift information and large statistics. We present cosmological results from the full 5000 sq. deg. of the Dark Energy Survey first three years of observations (Y3) combining those two-point functions, using for the first time a magnitude-limited lens sample (MagLim) of 11 million galaxies especially selected to optimize such combination, and 100 million background shapes. We consider two cosmological models, flat $\Lambda$CDM and $w$CDM. In $\Lambda$CDM we obtain for the matter density $\Omega_m = 0.320^{+0.041}_{-0.034}$ and for the clustering amplitude $S_8 = 0.778^{+0.037}_{-0.031}$, at 68\% C.L. The latter is only 1$\sigma$ smaller than the prediction in this model informed by measurements of the cosmic microwave background by the Planck satellite. In $w$CDM we find $\Omega_m = 0.32^{+0.044}_{-0.046}$, $S_8=0.777^{+0.049}_{-0.051}$, and dark energy equation of state $w=-1.031^{+0.218}_{-0.379}$. We find that including smaller scales while marginalizing over non-linear galaxy bias improves the constraining power in the $\Omega_m-S_8$ plane by $31\%$ and in the $\Omega_m-w$ plane by $41\%$ while yielding consistent cosmological parameters from those in the linear bias case. These results are combined with those from cosmic shear in a companion paper to present full DES-Y3 constraints from the three two-point functions (3x2pt).

We present a validation of the Dark Energy Survey Year 3 (DES Y3) $3\times2$-point analysis choices by testing them on Buzzard v2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak lensing analyses. We show that the Buzzard v2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak lensing statistics. We then show that our model for the $3\times2$-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the non-linear matter power spectrum, non-linear galaxy bias and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy-galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

This paper details the modeling pipeline and validates the baseline analysis choices of the DES Year 3 joint analysis of galaxy clustering and weak lensing (a so-called "3$\times$2pt" analysis). These analysis choices include the specific combination of cosmological probes, priors on cosmological and systematics parameters, model parameterizations for systematic effects and related approximations, and angular scales where the model assumptions are validated. We run a large number of simulated likelihood analyses using synthetic data vectors to test the robustness of our baseline analysis. We demonstrate that the DES Year 3 modeling pipeline, including the calibrated scale cuts, is sufficiently accurate relative to the constraining power of the DES Year 3 analyses. Our systematics mitigation strategy accounts for astrophysical systematics, such as galaxy bias, intrinsic alignments, source and lens magnification, baryonic effects, and source clustering, as well as for uncertainties in modeling the matter power spectrum, reduced shear, and estimator effects. We further demonstrate excellent agreement between two independently-developed modeling pipelines, and thus rule out any residual uncertainties due to the numerical implementation.

We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg$^2$. We perform an analysis combining three two-point correlation functions (3$\times$2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions. The analysis was designed to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results. We model the data within the flat $\Lambda$CDM and $w$CDM cosmological models. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude $S_8=0.776^{+0.017}_{-0.017}$ and matter density $\Omega_{\mathrm{m}} = 0.339^{+0.032}_{-0.031}$ in $\Lambda$CDM, mean with 68% confidence limits; $S_8=0.775^{+0.026}_{-0.024}$, $\Omega_{\mathrm{m}} = 0.352^{+0.035}_{-0.041}$, and dark energy equation-of-state parameter $w=-0.98^{+0.32}_{-0.20}$ in $w$CDM. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed $p=0.13$ to $0.48$. When combining DES 3$\times$2pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find $p=0.34$. Combining all of these data sets with Planck CMB lensing yields joint parameter constraints of $S_8 = 0.812^{+0.008}_{-0.008}$, $\Omega_{\mathrm{m}} = 0.306^{+0.004}_{-0.005}$, $h=0.680^{+0.004}_{-0.003}$, and $\sum m_{\nu}<0.13 \;\mathrm{eV\; (95\% \;CL)}$ in $\Lambda$CDM; $S_8 = 0.812^{+0.008}_{-0.008}$, $\Omega_{\mathrm{m}} = 0.302^{+0.006}_{-0.006}$, $h=0.687^{+0.006}_{-0.007}$, and $w=-1.031^{+0.030}_{-0.027}$ in $w$CDM. (abridged)

The debate surrounding fast magnetic energy dissipation by magnetic reconnection has remained a fundamental topic in the plasma universe, not only in the Earth's magnetosphere but in astrophysical objects such as pulsar magnetospheres and magnetars, for more than half a century. Recently, nonthermal particle acceleration and plasma heating during reconnection have been extensively studied, and it has been argued that rapid energy dissipation can occur for a collisionless "thin" current sheet, the thickness of which is of the order of the particle gyro-radius. However, it is an intriguing enigma as to how the fast energy dissipation can occur for a "thick" current sheet with thickness larger than the particle gyro-radius. Here we demonstrate, using a high-resolution particle-in-cell simulation for a pair plasma, that an explosive reconnection can emerge with the enhancement of the inertia resistivity due to the magnetization of the meandering particles by the reconnecting magnetic field and the shrinkage of the current sheet. In addition, regardless of the initial thickness of the current sheet, the time scale of the nonlinear explosive reconnection is tens of the Alfv\'{e}n transit time.

We study the $\delta$ Scuti -- red giant binary KIC9773821, the first double-pulsator binary of its kind. It was observed by \textit{Kepler} during its four-year mission. Our aims are to ascertain whether the system is bound, rather than a chance alignment, and to identify the evolutionary state of the red giant via asteroseismology. An extension of these aims is to determine a dynamical mass and an age prior for a $\delta$ Sct star, which may permit mode identification via further asteroseismic modelling. We determine spectroscopic parameters and radial velocities (RVs) for the red giant component using HERMES@Mercator spectroscopy. Light arrival-time delays from the $\delta$ Sct pulsations are used with the red-giant RVs to determine that the system is bound and to infer its orbital parameters, including the binary mass ratio. We use asteroseismology to model the individual frequencies of the red giant to give a mass of $2.10^{+0.20}_{-0.10}$ M$_{\odot}$ and an age of $1.08^{+0.06}_{-0.24}$ Gyr. We find that it is a helium-burning secondary clump star, confirm that it follows the standard $\nu_{\rm max}$ scaling relation, and confirm its observed period spacings match their theoretical counterparts in the modelling code MESA. Our results also constrain the mass and age of the $\delta$ Sct star. We leverage these constraints to construct $\delta$ Sct models in a reduced parameter space and identify four of its five pulsation modes.

Stress accumulation-relaxation meta-models of pulsar glitches make precise, microphysics-agnostic predictions of long-term glitch statistics, which can be falsified by existing and future timing data. Previous meta-models assume that glitches are triggered by an avalanche process, e.g. involving superfluid vortices, and that the probability density function (PDF) of the avalanche sizes is history-independent and specified exogenously. Here a recipe is proposed to generate the avalanche sizes endogenously in a history-dependent manner, by tracking the thresholds of occupied vortex pinning sites as a function of time. Vortices unpin spasmodically from sites with thresholds below a global, time-dependent stress and repin at sites with thresholds above the global stress, imbuing the system with long-term memory. The meta-model predicts PDFs, auto- and cross-correlations for glitch sizes and waiting times, which are provisionally inconsistent with current observations, unlike some previous meta-models (e.g. state dependent Poisson process), whose predictions are consistent. The theoretical implications are intriguing, albeit uncertain, because history-dependent avalanches embody faithfully the popular, idealized understanding in the literature of how vortex unpinning operates as a driven, stochastic process. The meta-model predicts aftershocks, which occur with larger than average sizes and longer than average waiting times after the largest, system-resetting glitches. This prediction will be tested, once more data are generated by the next generation of pulsar timing campaigns.

The corrugated horn is a high performance feed often used in radio telescopes. There has been a growing demand for wideband optics and corrugated horns in millimeter and submillimeter-wave receivers. It improves the observation efficiency and allows us to observe important emission lines such as CO in multiple excited states simultaneously. However, in the millimeter/submillimeter band, it has been challenging to create a conical corrugated horn with a fractional bandwidth of ~60% because the wavelength is very short, making it difficult to make narrow corrugations. In this study, we designed a conical corrugated horn with good return loss, low cross-polarization, and symmetric beam pattern in the 210-375GHz band (56% fractional bandwidth) by optimizing the dimensions of the corrugations. The corrugated horn was installed on the Osaka 1.85-m mm-submm telescope with the matched frequency-independent optics, and simultaneous observations of 12CO, 13CO, and C18O (J = 2-1, 3-2) were successfully made. In this paper, we describe the new design of the corrugated horn and report the performance evaluation results including the optics.

Context. Stellar spectral synthesis is essential for various applications, ranging from determining stellar parameters to comprehensive stellar variability calculations. New observational resources as well as advanced stellar atmosphere modelling, taking three dimensional (3D) effects from radiative magnetohydrodynamics calculations into account, require a more efficient radiative transfer. Aims. For accurate, fast and flexible calculations of opacity distribution functions (ODFs), stellar atmospheres and stellar spectra we developed an efficient code building on the well-established ATLAS9 code. The new code also paves the way for an easy and fast access to different elemental compositions in stellar calculations. Methods. For the generation of ODF tables we further developed the well-established DFSYNTHE code by implementing additional functionality, and a speed-up by employing a parallel computation scheme. In addition, the line lists used can be changed from Kurucz's recent lists. In particular, we implemented the VALD3 line list. Results. A new code, the Merged Parallelised Simplified ATLAS is presented. It combines the efficient generation of ODF, atmosphere modelling and spectral synthesis in local thermodynamic equilibrium, therefore being an all-in-one code. This all-in-one code provides more numerical functionality and is substantially faster compared to other available codes. The fully portable MPS-ATLAS code is validated against previous ATLAS9 calculations, the PHOENIX code calculations, and high quality observations.

The variable source OGLE-UCXB-01 is the first OGLE-discovered ultracompact X-ray binary (UCXB). The 12-year long-term OGLE optical photometry of this source shows a period of P= 12.8 min and a fast period decreasing rate Pdot= -9.2E-11 s s^-1. At a luminosity of L_X ~ 4E33 erg s^-1, its X-ray emission is also variable and correlated with the optical variability. To determine the nature of this variable source, specifically the masses and types of its binary components, we consider first an attractive possibility that the optical variation is due to the secondary's ellipsoidal variation and a strong gravitational wave emission drives the orbital decay. However, we can not find an allowable solution to the secondary that satisfies simultaneously the three constraints: an ultra-tight orbit, the bright absolute magnitude, and the large amplitude of the brightness variation. Moreover, the inferred mass transfer rate is too high. This scenario is therefore ruled out. We then find the system is fully consistent with an "intermediate polar" model, in which the optical and X-ray emission comes from a magnetized white dwarf (WD) accreting from a low-mass (<~ 0.7 M_sun) main-sequence secondary. The observed period decay is the accretion-driven spin-up of the WD. The WD spin period is 12.8 min and the orbital period is shorter than 10 hr. The method presented here can be applied to other UCXB candidates or impostors with time-domain data available only.

Structured jets are recently invoked to explain the complex emission of gamma ray bursts, such as GW 170817. Based on the accretion simulations, the jets are expected to have a structure that is more complex than a simple top-hat. Also, the structure of launching regions of blazar jets should influence their large scale evolution. This is recently revealed by the interactions of jet components in TXS 0506+056, where the jet is observed at a viewing angle close to zero. Observational studies have also shown an anti-correlation between the jet variability, measured e.g. by its minimum variability time scale, and the Lorentz factor, that spans several orders of magnitude and covers both blazars and GRBs samples. Motivated by those observational properties of black hole sources, we investigate the accretion inflow and outflow properties, by means of numerical GR MHD simulations. We perform axisymmetric calculations of the structure and evolution of central engine, composed of magnetized torus around Kerr black hole that is launching a non-uniform jet. We probe the jet energetics at different points along the line of sight, and we measure the jet time variability as localized in these specific regions. We quantify our results by computing the minimum variability timescales and power density spectra. We reproduce the MTS-$\Gamma$ correlation and we attribute it to the black hole spin as the main driving parameter of the engine. We also find that the PDS slope is not strongly affected by the black hole spin, while it differs for various viewing angles.

Content. Electron-cyclotron maser emission (ECME) is the favored mechanism for solar radio spikes and has been investigated extensively since the 1980s. Most studies relevant to solar spikes employ a loss-cone-type distribution of energetic electrons, generating waves mainly in the fundamental X/O mode (X1/O1), with a ratio of plasma oscillation frequency to electron gyrofrequency (${\omega}_ {pe}/{\Omega}_{ce}$) lower than 1. Despite the great progress made in this theory, one major problem is how the fundamental emissions pass through the second-harmonic absorption layer in the corona and escape. This is generally known as the escaping difficulty of the theory. Aims. We study the harmonic emissions generated by ECME driven by energetic electrons with the horseshoe distribution to solve the escaping difficulty of ECME for solar spikes. Methods. We performed a fully kinetic electromagnetic PIC simulation with ${\omega}_ {pe}/{\Omega}_{ce}$ = 0.1, corresponding to the strongly magnetized plasma conditions in the flare region, with energetic electrons characterized by the horseshoe distribution. We also varied the density ratio of energetic electrons to total electrons ($n_e/n_0$) in the simulation. Results. We obtain efficient amplification of waves in Z and X2 modes, with a relatively weak growth of O1 and X3. With a higher-density ratio, the X2 emission becomes more intense, and the rate of energy conversion from energetic electrons into X2 modes can reach $\sim$0.06% and 0.17%, with $n_e/n_0$= 5% and 10%, respectively. Conclusions. We find that the horseshoe-driven ECME can lead to an efficient excitation of X2 and X3 with a low value of ${\omega}_ {pe}/{\Omega}_{ce}$, providing novel means for resolving the escaping difficulty of ECME when applied to solar radio spikes. The simultaneous growth of X2 and X3 can be used to explain some harmonic structures observed in solar spikes.

The formation of Uranus' regular moons has been suggested to be linked to the origin of its enormous spin axial tilt (~98o). A giant impact between proto-Uranus and a 2-3 M_Earth impactor could lead to a large tilt and to the formation of a debris disc, where prograde and circular satellites are accreted. The most intriguing features of the current regular Uranian satellite system is that it possesses a positive trend in the mass-distance distribution and likely also in the bulk density, implying that viscous spreading of the debris disc after the giant impact plays a crucial role in shaping the architecture of the final system. In this paper, we investigate the formation of Uranus' satellites by combining results of SPH simulations for the giant impact, a 1D semi-analytic disc model for viscous spreading of the post-impact debris disc, and N-body simulations for the assembly of satellites from a disc of moonlets. Assuming the condensed rock (i.e., silicate) remains small and available to stick onto the relatively rapid growing condensed water-ice, we find that the best case in reproducing the observed mass and bulk composition of Uranus' satellite system is a pure-rocky impactor with 3 M_Earth colliding with the young Uranus with an impact parameter b = 0.75. Such an oblique collision could also naturally explain Uranus' large tilt and possibly, its low internal heat flux. The giant impact scenario can naturally explain the key features of Uranus and its regular moons. We therefore suggest that the Uranian satellite system formed as a result of an impact rather than from a circumplanetary disc. Objects beyond the water snow-line could be dominated by rocky objects similar to Pluto and Triton. Future missions to Uranus and its satellite system would further constrain the properties of Uranus and its moons and provide further insight on their formation processes.

We present new spectrophotometric and spectropolarimetric observations of Mrk 1239, one of the 8 prototypes that defines type-1 narrow-line Seyfert galaxies (NLS1s). Unlike the other typical NLS1s though, a high degree of polarization ($P\sim$5.6%) and red optical-IR ($g-W_4$ = 12.35) colors suggest that Mrk 1239 is more similar to type-2 active galactic nuclei like NGC 1068. Detailed analysis of spectral energy distribution in the UV-optical-IR yields two components from the nucleus: a direct and transmitted component that is heavily obscured ($E_{B-V} \approx 1.6$), and another indirect and scattered one with mild extinction ($E_{B-V} \sim$ 0.5). Such a two-light-paths scenario is also found in previous reports based on the X-ray data. Comparison of emission lines and the detection of He\,{\footnotesize I}*$\lambda$10830 BAL at [-3000,-1000] km s$^{-1}$ indicates that the obscuring clouds are at physical scale between the sublimation radius and that of the narrow emission line regions. The potential existence of powerful outflows is found as both the obscurer and scatterer are outflowing. Similar to many other type-2s, jet-like structure in the radio band is found in Mrk 1239, perpendicular to the polarization angle, suggesting polar scattering. We argue that Mrk 1239 is very probably a type-2 counterpart of NLS1s. The identification of 1 out of 8 prototype NLS1s as a type-2 counterpart implies that there can be a substantial amount of analogs of Mrk 1239 misidentified as type-1s in the optical band. Properties of these misidentified objects are going to be explored in our future works.

We present an updated version catalogue of gamma-ray source candidates, 12Y-MST, selected using the minimum spanning tree (MST) algorithm on the 12-year Fermi-LAT sky (Pass 8) at energies higher than 10 GeV. The high-energy sky at absolute Galactic latitudes above 20 degrees has been investigated using rather restrictive selection criteria, resulting in a total sample of 1664 photon clusters, or candidate sources. Of these, 230 are new detections, that is, candidate sources without any association in other gamma-ray catalogues. A large fraction of them have interesting counterparts, most likely blazars. We describe the main results on the catalogue selection and search of counterparts. We also present an additional sample of 224 candidate sources (12Y-MSTw), which are clusters that we extracted by applying weaker selection criteria: about 57% of them have not been reported in other catalogues.

We present a observational study of the dark matter fraction in 225 rotation supported star-forming galaxies at $z\approx 0.9$ having stellar mass range: $ 9.0 \leq log(M_* \ \mathrm{M_\odot}) \leq 11.0$ and star formation rate: $0.49 \leq log \left(SFR \ \mathrm{[M_{\odot}\ yr^{-1}]} \right) \leq 1.77$. This is a sub sample of KMOS redshift one spectroscopic survey (KROSS) previously studied by \citet{GS20}. The stellar masses ($M_*$) of these objects were previously estimated using mass-to-light ratios derived from fitting the spectral energy distribution of the galaxies. Star formation rates were derived from the H$_\alpha$ luminosities. The total gas masses ($M_{gas}$) are determined by scaling relations of molecular and atomic gas \citep[][respectively] {Tacconi2018, Lagos2011}. The dynamical masses ($M_{dyn}$) are directly derived from the rotation curves (RCs) at different scale lengths (effective radius: $R_e$, $\sim 2 \ R_e$ and $\sim 3 \ R_e$) and then the dark matter fractions ($f_{ DM }=1-M_{bar}/M_{dyn}$) at these radii are calculated. We report that at $z\sim 1$ only a small fraction ($\sim 5\%$) of our sample has a low ($< 20\%$) DM fraction within $\sim$ 2-3 $R_e$. The majority ($> 72\%$) of SFGs in our sample have dark matter dominated outer disks ($\sim 5-10$ kpc) in agreement with local SFGs. Moreover, we find a large scatter in the fraction of dark matter at a given stellar mass (or circular velocity) with respect to local SFGs, suggesting that galaxies at $z \sim 1$, a) span a wide range of stages in the formation of stellar disks, b) have diverse DM halo properties coupled with baryons.

Rapidly propagating fast magnetoacoustic wave trains guided by field-aligned plasma non-uniformities are confidently observed in the Sun's corona. Observations at large heights suggest that fast wave trains can travel long distances from the excitation locations. We study characteristic time signatures of fully developed, dispersive fast magnetoacoustic wave trains in field-aligned zero-$\beta$ plasma slabs in the linear regime. Fast wave trains are excited by a spatially localised impulsive driver and propagate along the waveguide as prescribed by the waveguide-caused dispersion. In slabs with steeper transverse density profiles, developed wave trains are shown to consist of three distinct phases: a long-period quasi-periodic phase with the oscillation period shortening with time, a multi-periodic (peloton) phase in which distinctly different periods co-exist, and a short-lived periodic Airy phase. The appearance of these phases is attributed to a non-monotonic dependence of the fast wave group speed on the parallel wavenumber due to the waveguide dispersion, and is shown to be different for axisymmetric (sausage) and non-axisymmetric (kink) modes. In wavelet analysis, this corresponds to the transition from the previously known tadpole shape to a new boomerang shape of the wave train spectrum, with two well-pronounced arms at shorter and longer periods. We describe a specific previously published radio observation of a coronal fast wave train, highly suggestive of a change of the wavelet spectrum from a tadpole to a boomerang, broadly consistent with our modelling. The applicability of these boomerang-shaped fast wave trains for probing the transverse structuring of the waveguiding coronal plasma is discussed.

We report on a search for pulsars at the positions of eight low-mass white dwarfs and one higher-mass white dwarf with the 100-m Effelsberg Radio Telescope. These systems have orbital parameters suggesting that their unseen companions are either massive white dwarfs or neutron stars. Our observations were performed at 1.36 GHz, reaching sensitivities of 0.1-0.2 mJy. We searched our data accounting for the possible acceleration and jerk of the pulsar signals due to orbital motion, but found no significant pulsar signals. Considering our result jointly with 20 non-detections of similar systems with the Greenbank Radio Telescope, we infer $f_{\rm NS}\leq 0.10$, for the fraction of NSs orbiting these white dwarfs. We discuss the sensitivity of this result to the underlying assumptions and conclude with a brief discussion on the prospects of targeted surveys for discovering millisecond pulsars.

We report the discovery of two planetary systems, namely G 264-012, an M4.0 dwarf with two terrestrial planets ($M_{\rm b}\sin{i} = 2.50^{+0.29}_{-0.30}$ M$_{\oplus}$ and $M_{\rm c}\sin{i} = 3.75^{+0.48}_{-0.47}$ M$_{\oplus}$), and Gl 393, a bright M2.0 dwarf with one terrestrial planet ($M_{\rm b}\sin{i} = 1.71 \pm 0.24$ M$_{\oplus}$). Although both stars were proposed to belong to young stellar kinematic groups, we estimate their ages to be older than about 700 Ma. The two planets around G 264-012 were discovered using only radial-velocity (RV) data from the CARMENES exoplanet survey, with estimated orbital periods of $2.30$ d and $8.05$ d, respectively. Photometric monitoring and analysis of activity indicators reveal a third signal present in the RV measurements, at about 100 d, caused by stellar rotation. The planet Gl 393 b was discovered in the RV data from the HARPS, CARMENES, and HIRES instruments. Its identification was only possible after modelling, with a Gaussian process (GP), the variability produced by the magnetic activity of the star. For the earliest observations, this variability produced a forest of peaks in the periodogram of the RVs at around the 34 d rotation period determined from {\em Kepler} data, which disappeared in the latest epochs. After correcting for them with this GP model, a significant signal showed at a period of $7.03$ d. No significant signals in any of our spectral activity indicators or contemporaneous photometry were found at any of the planetary periods. Given the orbital and stellar properties, the equilibrium temperatures of the three planets are all higher than that for Earth. Current planet formation theories suggest that these two systems represent a common type of architecture. This is consistent with formation following the core accretion paradigm.

Neutron stars in low mass X-ray binaries are hypothesised to emit continuous gravitational waves that may be detectable by ground-based observatories. The torque balance model predicts that a higher accretion rate produces larger-amplitude gravitational waves, hence low mass X-ray binaries with high X-ray flux are promising targets for gravitational wave searches. The detection of X-ray pulsations would identify the spin frequency of these neutron stars, and thereby improve the sensitivity of continuous gravitational-wave searches by reducing the volume of the search parameter space. We perform a semi-coherent search for pulsations in the two low mass X-ray binaries Scorpius X-1 and Cygnus X-2 using X-ray data from the \textit{ Rossi X-ray Timing Explorer} Proportional Counter Array. We find no clear evidence for pulsations, and obtain upper limits (at $90\%$ confidence) on the fractional pulse amplitude, with the most stringent being $0.034\%$ for Scorpius X-1 and $0.23\%$ for Cygnus X-2. These upper limits improve upon those of Vaughan et al. (1994) by factors of $\sim 8.2$ and $\sim 1.6$ respectively.

Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk is optically thick enough for stellar light, i.e., small-grain-to-gas mass ratio of $\gtrsim0.0001$. The mid-plane temperature varies as the waves propagate and hence gap and ring structures can be seen in both millimeter and infrared emission. The millimeter substructures can be observed even if the disk is fully optically thick since it is induced by the temperature variation, while density-induced substructures would disappear in the optically thick regime. The fractional separation between TWI-induced ring and gap is $\Delta r/r \sim$ 0.2-0.4 at $\sim$ 10-50 au, which is comparable to those found by ALMA. Due to the temperature variation, snow lines of volatile species move radially and multiple snow lines are observed even for a single species. The wave propagation velocity is as fast as $\sim$ 0.6 ${\rm au~yr^{-1}}$, which can be potentially detected with a multi-epoch observation with a time separation of a few years.

There is a longstanding discrepancy between the observed Galactic classical nova rate of $\sim 10$ yr$^{-1}$ and the predicted rate from Galactic models of $\sim 30$--50 yr$^{-1}$. One explanation for this discrepancy is that many novae are hidden by interstellar extinction, but the degree to which dust can obscure novae is poorly constrained. We use newly available all-sky three-dimensional dust maps to compare the brightness and spatial distribution of known novae to that predicted from relatively simple models in which novae trace Galactic stellar mass. We find that only half ($\sim 48$\%) of novae are expected to be easily detectable ($g \lesssim 15$) with current all-sky optical surveys such as the All-Sky Automated Survey for Supernovae (ASAS-SN). This fraction is much lower than previously estimated, showing that dust does substantially affect nova detection in the optical. By comparing complementary survey results from ASAS-SN, OGLE-IV, and the Palomar Gattini IR-survey in the context of our modeling, we find a tentative Galactic nova rate of $\sim 40$ yr$^{-1}$, though this could decrease to as low as $\sim 30$ yr$^{-1}$ depending on the assumed distribution of novae within the Galaxy. These preliminary estimates will be improved in future work through more sophisticated modeling of nova detection in ASAS-SN and other surveys.

Since neutrinos can escape from dense regions without being deflected, they are promising candidates to study the new physics at the sources that produce them. With the increasing development of more sensitive detectors in the coming years, we will infer several intrinsic properties from incident neutrinos. In particular, we centralise our study in those produced by thermal processes in short gamma-ray bursts (sGRBs) and their interactions within the central engine's anisotropic medium. On the one hand, we consider baryonic winds produced with a strong magnetic contribution, and on the other hand, we treat only neutrino-driven winds. First, we obtain the effective neutrino potential considering both baryonic density profiles around the central engine. Then, we get the three-flavour oscillation probabilities in this medium to finally calculate the expected neutrino ratios. We find a stronger angular dependence on the expected neutrino ratios, which, incidentally, contrast from the expected theoretical ratios without considering the winds' additional contribution. The joint analysis of this observable, together with the sGRB ejected jet angle, might lead to an effective mechanism to discriminate between the involved merger progenitors (black hole-neutron star or neutron star-neutron star), acting as an additional detection channel to gravitational waves.

Strange quark matter, which is composed of u, d, and s quarks, could be the true ground of matter. According to this hypothesis, compact stars may actually be strange quark stars, and there may even be stable strange quark dwarfs and strange quark planets. The detection of the binary neutron star merger event GW170817 provides us new clues on the equation of state of compact stars. In this study, the tidal deformability of strange quark planets and strange quark dwarfs are calculated. It is found that the tidal deformability of strange quark objects is smaller than that of normal matter counterparts. For a typical 0.6 M$_\odot$ compact star, the tidal deformability of a strange dwarf is about 1.4 times less than that of a normal white dwarf. The difference is even more significant between strange quark planets and normal matter planets. Additionally, if the strange quark planet is a bare one (i.e., not covered by a normal matter curst), the tidal deformability will be extremely small, which means bare strange quark planets will hardly be distorted by tidal forces. Our study clearly proves the effectiveness of identifying strange quark objects via searching for strange quark planets through gravitational-wave observations.

A cosmological zoom-in simulation which leads to a Milky Way like halo is started at redshift 7. The initial dark matter distribution is augmented with globular star clusters, composed of 4 M_sun star particles, placed in the sub-halos on three different sets of orbits: in the range of 10\% and 1\% of the sub-halo virial radii, and, at the exact center of the sub-halos. The three sets of star clusters end with median orbital radii of 30, 5 and 16 kpc, losing about 2, 50 and 15\% of their mass, respectively and have very substantially different stellar kinematics beyond about 3 half mass radii. The clusters started at larger sub-halo radii have a velocity dispersion that declines to 15-25\% of the central value in the outer 5-20 half mass radius region. The clusters started at sub-halo centers have a rise in velocity dispersion beyond 3-5 half mass radii, some reaching values of twice the central velocity dispersion at 10 half mass radii. The sub-halo centered clusters tend to have isotropic velocity distributions at large radii, whereas the isolated ones have a region of radial oriented velocity ellipsoid.

We report on the search for very-high-energy gamma-ray emission from the regions around three nearby supersonic pulsars (PSR B0355+54, PSR J0357+3205 and PSR J1740+1000) that exhibit long X-ray tails. To date there is no clear detection of TeV emission from any pulsar tail that is prominent in X-ray or radio. We provide upper limits on the TeV flux, and luminosity, and also compare these limits with other pulsar wind nebulae detected in X-rays and the tail emission model predictions. We find that at least one of the three tails is likely to be detected in observations that are a factor of 2-3 more sensitive. The analysis presented here also has implications for deriving the properties of pulsar tails, for those pulsars whose tails could be detected in TeV.

We revise the cosmological phenomenology of Macroscopic Dark Matter (MDM) candidates, also commonly dubbed as Macros. A possible signature of MDM is the capture of baryons from the cosmological plasma in the pre-recombination epoch, with the consequent injection of high-energy photons in the baryon-photon plasma. By keeping a phenomenological approach, we consider two broad classes of MDM in which Macros are composed either of ordinary matter or antimatter. In both scenarios, we also analyze the impact of a non-vanishing electric charge carried by Macros. We derive constraints on the Macro parameter space from three cosmological processes: the change in the baryon density between the end of the Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB) decoupling, the production of spectral distortions in the CMB and the kinetic coupling between charged MDM and baryons at the time of recombination. In the case of neutral Macros we find that the tightest constraints are set by the baryon density condition in most of the parameter space. For Macros composed of ordinary matter and with binding energy $I$, this leads to the following bound on the reduced cross-section: $\sigma_X/M_X \lesssim 6.8 \cdot 10^{-7} \left(I/\mathrm{MeV}\right)^{-1.56} \, \text{cm}^2 \, \text{g}^{-1}$. Charged Macros with surface potential $V_X$, instead, are mainly constrained by the tight coupling with baryons, resulting in $\sigma_X/M_X \lesssim 2 \cdot 10^{-11} \left(|V_X|/\mathrm{MeV}\right)^{-2} \text{cm}^2 \, \text{g}^{-1}$. Finally, we show that future CMB spectral distortions experiments, like PIXIE and SuperPIXIE, would have the sensitivity to probe larger regions of the parameter space: this would allow either for a possible evidence or for an improvement of the current bounds on Macros as dark matter candidates.

We present the results of a comprehensive, near-UV-to-near-IR Hubble Space Telescope WFC3 imaging study of the young planetary nebula (PN) NGC 6302, the archetype of the class of extreme bi-lobed, pinched-waist PNe that are rich in dust and molecular gas. The new WFC3 emission-line image suite clearly defines the dusty toroidal equatorial structure that bisects NGC 6302's polar lobes, and the fine structures (clumps, knots, and filaments) within the lobes. The most striking aspect of the new WFC3 image suite is the bright, S-shaped 1.64 micron [Fe II] emission that traces the southern interior of the east lobe rim and the northern interior of the west lobe rim, in point-symmetric fashion. We interpret this [Fe II] emitting region as a zone of shocks caused by ongoing, fast (~100 km/s), collimated, off-axis winds from NGC 6302's central star(s). The [Fe II] emission and a zone of dusty, N- and S-rich clumps near the nebular symmetry axis form wedge-shaped structures on opposite sides of the core, with boundaries marked by sharp azimuthal ionization gradients. Comparison of our new images with earlier HST/WFC3 imaging reveals that the object previously identified as NGC 6302's central star is a foreground field star. Shell-like inner lobe features may instead pinpoint the obscured central star's actual position within the nebula's dusty central torus. The juxtaposition of structures revealed in this HST/WFC3 imaging study of NGC 6302 presents a daunting challenge for models of the origin and evolution of bipolar PNe.

Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface by mass accretion from a companion star. Observation of X-ray bursts provides valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics of X-ray bursters. However, they were mainly focused on the burning layers above the NS surface. The effects of the EOS and the heating and cooling processes inside the NS are still unknown. In this study, we calculated a series of X-ray bursts using a general relativistic stellar-evolution code with several NS EOSs. We compared the X-ray burst models with the burst parameters of a clocked burster associated with GS 1826-24. We found a monotonic correlation between the NS radius and the light-curve profile. A larger radius shows a higher recurrence time and a large peak luminosity. In contrast, the dependence of light curves on the NS mass becomes more complicated, where the neutrino cooling suppress the efficiency of nuclear ignition. We also constrained the EOS and mass of GS~1826-24, i.e., stiffer EOSs, corresponding to larger NS radii, are unpreffered due to a too high peak luminosity. The EOS and the cooling and heating of NSs are important to discuss the theoretical and observational properties of X-ray bursts.

Presolar silicon carbide (SiC) grains in meteoritic samples can help constrain circumstellar condensation processes and conditions in C-rich stars and core-collapse supernovae. This study presents our findings on eight presolar SiC grains from AGB stars (four mainstream and one Y grain) and core-collapse supernovae (three X grains), chosen on the basis of {\mu}-Raman spectral features that were indicative of their having unusual non-3C polytypes and/or high degrees of crystal disorder. Analytical transmission electron microscopy (TEM), which provides elemental compositional and structural information, shows evidence for complex histories for the grains. Our TEM results confirm the presence of non-3C,2H crystal domains. Minor element heterogeneities and/or subgrains were observed in all grains analyzed for their compositions. The C/O ratios inferred for the parent stars varied from 0.98 to greater than or equal to 1.03. Our data show that SiC condensation can occur under a wide range of conditions, in which environmental factors other than temperature (e.g., pressure, gas composition, heterogeneous nucleation on pre-condensed phases) play a significant role. Based on previous {\mu}-Raman studies, about 10% of SiC grains may have infrared (IR) spectral features that are influenced by crystal defects, porosity, and/or subgrains. Future sub-diffraction limited IR measurements of complex SiC grains might shed further light on the relative contributions of each of these features to the shape and position of the characteristic IR 11-{\mu}m SiC feature and thus improve the interpretation of IR spectra of AGB stars like those that produced the presolar SiC grains.

Previous studies of the galaxy and galaxy cluster distribution in the local Universe found indications for a large extension of the Local Supercluster up to a radius of 190 h_70^-1 Mpc. We are using our large and highly complete CLASSIX survey of X-ray luminous galaxy clusters detected in the ROSAT All Sky Survey to trace the matter distribution in the local Universe and to explore the size of the flattened local density structure associated with the Local Supercluster. The Local Supercluster is oriented almost perpendicular to the Galactic plane. Since Galactic extinction increases towards the Galactic plane, objects are on average more easily visible perpendicular to the plane than close to it, also producing an apparent concentration of objects along the Local Supercluster. We can correct for this bias by a careful treatment of the survey selection function. We find a significant overdensity of clusters in a flattened structure along the Supergalactic plane with a thickness of about 50 Mpc and an extent of about 100 Mpc radius. Structures at a distance larger than 100 Mpc are not correlated to the Local Supercluster any more. The matter density contrast of the local superstructure to the surroundings is about a factor of 1.3 - 2.3. Within the Supergalactic plane the matter is concentrated mostly in two superclusters, the Perseus-Pisces Chain and Hydra-Centaurus supercluster. We have shown in our earlier work that the local Universe in a region with a radius of 100 - 170 Mpc has a lower density than the cosmic mean. For this reason, the Local Supercluster is not overdense with respect to the cosmic mean density. Therefore this local superstructure will not collapse as a whole in the future, but rather fragment.

The galaxy cluster system RXCJ0225.9-4154 with the two sub-clusters A3016 and A3017 is embedded in a large-scale structure filament with signatures of filamentary accretion. A Chandra observation of this system shows the two clusters as well as a filament of X-ray emission connecting them and a galaxy group therein. The main cluster, A3017, shows indications of shocks most probably from a recent interaction with cluster components along the filament axis as well as a cold front. The filament between the two clusters is likely to be heated by the accretion shocks of the clusters and we discuss two scenarios for the origin of this X-ray filament. The paper describes many properties of the components of the cluster merger system that are used to assist the interpretation of the observed configuration.

Tumbling asteroids belong to a small group of objects, whose angular velocity vector is unaligned with any of its principal axes of inertia. This leads to challenging efforts to model the trajectory of any spacecraft designed to orbit these bodies. In this work, we deepen a previous study on this topic, concerning the asteroid (99942) Apophis during its close encounter with the Earth in 2029. We analyze the orbital behaviour of a spacecraft orbiting the asteroid during this event, by including the effects of the changes of orientation of the spin axis of the asteroid, depending on two sets of initial conditions. The global dynamics of the spacecraft around the target are analyzed using three approaches, MEGNO, PMap, and Time-Series prediction. We confirm that no spacecraft with natural orbits could survive the high perturbations caused by the close encounter with our planet.

The Starobinsky inflation model is one of the simplest inflation models that is consistent with the cosmic microwave background observations. In order to explain dark matter of the universe, we consider a minimal extension of the Starobinsky inflation model with introducing the dark sector which communicates with the visible sector only via the gravitational interaction. In Starobinsky inflation model, a sizable amount of dark-sector particle may be produced by the inflaton decay. Thus, a scalar, a fermion or a vector boson in the dark sector may become dark matter. We pay particular attention to the case with dark non-Abelian gauge interaction to make a dark glueball a dark matter candidate. In the minimal setup, we show that it is difficult to explain the observed dark matter abundance without conflicting observational constraints on the coldness and the self-interaction of dark matter. We propose scenarios in which the dark glueball, as well as other dark-sector particles, from the inflaton decay become viable dark matter candidates. We also discuss possibilities to test such scenarios.

The evaluation of the electron-ion bremsstrahlung cross section exact to all orders in the Coulomb potential is computationally expensive due to the appearance of hypergeometric functions. Therefore tabulations are widely used. Here, we provide an approximate formula for the non-relativistic dipole process valid for all applicable relative velocities and photon energies. Its validity spans from the Born- to the classical regime and from soft-photon emission to the kinematic endpoint. The error remains below 3% (and widely below 1%) except at an isolated region of hard-photon emission at the quantum-to-classical crossover. We use the formula to obtain the thermally averaged emission spectrum and cooling function in a Maxwellian plasma and demonstrate that they are accurate to better than 2%.

We study the convective and absolute forms of azimuthal magnetorotational instability (AMRI) in a Taylor-Couette (TC) flow with an imposed azimuthal magnetic field. We show that the domain of the convective AMRI is wider than that of the absolute AMRI. Actually, it is the absolute instability which is the most relevant and important for magnetic TC flow experiments. The absolute AMRI, unlike the convective one, stays in the device, displaying a sustained growth that can be experimentally detected. We also study the global AMRI in a TC flow of finite height using DNS and find that its emerging butterfly-type structure -- a spatio-temporal variation in the form of upward and downward traveling waves -- is in a very good agreement with the linear stability analysis, which indicates the presence of two dominant absolute AMRI modes in the flow giving rise to this global butterfly pattern.

The DARWIN observatory is a proposed next-generation experiment to search for particle dark matter and other rare interactions. It will operate a 50 t liquid xenon detector, with 40 t in the time projection chamber (TPC). To inform the final detector design and technical choices, a series of technological questions must first be addressed. Here we describe a full-scale demonstrator in the vertical dimension, Xenoscope, with the main goal of achieving electron drift over a 2.6 m distance, which is the scale of the DARWIN TPC. We have designed and constructed the facility infrastructure, including the cryostat, cryogenic and purification systems, the xenon storage and recuperation system, as well as the slow control system. We have also designed a xenon purity monitor and the TPC, with the fabrication of the former nearly complete. In a first commissioning run of the facility without an inner detector, we demonstrated the nominal operational reach of Xenoscope and benchmarked the components of the cryogenic and slow control systems, demonstrating reliable and continuous operation of all subsystems over 40 days. The infrastructure is thus ready for the integration of the purity monitor, followed by the TPC. Further applications of the facility include R&D on the high voltage feedthrough for DARWIN, measurements of electron cloud diffusion, as well as measurements of optical properties of liquid xenon. In the future, Xenoscope will be available as a test platform for the DARWIN collaboration to characterise new detector technologies.

We present the first application of a hierarchical Markov Chain Monte Carlo (MCMC) follow-up on continuous gravitational-wave candidates from real-data searches. The follow-up uses an MCMC sampler to draw parameter-space points following the F-statistic. As outliers are narrowed down, coherence time increases, imposing more restrictive phase-evolution templates. We introduce a novel Bayes factor to compare results from different stages: The signal hypothesis is derived from first principles, while the noise hypothesis uses extreme value theory to derive a background model. The effectiveness of our proposal is evaluated on fake Gaussian data and applied to a set of 30 outliers produced by different continuous wave searches on O2 Advanced LIGO data. The results of our analysis suggest all but three outliers are inconsistent with an astrophysical origin under the standard continuous wave signal model. We successfully ascribe two of the surviving outliers to an instrumental artifact and a strong hardware injection present in the data. The behavior of the third outlier suggests an instrumental origin as well, but we could not relate it to any known instrumental cause.

Accurate line lists are important for the description of the spectroscopic nature of small molecules. While a line list for CN (an important molecule for chemistry and astrophysics) exists, no underlying energy spectroscopic model has been published, which is required to consider the sensitivity of transitions to a variation of the proton-to-electron mass ratio. Here we have developed a Duo energy spectroscopic model as well as a novel hybrid style line list for CN and its isotopologues, combining energy levels that are derived experimentally (Marvel), using the traditional/perturbative approach (Mollist), and the variational approach (from a Duo spectroscopic model using standard ExoMol methodology). The final Trihybrid ExoMol-style line list for 12C14N consists of 28,004 energy levels (6,864 experimental, 1,574 perturbative, the rest variational) and 2,285,103 transitions up to 60,000 cm-1 between the three lowest electronic states (X2Sigma+, A2Pi and B2Sigma+). The spectroscopic model created is used to evaluate CN as a molecular probe to constrain the variation of the proton-to-electron mass ratio; no overly promising sensitive transitions for extragalactic study were identified.

We propose a novel method for computing $p$-values based on nested sampling (NS) applied to the sampling space rather than the parameter space of the problem, in contrast to its usage in Bayesian computation. The computational cost of NS scales as $\log^2{1/p}$, which compares favorably to the $1/p$ scaling for Monte Carlo (MC) simulations. For significances greater than about $4\sigma$ in both a toy problem and a simplified resonance search, we show that NS requires orders of magnitude fewer simulations than ordinary MC estimates. This is particularly relevant for high-energy physics, which adopts a $5\sigma$ gold standard for discovery. We conclude with remarks on new connections between Bayesian and frequentist computation and possibilities for tuning NS implementations for still better performance in this setting.

Light axion fields, if they exist, can be sourced by neutron stars due to their coupling to nuclear matter, and play a role in binary neutron star mergers. We report on a search for such axions by analysing the gravitational waves from the binary neutron star inspiral GW170817. We find no evidence of axions in the sampled parameter space. The null result allows us to impose constraints on axions with masses below $10^{-11} {\rm eV}$ by excluding the ones with decay constants ranging from $1.6\times10^{16} {\rm GeV}$ to $10^{18} {\rm GeV}$ at $3\sigma$ confidence level. Our analysis provides the first constraints on axions from neutron star inspirals, and rules out a large region in parameter space that has not been probed by the existing experiments.

Gravitational waves emitted by black hole binary inspiral and mergers enable unprecedented strong-field tests of gravity, requiring accurate theoretical modelling of the expected signals in extensions of General Relativity.In this paper we model the gravitational wave emission of inspiralling binaries in scalar Gauss-Bonnet gravity theories. Going beyond the weak-coupling approximation, we derive the gravitational waveform to first post-Newtonian order beyond the quadrupole approximation and calculate new contributions from nonlinear curvature terms. We quantify the effect of these terms and provide ready-to-implement gravitational wave and scalar waveforms as well as the Fourier domain phase for quasi-circular binaries. We also perform a parameter space study, which indicates that the values of black hole scalar charges play a crucial role in the detectability of deviation from General Relativity. We also compare the scalar waveforms to numerical relativity simulations to assess the impact of the relativistic corrections to the scalar radiation. Our results provide important foundations for future precision tests of gravity.

The observed distributions of the source properties from gravitational-wave detections are biased due to the selection effects and detection criteria in the detections, analogous to the Malmquist bias. In this work, this observation bias is investigated through its fundamental statistical and physical origins. An efficient semi-analytical formulation for its estimation is derived which is as accurate as the standard method of numerical simulations, with only a millionth of the computational cost. Then, the estimated bias is used for model independent inferences on the binary black hole population. These inferences show additional structures, specifically two potential mass gaps in the joint mass distribution, which were not found via modelled inferences. Example ready-to-use scripts and some produced datasets for this method are shared in an online repository.

We consider k-essence, a scalar-tensor theory with first-order derivative self-interactions that can screen local scales from scalar fifth forces, while allowing for sizeable deviations from General Relativity on cosmological scales. We construct fully non-linear static stellar solutions that show the presence of this screening mechanism, and we use them as initial data for simulations of stellar oscillations and gravitational collapse in spherical symmetry. We find that for k-essence theories of relevance for cosmology, the screening mechanism works in the case of stellar oscillation, and suppresses the monopole scalar emission to undetectable levels. In collapsing stars, we find that the Cauchy problem, although locally well-posed, can lead to diverging characteristic speeds for the scalar field. By introducing a "fixing equation" in the spirit of J. Cayuso, N. Ortiz, and L. Lehner [Phys. Rev. D 96, 084043 (2017)], inspired in turn by dissipative relativistic hydrodynamics, we manage to evolve collapsing neutron stars past the divergence of the characteristic speeds. We show that in these systems, the screening mechanism is less efficient than for oscillating and static stars, because the collapsing star must shed away all of its scalar hair before forming a black hole. For k-essence theories of relevance for cosmology, the characteristic frequency of the resulting scalar monopole signal is too low for terrestrial detectors, but we conjecture that space-borne interferometers such as LISA might detect it if a supernova explodes in the Galaxy.

Our observed Universe has a very strong arrow of time rooted in its low entropy starting point. This low entropy start can be related to various "tuning puzzles" about the early state of the Universe. Here we explore the relationship between the arrow of time and the emergence of classical from quantum in the hopes of ultimately gaining insights into cosmological initial conditions. Our focus is on einselection, the process whereby interactions with an environment select preferred states for a quantum system. This process plays an essential role in the emergence of classical from quantum. Studies of einselection have so far been limited to cases that exhibit an arrow of time. Here we study the ability of equilibrium systems to exhibit einselection -- and investigate whether detailed balance prevents this -- motivated by the question of whether classicality requires an arrow of time. We present calculations in the adapted Caldeira-Leggett model which demonstrate that einselection can indeed take place in equilibrium systems, and show how this phenomenon is tied to histories which express an arrow of time, despite the global equilibrium. We discuss some interesting implications of our results for cosmology and cosmological initial conditions. We are intrigued and a bit surprised by the role the consistent histories formalism has ended up playing in our analysis.