Papers for Tuesday, Jul 27 2021

The IceCube Upgrade is the first step towards the next-generation neutrino observatory at the South Pole, IceCube-Gen2, and will be installed in the central region of the existing array. The Upgrade will consist of 693 newly developed, densely spaced optical sensors and 50 standalone calibration devices, which will enhance IceCube's capabilities both at low and high neutrino energies. Of the new sensors, 402 will be multi-PMT Digital Optical Modules (mDOMs). Consisting of 24 small photomultipliers arranged inside a pressure vessel, the mDOM features a large sensitive area distributed nearly homogeneously over the full solid angle. The use of multiple, individually read-out PMTs allows directional information to be obtained for the registered photons and enables the use of multiplicity triggering within a single module, e.g., for background suppression. The challenges driving the mDOM development included tight restrictions on module size, data-transfer rate, and power consumption as well as the harsh environment in the deep ice at the South Pole. In this contribution we present the final mDOM design that meets these challenges.

The detection of a PeV high-energy neutrino of astrophysical origin, observed by the IceCube Collaboration and correlated with a 3$\sigma$ significance with Fermi measurements to the gamma-ray blazar TXS 0506+056, further stimulated the discussion on the production channels of high-energy particles in blazars. Many models also consider a hadronic component that would not only contribute to the emission of electromagnetic radiation in blazars but also lead to the production of secondary high-energy neutrinos and gamma-rays. Relativistic and compact plasma structures, so-called plasmoids, have been discussed in such flares to be moving along the jet axis. The frequently used assumption in such models that diffusive transport can describe particles in jet plasmoids is investigated in the present contribution. While the transport in the stationary scenario is diffusive for most of the parameter space, a flaring scenario is always accompanied by a non-diffusive phase in the beginning. In this paper, we present those conditions that determine the time scale to reach the diffusion phase as a function of the model parameters in the jet. We show that the type of the charged-particle transport, diffusive or ballistic, has a large influence on many observables, including the spectral energy distribution of blazars.

Planet formation around one component of a tight, eccentric binary system such as $\gamma$ Cephei (with semimajor axis around 20 AU) is theoretically challenging because of destructive high-velocity collisions between planetesimals. Despite this fragmentation barrier, planets are known to exist in such orbital configurations. Here we present a novel numerical framework for carrying out multi-annulus coagulation-fragmentation calculations of planetesimal growth, which fully accounts for the specifics of planetesimal dynamics in binaries, details of planetesimal collision outcomes, and the radial transport of solids in the disk due to the gas drag-driven inspiral. Our dynamical inputs properly incorporate the gravitational effects of both the eccentric stellar companion and the massive non-axisymmetric protoplanetary disk in which planetesimals reside, as well as gas drag. We identify a set of disk parameters that lead to successful planetesimal growth in systems such as $\gamma$ Cephei or $\alpha$ Centauri starting from $1-10$ km size objects. We identify the apsidal alignment of a protoplanetary disk with the binary orbit as one of the critical conditions for successful planetesimal growth: It naturally leads to the emergence of a dynamically quiet location in the disk (as long as the disk eccentricity is of order several percent), where favorable conditions for planetesimal growth exist. Accounting for the gravitational effect of a protoplanetary disk plays a key role in arriving at this conclusion, in agreement with our previous results. These findings lend support to the streaming instability as the mechanism of planetesimal formation. They provide important insights for theories of planet formation around both binary and single stars, as well as for the hydrodynamic simulations of protoplanetary disks in binaries (for which we identify a set of key diagnostics to verify).

In 2019, Sgr A* -- the supermassive black hole in the Galactic Center -- underwent unprecedented flaring activity, brightening by up to a factor of 100 compared to quiescent values. Here we report ALMA observations of Sgr A*'s continuum variability at 1.3 mm (230 GHz) -- a tracer of the accretion rate -- conducted one month after the brightest detected near infrared (NIR) and in the middle of the flaring activity of 2019. We develop an innovative light curve extraction technique which (together with ALMA's excellent sensitivity) allows us to obtain the light curves which are simultaneously of high time resolution (2 seconds) and high signal-to-noise ratio (~ 500). We construct an accurate intrinsic structure function of the Sgr A* submm variability, improving on previous studies by about two orders of magnitude in timescale and one order of magnitude in sensitivity. We compare the June 2019 variability behavior with that of 2001-2017, and suggest that the most likely cause of the bright NIR flares is magnetic reconnection.

We present the orbital solution of a peculiar double-lined spectroscopic and eclipsing binary system, 2M17091769+3127589. This solution was obtained by a simultaneous fit of both APOGEE radial velocities and TESS and ASAS-SN light curves to determine masses and radii. This system consists of an $M=0.256^{+0.010}_{-0.006}$ $M_\odot$, $R=3.961^{+0.049}_{-0.032}$ $R_{\odot}$ red giant and a hotter $M=1.518 ^{+0.057}_{-0.031}$ $M_\odot$, $R=2.608^{+0.034}_{-0.321}$ $R_{\odot}$ subgiant. Modelling with the MESA evolutionary codes indicates that the system likely formed 5.26 Gyrs ago, with a $M=1.2$ $M_\odot$ primary that is now the system's red giant and a $M=1.11$ $M_\odot$ secondary that is now a more massive subgiant. Due to Roche-lobe overflow as the primary ascends the red giant branch, the more evolved "primary" (i.e., originally the more massive star of the pair) is now only one-sixth as massive as the "secondary". Such a difference between the initial and the current mass ratio is one of the most extreme detected so far. Evolutionary modelling suggests the system is still engaged in mass transfer, at a rate of $\dot{M} \sim 10^{-9}$ $M_\odot$ yr$^{-1}$, and it provides an example of a less evolved precursor to some of the systems that consist of white dwarfs and blue stragglers.

We present the results of a study aimed at exploring, by means of N-body simulations, the evolution of rotating multi-mass star clusters during the violent relaxation phase, in the presence of a weak external tidal field. We study the implications of the initial rotation and the presence of a mass spectrum for the violent relaxation dynamics and the final properties of the equilibria emerging at the end of this stage. Our simulations show a clear manifestation of the evolution towards spatial mass segregation and evolution towards energy equipartition during and at the end of the violent relaxation phase. We study the final rotational kinematics and show that massive stars tend to rotate more rapidly than low-mass stars around the axis of cluster rotation. Our analysis also reveals that during the violent relaxation phase, massive stars tend to preferentially segregate into orbits with angular momentum aligned with the cluster's angular momentum, an effect previously found in the context of the long-term evolution of star clusters driven by two-body relaxation.

The increasingly large amount of cosmological data coming from ground-based and space-borne telescopes requires highly efficient and fast enough data analysis techniques to maximise the scientific exploitation. In this work, we explore the capabilities of supervised machine learning algorithms to learn the properties of the large-scale structure of the Universe, aiming at constraining the matter density parameter, Omega m. We implement a new Artificial Neural Network for a regression data analysis, and train it on a large set of galaxy two-point correlation functions in standard cosmologies with different values of Omega m. The training set is constructed from log-normal mock catalogues which reproduce the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) galaxies. The presented statistical method requires no specific analytical model to construct the likelihood function, and runs with negligible computational cost, after training. We test this new Artificial Neural Network on real BOSS data, finding Omega m=0.309p/m0.008, which is remarkably consistent with standard analysis results.

We have performed a coordinated focused ion beam (FIB)-scanning and transmission electron microscopy (S/TEM), electron probe microanalysis (EMPA)-synchrotron X-ray fluorescence (SXRF) microprobe study to determine phase-specific microstructural characteristics and high-resolution in situ trace element concentrations of primary pyrrhotite, pentlandite, and associated metal grains from chondrules in CM2 and CR2 carbonaceous chondrites. This work is the first of its kind to link trace element chemical and microstructural observations in chondritic sulfides in an attempt to determine formation mechanisms and conditions of primary sulfides in these meteorite groups. SXRF microprobe analyses allowed the concentrations of the minor and trace elements, Co, Cu, Ge, Zn, and Se to be quantified, in addition to Fe and Ni, at a spatial resolution of 2 microns. The similarity between the CM and CR PPI sulfide trace element patterns provides evidence for a common formation mechanism for this type of sulfide grain in both meteorite groups. In addition, the SRM sulfide and metal have comparable trace element patterns that indicates a genetic relationship between the two, such as sulfidization of metal. Enrichments in Ni, Co, Cu, and Se are consistent with the chalcophile/siderophile behavior of these elements. The observed depletions in Ge suggest that it may have been lost by evaporation or else was never incorporated into the metal or sulfide precursor materials. The depletion in Zn may also be attributable to evaporation, but, being partially lithophile, may also have been preferentially incorporated into silicates during chondrule formation. Trace element concentrations support crystallization from an immiscible sulfide melt in chondrules for formation of the PPI grains and sulfidization of metal for the origin of the SRM grains.

The Draconid meteor shower shows strong bursts of activity at irregular intervals, with nearly no activity in intervening years. Five outbursts of the Draconid meteor shower were observed with specular meteor radars in Canada and Europe between 1999 and 2018. The outbursts generally lasted between 6 and 8 hours, and most were not fully visible at a single geographical site, emphasizing the need for observations at multiple longitudes for short-duration shower outbursts. There is at least a factor of two difference in the peak flux as measured on different radars; the initial trail radius effect is undercorrected for Draconid meteors, which are known to be fragile.

The cyano radical (CN), one of the first detected interstellar molecular species, is a key molecule in many astrochemical chains. Particularly, it is detected towards molecular cores, the birth places of the stars, and it is known that it is involved in the rich chemistry that takes place in these sites. At present there are not so many studies about the emission of this molecular species at small spatial scales towards massive young stellar objects. Thus, we present a high-angular resolution CN study towards a sample of massive protostars, with the aim of unveiling the spatial distribution at the small scale of the emission of this radical in relation to the star-forming processes. The interstellar CN has a strong emission line at the rest frequency 226874.764 MHz, thus, we search for observing projects in the ALMA database regarding high-mass star-forming regions observed at Band 6. A sample of ten high-mass star-forming regions were selected in base on that they present a clear emission of CN at the mentioned frequency. We found that the CN traces both molecular condensations and diffuse and extended gas surrounding them. In general, the molecular condensations traced by the maximums of the CN emission do not spatially coincide with the peaks of the continuum emission at 1.3 mm, which trace the molecular cores where the massive stars born. Based on the presence or lack of near-IR emission associated with such cores, we suggest that our sample is composed by sources at different stages of evolution. The CN is present at both, suggesting that this radical may be ubiquitous along the different star formation stages, and hence it may be involved in different chemical reactions occurring along the time in the formation of the stars. Additionally, other complex molecules were detected towards the continuum peaks of some of the analyzed cores.

In this work we discuss ongoing development of a hybrid fiber/copper data and timing infrastructure for the future IceCube-Gen2 detector. The IceCube Neutrino Observatory is a kilometer-scale detector operating with 86 strings of modules. These modules communicate utilizing a custom protocol to mitigate the signaling challenges of long distance copper cables. Moving past the limitations of a copper-based backbone will allow larger future IceCube detectors with extremely precise timing and a large margin of excess throughput to accommodate innovative future modules. To this end, the upcoming IceCube Upgrade offers an opportunity to deploy a pathfinder for the new fiber optic infrastructure, called the Fiber Test System. This design draws on experience from AMANDA and IceCube and incorporates recently matured technologies such as ruggedized fibers and White Rabbit timing to deliver robust and high-performance data and timing transfer.

On July 30th, 2019 IceCube detected a high-energy astrophysical muon neutrino candidate, IC-190730A, with a $67\%$ probability of astrophysical origin. The flat spectrum radio quasar (FSRQ) PKS 1502+106 is in the error circle of the neutrino. Motivated by this observation, we investigate whether the emission of IC-190730A from this source is plausible, considering the multi-wavelength (infrared/UV/optical/X-ray/gamma-ray) emission of PKS 1502+106 at the time of the neutrino arrival. We analyse UV/optical and X-ray data and collect additional observations from the literature to construct the multi-wavelength spectral energy distribution of PKS 1502+106. We perform leptohadronic modelling of the multi-wavelength emission of the source and determine the most plausible emission scenarios and the maximum expected accompanying neutrino flux. A model in which the multi-wavelength emission of PKS 1502+106 originates beyond the broad-line region and inside the dust torus is most consistent with the observations. In this scenario, PKS 1502+106 can have produced up to of order one muon neutrino with energy exceeding 100 TeV in the lifetime of IceCube. An appealing feature of this model is that the required proton luminosity is consistent with the average required proton luminosity if blazars power the observed ultra-high-energy-cosmic-ray flux and well below the source's Eddington luminosity. If such a model is ubiquitous among FSRQs, additional neutrinos can be expected from other bright sources with energy $\gtrsim 10$ PeV.

Measurements of neutrinos at and below 10 GeV provide unique constraints of neutrino oscillation parameters as well as probes of potential Non-Standard Interactions (NSI). The IceCube Neutrino Observatory's DeepCore array is designed to detect neutrinos down to GeV energies. IceCube has built the world's largest data set of neutrinos >10 GeV, making searches for NSI a computational challenge. This work describes the use of convolutional neural networks (CNNs) to improve the energy reconstruction resolution and speed of reconstructing O(10 GeV) neutrino events in IceCube. Compared to current likelihood-based methods which take seconds to minutes, the CNN is expected to provide approximately a factor of 2 improvement in energy resolution while reducing the reconstruction time per event to milliseconds, which is essential for processing large datasets.

The IceCube Observatory has collected over 577 billion cosmic-ray induced muon events in its final configuration from May 2011 to May 2020. We used this data set to provide an unprecedented statistically accurate map of the cosmic ray arrival direction distribution in the TeV-PeV energy range scale in the Southern Hemisphere. Such an increase in event statistics makes it possible to extend the sensitivity to anisotropies at higher cosmic ray energies and smaller angular scales. It will also facilitate a more detailed assessment of the observatory stability over both short- and long-time scales. This will enable us to study the time variability of the cosmic ray anisotropy on a yearly-base and over the entire data sample period covering most of the solar cycle 24. We present the preliminary results from the study with the extended event sample.

The solar gravitational lens (SGL) provides a factor of $10^{11}$ amplification for viewing distant point sources beyond our solar system. As such, it may be used for resolved imaging of extended sources, such as exoplanets, not possible otherwise. To use the SGL, a spacecraft carrying a modest telescope and a coronagraph must reach the SGLs focal region, that begins at $\sim$550 astronomical units (AU) from the Sun and is oriented outward along the line connecting the distant object and the Sun. No spacecraft has ever reached even a half of that distance; and to do so within a reasonable mission lifetime (e.g., less than 25 years) and affordable cost requires a new type of mission design, using solar sails and microsats ($<100$~kg). The payoff is high -- using the SGL is the only practical way we can ever get a high-resolution, multi-pixel image of an Earth-like exoplanet, one that we identify as potentially habitable. This paper describes a novel mission design starting with a rideshare launch from the Earth, spiraling in toward the Sun, and then flying around it to achieve solar system exit speeds of over $20$ AU/year. A new sailcraft design is used to make possible high area to mass ratio for the sailcraft. The mission design enables other fast solar system missions, starting with a proposed very low cost technology demonstration mission (TDM) to prove the functionality and operation of the microsat-solar sail design and then, building on the TDM, missions to explore distant regions of the solar system, and those to study Kuiper Belt objects (KBOs) and the recently discovered interstellar objects (ISOs) are also possible.

We study the initial conditions of a common envelope (CE) event resulting in a stellar merger. A merger's dynamics could be understood through its light curve, but no synthetic light curve has yet been created for the full evolution. Using the Smoothed Particle Hydrodynamics (SPH) code StarSmasher, we have created three-dimensional (3D) models of a 1.52 $M_\odot$ star that is a plausible donor in the V1309 Sco progenitor. The integrated total energy profiles of our 3D models match their initial one-dimensional (1D) models to within a 0.1 per cent difference in the top 0.1 $M_\odot$ of their envelopes. We have introduced a new method for obtaining radiative flux by linking intrinsically optically thick SPH particles to a single stellar envelope solution from a set of unique solutions. For the first time, we calculated our 3D models' effective temperatures to within a few per cent of the initial 1D models, and found a corresponding improvement in luminosity by a factor of $\gtrsim10^6$ compared to ray tracing. We let our highest resolution 3D model undergo Roche-lobe overflow with a 0.16 $M_\odot$ point-mass accretor ($P\simeq1.6$ days) and found a bolometric magnitude variability amplitude of $\sim0.3$ -- comparable to that of the V1309 Sco progenitor. Our 3D models are, in the top 0.1 $M_\odot$ of the envelope and in terms of total energy, the most accurate models so far of the V1309 Sco donor star. A dynamical simulation that uses the initial conditions we presented in this paper can be used to create the first ever synthetic CE evolution light curve.

We present a broadband map of polarized diffuse emission at 167-198 MHz developed from data from the Murchison Widefield Array (MWA). The map is designed to improve visibility simulation and precision calibration for 21 cm Epoch of Reionization (EoR) experiments. It covers a large swath - 11,000 sq. deg. - of the Southern Hemisphere sky in all four Stokes parameters and captures emission on angular scales of 1 to 9 degrees. The band-averaged diffuse structure is predominantly unpolarized but has significant linearly polarized structure near RA = 0 h. We evaluate the accuracy of the map by combining it with the GLEAM catalog and simulating an observation from the MWA, demonstrating that the accuracy of the short baselines (6.1-50 wavelengths) now approaches the accuracy of the longer baselines typically used for EoR calibration. We discuss how to use the map for visibility simulation for a variety of interferometric arrays. The map has potential to improve calibration accuracy for experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the forthcoming Square Kilometre Array (SKA) as well as the MWA.

We identify Local Group (LG) analogs in the IllustrisTNG cosmological simulation, and use these to study two mass estimators for the LG: one based on the timing argument (TA) and one based on the virial theorem (VT). Including updated measurements of the Milky Way-M31 tangential velocity and the cosmological constant, we show that the TA mass estimator slightly overestimates the true median LG-mass, though the ratio of the TA to the true mass is consistent at the approximate 90\% c.l. These are in broad agreement with previous results using dark matter-only simulations. We show that the VT estimator better estimates the true LG-mass, though there is a larger scatter in the virial mass to true mass ratio relative to the corresponding ratio for the TA. We attribute the broader scatter in the VT estimator to several factors, including the predominantly radial orbits for LG satellite galaxies, which differs from the VT assumption of isotropic orbits. With the systematic uncertainties we derive, the updated measurements of the LG mass at 90\% c.l. are $4.75_{-2.41}^{+2.22} \times 10^{12}$ M$_\odot$ from the TA and $2.0_{-1.5}^{+2.1} \times 10^{12}$ M$_\odot$ from the VT.

While there is evidence for the existence of dark matter, its properties have yet to be discovered. Simultaneously, the nature of high-energy astrophysical neutrinos detected by IceCube remains unresolved. If dark matter and neutrinos are coupled to each other, they may exhibit a non-zero elastic scattering cross section. Such an interaction between an isotropic extragalactic neutrino flux and dark matter would be concentrated in the Galactic Centre, where the dark matter column density is greatest. This scattering would attenuate the flux of high-energy neutrinos, which could be observed in IceCube. Using the seven-year Medium Energy Starting Events sample, we perform an unbinned likelihood analysis, searching for a signal based on a possible DM-neutrino interaction scenario. We search for a suppression of the high-energy astrophysical neutrino flux in the direction of the Galactic Centre, and compare these constraints to complementary low-energy information from large scale structure surveys and the cosmic microwave background.

The observed dark matter abundance in the Universe can be explained with non-thermal, heavy dark matter models. In order for dark matter to still be present today, its lifetime has to far exceed the age of the Universe. In these scenarios, dark matter decay can produce highly energetic neutrinos, along with other Standard Model particles. To date, the IceCube Neutrino Observatory is the world's largest neutrino telescope, located at the geographic South Pole. In 2013, the IceCube collaboration reported the first observation of high-energy astrophysical neutrinos. Since then, IceCube has collected a large amount of astrophysical neutrino data with energies up to tens of PeV, allowing us to probe the heavy dark matter models using neutrinos. We search the IceCube data for neutrinos from decaying dark matter in galaxy clusters and galaxies. The targeted dark matter masses range from 10 TeV to 10 PeV. In this contribution, we present the method and sensitivities of the analysis.

How do we bridge the gap between the Galactic and the extragalactic? By focusing on the topic of stellar dynamics and stellar populations of the Milky Way and its siblings this virtual meeting aimed at connecting both fields that each bring unique perspectives to understanding how disk galaxies form and evolve. As this meeting took place during a global pandemic, we also give our perspective on the challenges and best practises for running a virtual meeting.

Cosmological data from the late and the recent universe shows a puzzling $\sim 4.5\sigma$ tension in the rate of universe expansion considering the $\Lambda$CDM model. Along with several scenarios which were proposed to resolve the tension, understanding how a model might bias the results is very important. In this paper, we investigate the three most known non-parametric methods namely, smoothing method, genetic algorithm and Gaussian process for modeling a data set and reveal some aspects of such scenarios. Considering these three methods, we investigate the recent Hubble parameters data to reconstruct the rate of universe expansion and Pantheon sample to reconstruct the luminosity distance. In contrast to the similar studies in the literature, the $\chi^2$ distribution has been used as a criterion to select curves consistent with the data. Finally, we compute the $H_0$ from reconstructed curves and compare the results.

During the X-ray bursts of GS 1826$-$24, "clocked burster", the nuclear reaction flow that surges through the rapid-proton capture process path has to pass through the NiCu cycles before reaching the ZnGa cycles that moderate the further extent of hydrogen burning in the region above germanium and selenium isotopes. The $^{57}$Cu(p,$\gamma$)$^{58}$Zn reaction located in the NiCu cycles plays an important role in influencing the burst light curves as found by Cyburt et al. (2016). We deduce the $^{57}$Cu(p,$\gamma$)$^{58}$Zn reaction rate based on the experimentally determined important nuclear structure information, isobaric-multiplet-mass equation, and large-scale shell model calculations. With the isobaric-multiplet-mass equation, we propose a possible order of $1^+_1$ and $2^+_3$ dominant resonance states that the $2^+_3$ resonance state is higher than the $1^+_1$ state, and estimate the resonance energy of $1^+_2$ contributing resonance state. The new rate is up to a factor of five lower than the Forstner et al. (2001) rate recommended by JINA REACLIB v2.2. Using the one-dimensional implicit hydrodynamic code, KEPLER, to model the thermonuclear X-ray bursts of GS 1826$-$24 clocked burster, we find that the new $^{57}$Cu(p,$\gamma$)$^{58}$Zn reaction redistributes the reaction flow in the NiCu cycles and reduces the production of $^{58}$Zn, whereas the $^{59}$Cu(p,$\alpha$)$^{56}$Ni and $^{59}$Cu(p,$\gamma$)$^{60}$Zn reactions suppress the influence of the $^{57}$Cu(p,$\gamma$)$^{58}$Zn reaction and strongly diminish the impact of nuclear reaction flow that by-passes the important $^{56}$Ni waiting point induced by the $^{55}$Ni(p,$\gamma$)$^{56}$Cu reaction on burst light curve. The influence of the newly deduced $^{56}$Ni(p,$\gamma$)$^{57}$Cu is also discussed.

The energy and waiting time distributions are important properties to understand the physical mechanism of repeating fast radio bursts (FRBs). Recently, Five-hundred-meter Aperture Spherical radio Telescope (FAST) detected the largest sample of FRB 121102, containing 1652 bursts. The energy distribution at high-energy range ($>10^{38}$ erg) can be fitted with a single power-law function with an index of $-1.86$. However, the distribution at low-energy range deviates from the power-law function. The energy distributions of high-energy bursts at different epochs are inconsistent. We find the power-law index of $-1.70$ for early bursts and $-2.60$ for later bursts. For bursts observed in a single day, a linear repetition pattern is found. We use the Weibull function to fit the waiting time distribution. The shape parameter $k = 0.72^{+0.01}_{-0.02}$ and the event rate $r = 734.47^{+29.04}_{-27.58}$ day$ ^{-1} $ are derived. If the waiting times with $\delta_t < 28$ s are excluded, the burst behavior can be described by a Poisson process. The best-fitting values of $k$ are slightly different for low-energy and high-energy bursts. The event rates change significantly across the observing time, while the shape parameters $k$ vary slightly in different days.

Kepler's supernova remnant (SNR) which is produced by the most recent naked-eye supernova in our Galaxy is one of the best studied SNRs, but its gamma-ray detection has eluded us so far. Observations with modern imaging atmospheric Cherenkov telescopes (IACT) have enlarged the knowledge about nearby SNRs with ages younger than 500 years by establishing Cassiopeia A and Tycho's SNRs as very high energy (VHE) gamma-ray sources and setting a lower limit on the distance to Kepler's SNR. This SNR is significantly more distant than the other two and expected to be one of the faintest gamma-ray sources within reach of the IACT arrays of this generation. We report strong evidence for a VHE signal from Kepler's SNR based on deep observations of the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours, including 122 hours accumulated in 2017-2020. We further discuss implications of this result for cosmic-ray acceleration in young SNRs.

Extensive Air Showers generated by very low energy cosmic ray particles, because of the very steep cosmic ray energy spectrum, dominate the secondary particle flux measured by single detectors and small shower arrays. Showers observed by a number of detectors placed close together forming a small local array connected in extensive networks can be used to search for potentially interesting spatial correlations between individual showers that may shed new light on the nature of ultra-high energy cosmic rays. Quantitative interpretation of showers recorded by small arrays requires a different methodology to that used for normal large EAS arrays operating in the 'knee' area and above. To draw physical conclusions from such events, a suitable simulation tool is needed to determine, within a reasonable time, e.g., the expected registration rate for each configuration. In this paper we propose a small EAS generator, semi-analytical method to integrate cosmic-ray spectra at energies of interest for arbitrary detector configurations. For an efficient and fast summation over the mass composition of the primary cosmic rays, we have analysed different versions of the superposition model of heavy cosmic ray nuclei interactions and finally found a reliable description for a realistic, genuine wounded nucleon superposition model, which we have used in our generator. As an example of its use, a comparison of the results obtained from our generator with measurements of the secondary particle density spectrum at sea level is shown.

We derived a bi-dimensional calibration between the emission line ratios R23=([O II]3726+3729+[O II]4959+5007)/Hb, P=[([O II]4959+5007)/Hb]/R23 and the oxygen abundance relative to hydrogen (O/H) in the gas phase of Seyferts 1 and 2 nuclei. In view of this, emission-line intensity ratios for a sample of objects taken from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) measured by the MPA/JHU group and direct estimates of O/H based on Te-method, adapted for AGNs, are considered. We find no variation of R23 observed along the radii of AGNs which shows that this line ratio is a good oxygen abundance (O/H) indicator for the class of objects considered in this work. The derived O/H = f(R23, P) relation produces O/H values similar to estimations via Te-method in a wide range of metallicities [8.0 < 12+log(O/H) < 9.2]. Conversely to star-forming regions in the high metallicity regime, R23 shows a positive correlation trend with O/H in AGNs. This indicates that the hardness of ionizing radiation is not affected by the metallicities in these objects or Narrow Line Regions (NLRs) are not significantly modified by changes in the Spectral Energy Distribution due to metallicity variations.

The fundamental nature and extent of the coronal line region (CLR), which may serve as a vital tracer for Active Galactic Nucleus (AGN) activity, remain unresolved. Previous studies suggest that the CLR is produced by AGN-driven outflows and occupies a distinct region between the broad line region and the narrow line region, which places it tens to hundreds of parsecs from the galactic center. Here, we investigate 10 coronal line (CL; ionization potential $\ge$ 100 eV) emitting galaxies from the SDSS-IV MaNGA catalog with emission from one or more CLs detected at $\ge$ $5{\sigma}$ above the continuum in at least 10 spaxels - the largest such MaNGA catalog. We find that the CLR is far more extended, reaching out to 1.3 - 23 kpc from the galactic center. We cross-match our sample of 10 CL galaxies with the largest existing MaNGA AGN catalog and identify 7 in it; two of the remaining three are galaxy mergers and the final one is an AGN candidate. Further, we measure the average CLR electron temperatures to range between 12,331 K - 22,530 K, slightly above the typical threshold for pure AGN photoionization ($\sim$ 20,000 K) and indicative of shocks (e.g., merger-induced or from supernova remnants) in the CLR. We reason that ionizing photons emitted by the central continuum source (i.e. AGN photoionization) primarily generate the CLs, and that energetic shocks are an additional ionization mechanism that likely produce the most extended CLRs we measure.

Optical and near-UV continuum emissions in flares contribute substantially to flare energy budget. Two mechanisms play an important role for continuum emission in flares: hydrogen recombination after sudden ionization at chromospheric layers and transportation of the energy radiatively from the chromosphere to lower layers in the atmosphere, the so called back-warming. The aim of the paper is to disentangle between these two mechanisms for the excess of Balmer continuum observed in a flare. Methods. We combine the observations of Balmer continuum obtained with IRIS (spectra and SJIs 2832 A) and hard X-ray (HXR) emission detected by FERMI Gamma Burst Monitor (GBM) during a mini flare. Calibrated Balmer continuum is compared to non-LTE radiative transfer flare models and radiated energy is estimated. Assuming thick target HXR emission, we calculate the energy of non-thermal electrons detected by FERMI GBM and compare it to the radiated energy. The favorable argument of a relationship between the Balmer continuum excess and the HXR emission is that there is a good time coincidence between both of them. In addition, the shape of the maximum brightness in the 2832 SJIs, which is mainly due to this Balmer continuum excess, is similar to the FERMI/GBM light curve. The electron-beam flux estimated from FERMI/GBM is consistent with the beam flux required in non-LTE radiative transfer models to get the excess of Balmer continuum emission observed in the IRIS spectra. The low energy input by non thermal electrons above 20 keV is sufficient to produce the enhancement of Balmer continuum emission. This could be explained by the topology of the reconnection site. The reconnection starts in a tiny bald patch region which is transformed dynamically in a X-point current sheet. The size of the interacting region would be under the spatial resolution of the instrument.

Cosmic strings may have formed in the early universe due to the Kibble mechanism. While string networks are usually modeled as being of Nambu-Goto type, this description is understood to be a convenient approximation, which neglects the typically expected presence of additional degrees of freedom on the string worldsheet. Previous simulations of cosmic strings in expanding universes have established beyond doubt the existence of a significant amount of short-wavelength propagation modes (commonly called wiggles) on the strings, and a wiggly string extension of the canonical velocity-dependent one-scale model has been recently developed. Here we improve the physical interpretation of this model, by studying the possible asymptotic scaling solutions of this model, and in particular how they are affected by the expansion of the universe and the available energy loss or transfer mechanisms -- e.g., the production of loops and wiggles. In addition to the Nambu-Goto solution, to which the wiggly model reduces in the appropriate limit, we find that there are also solutions where the amount of wiggliness can grow as the network evolves or, for specific expansion rates, become a constant. Our results show that full scaling of the network, including the wiggliness, is much more likely in the matter era than in the radiation era, which is in agreement with numerical simulation results.

We present redshifts for 2753 low-redshift galaxies between $0.03 \lesssim z_{\rm spec}\lesssim0.5$ with 18 $\leq$ $r$ $\leq$ 22 obtained with Hectospec at the Multi-Mirror Telescope (MMT). The observations targeted the XMM-LSS, ELAIS-N1 and DEEP2-3 fields, each of which covers $\sim$ 1 deg$^2$. These fields are also part of the recently completed CFHT Large Area U-band Deep Survey (CLAUDS) and on-going Hyper Suprime-Cam deep fields surveys. The efficiency of our technique for selecting low-redshift galaxies is confirmed by the redshift distribution of our sources. In addition to redshifts, these high S/N spectra are used to measure ages, metallicities, and nuclear activity levels. In combination with the photometric catalogue in $u$, $g$, $r$, $i$, $z$, $y$ down to 27 AB mag, we are able to study the galaxy population down to stellar masses of $\sim$ 10$^8 M_\odot$ . This paper presents the observational strategy, the reduction procedure and properties of the galaxy sample.

The Spitzer Extended Deep Survey (SEDS) as a deep and wide mid-infrared (MIR) survey project provides a sample of 500000+ sources spreading 1.46 square degree and a depth of 26 AB mag (3$\sigma$). Combining with the previous available data, we build a PSF-matched multi-wavelength photometry catalog from u band to 8$\mu$m. We fit the SEDS galaxies spectral energy distributions by the local galaxy templates. The results show that the SEDS galaxy can be fitted well, indicating the high redshift galaxy ($z \sim 1$) shares the same templates with the local galaxies. This study would facilitate the further study of the galaxy luminosity and high redshift mass function.

We investigate the effects of ram pressure on the molecular ISM in the disk of the Coma cluster galaxy NGC 4921, via high resolution CO observations. We present 6" resolution CARMA CO(1-0) observations of the full disk, and 0.4" resolution ALMA CO(2-1) observations of the leading quadrant, where ram pressure is strongest. We find evidence for compression of the dense interstellar medium (ISM) on the leading side, spatially correlated with intense star formation activity in this zone. We also detect molecular gas along kiloparsec-scale filaments of dust extending into the otherwise gas stripped zone of the galaxy, seen in HST images. We find the filaments are connected kinematically as well as spatially to the main gas ridge located downstream, consistent with cloud decoupling inhibited by magnetic binding, and inconsistent with a simulated filament formed via simple ablation. Furthermore, we find several clouds of molecular gas $\sim 1-3$ kpc beyond the main ring of CO that have velocities which are blueshifted by up to 50 km s$^{-1}$ with respect to the rotation curve of the galaxy. These are some of the only clouds we detect that do not have any visible dust extinction associated with them, suggesting that they are located behind the galaxy disk midplane and are falling back towards the galaxy. Simulations have long predicted that some gas removed from the galaxy disk will fall back during ram pressure stripping. This may be the first clear observational evidence of gas re-accretion in a ram pressure stripped galaxy.

Molecular clouds are complex magnetized structures, with variations over a broad range of length scales. Ionization in dense, shielded clumps and cores of molecular clouds is thought to be caused by charged cosmic rays (CRs). These CRs can also contribute to heating the gas deep within molecular clouds, and their effect can be substantial in environments where CRs are abundant. CRs propagate predominantly by diffusion in media with disordered magnetic fields. The complex magnetic structures in molecular clouds therefore determine the propagation and spatial distribution of CRs within them, and hence regulate their local ionization and heating patterns. Optical and near-infrared (NIR) polarization of starlight through molecular clouds is often used to trace magnetic fields. The coefficients of CR diffusion in magnetized molecular cloud complexes can be inferred from the observed fluctuations in these optical/NIR starlight polarisations. Here, we present calculations of the expected CR heating patterns in the star-forming filaments of IC 5146, determined from optical/NIR observations. Our calculations show that local conditions give rise to substantial variation in CR propagation. This affects the local CR heating power. Such effects are expected to be severe in star-forming galaxies rich in CRs. The molecular clouds in these galaxies could evolve differently to those in galaxies where CRs are less abundant.

It was recently proposed that the intensity ratios of several extreme ultraviolet spectral lines from the Fe X ion can be used to measure the solar coronal magnetic field based on the magnetic-field-inducedtransition (MIT) theory. To verify the suitability of this method, we performed forward modelingwith a three-dimensional radiation magnetohydrodynamic model of a solar active region. Intensities of several spectral lines from Fe X were synthesized from the model. Based on the MIT theory, intensity ratios of the MIT line Fe X 257 A to several other Fe X lines were used to derive the magnetic field strengths, which were then compared with the field strengths in the model. We also developed a new method to simultaneously estimate the coronal density and temperature from the Fe X 174/175 and 184/345 A line ratios. Using these estimates, we demonstrated that the MIT technique can provide reasonably accurate measurements of the coronal magnetic field in both on-disk and off-limb solar observations. Our investigation suggests that a spectrometer that can simultaneously observe the Fe X 174, 175, 184, 257, and 345 A lines and allow an accurate radiometric calibration for these lines is highly desired to achieve reliable measurements of the coronal magnetic field. We have also evaluatedthe impact of the uncertainty in the Fe X 3p4 3d 4D5/2 and 4D7/2 energy difference on the magnetic field measurements.

We intend to use the impact of microlensing on the Fe III emission line blend along with a measure of its gravitational redshift to estimate the mass of the quasar's central supermassive black hole (SMBH). We fit the Fe III feature in multiple spectroscopic observations between 2008 and 2016 of the gravitationally lensed quasar Q 0957+561 with relatively high signal-to-noise ratios (at the adequate wavelength). Based on the statistics of microlensing magnifications, we used a Bayesian method to derive the size of its emitting region. The Fe III spectral feature appears systematically redshifted in all epochs of observation by a value of 17 angstroms on average. We find clear differences in the shape of the Fe III line blend between images A and B. Measuring the strength of those magnitude differences, we conclude that this blend may arise from a region of half-light radius of 15 lt-days, which is in good agreement with the accretion disk dimensions for this system. We obtain a mass for the central SMBH of (1.5 +/- 0.5) x 10^9 solar masses, consistent within uncertainties with previous mass estimates based on the virial theorem. The relatively small uncertainties in the mass determination (< 35%) make this method a compelling alternative to other existing techniques (e.g., the virial plus reverberation mapping based size) for measuring black hole masses. Combining the Fe III redshift-based method with the virial, we estimate a virial factor in the 1.2 to 1.7 range for this system.

The IceCube Neutrino Observatory, located at the geographic South Pole, is the world's largest neutrino telescope, instrumenting 1 km$^3$ of Antarctic ice with 5160 photosensors to detect Cherenkov light. For the IceCube Upgrade, to be deployed during the 2022-23 polar field season, and the enlarged detector IceCube-Gen2 several new optical sensor designs are under development. One of these optical sensors, the Wavelength-shifting Optical Module (WOM), uses wavelength-shifting and light-guiding techniques to measure Cherenkov photons in the UV range from 250 nm to 380 nm. In order to understand the potential gains from this new technology, a measurement of the scattering and absorption lengths of UV light was performed in the SPICEcore borehole at the South Pole during the winter seasons of 2018/2019 and 2019/2020. For this purpose, a calibration device with a UV light source and a detector using the wavelength shifting technology was developed. We present the design of the developed calibration device, its performance during the measurement campaigns, and the comparison of data to a Monte Carlo simulation.

Extremely red quasars (ERQs) in the redshift range of $2.0 - 3.4$ have extreme red colours of $i-W3\ge4.6$. Core ERQs have strong CIV emission lines with rest equivalent width of $\ge100$A. Many core ERQs also have CIV line profiles with peculiar boxy shapes which distinguish them from normal blue quasars. We show, using a combination of kernel density estimation and local outlier factor analyses on a space of the $i-W3$ colour, CIV rest equivalent width and line kurtosis, that core ERQs likely represent a separate population rather than a smooth transition between normal blue quasars and the quasars in the tail of the colour-REW distribution. We apply our analyses to find new criteria for selecting ERQs in this 3D parameter space. Our final selection produces $133$ quasars, which are three times more likely to have a visually verified CIV broad absorption line feature than the previous core ERQ sample.

Two important parameters inferred from the gravitational wave signals of binaries of precessing black holes are the spin tilt angles, i.e., the angles at which the black holes' spin axes are inclined with respect to the binary's orbital angular momentum. The LIGO-Virgo parameter estimation analyses currently provide spin tilts at a fiducial reference frequency, often the lowest frequency used in the data analysis. However, the most astrophysically interesting quantities are the spin tilts when the binary was formed, which can be significantly different from those at the reference frequency for strongly precessing binaries. The spin tilts at formally infinite separation are a good approximation to the tilts at formation in many formation channels and can be computed efficiently for binary black holes using precession-averaged evolution. Here, we present a new code for computing the tilts at infinity that combines the precession-averaged evolution with orbit-averaged evolution at high frequencies and illustrate its application to GW190521 and other binary black hole detections from O3a. We have empirically determined the transition frequency between the orbit-averaged and precession-averaged evolution to produce tilts at infinity with a given accuracy. We also have regularized the precession-averaged equations in order to obtain good accuracy for the very close-to-equal-mass binary parameters encountered in practice. This additionally allows us to investigate the singular equal-mass limit of the precession-averaged expressions, where we find an approximate scaling of $1/(1 - q)$ with the mass ratio $q$.

Interpretations of the very rich second solar spectrum of the MgH molecule face serious problems owing to the complete lack of any information about rates of collisions between the MgH and hydrogen atoms. This work seeks to begin the process of filling this lacuna by providing, for the first time, quantum excitation, depolarization, and polarization transfer collisional rates of the MgH ground state $X^2\Sigma$. To achieve the goals of this work, potential energy surfaces are calculated and then are included in the Schr\"odinger equation to obtain the probabilities of collisions and, thus, all collisional rates. Our rates are obtained for temperatures ranging from $T \!\!=$2000 K to $T \!\!=$15,000 K. Sophisticated genetic programming methods are adopted in order to fit all depolarization rates with useful analytical functions of two variables: the total molecular angular momentum and temperatures. We study the solar implications of our results, and we find that the $X^2\Sigma$ state of MgH is partially depolarized by isotropic collisions with neutral hydrogen in its ground state $^2S$. Our findings show the limits of applicability of the widely used approximation in which the lower-level polarization is neglected.

The hyperfine transitions of the ground-rotational state of the hydroxyl radical (OH) have emerged as a versatile tracer of the diffuse molecular interstellar medium. We present a novel automated Gaussian decomposition algorithm designed specifically for the analysis of the paired on-source and off-source optical depth and emission spectra of these transitions. In contrast to existing automated Gaussian decomposition algorithms, AMOEBA (Automated MOlecular Excitation Bayesian line-fitting Algorithm) employs a Bayesian approach to model selection, fitting all 4 optical depth and 4 emission spectra simultaneously. AMOEBA assumes that a given spectral feature can be described by a single centroid velocity and full width at half-maximum, with peak values in the individual optical depth and emission spectra then described uniquely by the column density in each of the four levels of the ground-rotational state, thus naturally including the real physical constraints on these parameters. Additionally, the Bayesian approach includes informed priors on individual parameters which the user can modify to suit different data sets. Here we describe AMOEBA and evaluate its validity and reliability in identifying and fitting synthetic spectra with known parameters.

Over 1500 DBZ or DZ white dwarfs (WDs) have been observed so far, and polluted atmospheres with metal elements have been found among these WDs. The surface heavy element abundances of known DBZ or DZ WDs show an evolutionary sequence. By using Modules for Experiments in Stellar Evolution, we create DB WDs, and simulate the element diffusion due to high gravitational fields and the metal-rich material accretion coming from the planet disrupted by the WD. In our models, the input parameters ($\alpha_{\rm MLT}$, $\alpha_{\rm th}$ and $Z$) have very weak effect on DB WD structures including interior temperatures, chemical profiles and convective zones.The mass-accretion rate and the effective temperature of DB WDs determine the abundances of heavy elements. The evolutionary sequence of Ca element for about 1500 observed DB or DBZ WDs cannot be explained by the model with a constant mass-accretion rate, but is consistent well with the model in which the mass-accretion rate decreases by one power law when $T_{\rm eff}>10$ kK and slightly increases by another power law when $T_{\rm eff}<10$ kK. The observed DB WD evolutionary sequence of heavy element abundances originates from WD cooling and the change of mass-accretion rate.

We study the image formation near point singularities (swallowtail and umbilics) in the simulated strongly lensed images of Hubble Ultra Deep Field (HUDF) by the Hubble Frontier Fields (HFF) clusters. In this work, we only consider nearly half of the brightest (a total of 5271) sources in the HUDF region. For every HFF cluster, we constructed 11 realizations of strongly lensed HUDF with an arbitrary translation of the cluster centre within the central region of HUDF and an arbitrary rotation. In each of these realizations, we visually identify the characteristic/exotic image formation corresponding to the different point singularities. We find that our current results are consistent with our earlier results based on different approaches. We also study time delay in these exotic image formations and compare it with typical five-image geometries. We find that the typical time delay in exotic image formations is an order of magnitude smaller than the typical time delay in a generic five-image geometry.

Several mechanisms for the transformation of blue star-forming to red quiescent galaxies have been proposed, and the green valley (GV) galaxies amid them are widely accepted in a transitional phase. Thus, comparing the morphological and environmental differences of the GV galaxies with early-type disks (ETDs; bulge dominated and having a disk) and late-type disks (LTDs; disk dominated) is suitable for distinguishing the corresponding quenching mechanisms. A large population of massive ($M_* \geqslant 10^{10}M_\odot$) GV galaxies at $0.5 \leqslant z \leqslant 1.5$ in 3D-HST/CANDELS is selected using extinction-corrected $(U-V)_{\rm rest}$ color. After eliminating any possible active galactic nucleus candidates and considering the "mass-matching", we finally construct two comparable samples of GV galaxies with either 319 ETD or 319 LTD galaxies. Compared to the LTD galaxies, it is found that the ETD galaxies possess higher concentration index and lower specific star formation rate, whereas the environments surrounding them are not different. This may suggest that the morphological quenching may dominate the star formation activity of massive GV galaxies rather than the environmental quenching. To quantify the correlation between the galaxy morphology and the star formation activity, we define a dimensionless morphology quenching efficiency $Q_{\rm mor}$ and find that $Q_{\rm mor}$ is not sensitive to the stellar mass and redshift. When the difference between the average star formation rate of ETD and LTD galaxies is about 0.7 $M_\odot \rm \;yr^{-1}$, the probability of $Q_{\rm mor}\gtrsim 0.2$ is higher than 90\%, which implies that the degree of morphological quenching in GV galaxies might be described by $Q_{\rm mor}\gtrsim 0.2$.

We performed a detailed spectral and timing analysis of a Seyfert 1 galaxy Mrk 509 using data from the Neil Gehrels Swift observatory that spanned over ~13 years between 2006 and 2019. To study the variability properties from the optical/UV to X-ray emission, we used a total of 275 pointed observations in this work. The average spectrum over the entire duration exhibits a strong soft X-ray excess above the power-law continuum. The soft X-ray excess is well described by two thermal components with temperatures of kT_BB1 ~120 eV and kT_BB2 ~460 eV. The warm thermal component is likely due to the presence of an optically thick and warm Comptonizing plasma in the inner accretion disk. The fractional variability amplitude is found to be decreasing with increasing wavelength, i.e. from the soft X-ray to UV/optical emission. However, the hard X-ray (2-8 keV) emission shows very low variability. The strength of the correlation within the UV and the optical bands (0.95-0.99) is found to be stronger than the correlation between the UV/Optical and X-ray bands (0.40-0.53). These results clearly suggest that the emitting regions of the X-ray and UV/optical emission are likely distinct or partly interacting. Having removed the slow variations in the light curves, we find that the lag spectrum is well described by the 4/3 rule for the standard Shakura-Sunyaev accretion disk when we omit X-ray lags. All these results suggest that the real disk is complex, and the UV emission is likely reprocessed in the accretion disk to give X-ray and optical emission.

In this paper, we present the analysis of incoherent non-thermal radio emission from a sample of hot magnetic stars, ranging from early-B to early-A spectral type. Spanning a wide range of stellar parameters and wind properties, these stars display a commonality in their radio emission which presents new challenges to the wind scenario as originally conceived. It was thought that relativistic electrons, responsible for the radio emission, originate in current sheets formed where the wind opens the magnetic field lines. However, the true mass-loss rates from the cooler stars are too small to explain the observed non-thermal broadband radio spectra. Instead, we suggest the existence of a radiation belt located inside the inner-magnetosphere, similar to that of Jupiter. Such a structure explains the overall indifference of the broadband radio emissions on wind mass-loss rates. Further, correlating the radio luminosities from a larger sample of magnetic stars with their stellar parameters, the combined roles of rotation and magnetic properties have been empirically determined. Finally, our sample of early-type magnetic stars suggests a scaling relationship between the non-thermal radio luminosity and the electric voltage induced by the magnetosphere's co-rotation, which appears to hold for a broader range of stellar types with dipole-dominated magnetospheres (like the cases of the planet Jupiter and the ultra-cool dwarf stars and brown dwarfs). We conclude that well-ordered and stable rotating magnetospheres share a common physical mechanism for supporting the generation of non-thermal electrons.

Our understanding of the unification of jetted AGNs has advanced greatly as the size of extragalactic sources increased. In the present paper, based on the large sample of radio sources, we compiled 680 blazars (279 FSRQs and 401 BL Lacs) from the 3FGL sample and 64 Seyfert galaxies (3 Narrow-line, 34 and 27 regular Seyfert 1 and Seyfert 2 respectively) from the INTEGRAL survey to statistically test the relationship between Seyfert galaxies and the blazar samples of FSRQs and BL Lacs. We compute the synchrotron (SS), Compton (CS) and inverse Compton (IC) continuous spectra from the low energy components of radio to X-ray, radio to gamma-ray and the high energy component of X-ray to gamma-ray bands, respectively. Results show from the distributions of the continuous spectra that Seyfert galaxies form the tail of the distributions, suggestive of similar underlying history and evolution. A two-sample Kolmogorov-Smirnov test (K.S test) of the continuous spectra showed that Seyfert galaxies differ from BL Lacs and FSRQs in the low energy components of the spectra, while there is no clear difference between them in the high energy component, which implies that high energy emissions in Seyfert galaxies, BL Lacs and FSRQs may be as a result of the same emission mechanism. There is a regular sequence of the distributions on SS/CS and IC/CS planes in each individual subsample. Linear regression analyses of our sample yield significant positive correlations (r greater than 0.60) between SS/CS and IC/CS data. This upturns into an anti-correlation (r greater than -0.60) in IC/SS data. These results are not only consistent with unified scheme for blazars but also show that Seyfert galaxies can be unified with the classical radio-loud AGNs counterparts.

We study correlation analysis for monopole components of stochastic gravitational wave backgrounds, including the maximally allowed polarization degrees. We show that, for typical detector networks, the correlation analysis can probe virtually five spectra: three for the intensities of the tensor, vector, and scalar modes and two for the chiral asymmetries of the tensor and vector modes. The chiral asymmetric signal for the vector modes has been left untouched so far. In this paper, we derive the overlap reduction function for this signal and thus complete the basic ingredients required for widely dealing with polarization degrees. We comprehensively analyze the geometrical properties of all the five overlap reduction functions. In particular, we point out the importance of reflection transformations for configuring preferable networks in the future.

We present ALMA CO ($J$=2--1) and 1.3 mm continuum observations of the high-velocity jet associated with the FIR 6b protostar located in the Orion Molecular Cloud-2. We detect a velocity gradient along the short axis of the jet in both the red- and blue-shifted components. The position-velocity diagrams along the short axis of the red-shifted jet show a typical characteristic of a rotating cylinder. We attribute the velocity gradient in the red-shifted component to rotation of the jet. The rotation velocity ($>20\,\ \rm{km s^{-1}}$) and specific angular momentum ($>10^{22}\, \rm{cm^{2}\, s^{-1}}$) of the jet around FIR 6b are the largest among all jets in which rotation has been observed. By combining disk wind theory with our observations, the jet launching radius is estimated to be in the range of $2.18-2.96$\,au. The rapid rotation, large specific angular momentum, and a launching radius far from the central protostar can be explained by a magnetohydrodynamic disk wind that contributes to the angular momentum transfer in the late stages of protostellar accretion.

We present optical spectroscopy together with ultraviolet, optical and near-infrared photometry of SN 2019hcc, which resides in a host galaxy at redshift 0.044, displaying a sub-solar metallicity. The supernova spectrum near peak epoch shows a `w' shape at around 4000 {\AA} which is usually associated with O II lines and is typical of Type I superluminous supernovae. SN 2019hcc post-peak spectra show a well-developed H alpha P-Cygni profile from 19 days past maximum and its light curve, in terms of its absolute peak luminosity and evolution, resembles that of a fast-declining Hydrogen-rich supernova (SN IIL). The object does not show any unambiguous sign of interaction as there is no evidence of narrow lines in the spectra or undulations in the light curve. Our tardis spectral modelling of the first spectrum shows that Carbon, Nitrogen and Oxygen (CNO) at 19000 K reproduce the `w' shape and suggests that a combination of non-thermally excited CNO and metal lines at 8000 K could reproduce the feature seen at 4000 {\AA}. The Bolometric light curve modelling reveals that SN 2019hcc could be fit with a magnetar model, showing a relatively strong magnetic field (B > 3 x 10^14 G), which matches the peak luminosity and rise time without powering up the light curve to superluminous luminosities. The high-energy photons produced by the magnetar would then be responsible for the detected O II lines. As a consequence, SN 2019hcc shows that a `w' shape profile at around 4000 {\AA}, usually attributed to O II, is not only shown in superluminous supernovae and hence it should not be treated as the sole evidence of the belonging to such a supernova type.

Context. We address the heating of the solar chromosphere and the related generation of plasma inflows and outflows. Aims. We attempt to detect variations in ion temperature and vertical plasma flows, which are driven by impulsively excited two-fluid Alfv\'en waves. We aim to investigate the possible contribution of these waves to solar chromosphere heating and plasma outflows. Methods. We performed numerical simulations of the generation and evolution of Alfv\'en waves with the use of the JOANNA code, which solves the two-fluid equations for ions+electrons and neutrals, coupled by collision terms. Results. We confirm that the damping of impulsively generated small-amplitude Alfv\'en waves slightly affects the temperature of the chromosphere and generates slow plasma flows. In contrast, the Alfv\'en waves generated by large-amplitude pulses increase the chromospheric plasma temperature more significantly and result in faster plasma outflows. The maximum heating occurs when the pulse is launched from the central photosphere, and the magnitude of the related plasma flows grows with the amplitude of the pulse. Conclusions. Large-amplitude two-fluid Alfv\'en waves can contribute significantly to the heating of the solar chromosphere and to the generation of plasma outflows.

This work studies the candidate double and multiple open clusters (OCs) in the galactic sector from l = 240o to l = 270o, which contains the Vela-Puppis star formation region. To do that, we have searched the most recent and complete catalogues of OCs by hand to get an extensive list of 22 groups of OCs involving 80 candidate members. Gaia EDR3 has been used to review some of the candidate OCs and look for new OCs near the candidate groups. Gaia data also permitted filtering out most of the field sources that are not member stars of the OCs. The plotting of combined colour-magnitude diagrams of candidate pairs has allowed, in several cases, endorsing or discarding their link. The most likely systems are formed by OCs less than 0.1 Gyr old, with only one eccentric OC in this respect. No probable system of older OCs has been found. Preliminary estimations of the fraction of known OCs that form part of groups (9.4 to 15%) support the hypothesis that the Galaxy and the Large Magellanic Cloud are similar in this respect. The results indicate that OCs are born in groups like stars are born in OCs.

Recent models for evolved Low Mass Stars (with $M \lesssim 3M_\odot$), undergoing the AGB phase assume that magnetic flux-tube buoyancy drives the formation of $^{13}$C reservoirs in He-rich layers. We illustrate their crucial properties, showing how the low abundance of $^{13}$C generated below the convective envelope hampers the formation of primary $^{14}$N and the ensuing synthesis of intermediate-mass nuclei, like $^{19}$F and $^{22}$Ne. In the mentioned models, their production is therefore of a purely secondary nature. Shortage of primary $^{22}$Ne has also important effects in reducing the neutron density. Another property concerns AGB winds, which are likely to preserve C-rich subcomponents, isolated by magnetic tension, even when the envelope composition is O-rich. Conditions for the formation of C-rich compounds are therefore found in stages earlier than previously envisaged. These issues, together with the uncertainties related to several nuclear physics quantities, are discussed in the light of the isotopic admixtures of s-process elements in presolar SiC grains of stellar origin, which provide important and precise constraints to the otherwise uncertain parameters. By comparing nucleosynthesis results with measured SiC data, it is argued that such a detailed series of constraints indicates the need for new measurements of weak interaction rates in ionized plasmas, as well as of neutron-capture cross sections, especially near the N = 50 and N = 82 neutron magic numbers. Nontheless, the peculiarity of our models allows us to achieve fits to the presolar grain data of a quality so far never obtained in previously published attempts.

Context. The origin of the heating of the solar atmosphere is still an unsolved problem. As the photosphere and chromosphere radiate more energy than the solar corona, it is challenging but important to reveal all the mechanisms that contribute to plasma heating there. Ion-neutral collisions could play an important role. Aims. We aim to investigate the impulsively generated two-fluid magnetoacoustic waves in the partially ionized solar chromosphere and to study the associated heating and plasma outflows, which higher up may result in nascent solar wind. Methods. To describe the plasma dynamics, we applied a two-fluid model in which ions+electrons and neutrals are treated as separate fluids. We solved the two-fluid equations numerically using the JOANNA code. Results. We show that magnetoacoustic waves triggered in the photosphere by localised velocity pulses can steepen into shocks which heat the chromosphere through ion-neutral collisions. Pulses of greater amplitude heat plasma more effectively and generate larger plasma outflows. Rising the altitude at which the pulse is launched results in opposite effects, mainly in local cooling of the chromosphere and slower plasma outflows. Conclusions. Even a solitary pulse results in a train of waves. These waves can transform into shock waves and release thermal energy, heating the chromosphere significantly. A pulse can drive vertical flows which higher up can result in the origin of the solar wind.

Open clusters provide clues to understand the evolution of Li7 at the surface of low-mass stars and its possible correlation with stellar rotation, which is a challenge for both stellar hydrodynamics and Galactic chemical evolution. We aim to quantify the efficiency of the transport processes for both angular momentum and chemicals that are required to explain simultaneously the observed behaviour of surface Li7 and rotation as well as the internal rotation profiles inferred from helio- and asteroseismology in F- and G-type main sequence stars. We apply the model for the transport of angular momentum and chemicals that we tailored in a previous work for solar-type stars to an extended range of initial masses and metallicities corresponding to F- an G-type stars in a sample of 20 Galactic open clusters. We evaluate its ability to explain the Li7, Be9, and rotation periods observations. Over the entire range of masses, metallicities, and ages explored, we reproduce the evolution of the surface rotation rates and predict, for the first time, the observed anti-correlation between the surface rotation rate and Li7 depletion as a consequence of the penetrative convection prescription. However, the ability of the model to reproduce the so-called Li7 dip centred around 6600K strongly depends on the adopted prescriptions for shear turbulence. It also requires a stellar mass dependence for the viscosity adopted for the transport of angular momentum, similar to the behaviour predicted for the generation and luminosity of internal gravity waves generated by stellar convective envelopes. We provide an efficient way to model G-type stars of different ages and metallicities successfully. However, the Li7 and Be9 dip constraints call for further hydrodynamical studies to better model turbulence in stars.

Context: Resistive Ohmic dissipation has been suggested as a mechanism for heating the solar chromosphere, but few studies have established this association. Aim: We aim to determine how Ohmic dissipation by electric currents can heat the solar chromosphere. Methods: We combine high-resolution spectroscopic Ca II data from the Dunn Solar Telescope and vector magnetic field observations from the Helioseismic and Magnetic Imager (HMI) to investigate thermal enhancements in a sunspot light bridge. The photospheric magnetic field from HMI was extrapolated to the corona using a non-force-free field technique that provided the three-dimensional distribution of electric currents, while an inversion of the chromospheric Ca II line with a local thermodynamic equilibrium and a nonlocal thermodynamic equilibrium spectral archive delivered the temperature stratifications from the photosphere to the chromosphere. Results: We find that the light bridge is a site of strong electric currents, of about 0.3 A/m^2 at the bottom boundary, which extend to about 0.7 Mm while decreasing monotonically with height. These currents produce a chromospheric temperature excess of about 600-800 K relative to the umbra. Only the light bridge, where relatively weak and highly inclined magnetic fields emerge over a duration of 13 hr, shows a spatial coincidence of thermal enhancements and electric currents. The temperature enhancements and the Cowling heating are primarily confined to a height range of 0.4-0.7 Mm above the light bridge. The corresponding increase in internal energy of 200 J/m^3 can be supplied by the heating in about 10 min. Conclusions: Our results provide direct evidence for currents heating the lower solar chromosphere through Ohmic dissipation.

We present results of variability study of a sample of 20 powerful blazars using Fermi/LAT (0.1--300 GeV) observations. We studied decade-long observations applying various analysis tools such as flux distribution, symmetry analysis, and RMS-flux relation. It was found that the $\gamma$-ray flux distribution closely resembles a log-normal probability distribution function and can be characterized by linear RMS-flux relation. The power spectral density analysis shows the statistical variability properties of the sources as studied are consistent with flicker noise, an indication of long-memory processes at work. Statistical analysis of the distribution of flux rise and decay rates in the light curves of the sources, aimed at distinguishing between particle acceleration and energy-dissipation timescales, counter-intuitively suggests that both kinds of rates follow a similar distribution and the derived mean variability timescales are on the order of a few weeks. The corresponding emission region size is used to constrain the location of $\gamma$-ray production sites in the sources to be a few parsecs. Additionally, using Lomb-Scargle periodogram and weighted wavelet z-transform methods and extensive Monte Carlo simulations, we detected year-timescale quasi-periodic oscillations in the sources S5 0716+714, Mrk 421, ON +325, PKS 1424-418, and PKS 2155-304. We also performed recurrence quantification analysis of the sources and directly measure the deterministic quantities, which suggest that the dynamical processes in blazars could be a combination of deterministic and stochastic processes, while some of the source light curves revealed significant deterministic content.

The field of deep learning has become increasingly important for particle physics experiments, yielding a multitude of advances, predominantly in event classification and reconstruction tasks. Many of these applications have been adopted from other domains. However, data in the field of physics are unique in the context of machine learning, insofar as their generation process and the laws and symmetries they abide by are usually well understood. Most commonly used deep learning architectures fail at utilizing this available information. In contrast, more traditional likelihood-based methods are capable of exploiting domain knowledge, but they are often limited by computational complexity. In this contribution, a hybrid approach is presented that utilizes generative neural networks to approximate the likelihood, which may then be used in a traditional maximum-likelihood setting. Domain knowledge, such as invariances and detector characteristics, can easily be incorporated in this approach. The hybrid approach is illustrated by the example of event reconstruction in IceCube.

In massive objects, such as galaxy clusters, the turbulent velocity dispersion, $\sigma_\mathrm{turb}$, is tightly correlated to both the object mass, $M$, and the thermal energy. Here, we investigate whether these scaling laws extend to lower-mass objects in dark-matter filaments. We perform a cosmological zoom-in simulation of a filament using an adaptive filtering technique for the resolved velocity component and a subgrid-scale model to account for the unresolved component. We then compute the mean turbulent and thermal energies for all halos in the zoom-in region and compare different definitions of halo averages. Averaging constrained by density and temperature thresholds is favored over averages solely based on virial spheres. We find no clear trend for the turbulent velocity dispersion versus halo mass, but significant correlation and a scaling law with exponent $\alpha\sim 0.5$ between the turbulent velocity dispersion and thermal energy that agrees with a nearly constant turbulent Mach number, similar to more massive objects. We conclude that the self-similar energetics proposed for galaxy clusters extends down to the CGM of individual galaxies.

Recent results from IceCube regarding TXS 0506+056 suggest the presence of neutrino flares that are not temporally coincident with a significant corresponding gamma ray flare. Such flares are particularly difficult to identify, as their presence must be inferred from the temporal distribution of neutrino data alone. Here we present the results of using a novel method to search for all such flares across the entire neutrino sky in 10 years of IceCube data, using both Gaussian and box-shaped flare hypotheses. Unlike for past searches, that looked for only the most significant neutrino flare in the data at a given direction, here we implement an algorithm to combine information from multiple flares associated with a single source candidate. This represents the most detailed description of the neutrino sky to date, providing the location and intensity of all neutrino cluster candidates in both space and time. These results can be used to further constrain potential populations of transient neutrino sources, serving as a complement to existing time-integrated and time-dependent methods.

As one of the main formation mechanisms of solar filament formation, the chromospheric evaporation-coronal condensation model has been confirmed by numerical simulations to explain the formation of filament threads very well in flux tubes with single dips. However, coronal magnetic extrapolations indicated that some magnetic field lines might possess more than one dip. It is expected that the formation process would be significantly different in this case compared to a single-dipped magnetic flux tube. In this paper, based on the evaporation-condensation model, we study filament thread formation in double-dipped magnetic flux tubes by numerical simulations. We find that only with particular combinations of magnetic configuration and heating, e.g., concentrated localized heating and a long magnetic flux tube with deep dips, can two threads form and persist in a double-dipped magnetic flux tube. Comparing our parametric survey with observations, we conclude that such magnetically connected threads due to multiple dips are more likely to exist in quiescent filaments than in active-region filaments. Moreover, we find that these threads are usually shorter than independently trapped threads, which might be one of the reasons why quiescent filaments have short threads. These characteristics of magnetically connected threads could also explain barbs and vertical threads in quiescent filaments.

An upgrade to the IceCube Neutrino Telescope is currently under construction. For this IceCube Upgrade, seven new strings will be deployed in the central region of the 86 string IceCube detector to enhance the capability to detect neutrinos in the GeV range. One of the main science objectives of the IceCube Upgrade is an improved calibration of the IceCube detector to reduce systematic uncertainties related to the optical properties of the ice. We have developed a novel optical camera and illumination system that will be part of 700 newly developed optical modules to be deployed with the IceCube Upgrade. A combination of transmission and reflection photographic measurements will be used to measure the optical properties of bulk ice between strings and refrozen ice in the drill hole, to determine module positions, and to survey the local ice environments surrounding the sensor module. In this contribution we present the production design, acceptance testing, and plan for post-deployment calibration measurements with the camera system.

The IceCube Neutrino Observatory is a cubic-kilometer scale neutrino detector embedded in the Antarctic ice of the South Pole. In the near future, the detector will be augmented by extensions, such as the IceCube Upgrade and the planned Gen2 detector. The sparseness of observed light in the detector for low-energy events, and the irregular detector geometry, have always been a challenge to the reconstruction of the detected neutrinos' parameters of interest. This challenge remains with the IceCube Upgrade, currently under construction, which introduces seven new detector strings with novel detector modules. The Upgrade modules will increase the detection rate of low-energy events and allow us to further constrain neutrino oscillation physics. However, the geometry of these modules render existing traditional reconstruction algorithms more difficult to use. We introduce a new reconstruction algorithm based on Graph Neural Networks, which we use to reconstruct neutrino events at much faster processing times than the traditional algorithms, while providing comparable resolution. We show that our algorithm is applicable not only to reconstructing data of the current IceCube detector, but also simulated events for next-generation extensions, such as the IceCube Upgrade.

We report the frequency analysis of a known roAp star, HD 86181 (TIC 469246567), with new inferences from TESS data. We derive the rotation frequency to be $\nu_{rot}$ = 0.48753 $\pm$ 0.00001d$^{-1}$. The pulsation frequency spectrum is rich, consisting of two doublets and one quintuplet, which we interpret to be oblique pulsation multiplets from consecutive, high-overtone dipole, quadrupole and dipole modes. The central frequency of the quintuplet is 232.7701d$^{-1}$ (2.694 mHz). The phases of the sidelobes, the pulsation phase modulation, and a spherical harmonic decomposition all show that the quadrupole mode is distorted. Following the oblique pulsator model, we calculate the rotation inclination, i, and magnetic obliquity, $\beta$, of this star, which provide detailed information about the pulsation geometry. The i and $\beta$ derived from the best fit of the pulsation amplitude and phase modulation to a theoretical model, including the magnetic field effect, slightly differ from those calculated for a pure quadrupole, indicating the contributions from l = 4, 6, 8, ... are small. Non-adiabatic models with different envelope convection conditions and physics configurations were considered for this star. It is shown that models with envelope convection almost fully suppressed can explain the excitation at the observed pulsation frequencies.

The recent NANOGrav finding of a common-spectrum process has invited interpretations as possible evidence of a primordial stochastic gravitational-wave background (SGWB) stronger than predicted by standard inflation+LCDM. Such an SGWB would contribute an extra radiation component to the background Universe which may affect its expansion history. As such, it may help alleviate the current Hubble tension, a novel connection between gravitational waves and cosmology. We demonstrate this by considering a cosmological model, the "standard inflation + stiff amplification" scenario, with two components added to the LCDM model: a stiff component (w=1) and the primordial SGWB. Previously, we showed that even for standard inflation, the SGWB may be detectable at the high frequencies probed by laser interferometers, if it is amplified by a possible early stiff era after reheating. Models that boost the SGWB enough to cause significant backreaction, however, must still preserve the well-measured radiation-matter equality, as precision cosmology demands. For that, we calculate the fully-coupled evolution of the SGWB and expansion history, sampling parameter space (tensor-to-scalar ratio, reheating temperature and temperature at stiff-to-radiation equality). We then perform a joint analysis of the NANOGrav results and latest upper bounds from Planck, big bang nucleosynthesis and Advanced LIGO-Virgo, to constrain the model. The resulting blue-tilted, stiff-amplified SGWB is still too small to explain the NANOGrav results. However, if someday, Advanced LIGO-Virgo detects the SGWB, our model can explain it within standard inflation (without requiring an initial blue tilt). Meanwhile, this model may bring current high-z measurements of the Hubble constant within 3.4 sigma of the low-z measurements by SH0ES (from 4.4 sigma) and within 2.6 sigma of those by H0LiCOW (from 3.1 sigma), reducing the tension.

The most accepted scenario for the evolution of massive galaxies across cosmic time predicts a regulation based on the interplay between AGN feedback, which injects large amounts of energy in the host environment, and galaxy mergers, being able to trigger massive star formation events and accretion onto the supermassive black holes. Interacting systems hosting AGN are useful laboratories to get key insights into both phenomena. In this context, we present the analysis of the optical spectral properties of IRAS 20210+1121 (I20210), a merging system at $z = 0.056$. According to X-ray data, this object comprises two interacting galaxies, each hosting an obscured AGN. The optical spectra confirm the presence of AGN features in both galaxies. In particular, we are able to provide a Seyfert classification for I20210 North. The spectrum of I20120 South shows broad blueshifted components associated with the most intense emission lines that indicate the presence of an ionized outflow, for which we derive a maximum velocity of $\sim$2000 km s$^{-1}$, an extension of $\sim$2 kpc and a mass rate of $\sim$0.6 M$_\odot$ yr$^{-1}$. We also report the existence of an ionized nebular component with $v \sim 1000$ km s$^{-1}$ at $\sim$6.5 kpc Southwards of I20210 South, that can be interpreted as disrupted gas ejected from the host galaxy by the action of the outflow. I20120 therefore exhibits a double obscured AGN, with one of them showing evidence of ongoing events for AGN-powered outflows. Future spatially-resolved spectroscopy will allow to accurately map on the gas kinematics in this AGN pair and evaluate the impact of the outflow on both the interstellar medium and galaxy environment.

Stellar tidal disruption events (TDEs) are typically discovered by transient emission due to accretion or shocks of the stellar debris. Yet this luminous flare can be reprocessed by gas or dust that inhabits a galactic nucleus, resulting in multiple reverberation signals. Nuclear dust heated by the TDE will lead to an echo at infrared wavelengths (1-10 $\mu$m) and transient coronal lines in optical spectra of TDEs trace reverberation by gas that orbits the black hole. Both of these signal have been detected, here we review this rapidly developing field. We also review the results that have been extracted from TDEs with high-quality X-ray light curves: quasi periodic oscillations (QPOs), reverberation lags of fluorescence lines, and cross-correlations with emission at other wavelengths. The observational techniques that are covered in this review probe the emission from TDEs over a wide range of scales: from light years to the innermost parts of the newly formed accretion disk. They provide insights into important properties of TDEs such as their bolometric output and the geometry of the accretion flow. While reverberation signals are not detected for every TDE, we anticipate they will become more commonplace when the next generation of X-ray and infrared instruments become operational.

Many experiments have confirmed the spectral hardening at a few hundred GV of cosmic-ray (CR) nuclei spectra, and 3 general different origins have been proposed: the primary source acceleration, the propagation, and the superposition of different kinds of sources. Here we report some new findings from the AMS-02 nuclei spectra of B and its dominating parents species (C, N, O, Ne, Mg, and Si): the nuclei spectral hardening in a few hundred GV should have hybrid origins. Besides the propagation origin, the superposition of different kinds of sources are also needed for different kinds of the CR primary nuclei species. All these results can be further confirmed by more precise CR nuclei spectra data in high rigidity regions (like that from DAMPE), and could provide us an opportunity to improve the current CR models.

In this paper, we report on the observational performance of the Swift Ultra-violet/Optical Telescope (UVOT) in response to the Gravitational Wave alerts announced by the Advanced Laser Interferometer Gravitational Wave Observatory and the Advanced Virgo detector during the O3 period. We provide the observational strategy for follow-up of GW alerts and provide an overview of the processing and analysis of candidate optical/UV sources. For the O3 period, we also provide a statistical overview and report on serendipitous sources discovered by Swift/UVOT. Swift followed 18 gravitational-wave candidate alerts, with UVOT observing a total of 424 deg^2. We found 27 sources that changed in magnitude at the 3 sigma level compared with archival u or g-band catalogued values. Swift/UVOT also followed up a further 13 sources reported by other facilities during the O3 period. Using catalogue information, we divided these 40 sources into five initial classifications: 11 candidate active galactic nuclei (AGN)/quasars, 3 Cataclysmic Variables (CVs), 9 supernovae, 11 unidentified sources that had archival photometry and 6 uncatalogued sources for which no archival photometry was available. We have no strong evidence to identify any of these transients as counterparts to the GW events. The 17 unclassified sources are likely a mix of AGN and a class of fast-evolving transient, and one source may be a CV.

We report the discovery of J0121+0025, an extremely luminous and young star-forming galaxy (M_UV = -24.11, log[L_Lya / erg s^-1] = 43.8) at z = 3.244 showing copious Lyman continuum (LyC) leakage (f_esc,abs ~ 40%). High SNR rest-frame UV spectroscopy with the Gran Telescopio Canarias reveals a high significance (7.9 sigma) emission below the Lyman limit (< 912A), with a flux density level f_900A = 0.78 +/- 0.10 uJy, and strong P-Cygni in wind lines of OVI 1033A, NV 1240A and CIV 1550A that are indicative of a young age of the starburst (<10 Myr). The spectrum is rich in stellar photospheric features, for which a significant contribution of an AGN at these wavelengths is ruled out. Low-ionization ISM absorption lines are also detected, but are weak (EW0 ~ 1A) and show large residual intensities, suggesting a clumpy geometry of the gas with a non-unity covering fraction or a highly ionized ISM. The contribution of a foreground and AGN contamination to the LyC signal is unlikely. Deep optical to Spitzer/IRAC 4.5um imaging show that the spectral energy distribution of J0121+0025 is dominated by the emission of the young starburst, with log(M*/Msun) = 9.9 +/- 0.1 and SFR = 981 +/- 232 Msun yr^-1. J0121+0025 is the most powerful LyC emitter known among the star-forming galaxy population. The discovery of such luminous and young starburst leaking LyC radiation suggests that a significant fraction of LyC photons can escape in sources with a wide range of UV luminosities and are not restricted to the faintest ones as previously thought. These findings might shed further light on the role of luminous starbursts to the cosmic reionization.

We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Criterion (BIC) and sampling joint posterior parameter distributions using MCMC methods. For analysis of quasar lightcurves, our implementation of ROA modelling can inter-calibrate lightcurve data from different telescopes, estimate the shape and thus the power-density spectrum of the lightcurve, and measure time delays among lightcurves at different wavelengths or from different images of a lensed quasar. Our noise model implements a robust treatment of outliers and error-bar adjustments to account for additional variance or poorly-quantified uncertainties. Tests with simulated data validate the parameter uncertainty estimates. We compare ROA delay measurements with results from cross-correlation and from JAVELIN, which models lightcurves with a prior on the power-density spectrum. We analyse published COSMOGRAIL lightcurves of multi-lensed quasar lightcurves and present the resulting measurements of the inter-image time delays and detection of microlensing effects.

Using numerical simulations of helical inflationary magnetogenesis in a low reheating temperature scenario, we show that the magnetic energy spectrum is strongly peaked at a particular wavenumber that depends on the reheating temperature. Gravitational waves (GWs) are produced at frequencies between 3 nHz and 50 mHz for reheating temperatures between 150 MeV and 3x10^5 GeV, respectively. At and below the peak frequency, the stress spectrum is always found to be that of white noise. This implies a linear increase of GW energy per logarithmic wavenumber interval, instead of a cubic one, as previously thought. Both in the helical and nonhelical cases, the GW spectrum is followed by a sharp drop for frequencies above the respective peak frequency. In this magnetogenesis scenario, the presence of a helical term extends the peak of the GW spectrum and therefore also the position of the aforementioned drop toward larger frequencies compared to the case without helicity. This might make a difference in it being detectable with space interferometers. The efficiency of GW production is found to be almost the same as in the nonhelical case, and independent of the reheating temperature, provided the electromagnetic energy at the end of reheating is fixed to be a certain fraction of the radiation energy density. Also, contrary to the case without helicity, the electric energy is now less than the magnetic energy during reheating. The fractional circular polarization is found to be nearly hundred per cent in a certain range below the peak frequency range.

The Global Meteor Network (GMN) utilizes highly sensitive low-cost CMOS video cameras which run open-source meteor detection software on Raspberry Pi computers. Currently, over 450 GMN cameras in 30 countries are deployed. The main goal of the network is to provide long-term characterization of the radiants, flux, and size distribution of annual meteor showers and outbursts in the optical meteor mass range. The rapid 24-hour publication cycle the orbital data will enhance the public situational awareness of the near-Earth meteoroid environment. The GMN also aims to increase the number of instrumentally observed meteorite falls and the transparency of data reduction methods. A novel astrometry calibration method is presented which allows decoupling of the camera pointing from the distortion, and is used for frequent pointing calibrations through the night. Using wide-field cameras ($88^{\circ}\times48^{\circ}$) with a limiting stellar magnitude of $+6.0 \pm 0.5$ at 25 frames per second, over 220,000 precise meteoroid orbits were collected since December 2018 until June 2021. The median radiant precision of all computed trajectories is $0.47^{\circ}$, $0.32^{\circ}$ for $\sim20\%$ of meteors which were observed from 4+ stations, a precision sufficient to measure physical dispersions of meteor showers. All non-daytime annual established meteor showers were observed during that time, including five outbursts. An analysis of a meteorite-dropping fireball is presented which showed visible wake, fragmentation details, and several discernible fragments. It had spatial trajectory fit errors of only \SI{\sim40}{\metre}, which translated into the estimated radiant and velocity errors of 3 arc minutes and tens of meters per second.

DBSP_DRP is a python package that provides fully automated data reduction of data taken by the Double Spectrograph (DBSP) at the 200-inch Hale Telescope at Palomar Observatory (Oke & Gunn, 1982). The underlying data reduction functionality to extract 1D spectra, perform flux calibration and correction for atmospheric absorption, and coadd spectra together is provided by PypeIt (Prochaska et al., 2020). The new functionality that DBSP_DRP brings is in orchestrating the complex data reduction process by making smart decisions so that no user input is required after verifying the correctness of the metadata in the raw FITS files in a table-like GUI. Though the primary function of DBSP_DRP is to automatically reduce an entire night of data without user input, it has the flexibility for astronomers to fine-tune the data reduction with GUIs for manually identifying the faintest objects, as well as exposing the full set of PypeIt parameters to be tweaked for users with particular science needs. DBSP_DRP also handles some of the occasional quirks specific to DBSP, such as swapping FITS header cards, adding (an) extra null byte/s to FITS files making them not conform to the FITS specification, and not writing the coordinates of the observation to file. Additionally, DBSP_DRP contains a quicklook script for making real-time decisions during an observing run, and can open a GUI displaying a minimally reduced exposure in under 15 seconds. Docker containers are available for ease of deploying DBSP_DRP in its quicklook configuration (without some large atmospheric model files) or in its full configuration.

A simulation tool which utilizes parallel processing is developed to describe molecular kinetics in 2D, single-and multi-component atmospheres on Callisto. This expands on our previous study on the role of collisions in 1D atmospheres on Callisto composed of radiolytic products (Carberry Mogan et al., 2020) by implementing a temperature gradient from noon to midnight across Callisto's surface and introducing sublimated water vapor. We compare single-species, ballistic and collisional O2, H2 and H2O atmospheres, as well as an O2+H2O atmosphere to 3-species atmospheres which contain H2 in varying amounts. Because the H2O vapor pressure is extremely sensitive to the surface temperatures, the density drops several order of magnitude with increasing distance from the subsolar point, and the flow transitions from collisional to ballistic accordingly. In an O2+H2O atmosphere the local temperatures are determined by H2O near the subsolar point and transition with increasing distance from the subsolar point to being determined by O2 When radiolytically produced H2 is not negligible in O2+H2O+H2 atmospheres, this much lighter molecule, with a scale height roughly an order of magnitude larger than that for the heavier species, can cool the local temperatures via collisions. In addition, if the H2 component is dense enough, particles originating on the day-side and precipitating into the night-side atmosphere deposit energy via collisions, which in turn heats the local atmosphere relative to the surface temperature. Finally, we discuss the potential implications of this study on the presence of H2 in Callisto's atmosphere and how the simulated densities correlate with expected detection thresholds at flyby altitudes of the proposed JUpiter ICy moons Explorer (JUICE) spacecraft.

Molecular kinetic simulations are typically used to accurately describe the tenuous regions of the upper atmospheres on planetary bodies. These simulations track the motion of particles representing real atmospheric atoms and/or molecules subject to collisions, the object's gravity, and external influences. Because particles can end up in very large ballistic orbits, upper boundary conditions (UBC) are typically used to limit the domain size thereby reducing the time for the atmosphere to reach steady-state. In the absence of a clear altitude at which all molecules are removed, such as a Hill sphere, an often used condition is to choose an altitude at which collisions become infrequent so that particles on escape trajectories are removed. The remainder are then either specularly reflected back into the simulation domain or their ballistic trajectories are calculated analytically or explicitly tracked so they eventually re-enter the domain. Here we examine the effect of the choice of the UBC on the escape rate and the structure of the atmosphere near the nominal exobase in the convenient and frequently used 1D spherically symmetric approximation. Using Callisto as the example body, we show that the commonly used specular reflection UBC can lead to significant uncertainties when simulating a species with a lifetime comparable to or longer than a dynamical time scale, such as an overestimation of escape rates and an inflated exosphere. Therefore, although specular reflection is convenient, the molecular lifetimes and body's dynamical time scales need to be considered even when implementing the convenient 1D spherically symmetric simulations in order to accurately estimate the escape rate and the density and temperature structure in the transition regime.

The third iteration of the International Celestial Reference Frame (ICRF3) is made up of 4536 quasars observed at S/X bands using Very Long baseline Interferometry (VLBI). These sources are high redshift quasars, typically between $1<z<2$, that are believed to host active galactic nuclei (AGN) at their centers. The position of compact radio sources can be determined better than sources with large amounts of extended radio structure. Here we report information on a series of 20 observations from January 2017 through December 2017 which were designed for precise astrometry and to monitor the structure of sources included in the ICRF3. We targeted 3627 sources over the one year campaign and found the median flux density of 2697 detected sources at S-band is 0.13 Jy, and the flux density of 3209 sources detected at X-band is 0.09 Jy. We find that $70\%$ of detected sources in our campaign are considered compact at X-band and ideal for use in the ICRF and $89\%$ of the 2615 sources detected at both frequencies have a flat spectral index, $\alpha>0.5$

Spin is a fundamental property of Primordial Black Holes, and it can substantially affect Primordial Black Holes (PBHs) evaporation rate. We consider the non-negligible spin of PBHs and study the bounds on the fraction of PBHs that accounts for dark matter as a function of mass and spin. Spinning Primordial Black Holes (SPBHs) can evaporate by Hawking radiation and inject energy into the intergalactic medium (IGM). The injection of the background radiation into IGM can significantly affect the thermal and ionization history of gas during cosmic dawn. Recently, the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) low-band observation reported a 21 cm absorption signal in the redshift range 15-20. The presence of Hawking radiation by evaporating PBHs can modify the 21 cm absorption profile. Considering the EDGES signal into account, we constrain the fraction of SPBHs that accounts for dark matter as a function of mass and spin. We also compare these bound with other available observations.

Young brown dwarfs are analogs to giant exoplanets, as they share effective temperatures, near-infrared colors, and surface gravities. Thus, the detailed characterization of young brown dwarfs might shed light on the study of giant exoplanets, that we are currently unable to observe with sufficient signal-to-noise to allow precise characterization of their atmospheres. 2MASS J22081363+2921215 is a young L3 brown dwarf, member of the beta-Pictoris young moving group (23 +/-3 Myr), that shares its effective temperature and mass with the beta Pictoris b giant exoplanet. We performed a ~2.5 hr spectro-photometric J-band monitoring of 2MASS J22081363+2921215 with the MOSFIRE multi-object spectrograph, installed at the Keck I telescope. We measured a minimum variability amplitude of 3.22 +/- 0.42 % for its J-band light curve. The ratio between the maximum and the minimum flux spectra of 2MASS J22081363+2921215 shows a weak wavelength dependence and a potential enhanced variability amplitude in the alkali lines. Further analysis suggests that the variability amplitude on the alkali lines is higher than the overall variability amplitude (4.5-11 %, depending on the lines). The variability amplitudes in these lines are lower if we degrade the resolution of the original MOSFIRE spectra to R~100, which explains why this potential enhanced variability in the alkali lines has not been found previously in HST/WFC3 light curves. Using radiative-transfer models, we obtained the different cloud layers that might be introducing the spectro-photometric variability we observe for 2MASS J22081363+2921215, that further support the measured enhanced variability amplitude in the alkali lines. We provide an artistic recreation of the vertical cloud structure of this beta-Pictoris b analog.

Long period comet C/2021 A1 (Leonard) will approach Venus to within 0.029 au on 2021 December 18 and may subsequently graze the planet with its dust trail less than two days later. We observed C/2021 A1 with the Lowell Discovery Telescope on 2021 January 13 and March 3, as well as with the Palomar Hale Telescope on 2021 March 20, while the comet was inbound at heliocentric distances of r=4.97 au, 4.46 au, and 4.28 au, respectively. Tail morphology suggests that the dust is optically dominated by ~0.1-1 mm radius grains produced in the prior year. Neither narrowband imaging photometry nor spectrophotometry reveal any definitive gas emission, placing 3-sigma upper bounds on CN production of <1e23 molec/s at both of the latter two epochs. Trajectory analysis indicates that large (>1 mm) grains ejected at extremely large heliocentric distances (r>30 au) are most strongly favored to reach Venus. The flux of such meteors on Venus, and thus their potential direct or indirect observability, is highly uncertain as the comet's dust production history is poorly constrained at these distances, but will likely fall well below the meteor flux from comet C/2013 A1 (Siding Spring)'s closer encounter to Mars in 2014, and thus poses negligible risk to any spacecraft in orbit around Venus. Dust produced in previous apparitions will not likely contribute substantially to the meteor flux, nor will dust from any future activity apart from an unlikely high speed (>0.5 km/s) dust outburst prior to the comet reaching r~2 au in 2021 September.

Cosmology relies on a coarse-grained description of the universe, assumed to be valid on large length scales. However, the nonlinearity of general relativity makes coarse-graining extremely difficult. We here address this problem by extending the Mori-Zwanzig projection operator formalism, a highly successful coarse-graining method from statistical mechanics, towards general relativity. Using the Buchert equations, we derive a new dynamic equation for the Hubble parameter which captures the effects of averaging through a memory function. This gives an empirical prediction for the cosmic jerk.

Matrix theory is a proposed non-perturbative definition of superstring theory in which space is emergent. We begin a study of cosmology in the context of matrix theory. Specifically, we show that matrix theory can lead to an emergent non-singular cosmology which, at late times, can be described by an expanding phase of Standard Big Bang cosmology. The horizon problem of Standard Big Bang cosmology is automatically solved. We show that thermal fluctuations in the emergent phase source an approximately scale-invariant spectrum of cosmological perturbations and a scale-invariant spectrum of gravitational waves. Hence, it appears that matrix theory can lead to a successful scenario for the origin of perturbations responsible for the currently observed structure in the universe while providing a consistent UV-complete description.

According to the model ($\Lambda$CDM), based on deep cosmological observations, the current universe is constituted of 5$\%$ baryonic matter and 25 $\%$ non-baryonic cold dark matter (of speculative origin). These include quanta of scalar filed like dilaton($\phi$) of scale symmetry origin and quanta of pseudoscalar field of extra standard model symmetry ( Peccei-Quinn) origin, like axion ($\phi'$). These fields couple to di-photons through dim-5 operators. In magnetized medium, they in principle can interact with the three degrees of freedom (two transverse ($A_{\parallel,\perp}$) and one longitudinal ($A_{L}$)) of photon($\gamma$) as long as the total spin is conserved. Because of intrinsic spin being zero, both $\phi$ and $\phi'$ could in principle have interacted with $A_{L}$ (having $s_{z}=0$). However, out of $\phi$ and $\phi'$ only one interacts with $A_{L}$. Furthermore, the ambient external magnetic field and media, breaks the intrinsic Lorentz symmetry of the system. Invoking Charge conjugation, Parity and Time reversal symmetries, we analyse the mixing dynamics of $\phi\gamma$ and $\phi'\gamma$ systems and the structural {\it difference} of their mixing matrices. The electromagnetic signals (EMS) due to $\phi\gamma$ and $\phi'\gamma$ interactions as a result would be {\it different}. We conclude by commenting on the possibility of detecting this {\it difference} -- in the EMS -- using the existing space-borne detectors.

Simulation results from a global magnetohydrodynamic model of the solar corona and solar wind are compared with Parker Solar Probe (PSP) observations during its first five orbits. The fully three-dimensional model is based on Reynolds-averaged mean-flow equations coupled with turbulence transport equations. The model includes the effects of electron heat conduction, Coulomb collisions, turbulent Reynolds stresses, and heating of protons and electrons via a turbulent cascade. Turbulence transport equations for average turbulence energy, cross helicity, and correlation length are solved concurrently with the mean-flow equations. Boundary conditions at the coronal base are specified using solar synoptic magnetograms. Plasma, magnetic field, and turbulence parameters are calculated along the PSP trajectory. Data from the first five orbits are aggregated to obtain trends as a function of heliocentric distance. Comparison of simulation results with PSP data shows good agreement, especially for mean-flow parameters. Synthetic distributions of magnetic fluctuations are generated, constrained by the local rms turbulence amplitude given by the model. Properties of this computed turbulence are compared with PSP observations.

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduce the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass $ 2\times10^{-14}\lesssim \mu/\textrm{eV} \lesssim2\times10^{-13}$ would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to $\approx60\%$ more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate mass black hole mergers to $\approx0.08$\,Gpc$^{-3}$\,yr$^{-1}$ in the local Universe. These results imply that a select range of ultralight boson mass can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario.

We derive spectral lineshapes of the expected signal for a haloscope experiment searching for axionlike dark matter. The knowledge of these lineshapes is needed to optimize the experimental design and data analysis procedure. We extend the previously known results for the axion-photon and axion-gluon couplings to the case of gradient (axion-fermion) coupling. A unique feature of the gradient interaction is its dependence not only on magnitudes but also on directions of velocities of galactic halo particles, which leads to directional sensitivity of the corresponding haloscope. We also discuss the daily and annual modulations of the gradient signal caused by the Earth's rotational and orbital motions. In the case of detection, these periodic modulations will be an important confirmation that the signal is sourced by axionlike particles in the halo of our galaxy.

In an endeavor to explain the light neutrino masses and dark matter (DM) simultaneously, we study a gauged $U(1)_{\rm B-L}$ extension of the standard model (SM). The neutrino masses are generated through a variant of type-II seesaw mechanism in which one of the scalar triplets has a mass in a scale that is accessible at the present generation colliders. Three right chiral fermions $\chi_{iR}$($i=e,\mu,\tau$) with $\rm B-L$ charges -4, -4, +5 are invoked to cancel the $\rm B-L$ gauge anomalies and the lightest one among these three fermions becomes a viable DM candidate as their stability is guaranteed by a remnant $\mathcal Z_2$ symmetry to which $U(1)_{\rm B-L}$ gauge symmetry gets spontaneously broken. Interestingly in this scenario, the neutrino mass and the co-annihilation of DM are interlinked through the breaking of $U(1)_{\rm B-L}$ symmetry. Apart from giving rise to the observed neutrino mass and dark matter abundance, the model also predicts exciting signals at the colliders especially regarding the discovery of the triplet scalar in presence of the $\rm B-L$ gauge boson. We see a $(34-54)\%$ enhancement in the production of the TeV scale doubly charged scalar in presence of the $Z_{\rm BL}$ gauge boson in a mass range $2.5$ TeV to $4.4$ TeV. We discuss all the relevant constraints on model parameters from observed DM abundance and null detection of DM at direct and indirect search experiments as well as the constraints on the $\rm B-L$ gauge boson from recent colliders.

The detection of an intermediate-mass black hole population ($10^2-10^6\ M_\odot$) will provide clues to their formation environments (e.g., disks of active galactic nuclei, globular clusters) and illuminate a potential pathway to produce supermassive black holes. Ground-based gravitational-wave detectors are sensitive to a subset of such mergers and have been used to detect one $142^{+28}_{-16}\ M_\odot$ intermediate-mass black hole formation event. However, ground-based detector data contain numerous incoherent short duration noise transients that can mimic the gravitational-wave signals from merging intermediate-mass black holes, limiting the sensitivity of searches. Here we search for binary black hole mergers using a Bayesian-inspired ranking statistic which measures the coherence or incoherence of triggers in multiple-detector data. We use this statistic to identify candidate events with lab-frame total masses $\gtrsim55\ M_\odot$ using data from LIGO's second observing run. Our analysis does not yield evidence for new intermediate-mass black holes. However, we find support for some stellar-mass binary black holes not reported in the first LIGO--Virgo gravitational-wave transient catalog, GWTC-1.

In models of minicharged dark matter associated with a hidden $U(1)$ symmetry, astrophysical black holes may acquire a "dark" charge, in such a way that the inspiral dynamics of binary black holes can be formally described by an Einstein-Maxwell theory. Charges enter the gravitational wave signal predominantly through a dipole term, but their effect is known to effectively first post-Newtonian order in the phase, which enables measuring the size of the charge-to-mass ratios, $|q_i/m_i|$, $i = 1,2$, of the individual black holes in a binary. We set up a Bayesian analysis to discover, or constrain, dark charges on binary black holes. After testing our framework in simulations, we apply it to selected binary black hole signals from the second Gravitational Wave Transient Catalog (GWTC-2), namely those with low masses so that most of the signal-to-noise ratio is in the inspiral regime. We find no evidence for charges on the black holes, and place typical 1-$\sigma$ bounds on the charge-to-mass ratios of $|q_i/m_i| \lesssim 0.2 - 0.3$.

Phenomenological implications of the Mimetic Tensor-Vector-Scalar theory (MiTeVeS) are studied. The theory is an extension of the vector field model of mimetic dark matter, where a scalar field is also incorporated, and it is known to be free from ghost instability. In the absence of interactions between the scalar field and the vector field, the obtained cosmological solution corresponds to the General theory of Relativity (GR) with a minimally-coupled scalar field. However, including an interaction term between the scalar field and the vector field yields interesting dynamics. There is a shift symmetry for the scalar field with a flat potential, and the conserved Noether current, which is associated with the symmetry, behaves as a dark matter component. Consequently, the solution contains a cosmological constant, dark matter and a stiff matter fluid. Breaking the shift symmetry with a non-flat potential gives a natural interaction between dark energy and dark matter.

The IceCube Neutrino Observatory at the South Pole has tremendous emotional appeal -- the extreme Antarctic environment coupled with the aura of a pioneering experiment that explores the universe in a new way. However, as with most cutting-edge experiments, it is still challenging to translate the exotic, demanding science into accessible language. We present three examples of recent successful education, outreach, and communication activities that demonstrate how we leverage efforts and sustain connections to produce engaging results. First we describe our participation in the PolarTREC program, which pairs researchers with educators to provide deployments in the Antarctic, and how we have sustained relationships with these educators to produce high-quality experiences to reach target audiences even during a pandemic. We then focus on three activities from the past year: a summer enrichment program for high school students that was also modified for a 10-week IceCube after school program, a virtual visit to the South Pole for the ScienceWriters 2020 conference, and a series of short videos in English and Spanish suitable for all ages that explain traveling, living, and working at the South Pole.

The Multimessenger Diversity Network (MDN), formed in 2018, extends the basic principle of multimessenger astronomy -- that working collaboratively with different approaches enhances understanding and enables previously impossible discoveries -- to equity, diversity, and inclusion (EDI) in science research collaborations. With support from the National Science Foundation INCLUDES program, the MDN focuses on increasing EDI by sharing knowledge, experiences, training, and resources among representatives from multimessenger science collaborations. Representatives to the MDN become engagement leads in their collaboration, extending the reach of the community of practice. An overview of the MDN structure, lessons learned, and how to join are presented.

We introduce a rigorous and general framework to study systematically self-gravitating elastic materials within general relativity, and apply it to investigate the existence and viability, including radial stability, of spherically symmetric elastic stars. We present the mass-radius ($M-R$) diagram for various families of models, showing that elasticity contributes to increase the maximum mass and the compactness up to a ${\cal O}(10\%)$ factor, thus supporting compact stars with mass well above two solar masses. Some of these elastic stars can reach compactness as high as $GM/(c^2R)\approx 0.35$ while remaining stable under radial perturbations and satisfying all energy conditions and subluminal wave propagation, thus being physically viable models of stars with a light ring. We provide numerical evidence that radial instability occurs for central densities larger than that corresponding to the maximum mass, as in the perfect-fluid case. Elasticity may be a key ingredient to build consistent models of exotic ultracompact objects and black-hole mimickers, and can also be relevant for a more accurate modelling of the interior of neutron stars.

We report the analysis of data taken during a pilot run in 2018 to study the feasibility of nuclear fragmentation measurements with the NA61/SHINE experiment at the CERN SPS. These nuclear reactions are important for the interpretation of secondary cosmic-ray nuclei production (Li, Be, and B) in the Galaxy. The pilot data were taken with $^{12}\text{C}$ projectiles at a beam momentum of 13.5 A GeV/c and two fixed targets, polyethylene (C$_2$H$_4$) and graphite. The specific focus here is the measurement of total Boron ($^{10}\text{B}$ and $^{11}\text{B}$) production cross section in $\text{C+p}$ interactions at 13.5 A GeV/c. The cosmic-ray nucleus $^{11}\text{C}$ is termed a `Ghost nucleus' on account of its short lifetime compared to the usual cosmic-ray diffusion time in the Galaxy and it ultimately decays to Boron as, $^{11}\text{C} \to ^{11}\text{B} + \beta^+$. Therefore, precise knowledge of the production cross section of $^{11}\text{C}$ is very relevant for the understanding of Boron production in the Galaxy. We present a preliminary measurement of the fragmentation cross section of $\text{C+p}\to ^{11}\text{C}$, which, together with our previously reported B-production cross section, provides a new constraint on Boron production in the Galaxy in the high-energy range relevant for modern space based cosmic-ray experiments like AMS-02.

Solar eclipses provide an excellent opportunity to study the effects of a sudden localized change in photoionization flux in the Earth's ionosphere and its consequent repercussion in the Geomagnetic field. We have focused on a subset of the data available from the North American 2017 eclipse in order to study VTEC measurements from GNSS data and geomagnetic field estimations from INTERMAGNET observatories near the eclipse path. Our simultaneous analysis of both datasets allowed us to quantify the ionosphere and magnetic field reaction to the eclipse event with which allowed us to compare how differently these take place in time. We found that studying the behaviour of VTEC differences with respect to reference values provides better insight of the actual eclipse effect and were able to characterize the dependence of parameters such as time delay of maximum depletion and recovery phase. We were also able to test models that link the ionospheric variations in a quantitative manner. Total electron content depletion measured from GNSS were fed into an approximation of Ashour-Chapman model at the locations of geomagnetic observatories and its predictions match the behaviour of magnetic field components in time and magnitude strikingly accurately.