Papers for Wednesday, Sep 15 2021

We examine the properties of spiral shocks from a steady, adiabatic, non-axisymmetric accretion disk around a compact star in binary. We first time incorporate all the possible influences from binary through adopting the Roche potential and Coriolis forces in the basic conservation equations. In this paper, we assume the spiral shocks to be point-wise self-similar, and the flow is in vertical hydrostatic equilibrium to simplify the study. We also investigate the mass outflow due to the shock compression and apply it to the accreting white dwarf in binary. We find that our model will be beneficial to overcome the ad hoc assumption of optically thick wind generally used in the studies of the progenitor of supernovae Ia.

The light curves of tidally-locked hot Jupiters transiting fast-rotating, early type stars are a rich source of information about both planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $\Psi$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $\Psi$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $\mu$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $\Psi=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we can separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $\Psi$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally-resolved JWST phase curve data.

If primordial black holes (PBHs) seeded the supermassive black holes (SMBHs) at the centers of high-redshift quasars, then the gas surrounding these black holes may reveal nucleosynthetic clues to their primordial origins. We present predictions of altered primordial abundances around PBHs massive enough to seed SMBHs at z~6-7.5. We find that if PBHs with initial masses of ~10^5 M$_{\odot}$ are responsible for such SMBHs, they may produce primordial Deuterium and Helium fractions enhanced by >~ 10%, and Lithium abundance depleted by >~10%, at distances of up to ~ a comoving kiloparsec away from the black hole after decoupling. We estimate that ~ 10^8 M$_{\odot}$ of gas is enhanced (or depleted) by at least one percent. Evidence of these modified primordial Deuterium, Helium, and Lithium abundances could still be present if this circum-PBH gas remains unaccreted by the SMBH and in or near the host galaxies of high-redshift quasars. Measuring the abundance anomalies will be challenging, but could offer a novel way to reveal the primordial origin of such SMBH seeds.

We present the accretion disk size estimates for a sample of 19 active galactic nuclei (AGN) using the optical $g$, $r$, and $i$ band light curves obtained from the Zwicky Transient Facility (ZTF) survey. All the AGN have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multi-band light curves are cross-correlated, and the reverberation lag is estimated using the Interpolated Cross-Correlation Function (ICCF) method and the Bayesian method using the {\sc javelin} code. As expected from the disk reprocessing arguments, the $g-r$ band lags are shorter than the $g-i$ band lags for this sample. The interband lags for all, but 5 sources, are larger than the sizes predicted from the standard Shakura Sunyaev (SS) analytical model. We fit the light curves directly using a thin disk model implemented through the {\sc javelin} code to get the accretion disk sizes. The disk sizes obtained using this model are on an average 3.9 times larger than the prediction based on the SS disk model. We find a weak correlation between the disk sizes and the known physical parameters, namely, the luminosity and the SMBH mass. In the near future, a large sample of AGN covering a range of luminosity and SMBH mass from large photometric surveys would be helpful in a better understanding of the structure and physics of the accretion disk.

We present approaches to quickly simulate systematics affecting CMB observations, including the effects of the scanning strategy. Using summary properties of the scan we capture features of full time ordered data (TOD) simulations, allowing maps and power spectra to be generated at much improved speed for a number of systematics - the cases we present experienced speed ups of 3-4 orders of magnitude when implementing the map-based approaches. We demonstrate the effectiveness of the approaches at capturing the salient features of the scan by directly comparing to full TOD simulations - seeing agreement at sub-percent levels of accuracy. We simulate the effects of differential gain, pointing, and ellipticity to show the effectiveness of the approaches, but note that one could extend these techniques to other systematics. We finally show how to apply these fast map-based simulations of systematic effects to a full focal plane showing their ability to incorporate thousands of detectors as seen in modern CMB experiments.

We develop and implement an inclination-dependent attenuation prescription for spectral energy distribution (SED) fitting and study its impact on derived star-formation histories. We apply our prescription within the SED fitting code Lightning to a clean sample of 82, z = 0.21-1.35 disk-dominated galaxies in the Great Observatories Origins Deep Survey (GOODS) North and South fields. To compare our inclination-dependent attenuation prescription with more traditional fitting prescriptions, we also fit the SEDs with the inclination-independent Calzetti et al. (2000) attenuation curve. From this comparison, we find that fits to a subset of 58, z < 0.7 galaxies in our sample, utilizing the Calzetti et al. (2000) prescription, recover similar trends with inclination as the inclination-dependent fits for the FUV-band attenuation and recent star-formation rates. However, we find a difference between prescriptions in the optical attenuation (AV) that is strongly correlated with inclination (p-value < 10^-10). For more face-on galaxies, with i < 50 deg, (edge-on, i = 90 deg), the average derived AV is 0.30 +/- 0.10 magnitudes lower (0.55 +/- 0.15 magnitudes higher) for the inclination-dependent model compared to traditional methods. Further, the ratio of stellar masses between prescriptions also has a significant (p-value < 10^-2) trend with inclination. For i = 0-65 deg, stellar masses are systematically consistent between fits, with log(Mstar_inc/Mstar_Calzetti) = -0.05 +/- 0.03 dex, and scatter of 0.11 dex. However, for i = 80-90 deg, derived stellar masses are lower for the Calzetti et al. (2000) fits by an average factor of 0.17 +/- 0.02 dex, and scatter of 0.13 dex. Therefore, these results suggest that SED fitting assuming the Calzetti et al. (2000) attenuation law potentially underestimates stellar masses in highly inclined disk-dominated galaxies.

As the Canadian Hydrogen Intensity Mapping Experiment (CHIME) has become the leading instrument for detecting Fast Radio Bursts (FRBs), CHIME/FRB Outriggers will use very-long-baseline interferometry (VLBI) to localize FRBs with milliarcsecond precision. The CHIME site uses a passive hydrogen maser frequency standard in order to minimize localization errors due to clock delay. However, not all outrigger stations will have access to a maser. This report presents techniques used to evaluate clocks for use at outrigger sites without a maser. More importantly, the resulting algorithm provides calibration methods for clocks that do not initially meet the stability requirements for VLBI, thus allowing CHIME/FRB Outriggers to remain true to the goal of having milliarcsecond precision.

Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. High-resolution observations of the CO emission around AGB stars now have the brightness temperature sensitivity to detect possible weak CO maser emission. We used high angular resolution observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the small-scale structure of CO $J=3-2$ emission around the oxygen-rich AGB star W Hya. We find CO maser emission amplifying the stellar continuum with an optical depth $\tau\approx-0.55$. The maser predominantly amplifies the limb of the star because CO $J=3-2$ absorption from the extended stellar atmosphere is strongest towards the centre of the star. The CO maser velocity corresponds to a previously observed variable component of high-frequency H$_2$O masers and with the OH maser that was identified as the amplified stellar image. This implies that the maser originates beyond the acceleration region and constrains the velocity profile since we find the population inversion primarily in the inner circumstellar envelope. We find that inversion can be explained by the radiation field at 4.6 $\mu$m and that the existence of CO maser emission is consistent with the estimated mass-loss rates for W Hya. However, the pumping mechanism requires a complex interplay between absorption and emission lines in the extended atmosphere. Excess from dust in the circumstellar envelope of W Hya is not sufficient to contribute significantly to the required radiation field at 4.6 $\mu$m. The interplay between molecular lines that cause the pumping can be constrained by future multi-level CO observations.

Impacts between planetary-sized bodies can explain the origin of satellites orbiting large ($R>500$~km) trans-Neptunian objects. Their water rich composition, along with the complex phase diagram of water, make it important to accurately model the wide range of thermodynamic conditions material experiences during an impact event and in the debris disk. Since differences in the thermodynamics may influence the system dynamics, we seek to evaluate how the choice of an equation of state (EOS) alters the system's evolution. Specifically, we compare two EOSs that are constructed by different approaches: either by a simplified analytic description (Tillotson), or by interpolation of tabulated data (Sesame). Approximately $50$ pairs of Smoothed Particle Hydrodynamics impact simulations were performed, with similar initial conditions but different EOSs, in the parameter space in which the Pluto-Charon binary is thought to form (slow impacts between Pluto-size, water rich bodies). Generally, we show that impact outcomes (e.g., circumplanetary debris disk) are consistent between EOSs. Some differences arise, importantly in the production of satellitesimals (large intact clumps) that form in the post-impact debris disk. When utilizing an analytic EOS, the emergence of satellitesimals is highly certain, while when using the tabulated EOS it is less common. This is because for the typical densities and energies experienced in these impacts, the analytic EOS predicts very low pressure values, leading to particles artificially aggregating by a tensile instability.

We report seven new eclipse timings for the novalike variable 1RXS J064434.5+334451. An analysis of our data, along with all previously available timings (36 published and 16 unpublished), yields a best-fitting linear ephemeris of BJD$_\mathrm{ecl} = 2,453,403.7611(2) + 0.269~374~43(2)~\mathrm{E}$. We find a somewhat improved fit with a quadratic ephemeris given by: BJD$_\mathrm{ecl} = 2,453,403.7598 + 0.269~374~87~\mathrm{E} - 2.0\times10^{-11}~\mathrm{E}^2$, which suggests that the orbital period may be decreasing at a rate given by $\dot P \simeq -1.5\times10^{-10}$.

A fundamental puzzle of our solar system's formation is understanding why the terrestrial bodies including the planets,comets,and asteroids are depleted in $^{16}$O compared to the Sun. The most favored mechanism,the selective photodissociation of CO gas to produce $^{16}$O depleted water,requires finely tuned mixing timescales to transport $^{16}$O depleted water from the cold outer solar system to exchange isotopically with dust grains to produce the $^{16}$O depleted planetary bodies observed today. Here we show that energetic particle irradiation of SiO$_2$ (and Al$_2$O$_3$) makes them susceptible to anomalous isotope exchange with H$_2$O ice at temperatures as low as 10 K. The observed magnitude of the anomalous isotope exchange (D$^{17}$O) is sufficient to generate the $^{16}$O depletion characteristic of the terrestrial bodies in the solar system. We calculated the cosmic-ray exposure times needed to produce the observed $^{16}$O depletions in silicate (SiO2) dust in the interstellar medium and early solar system and find that radiation damage induced oxygen isotope exchange could have rapidly (~10-100 yrs) depleted dust grains of $^{16}$O during the Sun's T-Tauri phase. Our model explains whythe oldest and most refractory minerals found in the solar system, the anhydrous Calcium with Aluminum Inclusions (CAIs),are generally $^{16}$O enriched compared to chondrules and the bulk terrestrial solids and provides a mechanism for producing $^{16}$O depleted grains very early in the solar system's history. Our findings have broad implications for the distribution of oxygen isotopes in the solar system, the interstellar medium, the formation of the planets and its building blocks as well as the nature of mass-independent isotope effects.

Compilation of papers presented by the VERITAS Collaboration at the 37th International Cosmic Ray Conference (ICRC), held July 12 through July 23, 2021 (online) in Berlin, Germany.

We present the design and lab performance of a prototype lenslet-slicer hybrid integral field spectrograph (IFS), validating the concept for use in future instruments like SCALES/PSI-Red. By imaging extrasolar planets with IFS, it is possible to measure their chemical compositions, temperatures and masses. Many exoplanet-focused instruments use a lenslet IFS to make datacubes with spatial and spectral information used to extract spectral information of imaged exoplanets. Lenslet IFS architecture results in very short spectra and thus low spectral resolution. Slicer IFSs can obtain higher spectral resolution but at the cost of increased optical aberrations that propagate through the down-stream spectrograph and degrade the spatial information we can extract. We have designed a lenslet/slicer hybrid that combines the minimal aberrations of the lenslet IFS with the high spectral resolution of the slicer IFS. The slicer output f/\# matches the lenslet f/\# requiring only additional gratings.

The Cherenkov Telescope Array (CTA) is the next-generation ground-based observatory for very-high-energy gamma-ray astronomy. An innovative 9.7 m aperture, dual-mirror Schwarzschild-Couder Telescope (SCT) design is a candidate design for CTA Medium-Sized Telescopes. A prototype SCT (pSCT) has been constructed at the Fred Lawrence Whipple Observatory in Arizona, USA. Its camera is currently partially instrumented with 1600 pixels covering a field of view of 2.7 degrees square. The small plate scale of the optical system allows densely packed silicon photomultipliers to be used, which combined with high-density trigger and waveform readout electronics enable the high-resolution camera. The camera's electronics are capable of imaging air shower development at a rate of one billion samples per second. We describe the commissioning and performance of the pSCT camera, including trigger and waveform readout performance, calibration, and absolute GPS time stamping. We also present the upgrade to the camera, which is currently underway. The upgrade will fully populate the focal plane, increasing the field of view to 8 degree diameter, and lower the front-end electronics noise, enabling a lower trigger threshold and improved reconstruction and background rejection.

The VHE component from at least two GRBs, i.e., GRB180720B and GRB190114C, has been detected in the afterglow phase. We systematically analyzed 199 GRBs detected by Fermi-LAT during 2008-2019. If an additional high-energy component exists in the afterglows of Fermi-LAT GRBs, the best-fit spectral model could be a broken power-law (BPL) model with an upturn above a break energy. We compare the afterglow spectra using PL and BPL representations. Out of the 30 GRBs with >10GeV photons that arrived after T90, 25 GRBs are tentatively or significantly detected at 0.1-200 GeV after 2*T90. The spectrum of GRB131231A shows an upturn above a break of 1.6+-0.8~GeV, supporting the BPL model. For GRB131231A, we performed a modeling of its X-ray and gamma-ray spectra, and found that the SSC model can explain the upturn with acceptable parameter values. In the cases of GRBs 190114C, 171210A, 150902A, 130907A, 130427A, and 090902B, the improvement of the BPL fit compared to the PL fit is tentative or marginal. There is no conclusive evidence that an additional higher energy component commonly exists in Fermi-LAT GRB afterglows, except for a group of Fermi-LAT GRBs mentioned above. Such an additional high-energy component may be explained by the synchrotron self-Compton mechanism. Current and future VHE observations will provide important constraints on the issue.

AstronomicAL is a human-in-the-loop interactive labelling and training dashboard that allows users to create reliable datasets and robust classifiers using active learning. This technique prioritises data that offer high information gain, leading to improved performance using substantially less data. The system allows users to visualise and integrate data from different sources and deal with incorrect or missing labels and imbalanced class sizes. AstronomicAL enables experts to visualise domain-specific plots and key information relating both to broader context and details of a point of interest drawn from a variety of data sources, ensuring reliable labels. In addition, AstronomicAL provides functionality to explore all aspects of the training process, including custom models and query strategies. This makes the software a tool for experimenting with both domain-specific classifications and more general-purpose machine learning strategies. We illustrate using the system with an astronomical dataset due to the field's immediate need; however, AstronomicAL has been designed for datasets from any discipline. Finally, by exporting a simple configuration file, entire layouts, models, and assigned labels can be shared with the community. This allows for complete transparency and ensures that the process of reproducing results is effortless

The recent measurement by LHAASO of gamma-ray emission extending up to 100s of TeV from multiple Galactic sources represents a major observational step forward in the search for the origin of the Galactic cosmic rays. The burning question is if this ultra-high-energy emission is associated to acceleration of protons and/or nuclei to PeV energies, or can be associated to PeV-electron accelerators. A strong Klein-Nishina suppression of inverse Compton emission at these energies is unavoidable, nevertheless we show here that inverse Compton emission can provide a natural explanation of the measured emission, and that an association to the established PeV-electron accelerating source class of pulsar wind nebulae is also rather natural. However, a clear distinction between different models requires taking into account multi-wavelength data, having good knowledge about the local environmental conditions and, in some cases, to perform multi-source modeling.

We report on the discovery of Swift J011511.0-725611, a rare Be X-ray binary system (BeXRB) with a White Dwarf (WD) compact object, in the Small Magellanic Cloud (SMC) by S-CUBED, a weekly X-ray/UV survey of the SMC by the Neil Gehrels Swift Observatory. Observations show an approximately 3 month outburst from Swift J011511.0-725611, the first detected by S-CUBED since it began in 2016 June. Swift J011511.0-725611 shows super-soft X-ray emission, indicative of a White Dwarf compact object, which is further strengthened by the presence of an 0.871 keV edge, commonly attributed to O viii K-edge in the WD atmosphere. Spectroscopy by SALT confirms the Be nature of the companion star, and long term light-curve by OGLE finds both the signature of a circumstellar disk in the system at outburst time, and the presence of a 17.4 day periodicity, likely the orbital period of the system. Swift J011511.0-725611 is suggested to be undergoing a Type-II outburst, similar to the previously reported SMC Be White Dwarf binary (BeWD), Swift J004427.3-734801. It is likely that the rarity of known BeWD is in part due to the difficulty in detecting such outbursts due to both their rarity, and their relative faintness compared to outbursts in Neutron Star BeXRBs.

Could we enjoy starry skies in our cities again? Arguably yes. The actual number of visible stars will depend, among other factors, on the spatial density of the overall city light emissions. In this paper it is shown that reasonably dark skies could be achieved in urban settings, even at the center of large metropolitan areas, if the light emissions are kept within admissible levels and direct glare from the light sources is avoided. These results may support the adoption of science-informed, democratic public decisions on the use of light in our municipalities, with the goal of recovering the possibility of contemplating the night sky everywhere in our planet.

We present a revised distance to the Cygnus Loop supernova remnant of $725\pm15$ pc based on Gaia Early Data Release 3 parallax measurements (EDR3) for several stars previously found to be located either inside or behind the supernova based on the presence of high-velocity absorption lines in their spectra. This revised distance estimate and error means the Cygnus Loop remnant now has an estimated distance uncertainty comparable to that of its $\simeq$18 pc radius.

We present the broadband spectral analysis of all the six hard, intermediate and soft state NuSTAR observations of the recently discovered transient black hole X-ray binary MAXI J1348-630 during its first outburst in 2019. We first model the data with a combination of a multi-colour disc and a relativistic blurred reflection, and, whenever needed, a distant reflection. We find that this simple model scheme is inadequate in explaining the spectra, resulting in a very high iron abundance. We, therefore, explore the possibility of reflection from a high-density disc. We use two different sets of models to describe the high-density disc reflection: relxill-based reflection models, and reflionx-based ones. The reflionx-based high-density disc reflection models bring down the iron abundance to around the solar value, while the density is found to be $10^{20.3-21.4} \rm cm^{-3}$. We also find evidence of a high-velocity outflow in the form of $\sim$7.3 keV absorption lines. The consistency between the best-fit parameters for different epochs and the statistical significance of the corresponding model indicates the existence of high-density disc reflection in MAXI J1348-630.

AQ Col (EC~05217-3914) is one of the first detected pulsating subdwarf B (sdB) stars and has been considered to be a single star. However, its periodic pulsation timing variations indicate that AQ Col may not be a single star. We present pulsation period variations observed over twenty-four years and derived orbital characteristics these would imply if these were a consequence of AQ Col being a pulsating hot subdwarf in a long-period binary. The derived orbital period is P = 486.0 days. In the sdB star binary evolution scenario, a Roche lobe overflow channel results in long period (450 < P < 1400 d) for sdB + Main Sequence (MS) binaries. However the derived orbital eccentricity of the system is 0.424, which is too large for a typical long period sdB+MS system. The Skymapper u - z vs. z - WISE W1 diagram is incompatible with sdB+MS binary systems, and suggests the system contains a white dwarf or other hot and faint object. The expected radial velocity amplitude of AQ Col due to this orbital motion is ~15 km/s. However, the radial velocity amplitude differences obtained from spectroscopy show that the amplitude could be more than ~300 km/s, which indicates the possibility that AQ Col also has a short period companion with orbital period of ~1 day. Therefore, the AQ Col system may be a triple star system. Because such systems have not yet been studied in detail, AQ Col may offer unique insight into the production of sdB stars and this system deserves continued time-series and spectroscopic monitoring.

The prime Kepler mission detected 34,032 transit-like signals, out of which 8,054 were identified as likely due to astrophysical planet transits or eclipsing binaries. We manually examined 306 of the remaining 25,978 detections, and found six plausible transiting or eclipsing objects, five of which are plausible planet candidates (PCs), and one stellar companion. One of our new PCs is a possible new second planet in the KOI 4302 system. Another new PC is a possible new planet around the KOI 4246, and when combined with a different possible planet rescued by the False Positive Working Group, we find that KOI 4246 may be a previously unrecognized three-planet system. \end{abstract}

The history of rivers on Mars is an important constraint on Martian climate evolution. The timing of relatively young, alluvial fan-forming rivers is especially important, as Mars' Amazonian atmosphere is thought to have been too thin to consistently support surface liquid water. Previous regional studies suggested that alluvial fans formed primarily between the Early Hesperian and the Early Amazonian. In this study, we describe how a combination of a global impact crater database, a global geologic map, a global alluvial fan database, and statistical models can be used to estimate the timing of alluvial fan formation across Mars. Using our global approach and improved statistical modeling, we find that alluvial fan formation likely persisted into the last ~2.5 Gyr, well into the Amazonian period. However, the data we analyzed was insufficient to place constraints on the duration of alluvial fan formation. Going forward, more crater data will enable tighter constraints on the parameters estimated in our models and thus further inform our understanding of Mars' climate evolution.

Axion as a coherently oscillating massive scalar field is known to behave as a zero-pressure irrotational fluid with characteristic quantum stress on a small scale. In relativistic perturbation theory, the case was most conveniently proved in the axion-comoving gauge up to fully nonlinear and exact order. Our basic assumption is that the Compton wavelength is smaller than the horizon scale. Here, we revisit the relativistic proof to the linear order in the other gauge conditions. The comoving gauge, the zero-shear gauge, and the uniform-curvature gauge give {\it the same} equation for density perturbation known in the non-relativistic quantum mechanical treatment in {\it all} scales. On the other hand, the quantum stress term is missing in the synchronous gauge, and inconsistency is found in the uniform-expansion gauge. In the absence of quantum stress, the simple density perturbation equations of the axion in the zero-shear gauge and the uniform-curvature gauge were {\it not} expected. Even in the zero-pressure fluid, the equations in the two gauges coincide with the one in the comoving gauge {\it only} in the sub-horizon scale. We clarify that our analysis is valid for scales larger than the Compton wavelength, which is negligible compared with the cosmological scale. For comparison, we review the non-relativistic quantum hydrodynamics and present the Schr\"odinger equation to first-order post-Newtonian expansion in the cosmology context.

Several satellite missions have uncovered a series of potential anomalies in the fluctuation spectrum of the cosmic microwave background temperature, including: (1) an unexpectedly low level of correlation at large angles, manifested via the angular correlation function, C(theta); and (2) missing power in the low multipole moments of the angular power spectrum, C_ell. Their origin is still debated, however, due to a persistent lack of clarity concerning the seeding of quantum fluctuations in the early Universe. A likely explanation for the first of these appears to be a cutoff, k_min=(3.14 +/- 0.36) x 10^{-4} Mpc^{-1}, in the primordial power spectrum, P(k). Our goal in this paper is twofold: (1) we examine whether the same k_min can also self-consistently explain the missing power at large angles, and (2) we confirm that the of this cutoff in P(k) does not adversely affect the remarkable consistency between the prediction of Planck-LCDM and the Planck measurements at ell > 30. We use the publicly available code CAMB to calculate the angular power spectrum, based on a line-of-sight approach. The code is modified slightly to include the additional parameter (i.e., k_min) characterizing the primordial power spectrum. In addition to this cutoff, the code optimizes all of the usual standard-model parameters. In fitting the angular power spectrum, we find an optimized cutoff, k_min = 2.04^{+1.4}_{-0.79} x 10^{-4} Mpc^{-1}, when using the whole range of ell's, and k_min=3.3^{+1.7}_{-1.3} x 10^{-4} Mpc^{-1}, when fitting only the range ell < 30, where the Sachs-Wolfe effect is dominant. These are fully consistent with the value inferred from C(theta), suggesting that both of these large-angle anomalies may be due to the same truncation in P(k).

The Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) spacecraft is the third NASA New Frontiers Program mission and arrived at the near-Earth asteroid (101955) Bennu in December 2018. Following completion of sample collection in October 2020, otherwise known as Touch-And-Go (TAG), the OSIRIS-REx spacecraft was set to verify its collected sample mass requirement (> 60g of material). The thoroughly tested Sample Mass Measurement (SMM) method was to be used for this verification. Imaging of the Touch-And-Go Sample Acquisition Mechanism (TAGSAM) was received shortly following the TAG event, intended to ensure mechanism health prior to moving forward with the SMM activity. These images displayed sample leakage, prompting discussion for alternative paths forward. Risk of continued sample loss and a desire to retain as much material as possible lead the team to pursue an accelerated sample stow schedule and forgo the planned SMM activity. Once the sample was safely stowed in the return capsule an alternative SMM method was proposed. The alternative SMM technique utilized reaction wheel momentum data from identical TAGSAM movements prior to and following the TAG event to estimate changes in spacecraft moment of inertia. Conservation of momentum was used to isolate the sample mass from this inertia change. Using this new method, the spacecraft team was able to successfully estimate collected sample mass to be 250.37 +/- 101 g.

The lithium abundance, A(Li), in stellar atmospheres suffers from various enhancement and depletion processes during the star's lifetime. While several studies have demonstrated that these processes are linked to the physics of stellar formation and evolution, the role that Galactic-scale events play in the galactic A(Li) evolution is not yet well understood. We aim to demonstrate that the observed A(Li) bi-modal distribution, in particular in the FGK-dwarf population, is not a statistical artefact and that the two populations connect through a region with a low number of stars. We also want to investigate the role that Galactic-scale events play in shaping the A(Li) distribution of stars in the thin disk. We use statistical techniques along with a Galactic chemical evolution model for A(Li) that includes most of the well-known $^7$Li production and depletion channels. We confirm that the FGK main-sequence stars belonging to the Milky Way's thin disk present a bi-modal A(Li) distribution. We demonstrate that this bi-modality can be generated by a particular Milky Way star formation history profile combined with the stellar evolution's $^7$Li depletion mechanisms. We show that A(Li) evolution can be used as an additional proxy for the star formation history of our Galaxy.

We report observations of a relatively long period of 3He-rich solar energetic particles (SEPs) measured by Solar Orbiter. The period consists of several well-resolved ion injections. The high-resolution STEREO-A imaging observations reveal that the injections coincide with EUV jets/brightenings near the east limb, not far from the nominal magnetic connection of Solar Orbiter. The jets originated in two adjacent, large, and complex active regions as observed by the Solar Dynamics Observatory when the regions rotated to the Earth's view. It appears that the sustained ion injections were related to the complex configuration of the sunspot group and the long period of 3He-rich SEPs to the longitudinal extent covered by the group during the analyzed time period.

Neutron stars are natural physical laboratories allowing us to study a plethora of phenomena in extreme conditions. In particular, these compact objects can have very strong magnetic fields with non-trivial origin and evolution. In many respects its magnetic field determines the appearance of a neutron star. Thus, understanding the field properties is important for interpretation of observational data. Complementing this, observations of diverse kinds of neutron stars enable us to probe parameters of electro-dynamical processes at scales unavailable in terrestrial laboratories. In this review we first briefly describe theoretical models of formation and evolution of magnetic field of neutron stars, paying special attention to field decay processes. Then we present important observational results related to field properties of different types of compact objects: magnetars, cooling neutron stars, radio pulsars, sources in binary systems. After that, we discuss which observations can shed light on obscure characteristics of neutron star magnetic fields and their behaviour. We end the review with a subjective list of open problems.

Context. While most (if not all) Type I Galactic globular clusters (GGCs) are characterised by spreads in the abundances of light chemical elements (e.g. Li, N, O, Na, Mg, Al), it is not yet well established whether similar spreads may exist in s-process elements as well. Aims. We investigated the possible difference in Ba abundance between the primordial (1P) and polluted (2P) stars in the Galactic globular cluster (GGC) 47 Tuc (NGC 104). For this, we obtained homogeneous abundances of Fe, Na, and Ba in a sample of 261 red giant branch (RGB) stars which is the largest sample used for Na and Ba abundance analysis in any GGC so far. Methods. Abundances of Na and Ba were determined using archival GIRAFFE/VLT spectra and 1D non-local thermodynamic equilibrium (NLTE) abundance analysis methodology. Results. Contrary to the finding of Gratton et al. (2013), we did not detect any significant Ba-Na correlation or 2P-1P Ba abundance difference in the sample of 261 RGB stars in 47 Tuc. This corroborates the result of D'Orazi et al. (2010) who found no statistically significant Ba-Na correlation in 110 RGB stars in this GGC. The average barium-to-iron ratio obtained in the sample of 261 RGB stars, $\langle{\rm Ba/Fe}_{\rm 1D~NLTE}\rangle = -0.01\pm0.06$, agrees well with those determined in Galactic field stars at this metallicity and may therefore represent the abundance of primordial proto-cluster gas that has not been altered during the subsequent chemical evolution of the cluster.

In solar filament formation mechanisms, magnetic reconnection between two sets of sheared arcades forms helical structures of the filament with numerous magnetic dips, and cooling and condensation of plasma trapped inside the helical structures supply mass to the filament. Although each of these processes, namely, magnetic reconnection and coronal condensation have been separately reported, observations that show the whole process of filament formation are rare. In this Letter, we present the formation of a sigmoid via reconnection between two sets of coronal loops, and the subsequent formation of a filament through cooling and condensation of plasma inside the newly formed sigmoid. On 2014 August 27, a set of loops in the active region 12151 reconnected with another set of loops that are located to the east. A longer twisted sigmoidal structure and a set of shorter lower-lying loops then formed. The observations coincide well with the tether-cutting model. The newly formed sigmoid remains stable and does not erupt as a coronal mass ejection. From the eastern endpoint, signatures of injection of material into the sigmoid (as brightenings) are detected, which closely outline the features of increasing emission measure at these locations. This may indicate the chromospheric evaporation caused by reconnection, supplying heated plasma into the sigmoid. In the sigmoid, thermal instability occurs, and rapid cooling and condensation of plasma take place, forming a filament. The condensations then flow bi-directionally to the filament endpoints. Our results provide a clear observational evidence of the filament formation via magnetic reconnection and coronal condensation.

We performed one-dimensional force-free magnetodynamic numerical simulations of the propagation of Alfven waves along magnetic field lines around a spinning black-hole-like object, the Banados--Teitelboim--Zanelli black string, to investigate the dynamic process of wave propagation and energy transport with Alfven waves. We considered axisymmetric and stationary magnetosphere and perturbed the background magnetosphere to obtain the linear wave equation for the Alfven wave mode. The numerical results show that the energy of Alfven waves monotonically increases as the waves propagate outwardly along the rotating curved magnetic field line around the ergosphere, where energy seems not to be conserved, in the case of energy extraction from the black string by the Blandford--Znajek mechanism. The apparent breakdown of energy conservation suggests the existence of an additional wave induced by the Alfven wave. Considering the additional fast magnetosonic wave induced by the Alfven wave, the energy conservation is recovered. Similar relativistic phenomena, such as the amplification of Alfven waves and induction of fast magnetosonic waves, are expected around a spinning black hole.

We analyze thermal emission spectra using the 2001 Mars Odyssey Thermal Emission Imaging System (THEMIS) and the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) to characterize grain size and mineralogical composition of dunes at Hargraves crater, Mars. Thermal inertia and bulk composition of the dunes were compared to inferred provenances from the thermal infrared response of surface constituent materials. We use a Markov Chain Monte Carlo (MCMC) technique to estimate the bulk amount of mineralogy contributed by each inferred provenance to the dune field composition. An average thermal inertia value of 238+/-17 Jm-2K-1s-0.5 was found for the dunes corresponding to a surface composed of an average effective grain size of ~391+/-172 um. This effective particle size suggests the presence of mostly medium sand-sized materials mixed with fine and coarse grain sands. The dunes are likely comprised of a weakly indurated surface mixed with unconsolidated materials. Compositional analysis specifies that the dunes are comprised of a mixture of feldspar, olivine, pyroxene, and relatively low bulk-silica content. Dune materials were likely derived from physical weathering, especially eolian erosion, predominantly from the crater ejecta unit at the crater, mixed with a small amount from the crater floor and crater rim and wall lithologies - indicating the dune materials were likely sourced locally.

Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$\alpha$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using measurements from $\sim$27,000 APOGEE stars (R=22,500, SNR$>$200), we build simple local linear models to predict a sample of elements (X = Si, O, Ca, Ti, Ni, Al, Mn, Cr) using (Fe, Mg) abundances alone, as fiducial tracers of supernovae production channels. Given [Fe/H] and [Mg/H], we predict these elements, [X/H], to about double the uncertainty of their measurements. The intrinsic dispersion, after subtracting measurement errors in quadrature is $\approx 0.015-0.04$~dex. The residuals of the prediction (measurement $-$ model) for each element demonstrate that each element has an individual link to birth properties at fixed (Fe, Mg). Residuals from primarily massive-star supernovae (i.e. Si, O, Al) partially correlate with guiding radius. Residuals from primarily supernovae Ia (i.e. Mn, Ni) partially correlate with age. A fraction of the intrinsic scatter that persists at fixed (Fe, Mg), however, after accounting for correlations, does not appear to further discriminate between birth properties that can be traced with present-day measurements. Presumably, this is because the residuals are also, in part, a measure of the typical (in)-homogeneity of the disk's stellar birth environments, previously inferred only using open-cluster systems. Our study implies at fixed birth radius and time, there is a median scatter of $\approx 0.01-0.015$ dex in elements generated in supernovae sources.

We have conducted mapping observations toward the n3 and n5 positions in the NGC\,2264-D cluster-forming region with the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. Observations with 10000 au scale beam reveal the chemical composition at the clump scale. The spatial distributions of the observed low upper-state-energy lines of CH$_{3}$OH are similar to those of CS and SO, and the HC$_{3}$N emission seems to be predominantly associated with clumps containing young stellar objects. The turbulent gas induced by the star formation activities produces large-scale shock regions in NGC\,2264-D, which are traced by the CH$_{3}$OH, CS and SO emissions. We derive the HC$_{3}$N, CH$_{3}$CN, and CH$_{3}$CHO abundances with respect to CH$_{3}$OH. Compared to the n5 field, the n3 field is farther (in projected apparent distance) from the neighboring NGC\,2264-C, yet the chemical composition in the n3 field tends to be similar to that of the protostellar candidate CMM3 in NGC\,2264-C. The HC$_{3}$N/CH$_{3}$OH ratios in the n3 field are higher than those in the n5 field. We find an anti-correlation between the HC$_{3}$N/CH$_{3}$OH ratio and their excitation temperatures. The low HC$_{3}$N/CH$_{3}$OH abundance ratio at the n5 field implies that the n5 field is an environment with more active star formation compared with the n3 field.

We revisit the geometrical meaning of statistical isotropy that is manifest in excursion sets of smooth random fields in two dimensions. Using the contour Minkowski tensor, $\W_1$, as our basic tool we first examine geometrical properties of single structures. For simple closed curves in two dimensions we show that $\W_1$ is proportional to the identity matrix if the curve has $m$-fold symmetry, with $m\ge 3$. Then we elaborate on how $\W_1$ maps any arbitrary shaped simple closed curve to an ellipse that is unique up to translations of its centroid. We also carry out a comparison of the shape parameters, $\alpha$ and $\beta$, defined using $\W_1$, with the filamentarity parameter defined using two scalar Minkowski functionals - area and contour length. We show that they contain complementary shape information, with $\W_1$ containing additional information of orientation of structures. Next, we apply our method to boundaries of excursion sets of random fields and examine what statistical isotropy means for the geometry of the excursion sets. Focusing on Gaussian isotropic fields, and using a semi-numerical approach we quantify the effect of finite sampling of the field on the geometry of the excursion sets. In doing so we obtain an analytic expression for $\alpha$ which takes into account the effect of finite sampling. Finally we derive an analytic expression for the ensemble expectation of $\W_1$ for Gaussian anisotropic random fields. Our results provide insights that are useful for designing tests of statistical isotropy using cosmological data.

We present the experimental phase function, degree of linear polarization (DLP), and linear depolarization (deltaL) curves of a set of forsterite samples representative of low-absorbing cosmic dust particles. The samples are prepared using state-of-the-art size-segregating techniques to obtain narrow size distributions spanning a broad range of the scattering size parameter domain. We conclude that the behavior of the phase function at the side- and back-scattering regions provides information on the size regime, the position and magnitude of the maximum of the DLP curve are strongly dependent on particle size, the negative polarization branch is mainly produced by particles with size parameters in the approx. 6 to 20 range, and the deltaL is strongly dependent on particle size at all measured phase angles except for the exact backward direction. From a direct comparison of the experimental data with computations for spherical particles, it becomes clear that the use of the spherical model for simulating the phase function and DLP curves of irregular dust produces dramatic errors in the retrieved composition and size of the scattering particles: The experimental phase functions are reproduced by assuming unrealistically high values of the imaginary part of the refractive index. The spherical model does not reproduce the bell-shaped DLP curve of dust particles with sizes in the resonance and/or geometric optics size domain. Thus, the use of the Mie model for analyzing polarimetric observations might prevent locating dust particles with sizes of the order of or larger than the wavelength of the incident light.

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Mini-Array (MA) project is an international collaboration led by the Italian National Institute for Astrophysics (INAF). ASTRI MA is composed of nine Cherenkov telescopes operating in the energy range 1-100 TeV, and it aims to study very high-energy gamma ray astrophysics and optical intensity interferometry of bright stars. ASTRI MA is currently under construction, and will be installed at the site of the Teide Observatory in Tenerife (Spain). The hardware and software system that is responsible of monitoring and controlling all the operations carried out at the ASTRI MA site is the Supervision Control and Data Acquisition (SCADA). The LOgging UnifieD (LOUD) subsystem is one of the main components of SCADA. It provides the service responsible for collecting, filtering, exposing and storing log events collected by all the array elements (telescopes, LIDAR, devices, etc.). In this paper, we present the LOUD architecture and the software stack explicitly designed for distributed computing environments exploiting Internet of Things technologies (IoT).

We present the current development of the Monitoring, Logging and Alarm subsystems in the framework of the Array Control and Data Acquisition System (ACADA) for the Cherenkov Telescope Array (CTA). The Monitoring System (MON) is the subsystem responsible for monitoring and logging the overall array (at each of the CTA sites) through the acquisition of monitoring and logging information from the array elements. The MON allows us to perform a systematic approach to fault detection and diagnosis supporting corrective and predictive maintenance to minimize the downtime of the system. We present a unified tool for monitoring data items from the telescopes and other devices deployed at the CTA array sites. Data are immediately available for the operator interface and quick-look quality checks and stored for later detailed inspection. The Array Alarm System (AAS) is the subsystem that provides the service that gathers, filters, exposes, and persists alarms raised by both the ACADA processes and the array elements supervised by the ACADA system. It collects alarms from the telescopes, the array calibration, the environmental monitoring instruments and the ACADA systems. The AAS sub-system also creates new alarms based on the analysis and correlation of the system software logs and the status of the system hardware providing the filter mechanisms for all the alarms. Data from the alarm system are then sent to the operator via the human-machine interface.

The remnant of SN 1987A is the best-studied object of its kind. The rich data-set of its thermal and non-thermal emission across the electromagnetic spectrum poses a unique testbed for the elaboration of particle-acceleration theory. We use 2D simulations of the progenitor's wind to obtain hydro-profiles for the medium around the supernova explosion. Various cones along prominent features of the ambient medium are then used in our time-dependent acceleration code RATPaC to model the evolution of the emission of SN 1987A and compare it to observational data. We solve for the transport of cosmic rays and the hydrodynamical flow, in the test-particle limit.The simulation code relies on 1D profiles but the large expansion speed of the young remnant renders lateral transport unimportant. We find that the increase in thermal X-ray emission predates the increase in the low-energy gamma-ray brightness by several years. The increase of the gamma-ray brightness at lower energies is followed by a smooth increase at the highest energies. The gamma-ray spectrum at the highest energies appears soft during the brightening but hardens as more material in the equatorial ring gets shocked. The X-ray and gamma-ray brightness remain almost constant once the SNR blast-wave passed the region of peak-density in the equatorial plane.

Light axion-like particles (ALPs) are expected to be abundantly produced in core-collapse supernovae (CCSNe), resulting in a $\sim$10-second long burst of ALPs. These particles subsequently undergo conversion into gamma-rays in external magnetic fields to produce a long gamma-ray burst (GRB) with a characteristic spectrum peaking in the 30--100-MeV energy range. At the same time, CCSNe are invoked as progenitors of {\it ordinary} long GRBs, rendering it relevant to conduct a comprehensive search for ALP spectral signatures using the observations of long GRB with the \textit{Fermi} Large Area Telescope (LAT). We perform a data-driven sensitivity analysis to determine CCSN distances for which a detection of an ALP signal is possible with the LAT's low-energy (LLE) technique which, in contrast to the standard LAT analysis, allows for a a larger effective area for energies down to 30~MeV. Assuming an ALP mass $m_a \lesssim 10^{-10}$~eV and ALP-photon coupling $g_{a\gamma} = 5.3\times 10^{-12}$ GeV$^{-1}$, values considered and deduced in ALP searches from SN1987A, we find that the distance limit ranges from $\sim\!0.5$ to $\sim\!10$~Mpc, depending on the sky location and the CCSN progenitor mass. Furthermore, we select a candidate sample of twenty-four GRBs and carry out a model comparison analysis in which we consider different GRB spectral models with and without an ALP signal component. We find that the inclusion of an ALP contribution does not result in any statistically significant improvement of the fits to the data. We discuss the statistical method used in our analysis and the underlying physical assumptions, the feasibility of setting upper limits on the ALP-photon coupling, and give an outlook on future telescopes in the context of ALP searches.

The Cherenkov Telescope Array (CTA), conceived as an array of tens of imaging atmospheric Cherenkov telescopes (IACTs), is an international project for a next-generation ground-based gamma-ray observatory, aiming to improve on the sensitivity of current-generation instruments a factor of five to ten and provide energy coverage from 20 GeV to more than 300 TeV. Arrays of IACTs probe the very-high-energy gamma-ray sky. Their working principle consists of the simultaneous observation of air showers initiated by the interaction of very-high-energy gamma rays and cosmic rays with the atmosphere. Cherenkov photons induced by a given shower are focused onto the camera plane of the telescopes in the array, producing a multi-stereoscopic record of the event. This image contains the longitudinal development of the air shower, together with its spatial, temporal, and calorimetric information. The properties of the originating very-high-energy particle (type, energy, and incoming direction) can be inferred from those images by reconstructing the full event using machine learning techniques. In this contribution, we present a purely deep-learning driven, full-event reconstruction of simulated, stereoscopic IACT events using CTLearn. CTLearn is a package that includes modules for loading and manipulating IACT data and for running deep learning models, using pixel-wise camera data as input.

A recent time-integrated analysis of a catalog of 110 candidate neutrino sources revealed a cumulative neutrino excess in the data collected by IceCube between April 6, 2008 and July 10, 2018. This excess, inconsistent with the background hypothesis in the Northern hemisphere at the $3.3~\sigma$ level, is associated with four sources: NGC 1068, TXS 0506+056, PKS 1424+240 and GB6 J1542+6129. This letter presents two time-dependent neutrino emission searches on the same data sample and catalog: a point-source search that looks for the most significant time-dependent source of the catalog by combining space, energy and time information of the events, and a population test based on binomial statistics that looks for a cumulative time-dependent neutrino excess from a subset of sources. Compared to previous time-dependent searches, these analyses enable a feature to possibly find multiple flares from a single direction with an unbinned maximum-likelihood method. M87 is found to be the most significant time-dependent source of this catalog at the level of $1.7~\sigma$ post-trial, and TXS 0506+056 is the only source for which two flares are reconstructed. The binomial test reports a cumulative time-dependent neutrino excess in the Northern hemisphere at the level of $3.0~\sigma$ associated with four sources: M87, TXS 0506+056, GB6 J1542+6129 and NGC 1068.

The role of non-ideal magnetohydrodynamics has been proven critical during the formation of the protoplanetary disk, particularly in regulating its size. We provide a simple model to predict the disk size under the interplay among the ambipolar diffusion, the Hall effect, and the Ohmic dissipation. The model predicts a small disk size of around 20 AU, that depends only sub-linearly on disk parameters, for a wide range of initial conditions of sub-Solar mass and moderate magnetization. It is able to explain phenomena manifested in existing numerical simulations, including the bimodal disk behavior under parallel and anti-parallel alignment between the rotation and magnetic field. In the parallel configuration, the disk size decreases and eventually disappears. In the anti-parallel configuration, and the disk has an outer partition (or pseudo-disk) that is flat, shrinking , and short-lived, as well as a inner partition that grows slowly with mass and is long-lived. Even with significant initial magnetization, the vertical field in the disk can only dominate at the early stage when the mass is low, and the toroidal field eventually dominates in all disks.

The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their model uncertainties and branching ratios. In this work, we show how one could add an additional and independent constraint to model selection by using the stochastic gravitational-wave background. In contrast to model selection analyses that have discriminating power only up to the gravitational-wave detector horizons (currently at redshifts $z\lesssim 1$ for LIGO-Virgo), the stochastic gravitational-wave background accounts for the redshift integration of all gravitational-wave signals in the Universe. As a working example, we consider the branching ratio results from a model selection study that includes potential contribution from astrophysical and primordial channels. We renormalize the relative contribution of each channel to the detected event rate to compute the total stochastic gravitational-wave background energy density. The predicted amplitude lies below the current observational upper limits of GWTC-2 by LIGO-Virgo, indicating that the results of the model selection analysis are not ruled out by current background limits. Furthermore, given the set of population models and inferred branching ratios, we find that, even though the predicted background will not be detectable by current generation gravitational-wave detectors, it will be accessible by third-generation detectors such as the Einstein Telescope and space-based detectors such as LISA.

Radioactive nuclei were present in the early Solar System, as inferred from analysis of meteorites. Many are produced in massive stars, either during their lives or their final explosions. In the first paper in this series (Brinkman et al. 2019), we focused on the production of $^{26}$Al in massive binaries. Here, we focus on the production of another two short-lived radioactive nuclei, $^{36}$Cl and $^{41}$Ca, and the comparison to the early Solar System data. We used the MESA stellar evolution code with an extended nuclear network and computed massive (10-80 M$ _{\odot} $), rotating (with initial velocities of 150 and 300 km/s) and non-rotating single stars at solar metallicity (Z=0.014) up to the onset of core collapse. We present the wind yields for the radioactive isotopes $^{26}$Al, $^{36}$Cl, and $^{41}$Ca, and the stable isotopes $^{19}$F and $^{20}$Ne. In relation to the stable isotopes, we find that only the most massive models, $\geq$ 60M$_{\odot}$ and $\geq$ 40M$_{\odot}$ give positive $^{19}$F and $^{20}$Ne yields, respectively, depending on the initial rotation rate. In relation to the radioactive isotopes, we find that the early Solar System abundances of $^{26}$Al and $^{41}$Ca can be matched with by models with initial masses $\geq$40M$_{\odot}$, while $^{36}$Cl is matched only by our most massive models, $\geq$60M$_{\odot}$. $^{60}$Fe is not significantly produced by any wind model, as required by the observations. Therefore, massive star winds are a favoured candidate for the origin of the very short-lived $^{26}$Al, $^{36}$Cl, and $^{41}$Ca in the early Solar System.

Knowledge of solar irradiance variability is critical to Earth's climate models and understanding the solar influence on Earth's climate. Direct solar irradiance measurements are only available since 1978. Reconstructions of past variability typically rely on sunspot data. These provide only indirect information on the facular and network regions, which are decisive contributors to irradiance variability on timescales of the solar cycle and longer. Our ultimate goal is to reconstruct past solar irradiance variations using historical full-disc Ca II K observations to describe the facular contribution independently of sunspot observations. Here, we develop the method and test it extensively by using modern CCD-based Ca II K observations and carry out initial tests on two photographic archives. We employ carefully reduced and calibrated Ca II K images from 13 datasets, such as those from the Meudon, Mt Wilson, and Rome observatories. We convert them to unsigned magnetograms and then use them as input to the adapted SATIRE model to reconstruct TSI variations over the period 1978-2019, for which direct irradiance measurements are available. The reconstructed TSI from the analysed Ca II K archives agrees well with direct TSI measurements and existing reconstructions. The model also returns good results on data taken with different bandpasses and images with low spatial resolution. Historical Ca II K archives suffer from numerous inconsistencies, but we show that these archives can still be used to reconstruct TSI with reasonable accuracy provided the observations are accurately processed. By using the unsigned magnetograms of the Sun reconstructed from high-quality Ca II K observations as input into the SATIRE model, we can reconstruct solar irradiance variations nearly as accurately as from directly recorded magnetograms.

Based on updated pulsation models for Classical Cepheids, computed for various assumptions about the metallicity and helium abundance, roughly representative of pulsators in the Small Magellanic Cloud ($Z$=$0.004$ and $Y$=$0.25$), Large Magellanic Cloud ($Z$=$0.008$ and $Y$=$0.25$), and M31 ($Z$=$0.03$ and $Y$=$0.28$), and self-consistent updated evolutionary predictions, we derived Period-Age and multi-band Period-Age-Color relations that also take into account variations in the Mass-Luminosity relation. These results, combined with those previously derived for Galactic Cepheids, were used to investigate the metallicity effect when using these variables as age indicators. In particular, we found that a variation in the metal abundance affects both the slope and the zero point of the above-mentioned relations. The new relations were applied to a sample of Gaia Early Data Release 3 Classical Cepheids. The retrieved distribution of the individual ages confirms that a brighter Mass-Luminosity relation produces older ages and that First Overtone pulsators are found to be concentrated towards older ages with respect to the Fundamental ones at a fixed Mass-Luminosity relation. Moreover, the inclusion of a metallicity term in the Period-Age and Period-Age-Color relations slightly modifies the predicted ages. In particular, the age distribution of the selected sample of Galactic Cepheids is found to be shifted towards slightly older values, when the F-mode canonical relations are considered, with respect to the case at a fixed solar chemical composition. A marginally opposite dependence can be found in the noncanonical F-mode and canonical FO-mode cases.

We discuss stationary and axisymmetric trans-magnetosonic outflows in the magnetosphere of a rotating black hole (BH). Ejected plasma from the plasma source located near the BH is accelerated far away to form a relativistic jet. In this study, the plasma acceleration efficiency and conversion of fluid energy from electromagnetic energy are considered by employing the trans-fast magnetosonic flow solution derived by Takahashi & Tomimatsu (2008). Considering the parameter dependence of magnetohydrodynamical flows, we search for the parameters of the trans-magnetosonic outflow solution to the recent M87 jet observations and obtain the angular velocity values of the magnetic field line and angular momentum of the outflow in the magnetized jet flow. Therefore, we estimate the locations of the outer light surface, Alfv\'en surface, and separation surface of the flow. We also discuss the electromagnetic energy flux from the rotating BH (i.e., the Blandford-Znajek process), which suggests that the energy extraction mechanism is effective for the M87 relativistic jet.

We present the public release of the Complete History of Interaction-Powered Supernovae (CHIPS) code, suited to model a variety of transients that arise from interaction with a dense circumstellar medium (CSM). Contrary to existing modellings which mostly attach the CSM by hand, CHIPS self-consistently simulates both the creation of the CSM from mass eruption of massive stars prior to core-collapse, and the subsequent supernova light curve. We demonstrate the performance of CHIPS by presenting examples of the density profiles of the CSM and the light curves. We show that the gross light curve properties of putative interaction-powered transients, such as Type IIn supernovae, rapidly evolving transients and recently discovered fast blue optical transients, can be comprehensively explained with the output of CHIPS.

The rate at which mass is lost during the Red Supergiant evolutionary stage may strongly influence how the star appears. Though there have been many studies discussing how RSGs appear in the mid and far infrared (IR) as a function of their mass-loss rate, to date there have been no such investigations at optical and near-IR wavelengths. In a preliminary study we construct model atmospheres for RSGs which include a wind, and use these models to compute synthetic spectra from the optical to the mid-infrared. The inclusion of a wind has two important effects. Firstly, higher mass-loss rates result in stronger absorption in the TiO bands, causing the star to appear as a later spectral type despite its effective temperature remaining constant. This explains the observed relation between spectral type, evolutionary stage and mid-IR excess, as well as the mismatch between temperatures derived from the optical and infrared. Secondly, the wind mimics many observed characteristics of a `MOLsphere', potentially providing an explanation for the extended molecular zone inferred to exist around nearby RSGs. Thirdly, we show that wind fluctuations can explain the spectral variability of Betelgeuse during its recent dimming, without the need for dust.

The atomic-to-molecular hydrogen (H/H2) transition has been extensively studied as it controls the fraction of gas in molecular state in an interstellar cloud. This fraction is linked to star-formation by the Schmidt-Kennicutt law. While theoretical estimates (Sternberg et al. 2014) of the column density of the H I layer have been proposed for static photodissociation regions (PDRs), Herschel and well-resolved ALMA observations have revealed dynamical effects in star forming regions, caused by the process of photo-evaporation. We extend the analytic study of the H/H2 transition to include the effects of the propagation of the ionization front, in particular in presence of photo-evaporation at the walls of blister H II regions, and find its consequences on the total atomic hydrogen column density at the surface of clouds in presence of a UV field, and on the properties of the H/H2 transition. We solve semi-analytically the differential equation giving the H2 column density profile taking into account H2 formation on grains, H2 photodissociation and the ionization front propagation dynamics modeled as an advection of the gas through the ionization front. Taking into account this advection reduces the width of the atomic region compared to static models. The atomic region may disappear if the ionization front velocity exceeds a certain value, leading the H/H2 transition and the ionization front to merge. For both dissociated and merged configurations, we provide analytical expressions to determine the total H I column density. Our results take into account the metallicity. Finally, we compared our results to observations of PDRs illuminated by O-stars, for which we concluded that the dynamical effects are strong, especially for low-excitation PDRs.

The KM3NeT research infrastructure is under construction in the Mediterranean Sea. KM3NeT will study atmospheric and astrophysical neutrinos with two multi-purpose neutrino detectors, ARCA and ORCA, primarily aimed at the GeV-PeV energy scale. Thanks to the multi-photomultiplier tube design of the digital optical modules, KM3NeT is capable of detecting the neutrino burst from a Galactic or near-Galactic core-collapse supernova. This potential is already exploitable with the first detection units deployed in the sea. This paper describes the real-time implementation of the supernova neutrino search, operating on the two KM3NeT detectors since the first months of 2019. A quasi-online astronomy analysis is introduced to study the time profile of the detected neutrinos for especially significant events. The mechanism of generation and distribution of alerts, as well as the integration into the SNEWS and SNEWS 2.0 global alert systems are described. The approach for the follow-up of external alerts with a search for a neutrino excess in the archival data is defined. Finally, an overview of the current detector capabilities and a report after the first two years of operation are given.

We present new spectral and photometric data of confirmed LBV star from the NGC4736 galaxy. The star NGC4736_1 (Mbol = -11.5 mag) showed noticeable spectral variability from 2015 to 2018, which was accompanied by a significant change in the brightness. We also have estimated possible initial mass of the object NGC4736_1 as ~130 Msun.

We investigate the dissolution process of young embedded star clusters with different primordial mass segregation levels using fractal distributions by means of N-body simulations. We combine several star clusters in virial and subvirial global states with Plummer and uniform density profiles to mimic the gas. The star clusters have masses of Mstars = 500 Mo which follow an initial mass function where the stars have maximum distances from the centre of r = 1.5 pc. The clusters are placed in clouds which at the same radius have masses of Mcloud = 2000 Mo, resulting in star formation efficiency of 0.2. We remove the background potential instantaneously at a very early phase, mimicking the most destructive scenario of gas expulsion. The evolution of the fraction of bound stellar mass is followed for a total of 16 Myr for simulations with stellar evolution and without. We compare our results with previous works using equal-mass particles where an analytical physical model was used to estimate the bound mass fraction after gas expulsion. We find that independent of the initial condition, the fraction of bound stellar mass can be well predicted just right after the gas expulsion, but tends to be lower at later stages, as these systems evolve due to the stronger two-body interactions resulting from the inclusion of a realistic initial mass function. This discrepancy is independent of the primordial mass segregation level.

Context. Colliding-wind binaries are massive stellar systems featuring strong, interacting winds. These binaries may be actual particle accelerators, making them variable gamma-ray sources due to changes in the wind collision region along the orbit. However, only two of these massive stellar binary systems have been identified as high-energy sources. The first and archetypical system of this class is Eta Carinae, a bright gamma-ray source with orbital variability peaking around its periastron passage. Aims. The origin of the high energy emission in Eta Carinae is still unclear, with both lepto-hadronic and hadronic scenarios being under discussion. Moreover, the gamma-ray emission seemed to differ between the two periastrons previously observed with the Fermi Large Area Telescope. Continuing observations might provide highly valuable information for the understanding of the emission mechanisms in this system. Methods. We have used almost 12 years of data from the Fermi Large Area Telescope. We studied both low and high energy components, searching for differences and similarities between both orbits, and made use of this large dataset to search for emission from nearby colliding-wind binaries. Results. We show how the energy component above 10 GeV of Eta Carinae peaks months before the 2014 periastron, while the 2020 periastron is the brightest to date. Additionally, upper limits are provided for the high-energy emission in other particle-accelerating colliding-wind systems. Conclusions. Current gamma-ray observations of Eta Carinae strongly suggest that the wind collision region of this system is perturbed from orbit to orbit affecting particle transport within the shock.

We present a new estimate for the binary fraction (the fraction of stars with a single companion) for M dwarfs using a log-normal fit to the orbital separation distribution. We use point estimates of the binary fraction (binary fractions over specific separation and companion mass ratio ranges) from four M dwarf surveys sampling distinct orbital radii to fit a log-normal function to the orbital separation distribution. This model, alongside the companion mass ratio distribution given by Reggiani & Meyer (2013), is used to calculate the frequency of companions over the ranges of mass ratio (q) and orbital separation (a) over which the referenced surveys were collectively sensitive - [0.60 $\leq$ q $\leq$ 1.00] and [0.00 $\leq$ a $\leq$ 10,000 AU]. This method was then extrapolated to calculate a binary fraction which encompasses the broader ranges of [0.10 $\leq$ q $\leq$ 1.00] and [0.00 $\leq$ a < $\infty$ AU]. Finally, the results of these calculations were compared to the binary fractions of other spectral types. The binary fraction over the constrained regions of [0.60 $\leq$ q $\leq$ 1.00] and [0.00 $\leq$ a $\leq$ 10,000 AU] was calculated to be $0.229 \pm 0.028$. This quantity was then extrapolated over the broader ranges of q (0.10 - 1.00) and a (0.00 - $\infty$ AU) and found to be $0.462^{+0.057}_{-0.052}$. We used a conversion factor to estimate the multiplicity fraction from the binary fraction and found the multiplicity fraction over the narrow region of [0.60 $\leq$ q $\leq$ 1.00] and [0.00 $\leq$ a $\leq$ 10,000 AU] to be $0.270 \pm 0.111$. Lastly, we estimate the multiplicity fractions of FGK, and A stars using the same method (taken over [0.60 $\leq$ q $\leq$ 1.00] and [0.00 $\leq$ a $\leq$ 10,000 AU]) and find that the multiplicity fractions of M, FGK, and A stars, when considered over common ranges of q and a, are more similar than generally assumed.

We introduce (H)DPGMM, a hierarchical Bayesian non-parametric method based on the Dirichlet Process Gaussian Mixture Model, designed to infer data-driven population properties of astrophysical objects without being committal to any specific physical model. We investigate the efficacy of our model on simulated datasets and demonstrate its capability to reconstruct correctly a variety of population models without the need of fine-tuning of the algorithm. We apply our method to the problem of inferring the black hole mass function given a set of gravitational wave observations from LIGO and Virgo, and find that the (H)DPGMM infers a binary black hole mass function that is consistent with previous estimates without the requirement of a theoretically motivated parametric model. Although the number of systems observed is still too small for a robust inference, (H)DPGMM confirms the presence of at least two distinct modes in the observed merging black holes mass function, hence suggesting in a model-independent fashion the presence of at least two classes of binary black hole systems.

Aims. This series of papers aims at building a new formalism specifically tailored to study the impact of turbulence on the global modes of oscillation in solar-like stars. This first paper aims at deriving a linear wave equation that directly and consistently contains the turbulence as an input to the model, and therefore naturally contains the information on the coupling between the turbulence and the modes, through the stochasticity of the equations. Methods. We use a Lagrangian stochastic model of turbulence based on Probability Density Function methods to describe the evolution of the properties of individual fluid particles through stochastic differential equations. We then transcribe these stochastic differential equations from a Lagrangian frame to an Eulerian frame, more adapted to the analysis of stellar oscillations. We combine this method with Smoothed Particle Hydrodynamics, where all the mean fields appearing in the Lagrangian stochastic model are estimated directly from the set of fluid particles themselves, through the use of a weighting kernel function allowing to filter the particles present in any given vicinity. The resulting stochastic differential equations on Eulerian variables are then linearised. Results. We obtain a stochastic, linear wave equation governing the time evolution of the relevant wave variables, while at the same time containing the effect of turbulence. The wave equation generalises the classical, unperturbed propagation of acoustic waves in a stratified medium to a form that, by construction, accounts for the impact of turbulence on the mode in a consistent way. The effect of turbulence consists in a non-homogeneous forcing term, responsible for the stochastic driving of the mode, and a stochastic perturbation to the homogeneous part of the wave equation, responsible for both the damping of the mode and the modal surface effects.

To complement our previous discovery of the young snake-like structure in the solar neighborhood and reveal the structure's full extent, we explore {\tt Gaia EDR3} data within and surround the "snake" territory. With the friends-of-friends algorithm, we identify 2694 and 9615 Snake member candidates from the two samples. Thirteen open clusters are embedded in these member candidates. By combining the spectroscopic data from multiple surveys, we investigate the comprehensive properties of the candidates and find that they actually belong to one sizable structure, since most of the components are well bridged in their spatial distributions, and follow a single stellar population with an age of $30-40$\,Myr and with solar metallicity. This sizable structure should be a hierarchically primordial structure, and probably formed from a filamentary giant molecular cloud with different formation history in the local regions. The whole structure is expanding. To further analyze the dynamics of the Snake, we divide the structure into five groups according to their tangential kinematics. We find that the groups are expanding at a coherent rate ($\kappa_X\sim3.0\,\times10^{-2}\,\rm km\,s^{-1}\,pc^ {-1}$) in along the length of the structure ($X$ direction). With over ten thousand member stars, the Snake is an ideal laboratory to study nearby coeval stellar formation, stellar physics, and environmental evolution over a large spatial extent.

It has been recently recognized that the observational relativistic effects, mainly arising from the light propagation in an inhomogeneous universe, induce the dipole asymmetry in the cross-correlation function of galaxies. In particular, the dipole asymmetry at small scales is shown to be dominated by the gravitational redshift effects. In this paper, we exploit a simple analytical description for the dipole asymmetry in the cross-correlation function valid at quasi-linear regime. In contrast to the previous model, a new prescription involves only one dimensional integrals, providing a faster way to reproduce the results obtained by Saga et al. (2020). Using the analytical model, we discuss the detectability of the dipole signal induced by the gravitational redshift effect from upcoming galaxy surveys. The gravitational redshift effect at small scales enhances the signal-to-noise ratio (S/N) of the dipole, and in most of the cases considered, the S/N is found to reach a maximum at $z\approx0.5$. We show that current and future surveys such as DESI and SKA provide an idealistic data set, giving a large S/N of $10\sim 20$. Two potential systematics arising from off-centered galaxies are also discussed (transverse Doppler effect and diminution of the gravitational redshift effect), and their impacts are found to be mitigated by a partial cancellation between two competitive effects. Thus, the detection of the dipole signal at small scales is directly linked to the gravitational redshift effect, and should provide an alternative route to test gravity.

Most young neutron stars belonging to the class of Central Compact Objects in supernova remnants (CCOs) do not have known periodicities. We investigated seven such CCOs to understand the common reasons for the absence of detected pulsations. Making use of XMM-Newton, Chandra, and NICER observations, we perform a systematic timing and spectral analysis to derive updated sensitivity limits for both periodic signals and multi-temperature spectral components that could be associated with radiation from hotspots on the neutron star surface. Based on these limits, we then investigated for each target the allowed viewing geometry that could explain the lack of pulsations. We estimate it is unlikely ($< 10^{-6}$) to attribute that we do not see pulsations to an unfavorable viewing geometry for five considered sources. Alternatively, the carbon atmosphere model, which assumes homogeneous temperature distribution on the surface, describes the spectra equally well and provides a reasonable interpretation for the absence of detected periodicities within current limits. The unusual properties of CCOs with respect to other young neutron stars could suggest a different evolutionary path, as that proposed for sources experiencing episodes of significant fallback accretion after the supernova event.

The recently reported Type II Gamma-ray Burst (GRB) 200826A challenges the collapsar models by questioning how they can generate a genuinely short duration of the event. This paper proposes that the burst can originate from the collapse of a Thorne-Zytkow-like Object (TZlO). The TZlO consists of a central neutron star (NS) with a dense white dwarf (WD) material envelope and a disk, which are formed as the aftermath of a WD-NS coalescence. We found the collapse of such a TZlO can naturally explain the short duration of GRB 200826A. Furthermore, the collapse can produce a magnetar as the central object, which provides additional energy injection via magnetic dipole radiation to the ejected WD materials, causing a bump-like feature in the optical band and a shallow decay of the X-ray band. The disk wind shell induced by the TZlO at a large radius also interacts with the ejected materials, which explains the ``supernova bump" observed at $\sim$ 28 days.

We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (<z>=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallicity. We find stellar metallicities, which mostly reflect iron abundances, scaling as $(Z_{Fe,\ast}/Z_{Fe,\odot})=-(0.81\pm0.01)+(0.32+0.03)\log(M_\ast/10^{10}M_\odot)$ across the mass range of $10^9\lesssim M_\ast/M_\odot\lesssim10^{11}$, being $\approx6\times$ lower than seen locally at the same masses. The instantaneous oxygen-to-iron ratio ($\alpha$-enhancement) inferred using the gas-phase oxygen MZRs, is on average found to be [O/Fe]$\approx0.47$, being higher than the local [O/Fe]$\approx0$. The observed changes in [O/Fe] and [Fe/H] are reproduced in simple flow-through gas-regulator models with steady star-formation histories (SFHs) that follow the evolving main sequence. Our models show that the [O/Fe] is determined almost entirely by the instantaneous specific star formation rate alone while being independent of the SFHs, mass, and the gas-regulation characteristics of the systems. We find that the locations of $\sim10^{10}M_\odot$ galaxies at z~2 in the [O/Fe]--metallicity planes are in remarkable agreement with the sequence of low-metallicity thick-disk stars in our Galaxy. This manifests a beautiful concordance between the results of Galactic archaeology and observations of high-redshift Milky Way progenitors. However, there remains a question of how and when the old metal-rich, low-$\alpha$/Fe stars seen in the bulge had formed by z~2 because such a stellar population is not seen in our data and difficult to explain in the context of our models.

In this study, the transit observations of HAT-P-16b and TrES-3b exoplanets were carried out at \c{C}\"U UZAYMER Observatory with a 50 cm Ritchey Chretien type telescope. By analyzing the light curves of both exoplanets, system parameters have been obtained with acceptable models. Some of these parameters for HAT-P-16 systems were found to be: M$_{P}$=4.172$\pm$0.163 M$_{J}$, R$_{P}$=1.309$\pm$0.111 R$_{J}$, a/R$_{*}$=7.1922$\pm$0.0017 and b=0.1003$\pm$0.1533 which are consistent with the values given in \cite{2010ApJ...720.1118B}. Also, for the TrES-3 system, the same parameters were calculated as M$_{P}$=1.959$\pm$0.111 M$_{J}$, R$_{P}$=1.320$\pm$0.169 R$_{J}$, a/R$_{*}$=6.0656$\pm$ 0.4899 and b=0.7892$\pm$0.0700. These values are also consistent with the values given in \cite{2007ApJ...663L..37O}. No significant Transit Timing Variations (TTV) was found as a result of the Lomb-Scargle (LS) periodogram for HAT-P-16b by using the transit times obtained from the transit light curves (False Alarm Probabilities -- FAP) = \%96). However, a significant periodicity was found at 32.38 days and a secondary periodicity was found at 41.07 days were determined as a result of the LS periodogram obtained from the TrES-3b data. We interpret that these periods with FAP of 3.41 \% and 1.88 \%, respectively, are close to the rotation period of the host star. Therefore, these periods could be due to the modulation of spot-induced light variations.

Sterile neutrinos can affect the evolution of the universe, and thus using the cosmological observations can search for sterile neutrinos. In this work, we use the cosmic microwave background (CMB) anisotropy data from the Planck 2018 release, combined with the latest baryon acoustic oscillation (BAO), type Ia supernova (SN), and Hubble constant ($H_0$) data, to constrain the cosmological models with considering sterile neutrinos. In order to test the influences of the properties of dark energy on the constraint results of searching for sterile neutrinos, in addition to the $\Lambda$ cold dark matter ($\Lambda$CDM) model, we also consider the $w$CDM model and the holographic dark energy (HDE) model. We find that sterile neutrinos cannot be detected when the $H_0$ local measurement is not included in the data combination. When the $H_0$ measurement is included in the joint constraints, it is found that $\Delta N_{\rm eff}>0$ is detected at about 2.7$\sigma$ level for the $\Lambda$CDM model and at about 1--1.7$\sigma$ level for the $w$CDM model. However, $m_{\nu,{\rm{sterile}}}^{\rm{eff}}$ still cannot be well constrained and only upper limits can be given. In addition, we find that the HDE model is definitely ruled out by the current data. We also discuss the issue of the Hubble tension, and we conclude that involving sterile neutrinos in the cosmological models cannot truly resolve the Hubble tension.

We present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos -- such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes, and show how it simplifies in the high energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schr\"odinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.

HD 141569 is a Herbig Ae/Be star that straddles the boundary between the transition disks and debris disks. It is a low dust mass disk that reveals numerous structural elements (e.g. gaps and rings) that may point to young planets. It also exhibits a reservoir of CO gas observed at both millimeter and IR wavelengths. Previous observations (Goto et al. 2006) reported a possible asymmetry in the CO gas emission. Herein the IR ro-vibrational emission lines are analyzed and modeled both spectroscopically and spectroastrometrically. We find emission features from both 12CO and 13CO isotopologues heated to a temperature of approximately 200 K in the radial extent of 13 to 60 au. We do not see evidenceof the previously reported asymmetry in CO emission, our results being consistent with a Keplerian, axisymmetric emitting region. This raises the question of whether the emission profile may be evolving in time, possibly as a result of an orbiting feature in the inner disk such as a planet.

In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) down to the faintest possible limits, spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of the GAMA-KiDS photometric catalogue: $m_r \lesssim 23$; more than a decade beyond the limits of the GAMA spectroscopic redshift sample via clustering-based redshift inference. This technique uses spatial cross-correlation statistics for a reference set with known redshifts (in our case, the main GAMA sample) to derive the redshift distribution for the target ensemble. For the calibration of the redshift distribution we use a simple parametrisation with an adaptive normalisation factor over the interval $0.005 < z < 0.48$ to derive the clustering redshift results. We find that the GLF has a relatively constant power-law slope $\alpha \approx -1.2$ for $-17 \lesssim M_r \lesssim -13$, and then appears to steepen sharply for $-13 \lesssim M_r \lesssim -10$. This upturn appears to be where Globular Clusters (GCs) take over to dominate the source counts as a function of luminosity. Thus we have mapped the GLF across the full range of the $z \sim 0$ field galaxy population from the most luminous galaxies down to the GC scale.

We introduce a semi-parametric model for the primary mass distribution of binary black holes (BBHs) observed with gravitational waves (GWs) that applies a cubic-spline perturbation to a power law. We apply this model to the 46 BBHs included in the second gravitational wave transient catalog (GWTC-2). The spline perturbation model recovers a consistent primary mass distribution with previous results, corroborating the existence of a peak at $35\,M_\odot$ ($>97\%$ credibility) found with the \textsc{Powerlaw+Peak} model. The peak could be the result pulsational pair-instability supernovae (PPISNe). The spline perturbation model finds potential signs of additional features in the primary mass distribution at lower masses similar to those previously reported by Tiwari and Fairhurst (2021). However, with fluctuations due to small number statistics, the simpler \textsc{Powerlaw+Peak} and \textsc{BrokenPowerlaw} models are both still perfectly consistent with observations. Our semi-parametric approach serves as a way to bridge the gap between parametric and non-parametric models to more accurately measure the BBH mass distribution. With larger catalogs we will be able to use this model to resolve possible additional features that could be used to perform cosmological measurements, and will build on our understanding of BBH formation, stellar evolution and nuclear astrophysics.

Adaptive optics systems are an essential technology for the modern astronomy for ground based telescopes. One of the most recent revolution in the field is the introduction of the pyramid wavefront sensor. The higher performance of this device is payed with increased complexity in the control. In this work we report about advances in the AO system control obtained with SOULat the Large Binocular Telescope. The first is an improved Tip/Tilt temporal control able to recover the nominal correction even in presence of high temporal frequency resonances. The second one is a modal gain optimization that has been successfully tested on sky for the first time. Pyramid wavefront sensors are the key technology for the first light AO systems of all ELTs and the reported advances can be relevant contributions for such systems.

The evolution of White Dwarfs (WDs) depends crucially on thermal processes. The plasma in their core can produce neutrinos which escape from the star, thus contributing to the energy loss. While in absence of a magnetic field the main cooling mechanism is plasmon decay at high temperature and photon surface emission at low temperature, a large magnetic field in the core hiding beneath the surface even of ordinary WDs, and undetectable to spectropolarimetric measurements, can potentially leave an imprint in the cooling. In this paper, we revisit the contribution to WD cooling stemming from neutrino pair synchrotron radiation and the effects of the magnetic field on plasmon decay. Our key finding is that even if observations limit the magnetic field strength at the stellar surface, strong magnetic fields in the interior of WDs -- with or without a surface magnetic field -- can be high enough to modify the cooling rate, which is sensitive to the magnetic field value due to neutrino pair synchrotron emission.

The significant neutrino flux at high rapidity at the LHC motivates dedicated forward detectors to study the properties of neutrinos at TeV energies. We investigate magnetic dipole interactions between the active neutrinos and new sterile states at emulsion and liquid argon experiments that could be located in a future Forward Physics Facility (FPF) downstream of the ATLAS interaction point. The up-scattering of neutrinos off electrons produces an electron recoil signature that can probe new regions of parameter space at the High Luminosity LHC (HL-LHC), particularly for liquid argon detectors due to low momentum thresholds. We also consider the decay of the sterile neutrino through the dipole operator, which leads to a photon that could be displaced from the production vertex. FPF detectors can test sterile neutrino states as heavy as 1 GeV produced through the dipole portal, highlighting the use of high energy LHC neutrinos as probes of new physics.

The problem of the gravitational radiation damping of neutron star fundamental ($f$) mode oscillations has received considerable attention. Many studies have looked at the stability of such oscillations in rapidly rotating stars, calculating the growth/decay rate of the mode amplitude. In this paper, we look at the relatively neglected problem of the radiation reaction on the spin of the star. We specialise greatly to the so-called Kelvin modes: the modes of oscillation of (initially) non-rotating incompressible stars. We find the unexpected result that the excitation of a mode of angular momentum $\delta J$ on an initially non-rotating star ends up radiating an angular momentum $2 \delta J$ to infinity, leaving the star itself with a bulk angular momentum of $-\delta J$. This result is interesting in itself, and also will have implications for the angular momentum budgets of spinning down neutron stars, should such modes be excited.

We improve the DHOST Genesis proposed in \cite{Ilyas:2020zcb}, such that the near scale invariant scalar power spectrum can be generated from the model itself, without involking extra mechanism like a string gas. Besides, the superluminality problem of scalar perturbation plagued in \cite{Ilyas:2020zcb} can be rescued by choosing proper DHOST action.

We probe the cosmological consequences of a recently proposed class of solutions to the cosmological constant problem. In these models, the universe undergoes a long period of inflation followed by a contraction and a bounce that sets the stage for the hot big bang era. A requirement of any successful early universe model is that it must reproduce the observed scale-invariant density perturbations at CMB scales. While these class of models involve a long period of inflation, the inflationary Hubble scale during their observationally relevant stages is at or below the current Hubble scale, rendering the de Sitter fluctuations too weak to seed the CMB anisotropies. We show that sufficiently strong perturbations can still be sourced thermally if the relaxion field serving as the inflaton interacts with a thermal bath, which can be generated and maintained by the same interaction. We present a simple model where the relaxion field is derivatively (i.e. technically naturally) coupled to a non-abelian gauge sector, which gets excited tachyonically and subsequently thermalizes due to its nonlinear self-interactions. This model explains both the smallness of the cosmological constant and the amplitude of CMB anisotropies.

We study an Abelian gauge extension of the standard model with fermion families having non-universal gauge charges. The gauge charges and scalar content are chosen in such an anomaly-free way that only the third generation fermions receive Dirac masses via renormalisable couplings with the Higgs boson. Incorporating additional vector like fermions and scalars with appropriate $U(1)$ charges can lead to radiative Dirac masses of first two generations with neutral fermions going in the loop being dark matter candidates. Focusing on radiative muon mass, we constrain the model from the requirement of satisfying muon mass, recently measured muon anomalous magnetic moment by the E989 experiment at Fermilab along with other experimental bounds including the large hadron collider (LHC) limits. The anomalous Higgs coupling to muon is constrained from the LHC measurements of Higgs to dimuon decay. The singlet fermion dark matter phenomenology is discussed showing the importance of both annihilation and coannihilation effects. Incorporating all bounds lead to a constrained parameter space which can be probed at different experiments.

We study isotropic and slowly-rotating stars made of dark energy adopting the extended Chaplygin equation-of-state. We compute the moment of inertia as a function of the mass of the stars, both for rotating and non-rotating objects. The solution for the non-diagonal metric component as a function of the radial coordinate for three different star masses is shown as well. We find that i) the moment of inertia increases with the mass of the star, ii) in the case of non-rotating objects the moment of inertia grows faster, and iii) the curve corresponding to rotation lies below the one corresponding to non-rotating stars.

The flow of an electrically conducting fluid in a thin disc under the action of an azimuthal Lorentz force is studied experimentally. At small forcing, the Lorentz force is balanced by either viscosity or inertia, yielding quasi-Keplerian velocity profiles. For very large current and moderate magnetic field, we observe a new regime, fully turbulent, which exhibits large fluctuations and a Keplerian mean rotation profile $\Omega\sim \frac{\sqrt{IB}}{r^{3/2}}$. In this turbulent regime, the dynamics is typical of thin layer turbulence, characterized by a direct cascade of energy towards the small scales and an inverse cascade to large scale. Finally, at very large magnetic field, this turbulent flow bifurcates to a quasi-bidimensional turbulent flow involving the formation of a large scale condensate in the horizontal plane. These results are well understood as resulting from an instability of the B\"odewadt-Hartmann layers at large Reynolds number and discussed in the framework of similar astrophysical flows.

These lecture notes aim to provide an introduction to dark matter from the perspective of astrophysics/cosmology. We start with a rapid overview of cosmology, including the evolution of the Universe, its thermal history and structure formation. Then we look at the observational evidence for dark matter, from observations of galaxies, galaxy clusters, the anisotropies in the cosmic microwave background radiation and large scale structure. To detect dark matter we need to know how it's distributed, in particular in the Milky Way, so next we overview relevant results from numerical simulations and observations. Finally, we conclude by looking at what astrophysical and cosmological observations can tell us about the nature of dark matter, focusing on two particular cases: warm and self-interacting dark matter.

A possible detection of sub-solar mass ultra-compact objects would lead to new perspectives on the existence of black holes that are not of astrophysical origin and/or pertain to formation scenarios of exotic ultra-compact objects. Both possibilities open new perspectives for better understanding of our universe. In this work, we investigate the significance of detection of sub-solar mass binaries with components mass in the range: $10^{-2} M_\odot$ up to 1$M_\odot$, within the expected sensitivity of the ground-based gravitational waves detectors of third-generation, viz., the Einstein Telescope (ET) and the Cosmic Explorer (CE). Assuming a minimum of amplitude signal-to-noise ratio for detection, viz., $\rho = 8$, we find that the maximum horizon distances for an ultra-compact binary system with components mass $10^{-2} \, M_\odot$ and 1$M_\odot$ are 40 Mpc and 1.89 Gpc, respectively, for ET, and 125 Mpc and 5.8 Gpc, respectively, for CE. Other cases are also presented in the text. We derive the merger rate, and discuss consequences on the abundances of primordial black hole (PBH), $f_{\rm PBH}$. Considering the entire mass range [$10^{-2}$ - 1]$M_\odot$, we find $f_{\rm PBH} < 0.70$ ($<$ $0.06$) for ET (CE), respectively.

Boosted dark matter can be produced by some mechanism and detecting it is a key to understand the nature of dark matter. We show that the semi-annihilation $\chi\chi \to \bar{\chi}\nu$ indicates signals distinctive from the other semi-annihilation and standard dark matter annihilation processes. Since the boosted dark matter produced by this semi-annihilation is regarded as a high energy neutrino, the total flux of the dark matter and the accompanying neutrino yields double peaks at the energy close to the dark matter mass. Both of the particles can be detectable at large volume neutrino detectors.

We investigate the cosmological aspects of the most general parity preserving Metric-Affine Gravity theory quadratic in torsion and non-metricity in the presence of a cosmological hyperfluid. The equations of motion are obtained by varying the action with respect to the metric and the independent affine connection. Subsequently, considering a Friedmann-Lema\^itre-Robertson-Walker background, we derive the most general form of the modified Friedmann equations for the full quadratic theory. We then focus on a characteristic sub-case involving only two quadratic contributions given in terms of torsion and non-metricity vectors. In this setup, studying the modified Friedmann equations along with the conservation laws of the perfect cosmological hyperfluid, we provide exact solutions both for purely dilation and for purely spin hypermomentum sources. We then discuss the physical consequences of our model and the prominent role of torsion and non-metricity in this cosmological setup.