Papers for Friday, Apr 30 2021

We characterize the extreme heartbeat star system MACHO 80.7443.1718 in the LMC using TESS photometry and spectroscopic observations from the Magellan Inamori Kyocera Echelle (MIKE) and SOAR Goodman spectographs. MACHO 80.7443.1718 was first identified as a heartbeat star system in the All-Sky Automated Survey for SuperNovae (ASAS-SN) with $P_{\rm orb}=32.836\pm0.008\,{\rm d}$. MACHO 80.7443.1718 is a young (${\sim}6$~Myr), massive binary, composed of a B0 Iae supergiant with $M_1 \simeq 35 M_\odot$ and an O9.5V secondary with $M_2 \simeq 16 M_\odot$ on an eccentric ($e=0.51\pm0.03$) orbit. In addition to having the largest variability amplitude amongst all known heartbeats stars, MACHO 80.7443.1718 is also one of the most massive heartbeat stars yet discovered. The B[e] supergiant has Balmer emission lines and permitted/forbidden metallic emission lines associated with a circumstellar disk. The disk rapidly dissipates at periastron which could indicate mass transfer to the secondary, but re-emerges immediately following periastron passage. MACHO 80.7443.1718 also shows tidally excited oscillations at the $N=25$ and $N=41$ orbital harmonics and has a rotational period of 4.4 d.

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be explained in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. In our search, we find that the signal from a first order phase transition is degenerate with that generated by Supermassive Black Hole Binary mergers. An interesting open question is how well gravitational wave observatories could separate such signals.

Several large stellar spectroscopic surveys are producing overwhelming amounts of data that can be used for determining stellar atmospheric parameters and chemical abundances. Nonetheless, the accuracy achieved in the derived astrophysical parameters is still insufficient, mainly because of the paucity of adequate calibrators, particularly in the metal-poor regime ([Fe/H] $\leq -$1.0). Here, we introduce the Titans metal-poor reference stars: a sample of 41 dwarf and subgiant stars with accurate parameters. Effective temperatures (Teff) were derived by fitting observed H$\alpha$ profiles with synthetic lines computed using 3D hydrodynamic NLTE models. Surface gravities (logg) were computed using evolutionary tracks and parallaxes from Gaia EDR3. The same methods recover the Teff values of the Gaia benchmark stars, which are mostly based on interferometric measurements, with a 1$\sigma$ dispersion of $\pm 50$ K. We assume this to be the accuracy of the H$\alpha$ profiles computed from 3D non-LTE models for metal-poor dwarfs and subgiants. We achieved an internal precision typically between 30-40 K, these errors dominated by instrumental effects. The final total uncertainty for the Teff values of the Titans are thus estimated to be of the order of $1\%$. The typical error for logg is $\leq$ 0.04 dex. In addition, we identified a few members of Gaia-Enceladus, of Sequoia, and of the Helmi stream in our sample. These stars can pave the way for the accurate chemical characterization of these Galactic substructures. Using the Titans as reference, large stellar surveys will be able to improve the internal calibration of their astrophysical parameters. Ultimately, this sample will help users of data from Gaia and large surveys in reaching their goal of redefining our understanding of stars, stellar systems, and the Milky Way.

We present an investigation of the interdependence of the optical-to-X-ray spectral slope ($\alpha_{\rm ox}$), the HeII equivalent-width (EW), and the monochromatic luminosity at 2500 Angstroms ($L_{2500}$). The values of $\alpha_{\rm ox}$ and HeII EW are indicators of the strength/shape of the quasar ionizing continuum, from the ultraviolet (UV; 1500--2500 Angstroms), through the extreme ultraviolet (EUV; 300--50 Angstroms), to the X-ray (2 keV) regime. For this investigation, we measure the HeII EW of 206 radio-quiet quasars devoid of broad absorption lines that have high-quality spectral observations of the UV and 2 keV X-rays. The sample spans wide redshift ($\approx$ 0.13--3.5) and luminosity (log$(L_{2500}$) $\approx$ 29.2--32.5 erg s$^{-1}$ Hz$^{-1}$) ranges. We recover the well-known $\alpha_{\rm ox}$--$L_{2500}$ and HeII EW--$L_{2500}$ anti-correlations, and we find a similarly strong correlation between $\alpha_{\rm ox}$ and HeII EW, and thus the overall spectral shape from the UV, through the EUV, to the X-ray regime is largely set by luminosity. A significant $\alpha_{\rm ox}$--HeII EW correlation remains after removing the contribution of $L_{2500}$ from each quantity, and thus the emission in the EUV and the X-rays are also directly tied. This set of relations is surprising, since the UV, EUV, and X-ray emission are expected to be formed in three physically distinct regions. Our results indicate the presence of a redshift-independent physical mechanism that couples the continuum emission from these three different regions, and thus controls the overall continuum shape from the UV to the X-ray regime.

The extensive timespan of modern radial velocity surveys have made the discovery of long-period substellar companions more common in recent years, however measuring the true masses of these objects remains challenging. Astrometry from the Gaia mission is expected to provide mass measurements for many of these long-period companions, but this data is not yet available. However, combining proper motion data from Gaia DR2 and the earlier Hipparcos mission makes it possible to measure true masses of substellar companions in favourable cases. In this work, we combine radial velocities with Hipparcos-Gaia astrometry to measure the true masses of two recently discovered long-period substellar companion candidates, HD 92987 B and HD 221420 b. In both cases, we find that the true masses are significantly higher than implied by radial velocities alone. A $2087 \pm 19$ m s$^{-1}$ astrometric signal reveals that HD 92987 B is not close to its $17$ $M_J$ minimum mass but is instead a $0.2562 \pm 0.0045$ $M_\odot$ star viewed at a near-polar orbital inclination, whereas the $22.9 \pm 2.2$ $M_J$ HD 221420 b can be plausibly interpreted as a high-mass "super-planet" or a low-mass brown dwarf. With semi-major axes of $\sim$10 AU both companions are interesting targets for direct imaging, and HD 221420 b in particular would be a benchmark metal-rich substellar object if it proves possible to directly detect. Our results demonstrate the power of Hipparcos-Gaia astrometry for studying long-period planet and brown dwarf candidates discovered from radial velocity surveys.

Astronomers require efficient automated detection and classification pipelines when conducting large-scale surveys of the (optical) sky for variable and transient sources. Such pipelines are fundamentally important, as they permit rapid follow-up and analysis of those detections most likely to be of scientific value. We therefore present a deep learning pipeline based on the convolutional neural network architecture called $\texttt{MeerCRAB}$. It is designed to filter out the so called 'bogus' detections from true astrophysical sources in the transient detection pipeline of the MeerLICHT telescope. Optical candidates are described using a variety of 2D images and numerical features extracted from those images. The relationship between the input images and the target classes is unclear, since the ground truth is poorly defined and often the subject of debate. This makes it difficult to determine which source of information should be used to train a classification algorithm. We therefore used two methods for labelling our data (i) thresholding and (ii) latent class model approaches. We deployed variants of $\texttt{MeerCRAB}$ that employed different network architectures trained using different combinations of input images and training set choices, based on classification labels provided by volunteers. The deepest network worked best with an accuracy of 99.5$\%$ and Matthews correlation coefficient (MCC) value of 0.989. The best model was integrated to the MeerLICHT transient vetting pipeline, enabling the accurate and efficient classification of detected transients that allows researchers to select the most promising candidates for their research goals.

The Gaia M-dwarf gap is a recently discovered feature in the colour-magnitude diagram that shows a deficiency of low-mass and low-metallicity stars at the lower end of the main sequence. We aim at performing theoretical stellar modelling at low metallicities using a fine mass step and a fine time step, looking specifically for the transition of models from partially to fully convective, since the convective kissing instability that occurs at this transition is believed to be the cause of the gap. Stellar evolution models with metallicities of Z = 0.01, Z = 0.001 and Z = 0.0001 are performed using MESA, with a mass step of 0.00025 M$_{\odot}$ and a time step of 50,000 years. The small time step produced models that experience loops in their evolutionary tracks in the Hertzsprung-Russell (HR) diagram. The fluctuations in effective temperature and luminosity correspond to repeated events in which the bottom of the convective envelope merges with the top of the convective core, transporting $^3$He from the core to the surface. In addition to the episodes of switching from partially to fully convective, several near-merger events that produced low amplitude fluctuations were also found. Low-metallicity models undergo the convective kissing instability for longer portions of their lifetime and with higher fluctuation amplitudes than models with higher metallicities. The small mass step used in the models revealed a discontinuity in the luminosity-mass relation at all three metallicities.

Modified Newtonian Dynamics (MOND) at low acceleration has been astonishingly powerful at explaining the flat rotation curve of galaxies and the relation between the baryonic content of the galaxies and their observed circular velocity, known as the Baryonic Tully-Fisher Relationship (BTFR). It is known that MOND fails at explaining the observed velocity dispersion of the ultra-faint dwarf galaxies (UFDs) with the justification that UFDs are more prone to tidal disruption in MOND compared to cold dark matter model. We show that: (i) the ratio of tidal to internal acceleration in UFDs is extremely low, (ii) there is no correlation between the deviation of UFDs from MOND's prediction as a function of tidal susceptibility, and (iii) recent constraints from Gaia proper motion analysis on the orbital parameters of the UFDs exacerbates the challenge to MOND. In particular, Gaia data indicates that Ursa Major I is experiencing a recent infall into the Milky Way's halo, and its inconsistency with MOND at 7-$\sigma$ level can not be attributed to being an early infall satellite. Moreover, the new data from Gaia DR2 shows Willman I to have the least eccentric orbit of all UFDs, and its deviation from MOND at 4-$\sigma$ level can not be attributed to a highly eccentric orbit as previously suggested. Finally, given that Tucana III is the only UFD observed to show tidal features, Reticulum II and Segue I are two other UFDs that potentially challenge MOND as they have comparable galactocentric distances to Tucana III while showing no tidal features. Whether wide binaries have inflated the velocity dispersion of the UFDs remains an open question to be addressed with future multi-epoch observations.

Multiple populations in globular clusters are usually explained by the formation of stars out of material with a chemical composition that is polluted to different degrees by the ejecta of short-lived, massive stars of various type. Among other things, these polluters differ by the amount of helium they spread in the surrounding medium. In this study we investigate whether the present-day photometric method used to infer the helium content of multiple populations indeed gives the true value or underestimates it by missing very He-rich, but rare stars. We focus on the specific case of NGC6752. We compute atmosphere models and synthetic spectra along isochrones produced for this cluster for a very broad range of He abundances covering the predictions of different pollution scenarios, including the extreme case of the fast-rotating massive star (FRMS) scenario. We calculate synthetic photometry in HST filters best suited to study the helium content. We subsequently build synthetic clusters with various distributions of stars. We finally determine the maximum helium mass fraction of these synthetic clusters using a method similar to that applied to observational data. We build toy models of clusters with various distributions of multiple populations and ensure that we are able to recover the input maximum Y. We then build synthetic clusters with the populations predicted by the FRMS scenario and find that while we slightly underestimate the maximum Y value, we are still able to detect stars much more He-rich than the current observed maximum Y. It is easier to determine the maximum Y on main sequence stars than on red giant branch stars, but qualitatively the results are unaffected by the sample choice. We show that in NGC6752 it is unlikely that stars more He-rich than the current observational limit of about 0.3 are present.

We report the detection for the first time in space of three new pure hydrocarbon cycles in TMC-1: c-C3HCCH (ethynyl cyclopropenylidene), c-C5H6 (cyclopentadiene) and c-C9H8 (indene). We derive a column density of 3.1e11 cm-2 for the former cycle and similar values, in the range (1-2)e13 cm-2, for the two latter molecules. This means that cyclopentadiene and indene, in spite of their large size, are exceptionally abundant, only a factor of five less abundant than the ubiquitous cyclic hydrocarbon c-C3H2. The high abundance found for these two hydrocarbon cycles, together with the high abundance previously found for the propargyl radical (CH2CCH) and other hydrocarbons like vinyl and allenyl acetylene (Agundez et al. 2021; Cernicharo et al. 2021a,b), start to allow us to quantify the abundant content of hydrocarbon rings in cold dark clouds and to identify the intermediate species that are probably behind the in situ bottom-up synthesis of aromatic cycles in these environments. While c-C3HCCH is most likely formed through the reaction between the radical CCH and c-C3H2, the high observed abundances of cyclopentadiene and indene are difficult to explain through currently proposed chemical mechanisms. Further studies are needed to identify how are five- and six-membered rings formed under the cold conditions of clouds like TMC-1.

We release RELIKE (Reionization Effective Likelihood), a fast and accurate effective likelihood code based on the latest Planck 2018 data that allows one constrain any model for reionization between $6 < z < 30$ using five constraints from the CMB reionization principal components (PC). We tested the code on two example models which showed excellent agreement with sampling the exact Planck likelihoods using either a simple Gaussian PC likelihood or its full kernel density estimate. This code enables a fast and consistent means for combining Planck constraints with other reionization data sets, such as kinetic Sunyaev-Zeldovich effects, line-intensity mapping, luminosity function, star formation history, quasar spectra, etc, where the redshift dependence of the ionization history is important. Since the PC technique tests any reionization history in the given range, we also derive model-independent constraints for the total Thomson optical depth $\tau_{\rm PC} = 0.0619^{+0.0056}_{-0.0068}$ and its $15\le z \le 30$ high redshift component $\tau_{\rm PC}(15, 30) < 0.020 $ (95\% C.L.). The upper limits on the high-redshift optical depth is a factor of $\sim3$ larger than those reported in the Planck 2018 cosmological parameter paper using the FlexKnot method and we validate our results with a direct analysis of a two-step model which permits this small high-$z$ component.

Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host more dynamically unstable planetary systems. Here we apply a sophisticated material sublimation model to white dwarf systems to show that these statistics are actually a natural result of the interplay of thermal and tidal forces, and show how they define the circumstellar regions where dusty debris disks can form. We demonstrate that these processes tend to prevent stability against both sublimative destruction and reaccretion into planetesimals for rocky materials until white dwarfs cool to below ~25,000-32,000 K, in agreement with the observed limit of ~27,000 K. For pure water ice, this critical temperature is less than 2,700 K (requiring a cooling age older the universe); this precludes pure water ice-rich debris disks forming through the accepted two-step mechanism. The critical temperature is size-dependent; more massive white dwarfs could potentially host dusty debris disks at warmer temperatures. Our model suggests that the location of the disks within the PG 0010+280, GD 56, GD 362, and PG 1541+651 systems are consistent with a forsterite-dominated olivine composition. We also find that very cool white dwarfs may simultaneously host multiple, independently formed dusty debris disks, consistent with observations of the LSPM J0207+3331 system.

Here, we present results on the intrinsic collision probabilities, $ P_I$, and range of collision speeds, $V_I$, as a function of the heliocentric distance, $r$, in the trans-Neptunian region. The collision speed is one of the parameters, that serves as a proxy to a collisional outcome e.g., complete disruption and scattering of fragments, or formation of crater, where both processes are directly related to the impact energy. We utilize an improved and de-biased model of the trans-Neptunian object (TNO) region from the "Outer Solar System Origins Survey" (OSSOS). It provides a well-defined orbital distribution model of TNOs, based on multiple opposition observations of more than 1000 bodies. In this work we compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using the \"{O}pik's approach, as revised and modified by Wetherill for non-circular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in ($r, V_I$) phase space. Although collisions are calculated to happen within $r\sim 20 - 200$~AU and $V_I \sim 0.1$~km/s to as high as $V_I\sim9$~km/s, most of the collisions are likely to happen at low relative velocities $V_I<1$~km/s and are dominated by the main classical belt.

Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section $\sigma/m$. We compare the shape of the main dark matter (DM) halo at redshift $z=0$ predicted by SIDM simulations (at $\sigma/m=0.1$, $1$, and $10$ cm$^2$ g$^{-1}$) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on $\sigma/m$. This transition could potentially be used to set limits on the SIDM cross-section in the MW.

A number of gamma-ray bursts (GRBs) exhibit the late simultaneous bumps in their optical and Xray afterglows around the jet break. Its origin is unclear. Based on the following two facts, we suggest that this feature may sound a transition of circum-burst environment from a free-wind medium to a homogeneous medium. (I) The late bump followed by a steep decay is strongly reminiscent of the afterglows of GRB 170817A, which is attributed to an off-axis observed external-forward shock (eFS) propagating in an interstellar medium. (II) Observations seem to feature a long shallow decay before the late optical bump, which is different from the afterglow of GRB 170817A. In this paper, we study the emission of an eFS propagating in a free-to-shocked wind for on/off-axis observers, where the mass density in the shocked-wind is almost constant. The late simultaneous bumps/plateaux in the optical and X-ray afterglows are really found around the jet break for high-viewing-angle observers. Moreover, there is a long plateau or shallow decay before the late bump in the theoretical light-curves, which is formed during the eFS propagating in the free-wind. For low-viewing-angle observers, the above bumps appear only in the situation that the structured jet has a low characteristic angle and the deceleration radius of the on-axis jet flow is at around or beyond the free-wind boundary. As examples, the X-ray and optical afterglows of GRBs 120326A, 120404A, and 100814A are fitted. We find that an off-axis observed eFS in a free-to-shocked wind can well explain the afterglows in these bursts.

One of the fundamental properties of an intermediate polar is the dynamical nature of the accretion flow as it encounters the white dwarf's magnetosphere. Many works have presumed a dichotomy between disk-fed accretion, in which the WD accretes from a Keplerian disk, and stream-fed accretion, in which the matter stream from the donor star directly impacts the WD's magnetosphere without forming a disk. However, there is also a third, poorly understood regime in which the accretion flow consists of a torus of diamagnetic blobs that encircles the WD. This mode of accretion is expected to exist at mass-transfer rates below those observed during disk-fed accretion, but above those observed during pure stream-fed accretion. We invoke the diamagnetic-blob regime to explain the exceptional TESS light curve of the intermediate polar TX Col, which transitioned into and out of states of enhanced accretion during Cycles 1 and 3. Power-spectral analysis reveals that the accretion was principally stream-fed. However, when the mass-transfer rate spiked, large-amplitude quasi-periodic oscillations (QPOs) abruptly appeared and dominated the light curve for weeks. The QPOs have two striking properties: they appear in a stream-fed geometry at elevated accretion rates, and they occur preferentially within a well-defined range of frequencies (~10-25 cycles per day). We propose that during episodes of enhanced accretion, a torus of diamagnetic blobs forms near the binary's circularization radius and that the QPOs are beats between the white dwarf's spin frequency and unstable blob orbits within the WD's magnetosphere. We discuss how such a torus could be a critical step in producing an accretion disk in a formerly diskless system.

A line-of-sight towards the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys towards the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hours of deep observations across 0.7--93 GHz. We report preliminary results from our survey for ETI beacons across 1--8 GHz with 7.0 and 11.2 hours of observations with Parkes and GBT, respectively. With our narrowband drifting signal search, we were able to place meaningful constraints on ETI transmitters across 1--4 GHz and 3.9--8 GHz with EIRP limits of $\geq$4$\times$10$^{18}$ W among 60 million stars and $\geq$5$\times$10$^{17}$ W among half a million stars, respectively. For the first time, we were able to constrain the existence of artificially dispersed transient signals across 3.9--8 GHz with EIRP $\geq$1$\times$10$^{14}$ W/Hz with a repetition period $\leq$4.3 hours. We also searched our 11.2 hours of deep observations of the GC and its surrounding region for Fast Radio Burst-like magnetars with the DM up to 5000 pc cm$^{-3}$ with maximum pulse widths up to 90 ms at 6 GHz. We detected several hundred transient bursts from SGR J1745$-$2900, but did not detect any new transient burst with the peak luminosity limit across our observed band of $\geq$10$^{31}$ erg s$^{-1}$ and burst-rate of $\geq$0.23 burst-hr$^{-1}$. These limits are comparable to bright transient emission seen from other Galactic radio-loud magnetars, constraining their presence at the GC.

We investigate the Hi envelope of the young, massive GMCs in the star-forming regions N48 and N49, which are located within the high column density Hi ridge between two kpc-scale supergiant shells, LMC 4 and LMC 5. New long-baseline Hi 21 cm line observations with the Australia Telescope Compact Array (ATCA) were combined with archival shorter baseline data and single dish data from the Parkes telescope, for a final synthesized beam size of 24.75" by 20.48", which corresponds to a spatial resolution of ~ 6 pc in the LMC. It is newly revealed that the Hi gas is highly filamentary, and that the molecular clumps are distributed along filamentary Hi features. In total 39 filamentary features are identified and their typical width is ~ 21 (8-49) [pc]. We propose a scenario in which the GMCs were formed via gravitational instabilities in atomic gas which was initially accumulated by the two shells and then further compressed by their collision. This suggests that GMC formation involves the filamentary nature of the atomic medium.

We seek to find the precursors of the Herbig Ae/Be stars in the solar vicinity within 500 pc from the Sun. We do this by creating an optically selected sample of intermediate mass T-Tauri stars (IMTT stars) here defined as stars of masses $1.5 M_{\odot}\leq M_* \leq 5 M_{\odot}$ and spectral type between F and K3, from literature. We use literature optical photometry (0.4-1.25$\mu$m) and distances determined from \textit{Gaia} DR2 parallax measurements together with Kurucz stellar model spectra to place the stars in a HR-diagram. With Siess evolutionary tracks we identify intermediate mass T-Tauri stars from literature and derive masses and ages. We use Spitzer spectra to classify the disks around the stars into Meeus Group I and Group II disks based on their [F$_{30}$/F$_{13.5}$] spectral index. We also examine the 10$\mu$m silicate dust grain emission and identify emission from Polycyclic Aromatic Hydrocarbons (PAH). From this we build a qualitative picture of the disks around the intermediate mass T-Tauri stars and compare this with available spatially resolved images at infrared and at sub-millimeter wavelengths to confirm our classification. We find 49 intermediate mass T-Tauri stars with infrared excess. The identified disks are similar to the older Herbig Ae/Be stars in disk geometries and silicate dust grain population. Spatially resolved images at infra-red and sub-mm wavelengths suggest gaps and spirals are also present around the younger precursors to the Herbig Ae/Be stars. Comparing the timescale of stellar evolution towards the main sequence and current models of protoplanetary disk evolution the similarity between Herbig Ae/Be stars and the intermediate mass T-Tauri stars points towards an evolution of Group I and Group II disks that are disconnected, and that they represent two different evolutionary paths.

(abridged) Context. The origin of hot exozodiacal dust and its connection with outer dust reservoirs remains unclear. Aims. We aim to explore the possible connection between hot exozodiacal dust and warm dust reservoirs (> 100 K) in asteroid belts. Methods. We use precision near-infrared interferometry with VLTI/PIONIER to search for resolved emission at H band around a selected sample of nearby stars. Results. Our observations reveal the presence of resolved near-infrared emission around 17 out of 52 stars, four of which are shown to be due to a previously unknown stellar companion. The 13 other H-band excesses are thought to originate from the thermal emission of hot dust grains. Taking into account earlier PIONIER observations, and after reevaluating the warm dust content of all our PIONIER targets through spectral energy distribution modeling, we find a detection rate of 17.1(+8.1)(-4.6)% for H-band excess around main sequence stars hosting warm dust belts, which is statistically compatible with the occurrence rate of 14.6(+4.3)(-2.8)% found around stars showing no signs of warm dust. After correcting for the sensitivity loss due to partly unresolved hot disks, under the assumption that they are arranged in a thin ring around their sublimation radius, we however find tentative evidence at the 3{\sigma} level that H-band excesses around stars with outer dust reservoirs (warm or cold) could be statistically larger than H-band excesses around stars with no detectable outer dust. Conclusions. Our observations do not suggest a direct connection between warm and hot dust populations, at the sensitivity level of the considered instruments, although they bring to light a possible correlation between the level of H-band excesses and the presence of outer dust reservoirs in general.

X-ray flux from the inner hot region around central compact object in a binary system illuminates the upper surface of an accretion disc and it behaves like a corona. This region can be photoionised by the illuminating radiation, thus can emit different emission lines. We study those line spectra in black hole X-ray binaries for different accretion flow parameters including its geometry. The varying range of model parameters captures maximum possible observational features. We also put light on the routinely observed Fe line emission properties based on different model parameters, ionization rate, and Fe abundances. We find that the Fe line equivalent width $W_{\rm E}$ decreases with increasing disc accretion rate and increases with the column density of the illuminated gas. Our estimated line properties are in agreement with observational signatures.

"Changing-look quasars" (CLQs) are active galactic nuclei (AGN) showing extreme variability that results in a transition from Type 1 to Type 2. The short timescales of these transitions present a challenge to the unified model of AGN and the physical processes causing these transitions remain poorly understood. CLQs also provide interesting samples for the study of AGN host galaxies since the central emission disappears almost entirely. Previous searches for CLQs have utilised photometric variability or SDSS classification changes to systematically identify CLQs, this approach may miss lower luminosity CLQs. In this paper, we aim to use spectroscopic data to asses if analysis difference spectra can be used to detect further changing look quasars missed by photometric searches. We search SDSS-II DR 7 repeat spectra for sources that exhibit either a disappearance or appearance of both broad line emission and accretion disk continuum emission by directly analysing the difference spectrum between two epochs of observation. From a sample of 24,782 objects with difference spectra, our search yielded six CLQs within the redshift range $0.1 \leq z \leq 0.3$, including four newly identified sources. Spectral analysis indicates that changes in accretion rate can explain the changing-look behaviour. While a change in dust extinction fits the changes in spectral shape, the time-scales of the changes observed are too short for obscuration from torus clouds. Using difference spectra was shown to be an effective and sensitive way to detect CLQs. We recover CLQs an order of magnitude lower in luminosities than those found by photometric searches and achieve higher completeness than spectroscopic searches relying on pipeline classification.

Results. We illustrate our profile-fitting technique and present the K\,{\sc i} velocity structure of the dense ISM along the paths to all targets. As a validation test of the dust map, we show comparisons between distances to several reconstructed clouds with recent distance assignments based on different techniques. Target star extinctions estimated by integration in the 3D map are compared with their K\,{\sc i} 7699 A absorptions and the degree of correlation is found comparable to the one between the same K\,{\sc i} line and the total hydrogen column for stars distributed over the sky that are part of a published high resolution survey. We show images of the updated dust distribution in a series of vertical planes in the Galactic longitude interval 150-182.5 deg and our estimated assignments of radial velocities to the opaque regions. Most clearly defined K\,{\sc i} absorptions may be assigned to a dense dust cloud between the Sun and the target star. It appeared relatively straightforward to find a velocity pattern consistent will all absorptions and ensuring coherence between adjacent lines of sight, at the exception of a few weak lines. We compare our results with recent determinations of velocities of several clouds and find good agreement. These results demonstrate that the extinction-K\,{\sc i} relationship is tight enough to allow linking the radial velocity of the K\,{\sc i} lines to the dust clouds seen in 3D, and that their combination may be a valuable tool in building a 3D kinetic structure of the dense ISM. We discuss limitations and perspectives for this technique.

Gamma-Ray Integrated Detectors (GRID) is a student project designed to use multiple gamma-ray detectors carried by nanosatellites (CubeSat), forming a full-time and all-sky gamma-ray detection network to monitor the transient gamma-ray sky in the multi-messenger astronomy era. A compact CubeSat gamma-ray detector has been designed and implemented for GRID, including its hardware and firmware. The detector employs four Gd2Al2Ga3O12 : Ce (GAGG:Ce) scintillators coupled with four silicon photomultiplier (SiPM) arrays to achieve a high detection efficiency of gamma rays between 10 keV and 2 MeV with low power and small dimensions. The first detector designed by the undergraduate student team onboard a commercial CubeSat was launched into a Sun-synchronous orbit on 29 October 2018. The detector has been in a normal observation state and accumulated data for approximately 1 month after on-orbit functional and performance tests in 2019.

Although the expansion of the Universe explicitly breaks the time-translation symmetry, cosmological predictions for the stochastic gravitational wave background (SGWB) are usually derived under the so-called stationary hypothesis. By dropping this assumption and keeping track of the time dependence of gravitational waves at all length scales, we derive the expected unequal-time (and equal-time) waveform of the SGWB generated by scaling sources, such as cosmic defects. For extinct and smooth enough sources, we show that all observable quantities are uniquely and analytically determined by the holomorphic Fourier transform of the anisotropic stress correlator. Both the strain power spectrum and the energy density parameter are shown to have an oscillatory fine structure, they significantly differ on large scales while running in phase opposition at large wavenumbers $k$. We then discuss scaling sources that are never extinct nor smooth and which generate a singular Fourier transform of the anisotropic stress correlator. For these, we find the appearance of interferences on top of the above-mentioned fine-structure as well as atypical behaviour at small scales. For instance, we expect the rescaled strain power spectrum $k^2 \mathcal{P}_h$ generated by long cosmic strings in the matter era to oscillate around a scale invariant plateau. These singular sources are also shown to produce orders of magnitude difference between the rescaled strain spectra and the energy density parameter suggesting that only the former should be used for making reliable observable predictions. Finally, we discuss how measuring such a fine structure in the SGWB could disambiguate the possible cosmological sources.

Transit photometry is perhaps the most successful method for detecting exoplanets to date. However, a substantial amount of signal processing is needed since the dip in the signal detected, an indication that there is a planet in transit, is minuscule compared to the overall background signal due mainly to its host star. In this paper, we put forth a doable and straightforward method to enhance the signal and reduce noise. We discuss how to achieve higher planetary signals by subtracting equal halves of the host star - a folded detection. This results in a light curve with a double peak-to-peak signal, 2R_p^2/R_s^2, compared to the usual transit. We derive an expression of the light curve and investigate the effect of two common noises: the white Gaussian background noise and the noise due to the occurrences of sunspots. We show that in both simulation and analytical expression, the folded transit reduces the effective noise by a factor of 1/sqrt(2). This reduction and the doubling of the signal enables: (1) less number of transit measurements to get a definitive transiting planet signal and (2) detection of smaller planetary radii with the usual transit with the same number of transit data. Furthermore, we show that in the presence of multiple sunspots, the estimation of planetary parameters is more accurate. While our calculations may be very simple, it covers the basic concept of planetary transits.

The Dual-Frequency synthetic aperture radar (DFSAR) system manifested on the Chandrayaan-2 spacecraft represents a significant step forward in radar exploration of solid solar system objects. It combines SAR at two wavelengths (L- and S-bands) and multiple resolutions with several polarimetric modes in one lightweight ($\sim$ 20 kg) package. The resulting data from DFSAR support calculation of the 2$\times$2 complex scattering matrix for each resolution cell, which enables lunar near surface characterization in terms of radar polarization properties at different wavelengths and incidence angles. In this paper, we report on the calibration and preliminary performance characterization of DFSAR data based on the analysis of a sample set of crater regions on the Moon. Our calibration analysis provided a means to compare on-orbit performance with pre-launch measurements and the results matched with the pre-launch expected values. Our initial results show that craters in both permanently shadowed regions (PSRs) and non-PSRs that are classified as Circular Polarization Ratio (CPR)-anomalous in previous S-band radar analyses appear anomalous at L-band also. We also observe that material evolution and physical properties at their interior and proximal ejecta are decoupled. For Byrgius C crater region, we compare our analysis of dual-frequency radar data with the predicted behaviours of theoretical scattering models. If crater age estimates are available, comparison of their radar polarization properties at multiple wavelengths similar to that of the three unnamed south polar crater regions shown in this study may provide new insights into how the rockiness of craters evolves with time.

We report the detection of a short-lived narrow quasi-periodic oscillation (QPO) at ~88 mHz in an Insight-HXMT observation during the soft state of the persistent black hole high mass X-ray binary Cygnus X-1. This QPO is significantly detected in all the three instruments of Insight-HXMT, so in the broad energy range 1-250 keV. The fractional RMS of the QPO does not show significant variations above 3 keV (~5%) while decreases at lower energy (~2%). We show that this QPO is different from the type-A, -B, and -C QPOs usually observed in black hole X-ray binaries. We compare QPOs at similar frequencies that have been previously detected in another persistent high mass X-ray binaries in the soft state, we speculate that such QPOs might relate to some local inhomogeneity rarely formed in the accretion flow of wind-fed accretion systems.

We revisit large spectroscopic data sets for field stars from the literature to derive the upper Li envelope in the high metallicity regime in our Galaxy. We take advantage of Gaia EDR3 data and state-of-the-art stellar models to precisely determine the position of the sample dwarf stars in the Hertzsprung-Russell diagram. The highest Li abundances are found in field metal-rich warm dwarfs from the GALAH survey, located on the hot side of the Li-dip. Their mean Li value agrees with what was recently derived for warm dwarfs in metal-rich clusters, pointing towards a continuous increase of Li up to super-solar metallicity. However, if only cool dwarfs are considered in GALAH, as done in the other literature surveys, it is found that the upper Li envelope decreases at super-solar metallicities, blurring the actual Li evolution picture. We confirm the suggestion that field and open cluster surveys that found opposite Li behaviour in the high metallicity regime do not sample the same types of stars: The first ones, with the exception of GALAH, miss warm dwarfs that can potentially preserve their original Li content. Although we can discard the bending of the Li upper envelope at high metallicity derived from the analysis of cool star samples, we still need to evaluate the effects of atomic diffusion on warm, metal-rich early-F and late-A type dwarfs before deriving the actual Li abundance at high metallicity.

We develop a reliable, fully automatic method for the detection of coronal holes, that provides consistent full-disk segmentation maps over the full solar cycle and can perform in real-time. We use a convolutional neural network to identify the boundaries of coronal holes from the seven EUV channels of the Atmospheric Imaging Assembly (AIA) as well as from line-of-sight magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). For our primary model (Coronal Hole RecOgnition Neural Network Over multi-Spectral-data; CHRONNOS) we use a progressively growing network approach that allows for efficient training, provides detailed segmentation maps and takes relations across the full solar-disk into account. We provide a thorough evaluation for performance, reliability and consistency by comparing the model results to an independent manually curated test set. Our model shows good agreement to the manual labels with an intersection-over-union (IoU) of 0.63. From the total of 261 coronal holes with an area $>1.5\cdot10^{10}$ km$^2$ identified during the time range 11/2010 - 12/2016, 98.1% were correctly detected by our model. The evaluation over almost the full solar cycle no. 24 shows that our model provides reliable coronal hole detections, independent of the level of solar activity. From the direct comparison over short time scales of days to weeks, we find that our model exceeds human performance in terms of consistency and reliability. In addition, we train our model to identify coronal holes from each channel separately and show that the neural network provides the best performance with the combined channel information, but that coronal hole segmentation maps can be also obtained solely from line-of-sight magnetograms.

We present the far ultraviolet (FUV) imaging of the nearest Jellyfish or Fireball galaxy IC3418/VCC 1217, in the Virgo cluster of galaxies, using Ultraviolet Imaging Telescope (UVIT) onboard the ASTROSAT satellite. The young star formation observed here in the 17 kpc long turbulent wake of IC3418, due to ram pressure stripping of cold gas surrounded by hot intra-cluster medium, is a unique laboratory that is unavailable in the Milkyway. We have tried to resolve star forming clumps, seen compact to GALEX UV images, using better resolution available with the UVIT and incorporated UV-optical images from Hubble Space Telescope archive. For the first time, we resolve the compact star forming clumps (fireballs) into sub-clumps and subsequently into a possibly dozen isolated stars. We speculate that many of them could be blue supergiant stars which are cousins of SDSS J122952.66+112227.8, the farthest star (~17 Mpc) we had found earlier surrounding one of these compact clumps. We found evidence of star formation rate (4 - 7.4 x 10^-4 M_sun per yr ) in these fireballs, estimated from UVIT flux densities, to be increasing with the distance from the parent galaxy. We propose a new dynamical model in which the stripped gas may be developing vortex street where the vortices grow to compact star forming clumps due to self-gravity. Gravity winning over turbulent force with time or length along the trail can explain the puzzling trend of higher star formation rate and bluer/younger stars observed in fireballs farther away from the parent galaxy.

Models of stellar structure and evolution can be constrained using accurate measurements of the parameters of eclipsing binary members of open clusters. Multiple binary stars provide the means to tighten the constraints and, in turn, to improve the precision and accuracy of the age estimate of the host cluster. In the previous two papers of this series, we have demonstrated the use of measurements of multiple eclipsing binaries in the old open cluster NGC6791 to set tighter constraints on the properties of stellar models than was previously possible, thereby improving both the accuracy and precision of the cluster age. We identify and measure the properties of a non-eclipsing cluster member, V56, in NGC\,6791 and demonstrate how this provides additional model constraints that support and strengthen our previous findings. We analyse multi-epoch spectra of V56 from FLAMES in conjunction with the existing photometry and measurements of eclipsing binaries in NGC6971. The parameters of the V56 components are found to be $M_{\rm p}=1.103\pm 0.008 M_{\odot}$ and $M_{\rm s}=0.974\pm 0.007 M_{\odot}$, $R_{\rm p}=1.764\pm0.099 R_{\odot}$ and $R_{\rm s}=1.045\pm0.057 R_{\odot}$, $T_{\rm eff,p}=5447\pm125$ K and $T_{\rm eff,s}=5552\pm125$ K, and surface [Fe/H]=$+0.29\pm0.06$ assuming that they have the same abundance. The derived properties strengthen our previous best estimate of the cluster age of $8.3\pm0.3$ Gyr and the mass of stars on the lower red giant branch (RGB), which is $M_{\rm RGB} = 1.15\pm0.02M_{\odot}$ for NGC6791. These numbers therefore continue to serve as verification points for other methods of age and mass measures, such as asteroseismology.

In a previous paper, we computed the energy density and the non-linear energy cascade rate for transverse kink waves using Elsasser variables. In this paper, we focus on the standing kink waves, which are impulsively excited in coronal loops by external perturbations. We present an analytical calculation to compute the damping time due to the non-linear development of the Kelvin-Helmholtz instability. The main result is that the damping time is inversely proportional to the oscillation amplitude. We compare the damping times from our formula with the results of numerical simulations and observations. In both cases we find a reasonably good match. The comparison with the simulations show that the non-linear damping dominates in the high amplitude regime, while the low amplitude regime shows damping by resonant absorption. In the comparison with the observations, we find a power law inversely proportional to the amplitude $\eta^{-1}$ as an outer envelope for our Monte Carlo data points.

The dramatic increase in sensitivity, spectral coverage and resolution of radio astronomical facilities in recent years has opened new possibilities for observation of chemical differentiation and isotopic fractionation in protostellar sources to shed light on their spatial and temporal evolution. In warm interstellar environments, methanol is an abundant species, hence spectral data for its isotopic forms are of special interest. In the present work, the millimeter-wave spectrum of the $^{13}$CH$_3$OD isotopologue has been investigated over the region from 150$-$510 GHz to provide a set of transition frequencies for potential astronomical application. The focus is on two types of prominent $^{13}$CH$_3$OD spectral groupings, namely the $a$-type $^qR$-branch multiplets and the $b$-type $Q$-branches. Line positions are reported for the $^qR(J)$ clusters for $J = 3$ to 10 for the $v_{\rm t} = 0$ and 1 torsional states, and for a number of $v_{\rm t} = 0$ and 1 $^rQ(J)$ or $^pQ(J)$ line series up to $J = 25$. The frequencies have been fitted to a multi-parameter torsion-rotation Hamiltonian, and upper level excitation energies have been calculated from the resulting molecular constants.

Thermal dissociation and recombination of molecular hydrogen, H_2, in the atmospheres of ultra-hot Jupiters (UHJs) has been shown to play an important role in global heat redistribution. This, in turn, significantly impacts their planetary emission, yet only limited investigations on the atmospheric effects have so far been conducted. Here we investigate the heat redistribution caused by this dissociation/recombination reaction, alongside feedback mechanisms between the atmospheric chemistry and radiative transfer, for a planetary and stellar configuration typical of UHJs. To do this, we have developed a time-dependent pseudo-2D model, including a treatment of time-independent equilibrium chemical effects. As a result of the reaction heat redistribution, we find temperature changes of up to $\sim$400 K in the atmosphere. When TiO and VO are additionally considered as opacity sources, these changes in temperature increase to over $\sim$800 K in some areas. This heat redistribution is found to significantly shift the region of peak atmospheric temperature, or hotspot, towards the evening terminator in both cases. The impact of varying the longitudinal wind speed on the reaction heat distribution is also investigated. When excluding TiO/VO, increased wind speeds are shown to increase the impact of the reaction heat redistribution up to a threshold wind speed. When including TiO/VO there is no apparent wind speed threshold, due to thermal stabilisation by these species. We also construct pseudo-2D phase curves from our model, and highlight both significant spectral flux damping and increased phase offset caused by the reaction heat redistribution.

Radio-loud quasars (RLQs) are known to produce excess X-ray emission, compared to radio-quiet quasars (RQQs) of the same luminosity, commonly attributed to jet-related emission. Recently, we found that the HeII EW and $\alpha_{\rm{ox}}$ in RQQs are strongly correlated, which suggests that their extreme-ultraviolet (EUV) and X-ray emission mechanisms are tightly related. Using 48 RLQs, we show that steep-spectrum radio quasars (SSRQs) and low radio-luminosity ($L_{\rm R}$) flat-spectrum radio quasars (FSRQs) follow the $\alpha_{\rm ox}$--HeII EW relation of RQQs. This suggests that the X-ray and EUV emission mechanisms in these types of RLQs is the same as in RQQs, and is not jet related. High-$L_{\rm R}$ FSRQs show excess X-ray emission given their HeII EW by a factor of $\approx$ 3.5, which suggests that only in this type of RLQ is the X-ray production likely jet related.

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

Inconsistencies regarding the nature of globular cluster multiple population radial distributions is a matter for concern given their role in testing or validating cluster dynamical evolution modeling. In this study, we present a re-analysis of eight globular cluster radial distributions using publicly available ground-based ugriz and UBVRI photometry; correcting for a systematic error identified in the literature. We detail the need for including and considering not only K-S probabilities but critical K-S statistic values as well when drawing conclusions from radial distributions, as well as the impact of sample incompleteness. Revised cumulative radial distributions are presented, and the literature of each cluster reviewed to provide a fuller picture of our results. We find that many multiple populations are not as segregated as once thought, and that there is a pressing need for better understanding of the spatial distributions of multiple populations in globular clusters.

There have recently been detections of radio emission from low-mass stars, some of which are indicative of star-planet interactions. Motivated by these exciting new results, here we present stellar wind models for the active planet-hosting M dwarf AU Mic. Our models incorporate the large-scale photospheric magnetic field map of the star, reconstructed using the Zeeman-Doppler Imaging method. We use our models to assess if planet-induced radio emission could be generated in the corona of AU Mic, through a mechanism analogous to the sub-Alfv\'enic Jupiter-Io interaction. In the case that AU Mic has a mass-loss rate of 27 times that of the Sun, we find that both planets b and c in the system can induce radio emission from 10 MHz to 3 GHz in the corona of the host star for the majority of their orbits, with peak flux densities of 10 mJy. Our predicted emission bears a striking similarity to that recently reported from GJ 1151 by Vedantham et al. (2020), which is indicative of being induced by a planet. Detection of such radio emission would allow us to place an upper limit on the mass-loss rate of the star.

Hybrid Morphology Radio Sources (HyMoRS) are a very rare and newly discovered subclass of radio galaxies that have mixed FR morphology i.e., these galaxies have FR-I structure on one side of the core and FR-II structure on the other side of the core. We systematically searched for HyMoRS using VLA Faint Images of the Radio Sky at Twenty-cm (FIRST) survey at 1400 MHz and identified forty-five confirmed HyMoRS and five candidates HyMoRS. Our finding significantly increased the known sample size of HyMoRS. HyMoRS may play an essential role in understanding the interaction of jets with the interstellar medium and a very debated topic of the FR dichotomy. We identified optical/IR counterparts for thirty-nine sources in our catalogue. In our sample of sources, five sources had Quasar-like behavior. We had estimated the spectral index and radio luminosity of HyMoR sources in our catalogue, when possible. We found that the source J1336+2329 ($\log L=26.93$ W Hz$^{-1}$sr$^{-1}$) was the most luminous and the source J1204+3801, a Quasar, was the farthest HyMoRS (with redshift $z$=1.28) in our sample. With the help of a large sample size of the newly discovered sources, various statistical properties were studied.

We use an up to date compilation of Tully-Fisher data to search for transitions in the evolution of the Tully-Fisher relation. Using a recently published data compilation, we find hints at $\approx 3\sigma$ level for a transition at a critical distance $D_c \simeq 17 Mpc$. The zero point (intercept) amplitude of the transition is $\Delta \log A_B \simeq 0.2 \pm 0.06$ while the slope remains practically unchanged. If the transition is interpreted as due to a gravitational strength transition, it would imply a shift of the effective gravitational constant to lower values for distances larger than $D_c\simeq 17 Mpc$ by $\frac{\Delta G}{G}=-0.1 \pm 0.03$. Such a shift is of the anticipated sign and magnitude but at somewhat lower distance (redshift) than the gravitational transition recently proposed to address the Hubble and growth tensions ($\frac{\Delta G}{G}\simeq -0.1$ at transition redshift $z_t\lesssim 0.01$ ($D_c\lesssim 40 Mpc$)).

In addition to its suite of narrow dense rings, Uranus is surrounded by an extremely complex system of dusty rings that were most clearly seen by the Voyager spacecraft after it flew past the planet. A new analysis of the highest resolution images of these dusty rings reveals that a number of them are less than 20 km wide. The extreme narrowness of these rings, along with the fact that most of them do not appear to fall close to known satellite resonances, should provide new insights into the forces responsible for sculpting the Uranian ring system.

The DES-CMASS sample (DMASS) is designed to optimally combine the weak lensing measurements from the Dark Energy Survey (DES) and redshift-space distortions (RSD) probed by the CMASS galaxy sample from the Baryonic Oscillation Spectroscopic Survey (BOSS). In this paper, we demonstrate the feasibility of adopting DMASS as the equivalent of BOSS CMASS for a joint analysis of DES and BOSS in the framework of modified gravity. We utilize the angular clustering of the DMASS galaxies, cosmic shear of the DES METACALIBRATION sources, and cross-correlation of the two as data vectors. By jointly fitting the combination of the data with the RSD measurements from the BOSS CMASS sample and Planck data, we obtain the constraints on modified gravity parameters $\mu_0 = -0.37^{+0.47}_{-0.45}$ and $\Sigma_0 = 0.078^{+0.078}_{-0.082}$. We do not detect any significant deviation from General Relativity. Our constraints of modified gravity measured with DMASS are tighter than those with the DES Year 1 redMaGiC galaxy sample with the same external data sets by $29\%$ for $\mu_0$ and $21\%$ for $\Sigma_0$, and comparable to the published results of the DES Year 1 modified gravity analysis despite this work using fewer external data sets. This improvement is mainly because the galaxy bias parameter is shared and more tightly constrained by both CMASS and DMASS, effectively breaking the degeneracy between the galaxy bias and other cosmological parameters. Such an approach to optimally combine photometric and spectroscopic surveys using a photometric sample equivalent to a spectroscopic sample can be applied to combining future surveys having a limited overlap such as DESI and LSST.

We present a mechanism for accelerated expansion of the universe in the generic case of negative-curvature compactifications of M-theory, with minimal ingredients. M-theory on a hyperbolic manifold with small closed geodesics supporting Casimir energy -- along with a single classical source (7-form flux) -- contains an immediate 3-term structure for volume stabilization at positive potential energy. Hyperbolic manifolds are well-studied mathematically, with an important rigidity property at fixed volume. They and their Dehn fillings to more general Einstein spaces exhibit explicit discrete parameters that yield small closed geodesics supporting Casimir energy. The off-shell effective potential derived by M. Douglas incorporates the warped product structure via the constraints of general relativity, screening negative energy. Analyzing the fields sourced by the localized Casimir energy and the available discrete choices of manifolds and fluxes, we find a regime where the net curvature, Casimir energy, and flux compete at large radius and stabilize the volume. Further metric and form field deformations are highly constrained by hyperbolic rigidity and warping effects, leading to calculations giving strong indications of a positive Hessian, and residual tadpoles are small. We test this via explicit back reacted solutions and perturbations in patches including the Dehn filling regions, initiate a neural network study of further aspects of the internal fields, and derive a Maldacena-Nunez style no-go theorem for Anti-de Sitter extrema. A simple generalization incorporating 4-form flux produces axion monodromy inflation. As a relatively simple de Sitter uplift of the large-N M2-brane theory, the construction applies to de Sitter holography as well as to cosmological modeling, and introduces new connections between mathematics and the physics of string/M theory compactifications.

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath, and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in 2D. We solve our island kinetic equation numerically for three different types of initial distribution: delta-distribution, Gaussian and power-law distribution. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde t^{-1}$, the number of islands decreases as $\tilde t^{-1}$, and the averaged area of islands grows as $\tilde t$, where $\tilde t$ is the time normalized to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied.

The Milky-Way magnetic field can trigger conversions between photons and axion-like particles (ALPs), leading to peculiar features on the observable photon spectra. Previous studies considered only the regular component of the magnetic field. However, observations consistently show the existence of an additional turbulent component, with a similar strength and correlated on a scale of a few 10$\,$pc. We investigate the impact of the turbulent magnetic field on the ALP-photon conversions, characterizing the effects numerically and analytically. We show that the turbulent magnetic field can change the conversion probability by up to a factor of two and may lead to observable irregularities in the observable photon spectra from different astrophysical sources.

The freeze-in production of Feebly Interacting Massive Particle (FIMP) dark matter in the early universe is an appealing alternative to the well-known - and constrained - Weakly Interacting Massive Particle (WIMP) paradigm. Although challenging, the phenomenology of FIMP dark matter has been receiving growing attention and is possible in a few scenarios. In this work, we contribute to this endeavor by considering a $Z^\prime$ portal to fermionic dark matter, with the $Z^\prime$ having both vector and axial couplings and a mass ranging from MeV up to PeV. We evaluate the bounds on both freeze-in and freeze-out from direct detection, atomic parity violation, leptonic anomalous magnetic moments, neutrino-electron scattering, collider, and beam dump experiments. We show that FIMPs can already be tested by most of these experiments in a complementary way, whereas WIMPs are especially viable in the $Z^\prime$ low mass regime, in addition to the $Z^\prime$ resonance region. We also discuss the role of the axial couplings of $Z^\prime$ in our results. We therefore hope to motivate specific realizations of this model in the context of FIMPs, as well as searches for these elusive dark matter candidates.

We discuss cosmological solutions for a diffeomorphism invariant gauge theory of the non-compact Lorentz group $SO(1,3)$. Besides the gauge bosons our model of pregeometry contains a vector field in the vector representation of $SO(1,3)$ and a scalar singlet. General relativity and variable gravity emerge as effective theories for large distances and times in Planck units. We propose an approximation to the effective action with up to two derivatives. For a suitable range of parameters the universe approaches for large times stable Minkowski space. For late cosmology the model predicts dynamical dark energy and provides for a candidate for dark matter. Early cosmology is characterized by an inflationary epoch. The beginning of the universe in the infinite past is great emptiness, corresponding to an ultraviolet fixed point with the associated quantum scale symmetry. The beginning universe is a vacuum state with vanishing expectation values and finite non-vanishing correlation functions for the fluctuations of all fields. There is no physical big bang singularity.

We tested four $f(Q,T)$ models in an extension of symmetric teleparallel gravity whose Friedmann equations reduce to $\Lambda CDM$ for certain parameters. Using low-redshift data we found that all our models were 2$\sigma$ consistent with $\Lambda CDM$ and the Hubble constant value were $\sigma$ consistent with the one of the SH0ES collaboration and at $12\sigma$ tension with the one of the Planck Collaboration. To see whether one of our models can challenge $\Lambda CDM$ at a background perspective, we computed the Bayesian evidence for our four models and $\Lambda CDM$. The concordance model was preferred over all our $f(Q,T)$ alternate models, showing a weak preference against models $f(Q,T) = -Q/G_N + bT$ and $-(Q+2\Lambda)/G_N + bT$ and a substantial preference against models $f(Q,T) = -(Q+2H_0c(Q/(6H_0^2))^{n+1})/G_N + bT$ and $f(Q,T) = -(Q+2H_0c(Q/(6H_0^2))^{n+1} + 2\Lambda)/G_N + bT$. Our models were successful to reproduce $\Lambda CDM$ acceleration from low redshifts.

Inspired by an interesting counterexample to the cosmic no-hair conjecture found in a supergravity-motivated model recently, we propose a multi-field extension, in which two scalar fields are allowed to non-minimally couple to two vector fields, respectively. This model is shown to admit an exact Bianchi type I power-law solution. Furthermore, stability analysis based on the dynamical system method is performed to show that this anisotropic solution is indeed stable and attractive if both scalar fields are canonical. Nevertheless, if one of the two scalar fields is phantom then the corresponding anisotropic power-law inflation turns unstable as expected.

The goal of this work is to find the simplest UV completion of Accidental Composite Dark Matter Models (ACDM) that can dynamically generate an asymmetry for the DM candidate, the lightest \textit{dark baryon} (DCb), and simultaneously annihilate the symmetric component. In this framework the DCb is a bound state of a confining $\text{SU}(N)_{\text{DC}}$ gauge group, and can interact weakly with the visible sector. The constituents of the DCb can possess non-trivial charges under the Standard Model gauge group. The generation of asymmetry for such candidate is a two-flavor variation of the \emph{out-of-equilibrium} decay of a heavy scalar, with mass $M_\phi\gtrsim 10^{15}$ GeV. Below the scale of the scalars, the models recover accidental stability, or long-livedness, of the DM candidate. The symmetric component is annihilated by residual confined interactions provided that the mass of the DCb $m_{\text{DCb}} \lesssim 75$ TeV. We implement the mechanism of asymmetry generation, or a variation of it, in all the original ACDM models, managing to generate the correct asymmetry for DCb of masses in this range. For some of the models found, the stability of the DM candidate is not spoiled even considering generic GUT completions or asymmetry generation mechanisms in the visible sector.

We construct the global phase portraits of inflationary dynamics in teleparallel gravity models with a scalar field nonminimally coupled to torsion scalar. The adopted set of variables can clearly distinguish between different asymptotic states as fixed points, including the kinetic and inflationary regimes. The key role in the description of inflation is played by the heteroclinic orbits which run from the asymptotic saddle points to the late time attractor point and are approximated by nonminimal slow roll conditions. To seek the asymptotic fixed points we outline a heuristic method in terms of the "effective potential" and "effective mass", which can be applied for any nonminimally coupled theories. As particular examples we study positive quadratic nonminimal couplings with quadratic and quartic potentials, and note how the portraits differ qualitatively from the known scalar-curvature counterparts. For quadratic models inflation can only occur at small nonminimal coupling to torsion, as for larger coupling the asymptotic de Sitter saddle point disappears from the physical phase space. Teleparallel models with quartic potentials are not viable for inflation at all, since for small nonminimal coupling the asymptotic saddle point exhibits weaker than exponential expansion, and for larger coupling disappears too.

White dwarfs, the most abundant stellar remnants, provide a promising means of probing dark matter interactions, complimentary to terrestrial searches. The scattering of dark matter from stellar constituents leads to gravitational capture, with important observational consequences. In particular, white dwarf heating occurs due to the energy transfer in the dark matter capture and thermalisation processes, and the subsequent annihilation of captured dark matter. We consider the capture of dark matter by scattering on either the ion or the degenerate electron component of white dwarfs. For ions, we account for the stellar structure, the star opacity, realistic nuclear form factors that go beyond the simple Helm approach, and finite temperature effects pertinent to sub-GeV dark matter. Electrons are treated as relativistic, degenerate targets, with Pauli blocking, finite temperature and multiple scattering effects all taken into account. We also estimate the dark matter evaporation rate. The dark matter-nucleon/electron scattering cross sections can be constrained by comparing the heating rate due to dark matter capture with observations of cold white dwarfs in dark matter-rich environments. We apply this technique to observations of old white dwarfs in the globular cluster Messier 4, which we assume to be located in a DM subhalo. For dark matter-nucleon scattering, we find that white dwarfs can probe the sub-GeV mass range inaccessible to direct detection searches, with the low mass reach limited only by evaporation, and can be competitive with direct detection in the $1-10^4$ GeV range. White dwarf limits on dark matter-electron scattering are found to outperform current electron recoil experiments over the full mass range considered, and extend well beyond the $\sim 10$ GeV mass regime where the sensitivity of electron recoil experiments is reduced.

This article discusses a dark energy cosmological model in the standard theory of gravity - general relativity with a broad scalar field as a source. Exact solutions of Einstein's field equations are derived by considering a particular form of deceleration parameter $q$, which shows a smooth transition from decelerated to accelerated phase in the evolution of the universe. The external datasets such as Hubble ($H(z)$) datasets, Supernovae (SN) datasets, and Baryonic Acoustic Oscillation (BAO) datasets are used for constraining the model par parameters appearing in the functional form of $q$. The transition redshift is obtained at $% z_{t}=0.67_{-0.36}^{+0.26}$ for the combined data set ($H(z)+SN+BAO$), where the model shows signature-flipping and is consistent with recent observations. Moreover, the present value of the deceleration parameter comes out to be $q_{0}=-0.50_{-0.11}^{+0.12}$ and the jerk parameter $% j_{0}=-0.98_{-0.02}^{+0.06}$ (close to 1) for the combined datasets, which is compatible as per Planck2018 results. The analysis also constrains the omega value i.e., $\Omega _{m_{0}}\leq 0.269$ for the smooth evolution of the scalar field EoS parameter. It is seen that energy density is higher for the effective energy density of the matter field than energy density in the presence of a scalar field. The evolution of the physical and geometrical parameters is discussed in some details with the model parameters' numerical constrained values. Moreover, we have performed the state-finder analysis to investigate the nature of dark energy.

A striking feature of our fundamentally indeterministic quantum universe is its quasiclassical realm -- the wide range of time place and scale in which the deterministic laws of classical physics hold. Our quasiclassical realmis an emergent feature of the fundamental theories of our universe's quantum state and dynamics. There are many types of quasiclassical realms our Universe could exhibit characterized by different variables, different levels of coarse-graining, different locations in spacetime, different classical physics, and different levels of classicality.We propose a measure of classicality for quasiclassical realms, We speculate on the observable consequences of different levels of classicality especially for information gathering and utilizing systems (IGUSes) such ourselves as observers of the Universe.

We study the possibility of low scale leptogenesis along with dark matter (DM) in the presence of primordial black holes (PBH). For a common setup to study both leptogenesis and DM we consider the minimal scotogenic model which also explains light neutrino mass at radiative level. While PBH in the mass range of $0.1-10^5$ g can, in principle, affect leptogenesis, the required initial PBH fraction usually leads to overproduction of DM whose thermal freeze-out occurs before PBH evaporation. PBH can lead to non-thermal source of leptogenesis as well as dilution of thermally generated lepton asymmetry via entropy injection, with the latter being dominant. The parameter space of scotogenic model which leads to overproduction of baryon or lepton asymmetry in standard cosmology can be made consistent in the presence of PBH with appropriate initial mass and energy fraction. On the other hand, for such PBH parameters, the DM is constrained to be in light mass regime where its freeze-out occurs after PBH evaporation.

We discuss a model where a mixed warm and hot keV neutrino dark matter rises naturally. We arrange active and sterile neutrinos in the same $SU(3)_L$ multiplet, with the lightest sterile neutrino being dark matter. The other two heavy sterile neutrinos, through their out-of-equilibrium decay, contribute both to the dilution of dark matter density and its population, after freeze-out. We show that this model features all ingredients to overcome the overproduction of keV neutrino dark matter, and explore the phenomenological implications for Big Bang Nucleosynthesis and the number of relativistic degrees of freedom.