Papers for Friday, Jul 23 2021

Context: Late on 2013 August 19, STEREO-A, STEREO-B, MESSENGER, Mars Odyssey, and the L1 spacecraft, spanning a longitudinal range of 222{\deg} in the ecliptic plane, observed an energetic particle flux increase. The widespread solar energetic particle (SEP) event was associated with a coronal mass ejection (CME) that came from a region located near the far-side central meridian from Earth's perspective. The CME erupted in two stages, and was accompanied by a late M-class flare observed as a post-eruptive arcade, persisting low-frequency (interplanetary) type II and groups of shock-accelerated type III radio bursts, all of them making this SEP event unusual. Aims: There are two main objectives of this study, disentangling the reasons for the different intensity-time profiles observed by the spacecraft, especially at MESSENGER and STEREO-A locations, longitudinally separated by only 15{\deg}, and unravelling the single solar source related with the widespread SEP event. Results: The solar source associated with the widespread SEP event is the shock driven by the CME, as the flare observed as a post-eruptive arcade is too late to explain the estimated particle onset. The different intensity-time profiles observed by STEREO-A, located at 0.97 au, and MESSENGER, at 0.33 au, can be interpreted as enhanced particle scattering beyond Mercury's orbit. The longitudinal extent of the shock does not explain by itself the wide spread of particles in the heliosphere. The particle increase observed at L1 may be attributed to cross-field diffusion transport, and this is also the case for STEREO-B, at least until the spacecraft is eventually magnetically connected to the shock when it reaches ~0.6 au.

We present the results of a spectroscopic monitoring program of the Pleiades region aimed at completing the census of spectroscopic binaries in the cluster, extending it to longer periods than previously reachable. We gathered 6104 spectra of 377 stars between 1981 and 2021, and merged our radial velocities with 1151 measurements from an independent survey by others started three years earlier. With the combined data spanning more than 43 yr we have determined orbits for some 30 new binary and multiple systems, more than doubling the number previously known in the Pleiades. The longest period is 36.5 yr. A dozen additional objects display long-term trends in their velocities, implying even longer periods. We examine the collection of orbital elements for cluster members, and find that the shape of the incompleteness-corrected distribution of periods (up to $10^4$ days) is similar to that of solar-type binaries in the field, while that of the eccentricities is different. The mass-ratio distribution is consistent with being flat. The binary frequency in the Pleiades for periods up to $10^4$ days is $25 \pm 3$%, after corrections for undetected binaries, which is nearly double that of the field up to the same period. The total binary frequency including known astrometric binaries is at least 57%. We estimate the internal radial velocity dispersion in the cluster to be $0.48 \pm 0.04$ km s$^{-1}$. We revisit the determination of the tidal circularization period, and confirm its value to be $7.2 \pm 1.0$ days, with an improved precision compared to an earlier estimate.

Cosmic-ray positrons have long been considered a powerful probe of dark matter annihilation. In particular, myriad studies of the unexpected rise in the positron fraction have debated its dark matter or pulsar origins. In this paper, we instead examine the potential for extremely precise positron measurements by AMS-02 to probe hard leptophilic dark matter candidates that do not have spectral features similar to the bulk of the observed positron excess. Utilizing a detailed cosmic-ray propagation model that includes a primary positron flux generated by Galactic pulsars in addition to a secondary component constrained by Helium and proton measurements, we produce a robust fit to the local positron flux and spectrum. We find no evidence for a spectral bump correlated with leptophilic dark matter, and set strong constraints on the dark matter annihilation cross-section that fall below the thermal annihilation cross-section for dark matter masses below 60 GeV and 380 GeV for annihilation into $\tau^+\tau^-$ and $e^+e^-$, respectively, in our default model.

Using Subaru Hyper Suprime-Cam (HSC) year 1 data, we perform the first $k$-cut cosmic shear analysis constraining both $\Lambda$CDM and $f(R)$ Hu-Sawicki modified gravity. To generate the $f(R)$ cosmic shear theory vector, we use the matter power spectrum emulator trained on COLA (COmoving Lagrangian Acceleration) simulations. The $k$-cut method is used to significantly down-weight sensitivity to small scale ($k > 1 \ h {\rm Mpc }^{-1}$) modes in the matter power spectrum where the emulator is less accurate, while simultaneously ensuring our results are robust to baryonic feedback model uncertainty. We have also developed a test to ensure that the effects of poorly modeled small scales are nulled as intended. For $\Lambda$CDM we find $S_8 = \sigma_8 (\Omega_m / 0.3) ^ {0.5} = 0.789 ^{+0.039}_{-0.022}$, while the constraints on the $f(R)$ modified gravity parameters are prior dominated. In the future, the $k$-cut method could be used to constrain a large number of theories of gravity where computational limitations make it infeasible to model the matter power spectrum down to extremely small scales.

We present the first high signal-to-noise broadband X-ray spectrum of the radio-quiet type-2 Seyfert ESO 033-G002, combining data from $XMM$-$Newton$ and $NuSTAR$. The nuclear X-ray spectrum is complex, showing evidence for both neutral and ionised absorption, as well as reflection from both the accretion disc and more distant material, but our broadband coverage allows us to disentangle all of these different components. The total neutral column during this epoch is $N_{\rm{H}} \sim 5-6 \times 10^{22}$ cm$^{-2}$, consistent with the optical classification of ESO 033-G002 as a type-2 Seyfert but not so large as to prevent us from robustly determining the properties of the innermost accretion flow. The ionised absorption - dominated by lines from Fe XXV and Fe XXVI - reveals a moderately rapid outflow ($v_{\rm{out}} \sim 5400$ km s$^{-1}$) which has a column comparable to the neutral absorption. We find the disc reflection from the innermost regions to be extreme, with a reflection fraction of $R_{\rm{frac}} \sim 5$. This requires strong gravitational lightbending and, in turn, both an extremely compact corona (within $\sim$2 $R_{\rm{G}}$ of the black hole) and a rapidly rotating black hole ($a^* > 0.96$). Despite this tight size constraint, with a temperature of $kT_{\rm{e}} = 40-70$ keV the X-ray corona in ESO 033-G002 appears similar to other AGN in terms of its placement in the compactness-temperature plane, consistent with sitting close to the limit determined by runaway pair production. Finally, combining X-ray spectroscopy, timing and updated optical spectroscopy, we also estimate the mass of the black hole to be $\log[M_{\rm{BH}} / M_{\odot}] \sim 7.0 - 7.5$.

Perhaps the most exciting promise of the Rubin Observatory Legacy Survey of Space and Time (LSST) is its capability to discover phenomena never before seen or predicted from theory: true astrophysical novelties, but the ability of LSST to make these discoveries will depend on the survey strategy. Evaluating candidate strategies for true novelties is a challenge both practically and conceptually: unlike traditional astrophysical tracers like supernovae or exoplanets, for anomalous objects the template signal is by definition unknown. We present our approach to solve this problem, by assessing survey completeness in a phase space defined by object color, flux (and their evolution), and considering the volume explored by integrating metrics within this space with the observation depth, survey footprint, and stellar density. With these metrics, we explore recent simulations of the Rubin LSST observing strategy across the entire observed footprint and in specific regions in the Local Volume: the Galactic Plane and Magellanic Clouds. Under our metrics, observing strategies with greater diversity of exposures and time gaps tend to be more sensitive to genuinely new phenomena, particularly over time-gap ranges left relatively unexplored by previous surveys. To assist the community, we have made all the tools developed publicly available. Extension of the scheme to include proper motions and the detection of associations or populations of interest, will be communicated in paper II of this series. This paper was written with the support of the Vera C. Rubin LSST Transients and Variable Stars and Stars, Milky Way, Local Volume Science Collaborations.

Ultraviolet spectra of protoplanetary disks trace distributions of warm gas at radii where rocky planets form. We combine HST-COS observations of H2 and CO emission from 12 classical T Tauri stars to more extensively map inner disk surface layers, where gas temperature distributions allow radially stratified fluorescence from the two species. We calculate empirical emitting radii for each species under the assumption that the line widths are entirely set by Keplerian broadening, demonstrating that the CO fluorescence originates further from the stars (r ~ 20 AU) than the H2 (r ~ 0.8 AU). This is supported by 2-D radiative transfer models, which show that the peak and outer radii of the CO flux distributions generally extend further into the outer disk than the H2. These results also indicate that additional sources of LyA photons remain unaccounted for, requiring more complex models to fully reproduce the molecular gas emission. As a first step, we confirm that the morphologies of the UV-CO bands and LyA radiation fields are significantly correlated and discover that both trace the degree of dust disk evolution. The UV tracers appear to follow the same sequence of disk evolution as forbidden line emission from jets and winds, as the observed LyA profiles transition between dominant red wing and dominant blue wing shapes when the high-velocity optical emission disappears. Our results suggest a scenario where UV radiation fields, disk winds and jets, and molecular gas evolve in harmony with the dust disks throughout their lifetimes.

We investigate the application of Hybrid Effective Field Theory (HEFT) -- which combines a Lagrangian bias expansion with subsequent particle dynamics from $N$-body simulations -- to the modeling of $k$-Nearest Neighbor Cumulative Distribution Functions ($k{\rm NN}$-${\rm CDF}$s) of biased tracers of the cosmological matter field. The $k{\rm NN}$-${\rm CDF}$s are sensitive to all higher order connected $N$-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the $k{\rm NN}$-${\rm CDF}$s of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how $k{\rm NN}$-${\rm CDF}$ statistics of a set of biased tracers, such as halos or galaxies, of the cosmological matter field can be modeled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of halos, we show that both the two-point clustering \textit{and} the $k{\rm NN}$-${\rm CDF}$s can be well-fit on quasi-linear scales ($\gtrsim 20 h^{-1}{\rm Mpc}$) by the second-order HEFT formalism with the \textit{same values} of the bias parameters, implying that joint modeling of the two is possible. Finally, using a Fisher matrix analysis, we show that including $k{\rm NN}$-${\rm CDF}$ measurements over the range of allowed scales in the HEFT framework can improve the constraints on $\sigma_8$ by roughly a factor of $3$, compared to the case where only two-point measurements are considered. Combining the statistical power of $k{\rm NN}$ measurements with the modeling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

We present the discovery of rare active galactic nuclei (AGNs) in nearby (z < 0.05) compact elliptical galaxies (cEs) located in isolated environments. Using spectroscopic data from the Sloan Digital Sky Survey (SDSS) Data Release 12, four AGNs were identified based on the optical emission-line diagnostic diagram. SDSS optical spectra of AGNs show the presence of distinct narrow-line emissions. Utilizing the black hole (BH) mass-stellar velocity dispersion scaling relation and the correlation between the narrow L([OIII])/L(H\b{eta}) line ratio and the width of the broad H{\alpha} emission line, we estimated the BH masses of the cEs to be in the range of 7 x 10^5 - 8 x 10^7 Msun. The observed surface brightness profiles of the cEs were fitted with a double Sersic function using the Dark Energy Camera Legacy Survey r-band imaging data. Assuming the inner component as the bulge, the K-band bulge luminosity was also estimated from the corresponding Two Micron All Sky Survey images. We found that our cEs follow the observed BH mass-stellar velocity dispersion and BH mass-bulge luminosity scaling relations, albeit there was a large uncertainty in the derived BH mass of one cE. In view of the observational properties of BHs and those of the stellar populations of cEs, we discuss the proposition that cEs in isolated environments are bona fide low-mass early-type galaxies (i.e., a nature origin).

Despite the remarkable success of the $\Lambda$Cold Dark Matter ($\Lambda$CDM) cosmological model, a growing discrepancy has emerged (currently measured at the level of $\sim 4-6 \sigma$) between the value of the Hubble constant $H_0$ measured using the local distance ladder and the value inferred using the cosmic microwave background and galaxy surveys. While a vast array of $\Lambda$CDM extensions have been proposed to explain these discordant observations, understanding the (relative) success of these models in resolving the tension has proven difficult -- this is a direct consequence of the fact that each model has been subjected to differing, and typically incomplete, compilations of cosmological data. In this review, we attempt to make a systematic comparison of sixteen different models which have been proposed to resolve the $H_0$ tension (spanning both early- and late-Universe solutions), and quantify the relative success of each using a series of metrics and a vast array of data combinations. Owing to the timely appearance of this article, we refer to this contest as the ''$H_0$ Olympics''; the goal being to identify which of the proposed solutions, and more broadly which underlying mechanisms, are most likely to be responsible for explaining the observed discrepancy (should unaccounted for systematics not be the culprit). This work also establishes a foundation of tests which will allow the success of novel proposals to be meaningful ''benchmarked''.

The evolution of the universe prior to Big Bang Nucleosynthesis could have gone through a phase of early matter domination (EMD) which enhanced the growth of small-scale dark matter structure. If EMD was long enough, self-gravitating objects formed prior to reheating. We study the evolution of these dense early halos (EHs) through reheating. At the end of EMD, EHs undergo rapid expansion and eventually eject their matter. We find that this process washes out structure on scales much larger than naively expected from the size of the original halos. We compute the density profiles of the EH remnants and use them to construct late-time power spectra that include these non-linear effects. EH dynamics limits the maximum enhancement that can be generated by EMD in a way that is independent of the dark matter microphysics. We evolve an extrapolated $\Lambda$CDM power spectrum to estimate the properties of microhalos that would form after matter-radiation equality. Surprisingly, cosmologies with a short period of EMD lead to an earlier onset of microhalo formation compared to those with a long period of EMD. In either case, dark matter structure formation begins much earlier than in the standard cosmology, with most DM bound in microhalos.

Strongly lensed quasars can provide measurements of the Hubble constant ($H_{0}$) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand on the previous analysis of three time-delay lenses with AO imaging (RXJ1131-1231, HE0435-1223, and PG1115+080), and perform a joint analysis of J0924+0219 by using AO imaging from the Keck Telescope, obtained as part of the SHARP (Strong lensing at High Angular Resolution Program) AO effort, with HST imaging to constrain the mass distribution of the lensing galaxy. Under the assumption of a flat $\Lambda$CDM model with fixed $\Omega_{\rm m}=0.3$, we show that by marginalizing over two different kinds of mass models (power-law and composite models) and their transformed mass profiles via a mass-sheet transformation, we obtain $\Delta t_{\rm BA}h\hat{\sigma}_{v}^{-2}=6.89\substack{+0.8\\-0.7}$ days, $\Delta t_{\rm CA}h\hat{\sigma}_{v}^{-2}=10.7\substack{+1.6\\-1.2}$ days, and $\Delta t_{\rm DA}h\hat{\sigma}_{v}^{-2}=7.70\substack{+1.0\\-0.9}$ days, where $h=H_{0}/100~\rm km\,s^{-1}\,Mpc^{-1}$ is the dimensionless Hubble constant and $\hat{\sigma}_{v}=\sigma^{\rm ob}_{v}/(280~\rm km\,s^{-1})$ is the scaled dimensionless velocity dispersion. Future measurements of time delays with 10% uncertainty and velocity dispersion with 5% uncertainty would yield a $H_0$ constraint of $\sim15$% precision.

We report the possible discovery of a new stellar system (YMCA-1), identified during a search for small scale overdensities in the photometric data of the YMCA survey. The object's projected position lies on the periphery of the Large Magellanic Cloud about $13^\circ$ apart from its center. The most likely interpretation of its color-magnitude diagram, as well as of its integrated properties, is that YMCA-1 may be an old and remote star cluster of the Milky Way at a distance of 100 kpc from the Galactic center. If this scenario could be confirmed, then the cluster would be significantly fainter and more compact than most of the known star clusters residing in the extreme outskirts of the Galactic halo, but quite similar to Laevens~3. However, much deeper photometry is needed to firmly establish the actual nature of the cluster and the distance to the system.

An asteroid pair can be described as two asteroids with highly similar heliocentric orbits that are genetically related but not gravitationally bound. They can be produced by asteroid collisions or rotational fission. Although over 200 asteroid pairs are known, many more are remaining to be identified, especially among the newly discovered asteroids. The purpose of our work is to find new asteroid pairs in the inner part of the main belt with a new pipeline for asteroid pair search, as well as to validate the pipeline on a sample of known asteroid pairs. Initially, we select pair candidates in the five-dimensional space of osculating orbital elements. Then the candidates are confirmed using numerical modeling with the backtrack integration of their orbits including the perturbations from the largest main-belt asteroids, as well as the influence of the non-gravitational Yarkovsky effect. We performed a survey of the inner part of the main belt and found 10 new probable asteroid pairs. Their estimated formation ages lie between 30 and 400 kyr. In addition, our pipeline was tested on a sample of 17 known pairs, and our age estimates agreed with the ones indicated in literature in most of the cases.

Studying materials released from Jupiter-family comets (JFCs) as seen in their inner com\ae, the envelope of gas and dust forming as the comet approaches the Sun, provides an improved understanding of their origin and evolutionary history. As part of a coordinated, multi-wavelength observing campaign, we observed comet 45P/Honda-Mrkos-Pajdu\v{s}\'{a}kov\'{a} during its close approach to Earth in February 2017. Narrowband observations were taken using the Bok 90'' telescope at KPNO on February 16 and 17 UT, revealing gas and dust structures. We observed different jet directions for different volatile species, implying source region heterogeneity, consistent with other ground-based as well as in situ observations of comet nuclei. A repeating feature visible in CN and C$_2$ images on February 16, was recovered on February 17 with an interval of $7.6\pm0.1$ hours, consistent with the rotation period of the comet derived from Arecibo Observatory radar observations. The repeating feature's projected gas velocity away from the nucleus is 0.8 km/s, with an expansion velocity as 0.5 km/s. The amount of CN material released in one cycle has a lower limit of 11 kg, depending on composition, a quantity small enough to be produced by repeated exposure of nucleus ices to sunlight. This repeating CN jet forming within 400 km of the nucleus may be typical of inner coma behavior in JFCs. Similar repeating CN features could exist and be common in other observed comets, but obscured by other processes and daughter product species as viewed from distances further than the scale length of CN molecules.

Stellar bow shocks are observed in a variety of interstellar environments and are shaped by the conditions of gas in the interstellar medium (ISM). In situ measurements of turbulent density fluctuations near stellar bow shocks are only achievable with a few observational probes, including H$\alpha$ emitting bow shocks and the Voyager Interstellar Mission (VIM). In this paper, we examine density variations around the Guitar Nebula, an H$\alpha$ bow shock associated with PSR B2224$+$65, in tandem with density variations probed by VIM near the boundary of the solar wind and ISM. High-resolution Hubble Space Telescope observations of the Guitar Nebula taken between 1994 and 2006 trace density variations over scales from 100s to 1000s of au, while VIM density measurements made with the Voyager 1 Plasma Wave System constrain variations from 1000s of meters to 10s of au. The power spectrum of density fluctuations constrains the amplitude of the turbulence wavenumber spectrum near the Guitar Nebula to ${\rm log}_{10}C_{\rm n}^2 = -0.8\pm0.2$ m$^{-20/3}$ and for the very local ISM probed by Voyager ${\rm log}_{10}C_{\rm n}^2 = -1.57\pm0.02$ m$^{-20/3}$. Spectral amplitudes obtained from multi-epoch observations of four other H$\alpha$ bow shocks also show significant enhancements in $C_{\rm n}^2$ from values that are considered typical for the diffuse, warm ionized medium, suggesting that density fluctuations near these bow shocks may be amplified by shock interactions with the surrounding medium, or by selection effects that favor H$\alpha$ emission from bow shocks embedded in denser media.

The large scale structure of the Solar System has been shaped by a transient dynamical instability that may have been triggered by the interaction of the giants planets with a massive primordial disk of icy debris. In this work, we investigate the conditions under which this primordial disk could have coalesced into planets using analytic and numerical calculations. In particular, we perform numerical simulations of the Solar System's early dynamical evolution that account for the viscous stirring and collisional damping within the disk. We demonstrate that if collisional damping would have been sufficient to maintain a temperate velocity dispersion, Earth mass trans-Neptunian planets could have emerged within a timescale of 10 Myr. Therefore, our results favor a scenario wherein the dynamical instability of the outer Solar System began immediately upon the dissipation of the gaseous nebula to avoid the overproduction of Earth mass planets in the outer Solar System.

We investigate equilibrium chemistry between a metal-core, a silicate-mantle, and a hydrogen-rich atmosphere (reactive core model) using 18 independent reactions among 25 phase components for sub-Neptune-like exoplanets. We find hydrogen and oxygen typically comprise 1-2% and ~10% by weight of the metal-core, respectively, leading to under-dense cores and thereby offering a possible alternative explanation for the densities of the Trappist-1 planets. In addition, hydrogen occurs at about 0.1% by mass in the silicate mantle, setting a maximum limit to the hydrogen-budget for out-gassing by future super-Earths. The total hydrogen-budget of most sub-Neptunes can be, to first order, well estimated from their atmospheres alone, as more than ~93% of all H resides in their atmospheres. However, reactions with the magma ocean produce significant amounts of SiO and H_2O in the atmospheres which increase the mean molecular weight averaged over the whole atmosphere, by about a factor of two, to ~4 amu. We also investigated the case where metal is excluded from the equilibrium chemistry (unreactive core model). In this case, we find most noticeably that, as the hydrogen mass fraction is reduced from 2% to 1%, the atmosphere becomes water dominated and large fractions of H are absorbed by the magma. As water dominated atmospheres appear inconsistent with observations, we conclude that either the unreactive core model does not apply to sub-Neptunes and that their evolution is better described by a reactive core, or that in-gassing of hydrogen into the mantle is much less efficient than permitted by equilibrium chemistry.

The optical sensors of the IceCube Neutrino Observatory are attached on vertical strings of cables. They were frozen into the ice in the deployment holes made by hot water drill. This hole ice, to the best of our knowledge, consists of a bubbly central column, with the remainder of the re-frozen volume being optically clear. The bubbly ice often blocks one or several of the calibration LEDs in every optical sensor and significantly distorts the angular profile of the calibration light pulses. It also affects the sensors' response to in-coming photons at different locations and directions. We present our modeling of the hole ice optical properties as well as optical sensor location and orientation within the hole ice. The shadowing effects of cable string and possible optical sensor tilt away from the nominal vertical alignment are also discussed.

The Chinese CubeSat Mission, Gamma Ray Integrated Detectors (GRID), recently detected its first gamma-ray burst, GRB 210121A, which was jointly observed by the Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM). This burst is confirmed by several other missions, including Fermi and Insight-HXMT. We combined multi-mission observational data and performed a comprehensive analysis of the burst's temporal and spectral properties. Our results show that the burst is special in its high peak energy, thermal-like low energy indices, and large fluence. By putting it to the $E_{\rm p}$-$E_{\rm\gamma, iso}$ relation diagram with assumed distance, we found this burst can be constrained at the redshift range of [0.3,3.0]. The thermal spectral component is also confirmed by the direct fit of the physical models to the observed spectra. Interestingly, the physical photosphere model also constrained a redshift of $z\sim$ 0.3 for this burst, which help us to identify a host galaxy candidate at such a distance within the location error box. Assuming the host galaxy is real, we found the burst can be best explained by the photosphere emission of a typical fireball with a radius of $r_0\sim 6.7\times 10^7$ cm.

Altair is the fastest rotating star at less than 10 parsecs from the Sun. Its precise modelling is a landmark for our understanding of stellar evolution with fast rotation, and all observational constraints are most welcome to better determine the fundamental parameters of this star. We wish to improve the seismic spectrum of Altair and confirm the $\delta$-Scuti nature of this star. We used the photometric data collected by the Microvariability and Oscillations of STars (MOST) satellite in the form of a series of Fabry images to derive Altair light curves at four epochs, namely in 2007, 2011, 2012, and 2013. We first confirm the presence of $\delta$-Scuti oscillations in the light curves of Altair. We extend the precision of some eigenfrequencies and add new ones to the spectrum of Altair, which now has 15 detected eigenmodes. The rotation period, which is expected at $\sim$7h46min from models reproducing interferometric data, seems to appear in the 2012 data set, but it still needs confirmation. Finally, Altair modal oscillations show noticeable amplitude variations on a timescale of 10 to 15 days, which may be the signature of a coupling between oscillations and thermal convection in the layer where the kappa-mechanism is operating.The Altair oscillation spectrum does not contain a large number of excited eigenmodes, which is similar to the fast rotating star HD220811. This supports the idea that fast rotation hinders the excitation of eigenmodes as already pointed out by theoretical investigations.

We present a new instability driven by a combination of coagulation and radial drift of dust particles. We refer to this instability as ``coagulation instability" and regard it as a promising mechanism to concentrate dust particles and assist planetesimal formation in the very early stages of disk evolution. Because of dust-density dependence of collisional coagulation efficiency, dust particles efficiently (inefficiently) grow in a region of positive (negative) dust density perturbations, which lead to a small radial variation of dust sizes and as a result radial velocity perturbations. The resultant velocity perturbations lead to dust concentration and amplify dust density perturbations. This positive feedback makes a disk unstable. The growth timescale of coagulation instability is a few tens of orbital periods even when dust-to-gas mass ratio is of the order of $10^{-3}$. In a protoplanetary disk, radial drift and coagulation of dust particles tend to result in dust depletion. The present instability locally concentrates dust particles even in such a dust-depleted region. The resulting concentration provides preferable sites for dust-gas instabilities to develop, which leads to further concentration. Dust diffusion and aerodynamical feedback tend to stabilize short-wavelength modes, but do not completely suppress the growth of coagulation instability. Therefore, coagulation instability is expected to play an important role in setting up the next stage for other instabilities to further develop toward planetesimal formation, such as streaming instability or secular gravitational instability.

Magnetic monopoles are hypothetical particles that carry magnetic charge. Depending on their velocity, different light production mechanisms exist to facilitate detection. In this work, a previously unused light production mechanism, luminescence of ice, is introduced. This light production mechanism is nearly independent of the velocity of the incident magnetic monopole and becomes the only viable light production mechanism in the low relativistic regime (0.1-0.55c). An analysis in the low relativistic regime searching for magnetic monopoles in seven years of IceCube data is presented. While no magnetic monopole detection can be claimed, a new flux limit in the low relativistic regime is presented, superseding the previous best flux limit by 2 orders of magnitude.

4C+28.07 is a $\gamma$-ray Flat Spectrum Radio Quasar (FSRQ) type source. It is often monitored at different frequencies, though long-term multi-wavelength data of this source have not been modelled in detail before. We have analyzed $\sim 12$ years (Aug, 2008 - May, 2020) of Fermi-LAT data with a binning of 10 days time scale and observed three distinctive flaring states. Each flaring state consists of different phases of activity, namely, pre-flare, flare \& post-flare regions. $\gamma$-ray spectral analysis of these different activity phases has been performed and the best fit model for its spectra is found to be a Log-parabola model. We have also studied the correlation of simultaneous $\gamma$-ray light curves with the optical \& radio counterparts in these flaring states and report the DCF with a 95\% significance level. A large time delay is found between radio and gamma-ray data for two flares, indicating two zones of emission. We have fitted the multi-wavelength data with a two-zone leptonic model. In our two-zone leptonic model, the maximum required power in the jet is 9.64 $\times$ 10$^{46}$ erg sec$^{-1}$, which is lower than its Eddington luminosity $2.29\times 10^{47}$ erg sec$^{-1}$.

The solar gravitational moments $J_{2n}$ are important astronomical quantities whose precise determination is relevant for solar physics, gravitational theory and high precision astrometry and celestial mechanics. Accordingly, we propose in the present work to calculate new values of $J_{2n}$ (for $n$=1,2,3,4 and 5) using recent two-dimensional rotation rates inferred from the high resolution SDO/HMI helioseismic data spanning the whole solar activity cycle 24. To this aim, a general integral equation relating $J_{2n}$ to the solar internal density and rotation is derived from the structure equations governing the equilibrium of slowly rotating stars. For comparison purpose, the calculations are also performed using rotation rates obtained from a recently improved analysis of SoHO/MDI heliseismic data for solar cycle 23. In agreement with earlier findings, the results confirmed the sensitivity of high order moments ($n>1$) to the radial and latitudinal distribution of rotation in the convective zone. The computed value of the quadrupole moment $J_{2}$ ($n=1$) is in accordance with recent measurements of the precession of Mercury's perihelion deduced from high precision ranging data of the MESSENGER spacecraft. The theoretical estimate of the related solar oblateness $\Delta_{\odot}$ is consistent with the most accurate space-based determinations, particularly the one from RHESSI/SAS.

HCN and its isomer HNC play an important role in molecular cloud chemistry and the formation of more complex molecules. We investigate here the impact of protostellar shocks on the HCN and HNC abundances from high-sensitivity IRAM 30m observations of the prototypical shock region L1157-B1 and the envelope of the associated Class 0 protostar, as a proxy for the pre-shock gas. The isotopologues H$^{12}$CN, HN$^{12}$C, H$^{13}$CN, HN$^{13}$C, HC$^{15}$N, H$^{15}$NC, DCN and DNC were all detected towards both regions. Abundances and excitation conditions were obtained from radiative transfer analysis of molecular line emission under the assumption of Local Thermodynamical Equilibrium. In the pre-shock gas, the abundances of the HCN and HNC isotopologues are similar to those encountered in dark clouds, with a HCN/HNC abundance ratio $\approx 1$ for all isotopologues. A strong D-enrichment (D/H$\approx 0.06$) is measured in the pre-shock gas. There is no evidence of $^{15}$N fractionation neither in the quiescent nor in the shocked gas. At the passage of the shock, the HCN and HNC abundances increase in the gas phase in different manners so that the HCN/HNC relative abundance ratio increases by a factor 20. The gas-grain chemical and shock model UCLCHEM allows us to reproduce the observed trends for a C-type shock with pre-shock density $n$(H)= $10^5$cm$^{-3}$ and shock velocity $V_s= 40$km/s. We conclude that the HCN/HNC variations across the shock are mainly caused by the sputtering of the grain mantle material in relation with the history of the grain ices.

Slow magnetoacoustic waves observed in the solar corona are used as seismological probes of plasma parameters. It has been shown that dispersion properties of such waves can vary significantly under the influence of the wave-induced thermal misbalance. In the current research, we study the effect of misbalance on waves inside the magnetic-flux tube under the second-order thin-flux-tube approximation. Using the parameters of active region fan coronal loops, we calculated wave properties such as the phase speed and decrement. It is shown that neglecting thermal misbalance may be the reason for the substantial divergence between seismological and spectrometric estimations of plasma parameters. We also show that the frequency dependence of phase speed is affected by two features, namely the geometric dispersion and the dispersion caused by the thermal misbalance. In contrast to the phase speed, the wave decrement primarily is affected by the thermal misbalance only. The dependencies of the phase speed and decrement of the slow wave on magnetic field and tube cross-section are also analyzed.

General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on 'conserved' variables which then need to be converted back into the fundamental ('primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the evolution of a BH-accrection disk system, show that RePrimAnd is able to support highly magnetized, low density environments, even for magnetizations as high as $10^4$, for which the previously used recovery scheme fails.

The generation of large-scale magnetic fields ($\overline{\mathbf{B}}$) in astrophysical systems is driven by the mean turbulent electromotive force. This can depend non-locally on $\overline{\mathbf{B}}$ through a convolution kernel $K_{ij}$. In a new approach to find $K_{ij}$, we directly fit the data from a galactic dynamo simulation using singular value decomposition. We calculate the usual turbulent transport coefficients as moments of $K_{ij}$, show the importance of including non-locality over eddy length scales to fully capture their amplitudes and that their higher order corrections are small.

The remarkable results by the Event Horizon Telescope collaboration concerning the emission from M87* and, more recently, its polarization properties, require an increasingly accurate modeling of the plasma flows around the accreting black hole. Radiatively inefficient sources such as M87* and Sgr A* are typically modeled with the SANE (standard and normal evolution) paradigm, if the accretion dynamics is smooth, or with the MAD (magnetically arrested disk) paradigm, if the black hole's magnetosphere reacts by halting the accretion sporadically, resulting in a highly dynamical process. While the recent polarization studies seem to favor MAD models, this may not be true for all sources, and SANE accretion surely still deserves attention. In this work, we investigate the possibility of reaching the typical degree of magnetization and other accretion properties expected for SANE disks by resorting to the mean-field dynamo process in axisymmetric GRMHD simulations, which are supposed to mimic the amplifying action of an unresolved magnetorotational instability-driven turbulence. We show that it is possible to reproduce the main diagnostics present in the literature by starting from very unfavorable initial configurations, such as a purely toroidal magnetic field with negligible magnetization.

We describe the special relativistic extension of the CRONOS code, which has been used for studies of gamma-ray binaries in recent years. The code was designed to be easily adaptable, allowing the user to easily change existing functionalities or introduce new modules tailored to the problem at hand. Numerically, the equations are treated using a finite-volume Godunov scheme on rectangular grids, which currently support Cartesian, spherical, and cylindrical coordinates. The employed reconstruction technique, the approximate Riemann solver and the equation of state can be chosen dynamically by the user. Further, the code was designed with stability and robustness in mind, detecting and mitigating possible failures early on. We demonstrate the code's capabilities on an extensive set of validation problems.

We have observed the Class I protostar TMC-1A in the Taurus molecular cloud using the Submillimeter Array (SMA) and the Atacama Large Millimeter/submillimeter Array (ALMA) in the linearly polarized 1.3 mm continuum emission at angular resolutions of ~3" and ~0.3", respectively. The ALMA observations also include CO, 13CO, and C18O J=2-1 spectral lines. The SMA observations trace magnetic fields on the 1000-au scale, the directions of which are neither parallel nor perpendicular to the outflow direction. Applying the Davis-Chandrasekhar-Fermi method to the SMA polarization angle dispersion, we estimate a field strength in the TMC-1A envelope of 1-5 mG. It is consistent with the field strength needed to reduce the radial infall velocity to the observed value, which is substantially less than the local} free-fall velocity. The ALMA polarization observations consist of two distinct components -- a central component and a north/south component. The central component shows polarization directions in the disk minor axis to be azimuthal, suggesting dust self-scattering in the TMC-1A disk. The north/south component is located along the outflow axis and the polarization directions are aligned with the outflow direction. We discuss possible origins of this polarization structure, including grain alignment by a toroidal magnetic field and mechanical alignment by the gaseous outflow. In addition, we discover a spiral-like residual in the total intensity (Stokes I) for the first time. The C18O emission suggests that material in the spiral-like structure is infalling at a speed that is 20% of the local Keplerian speed.

Previous works on the divergence of first-order mean-motion resonances (MMRs) have studied in detail the extent of the pericentric and apocentric libration zones of adjacent first-order MMRs, highlighting possible bridges between them in the low-eccentricity circular restricted three-body problem. Here, we describe the previous results in the context of periodic orbits and show that the so-called circular family of periodic orbits is the path that can drive the passage between neighbouring resonances under dissipative effects. We illustrate that the circular family can bridge first and higher order resonances while its gaps at first-order MMRs can serve as boundaries that stop transitions between resonances. In particular, for the Sun-asteroid-Jupiter problem, we show that, during the migration of Jupiter in the protoplanetary disc, a system initially evolving below the apocentric branch of a first-order MMR follows the circular family and can either be captured into the pericentric branch of an adjacent first-order MMR if the orbital migration is rapid or in a higher order MMR in case of slow migration. Radial transport via the circular family can be extended to many small body and planetary system configurations undergoing dissipative effects (e.g., tidal dissipation, solar mass-loss and gas drag).

Improved measurement of the Cosmic Microwave Background polarization from Planck allows a detailed study of reionization beyond the average optical depth. The lower value of the optical depth disfavours an early onset and an early completion of reionization in favour of a redsfhit range where different astrophysical probes provide sensible information on the sources of reionization and the status of the intergalactic medium. In this work we extend our previous study in which we constrained reionization by combining three different probes - CMB, UV luminosity density and neutral hydrogen fraction data - in both treatment and data: we first allow variation in the UV source term varying the product of the efficiency of conversion of UV luminosity into ionizing photons and the escape fraction together with the reionization and cosmological parameters, and then we investigate the impact of a less conservative cut for the UV luminosity function. We find that the estimate for the efficiency is consistent within 95% C.L. with the fixed value we considered in our previous results and is mostly constrained by the QHII data. We find that allowing the efficiency to vary does not affect significantly our results for the average optical depth for monotonic reionization histories, recovering $\tau=0.0519_{-0.0008}^{+0.0010}$ at 68% CL , consistent with our previous studies. Using a less conservative cut for the UV luminosity function, we find $\tau=0.0541_{-0.0016}^{+0.0013}$ at 68% CL, due to the faint end of the luminosity function in the data we use, that also prefers a larger contribution from higher redshifts.

The TOPG\"ot project studies a sample of 86 high-mass star-forming regions in different evolutionary stages from starless cores to ultra compact HII regions. The aim of the survey is to analyze different molecular species in a statistically significant sample to study the chemical evolution in high-mass star-forming regions, and identify chemical tracers of the different phases. The sources have been observed with the IRAM 30m telescope in different spectral windows at 1, 2, and 3 mm. In this first paper, we present the sample and analyze the spectral energy distributions (SEDs) of the TOPG\"ot sources to derive physical parameters. We use the MADCUBA software to analyze the emission of methyl cyanide (CH$_3$CN), a well-known tracer of high-mass star formation. The emission of the $\rm{CH_3CN(5_{K}-4_{K})}$ K-transitions has been detected towards 73 sources (85% of the sample), with 12 non-detections and one source not observed in the frequency range of $\rm{CH_3CN(5_{K}-4_{K})}$. The emission of CH$_3$CN has been detected towards all evolutionary stages, with the mean abundances showing a clear increase of an order of magnitude from high-mass starless-cores to later evolutionary stages. We found a conservative abundance upper limit for high-mass starless cores of $X_{\rm CH_3CN}<4.0\times10^{-11}$, and a range in abundance of $4.0\times10^{-11}<X_{\rm CH_3CN}<7.0\times10^{-11}$ for those sources that are likely high-mass starless cores or very early high-mass protostellar objects. In fact, in this range of abundance we have identified five sources previously not classified as being in a very early evolutionary stage. The abundance of $\rm{CH_3CN}$ can thus be used to identify high-mass star-forming regions in early phases of star-formation.

The spectral energy distribution (SED) in the millimetre (mm) to centimetre (cm) range is a useful tool for characterising the dust in protostellar envelopes as well as free-free emission from the protostar and outflow. While many studies have been carried out towards low- and high-mass protostars, little exists so far about solar-type protostars in high-mass star-forming regions, which are likely to be representatives of the conditions where the Solar System was born. We focus here on the OMC-2/3 solar-type protostars, which are bounded by nearby HII regions and which are, therefore, potentially affected by the high-UV illumination. We aim to understand whether the small-scale structure ($\leq$1000 au) and the evolutionary status of these solar-type protostars are affected by the nearby HII regions, as is the case for the large-scale ($\leq$10$^4$ au) gas chemical composition. We used ALMA in the 1.3 mm band (246.2 GHz) to image the continuum of 16 OMC-2/3 solar-type protostars, with an angular resolution of 0.25$''$ (100 au). We completed our data with archival data from the VANDAM survey of Orion Protostars at 333 and 32.9 GHz, respectively, to construct the dust SED, extract several dust parameters and to assess whether free-free emission is contaminating their dust SED in the cm range. From the mm to cm range dust SED, we found low dust emissivity spectral indexes ($\beta < 1$) for the majority of our source sample and free-free emission towards only 5 of the 16 sample sources. We were also able to confirm or correct the evolutionary status of the source sample. Finally, we did not find any dependence of the source dust parameters on their location in the OMC-2/3 filament. Our results show that the small-scale dust properties of the OMC-2/3 protostars are not affected by the high- UV illumination from the nearby HII regions.

We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring.

Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize $\sim$1700 variables from a sample of $\sim$5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability ($Linear, Curved, Periodic$) and stochastic variability ($Burst, Drop, Irregular$). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as $\Delta$W1 or $\Delta$W2 $>1$ mag followed by inspection of candidates) of $\sim$1000 years in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction ($\sim55\%$) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process.

We study a simple model describing thermal radiation spectra from magnetized neutron stars. The model assumes that a star is nearly spherical and isothermal inside and possesses dipole magnetic fields ($B \lesssim 10^{14}$ G) near the surface, which make the surface temperature distribution non--uniform. We assume further that any surface element emits a blackbody (BB) spectrum with a local effective temperature. We show that such thermal spectra (including phase--resolved) are accurately approximated by simple equivalent two--BB (2BB) models. We introduce and study phase--space maps of 2BB parameters and show that these maps can be useful for interpreting neutron star observations, in which 2BB spectral fits have been done.

We have determined the apsidal motion rate of 27 double-lined eclipsing binaries with precise physical parameters. The obtained values, corrected for their relativistic contribution, yield precise empirical parameters of the internal stellar density concentration. The comparison of these results with the predictions based on new theoretical models shows very good agreement. Small deviations are identified but remain within the observational uncertainties and the path for a refined comparison is indicated.

Contrary to expectations from scenarios of black hole growth driven by galaxy interactions and mergers, dual active galactic nuclei (AGN) with kiloparsec separations are rarely observed and are very difficult to identify, in particular at high redshifts (i.e. z>2). Focussing on the recently discovered dual AGN system LBQS 0302-0019 at z=3.29, we seek to identify further group members in its environment and to understand their formation history through deep high-angular-resolution imaging. We present deep Hubble Space Telescope (HST) Wide-field Camera 3 near-infrared imaging of LBQS 0302-0019. In combination with ground-based VLT/HAWK-I imaging, we infer accurate sizes, colours, ages, and stellar masses of companion galaxies. We clearly detect four companion objects close to LBQS 0302-0019 that also have faint signatures in the ground-based images. We constrain light-weighted ages and masses for the two most prominent companions, Jil1 and Jil2, to $t_\star=252_{-109}^{+222}$Myr with $\log(M_\star/[\mathrm{M}_\odot])= 11.2_{-0.1}^{+0.3}$ and $t_{\star}=19_{-14}^{+74}$Myr with $\log(M_\star/[\mathrm{M}_\odot])= 9.4_{-0.4}^{+0.9}$, respectively. The HST data also show that the obscured AGN, previously identified by strong nebular HeII emission, is associated with the young massive companion Jil2. Because very massive stars of the starburst cannot be solely responsible for the HeII emission, we strengthen our initial conclusion that Jil2 has been hosting an AGN. If the young starburst of Jil2 had been accompanied by sustained black hole growth, Jil2 may have contributed HeII-ionising flux to create the large HeII Ly$\alpha$ proximity zone around LBQS 0302-0019. Hence, the duration of the current luminous AGN episode of LBQS 0302-0019 may have been overestimated.

The IceCube neutrino observatory--installed in the Antarctic ice--is the largest neutrino telescope to date. It consists of 5,160 photomultiplier-tubes spread among 86 vertical strings making a total detector volume of more than a cubic kilometer. IceCube detects neutrinos via Cherenkov light emitted by charged relativistic particles produced when a neutrino interacts in or near the detector. The detector is particularly sensitive to high-energy neutrinos of due to its size and photosensor spacing. In this analysis we search for dark matter that annihilates into a metastable mediator that subsequently decays into Standard Model particles. These models yield an enhanced high-energy neutrino flux from dark matter annihilation inside the Sun compared to models without a mediator. Neutrino signals that are produced directly inside the Sun are strongly attenuated at higher energies due to interactions with the solar plasma. In the models considered here, the mediator can escape the Sun before producing any neutrinos, thereby avoiding attenuation. We present the results of an analysis of six years of IceCube data looking for dark matter in the Sun. We consider mediator lifetimes between 1 ms to 10 s and dark matter masses between 200 GeV and 75 TeV.

We report three AstroSat observations of BL Lacertae object OJ 287. The three observations caught it in very different flux states that are connected to different broadband spectral states. These observations trace the source spectral evolution from the end-phase of activity driven by a new, additional HBL like emission component in 2017 to its complete disappearance in 2018 and re-emergence in 2020. The 2017 observation shows a comparatively flatter optical-UV and X-ray spectrum. Supplementing it with the simultaneous NuSTAR monitoring indicates a hardening at the high-energy-end. The 2018 observation shows a harder X-ray spectrum and a sharp decline or cutoff in the optical-UV spectrum, revealed thanks to the Far-UV data from AstroSat. The brightest of all, the 2020 observation shows a hardened optical-UV spectrum and an extremely soft X-ray spectrum, constraining the low-energy peak of spectral energy distribution at UV energies -- a characteristic of HBL blazars. The contemporaneous MeV-GeV spectra from LAT show the well-known OJ 287 spectrum during 2018 but a flatter spectrum during 2017 and a hardening above ~1 GeV during 2020. Modeling broadband SEDs show that 2018 emission can be reproduced with a one-zone leptonic model while 2017 and 2020 observations need a two-zone model, with the additional zone emitting an HBL radiation.

IceCube is a cubic-kilometer scale neutrino telescope located at the geographic South Pole. The detector utilizes the extremely transparent Antarctic ice as a medium for detecting Cherenkov radiation from neutrino interactions. While the optical properties of the glacial ice are generally well modeled and understood, the uncertainties which remain are still the dominant source of systematic uncertainties for many IceCube analyses. A camera and LED system is being built for the IceCube Upgrade that will enable the observation of optical properties throughout the Upgrade array. The SPICEcore hole, a 1.7 km deep ice-core hole located near the IceCube detector, has given the opportunity to test the performance of the camera system ahead of the Upgrade construction. In this contribution, we present the results of the camera and LED system deployment during the 2019/2020 austral summer season as part of a SPICEcore luminescence logger system.

The gamma-ray bursts (GRBs), as astronomical phenomena of cosmological origin, may be lensed during their propagation. The detection of lensed pulses with millisecond-to-second time delays is a clue to the existence of intermediate-mass black holes, which has not been directly observed, yet. The prompt emissions of lensed GRBs are generally considered to have repeated pulses with identical light curves and spectra but different fluxes and slightly offset positions. In this paper, we chose to search within individual GRB to exclude the effect of spatial location. We report our finding of an interesting case GRB 200716C in a systematic search for lensed GRBs with the data of the Fermi/GBM, Swift/BAT and HXMT/HE. The autocorrelation of GRB 200716C light curves in different energy bands satisfy the selection criteria of a lensing candidate. Using the Bayesian inference to quantitatively compare models also results in a Bayesian factor more inclined to the lens model. Besides, the spectral analysis of the two pulses show that they have similar spectral properties and both follow the empirical relation of short GRBs. All these facts suggest that although the GRB 200718C is a long GRB in duration, it may actually be a short GRB undergoing gravitational lensing. If interpreting as a lensed GRB, the redshifted mass $(1+z_{L})M_L$ of the lensing object is estimated to be $4.72^{+2.31}_{-1.06}$ $\times$ $10^5$ $M_{\odot}$ (90$\%$ credibility), indicating that it may be an intermediate-mass black hole.

The non-Gaussianity of inflationary perturbations, as encoded in the bispectrum (or 3-point correlator), has become an important additional way of distinguishing between inflation models, going beyond the linear Gaussian perturbation quantities of the power spectrum. This habilitation thesis provides a review of my work on both the theoretical and the observational aspects of these non-Gaussianities. In the first part a formalism is described, called the long-wavelength formalism, that provides a way to compute the non-Gaussianities in multiple-field inflation. Applications of this formalism to various classes of models, as well as its extensions, are also treated. In the second part an estimator is described, called the binned bispectrum estimator, that allows the extraction of information about non-Gaussianities from data of the cosmic microwave background radiation (CMB). It was in particular one of the three estimators applied to the data of the Planck satellite to provide the currently best constraints on primordial non-Gaussianity. Various extensions of the estimator and results obtained are also discussed.

In this paper, we study how to directly measure the effect of peculiar velocities in the observed angular power spectra. We do this by constructing a new anti-symmetric estimator of Large Scale Structure using different dark matter tracers. We show that the Doppler term is the major component of our estimator and we show that we can measure it with a signal-to-noise ratio up to $\sim 50$ using a futuristic SKAO HI galaxy survey. We demonstrate the utility of this estimator by using it to provide constraints on the Euler equation.

Neutral Hydrogen Intensity Mapping (HI IM) surveys will be a powerful new probe of cosmology. However, strong astrophysical foregrounds contaminate the signal and their coupling with instrumental systematics further increases the data cleaning complexity. In this work, we simulate a realistic single-dish HI IM survey of a $5000$~deg$^2$ patch in the $950 - 1400$ MHz range, with both the MID telescope of the SKA Observatory (SKAO) and MeerKAT, its precursor. We include a state-of-the-art HI simulations and explore different foreground models and instrumental effects such as non-homogeneous thermal noise and beam side-lobes. We perform the first Blind Foreground Subtraction Challenge for HI IM on these synthetic data-cubes, aiming to characterise the performance of available foreground cleaning methods with no prior knowledge of the sky components and noise level. Nine foreground cleaning pipelines joined the Challenge, based on statistical source separation algorithms, blind polynomial fitting, and an astrophysical-informed parametric fit to foregrounds. We devise metrics to compare the pipeline performances quantitatively. In general, they can recover the input maps' 2-point statistics within 20 per cent in the range of scales least affected by the telescope beam. However, spurious artefacts appear in the cleaned maps due to interactions between the foreground structure and the beam side-lobes. We conclude that it is fundamental to develop accurate beam deconvolution algorithms and test data post-processing steps carefully before cleaning. This study was performed as part of SKAO preparatory work by the HI IM Focus Group of the SKA Cosmology Science Working Group.

We compute the real space galaxy power spectrum, including the leading order effects of General Relativity and primordial non-Gaussianity from the $f_{\mathrm{NL}}$ and $g_{\mathrm{NL}}$ parameters. Such contributions come from the one-loop matter power spectrum terms dominant at large scales, and from the factors of the non-linear bias parameter $b_{\mathrm{NL}}$ (akin to the Newtonian $b_{\phi}$). We use our modelling to assess the ability of Stage-IV surveys to constrain primordial non-Gaussianity. In addition, we show how this non-linear bias parameter can effectively renormalize diverging relativistic contributions at large scales.

An observer in relative motion to the Cosmic Microwave Background (CMB) rest frame is sensitive to both aberration and Doppler effects. Both effects introduce similar but non-identical off-diagonal couplings in the spherical harmonic coefficients. The CMB temperature dipole may have additional contributions from an intrinsic component, which in turn produces different aberration and Doppler couplings. Moreover, the standard conversion from intensity measurements into temperature also introduces spurious Doppler-like couplings. In order to learn about the intrinsic dipole it is therefore important to measure both aberration and Doppler couplings in an independent manner while also removing the spurious contributions from unit conversion, which are degenerate with the dipole. Here we present a pipeline to measure the Doppler and aberration signal independently from each other and from the dipole itself. We also consider realistic beaming, noise and mask effects. Our pipeline results in independent and unbiased estimators which have uncertainties only ~20% larger than the simple theoretical expectations. We discuss the achievable precision in each measurement for Planck 2018, and also forecast them for future ground-based experiments with the Simons Observatory and CMB-S4. An alternative pipeline is presented in order to cross-check results and improve robustness.

We study momentum space dispersion formulas in general QFTs and their applications for CFT correlation functions. We show, using two independent methods, that QFT dispersion formulas can be written in terms of causal commutators. The first derivation uses analyticity properties of retarded correlators in momentum space. The second derivation uses the largest time equation and the defining properties of the time-ordered product. At four points we show that the momentum space QFT dispersion formula depends on the same causal double-commutators as the CFT dispersion formula. At $n$-points, the QFT dispersion formula depends on a sum of nested advanced commutators. For CFT four-point functions, we show that the momentum space dispersion formula is equivalent to the CFT dispersion formula, up to possible semi-local terms. We also show that the Polyakov-Regge expansions associated to the momentum space and CFT dispersion formulas are related by a Fourier transform. In the process, we prove that the momentum space conformal blocks of the causal double-commutator are equal to cut Witten diagrams. Finally, by combining the momentum space dispersion formulas with the AdS Cutkosky rules, we find a complete, bulk unitarity method for AdS/CFT correlators in momentum space.

We study the analytic properties of tree-level wavefunction coefficients in quasi-de Sitter space. We focus on theories which spontaneously break dS boost symmetries and can produce significant non-Gaussianities. The corresponding inflationary correlators are (approximately) scale invariant, but are not invariant under the full conformal group. We derive cutting rules and dispersion formulas for the late-time wavefunction coefficients by using factorization and analyticity properties of the dS bulk-to-bulk propagator. This gives a unitarity method which is valid at tree-level for general $n$-point functions and for fields of arbitrary mass. Using the cutting rules and dispersion formulas, we are able to compute $n$-point functions by gluing together lower-point functions. As an application, we study general four-point, scalar exchange diagrams in the EFT of inflation. We show that exchange diagrams constructed from boost-breaking interactions can be written as a finite sum over residues. Finally, we explain how the dS identities used in this work are related by analytic continuation to analogous identities in Anti-de Sitter space.

We study the impact of gravitational waves originating from a first order phase transition on structure formation. To do so, we perform a second order perturbation analysis in the $1+3$ covariant framework and derive a wave equation in which second order, adiabatic density perturbations of the photon-baryon fluid are sourced by the gravitational wave energy density during radiation domination and on sub-horizon scales. The scale on which such waves affect the energy density perturbation spectrum is found to be proportional to the horizon size at the time of the phase transition times its inverse duration. Consequently, structure of the size of galaxies and bigger can only be affected in this way by relatively late phase transitions at $\ge 10^{6}\,\text{s}$. Using cosmic variance as a bound we derive limits on the strength $\alpha$ and the relative duration $(\beta/H_*)^{-1}$ of phase transitions as functions of the time of their occurrence which results in a new exclusion region for the energy density in gravitational waves today. We find that the cosmic variance bound forbids only relative long lasting phase transitions, e.g. $\beta/H_*\lesssim 6.8$ for $t_*\approx 5\times10^{11}\,\text{s}$, which exhibit a substantial amount of supercooling $\alpha>20$ to affect the matter power spectrum.

Within the context of Magnetohydrodynamics (MHD), the properties of a parallel shock do not depend on the field strength, as the field and the fluid are disconnected for such a geometry. However, in the collisionless case, the field can sustain a stable anisotropy in the downstream, triggering a departure from the expected MHD behavior. In a recent work [A. Bret and R. Narayan, J. Plasma Phys. \textbf{84}, 905840604 (2018)], a theoretical model was presented allowing to derive the density ratio of a non-relativistic parallel collisionless shock in an electron/positron plasma, as a function of the field. Here we derive the entropy, pressure and temperature jumps stemming from this model. It is found to offer a transition between a 3D and a 1D downstream for the jumps in density, entropy, parallel temperature and parallel pressure.

We report an investigation of X-ray induced desorption of neutrals, cations and anions from CO ice. The desorption of neutral CO, by far the most abundant, is quantified and discussed within the context of its application to astrochemistry. The desorption of many different cations, including large cations up to the mass limit of the spectrometer, are observed. In contrast, the only desorbing anions detected are O$^-$ and C$^-$. The desorption mechanisms of all these species are discussed with the aid of their photodesorption spectrum. The evolution of the X-ray absorption spectrum shows significant chemical modifications of the ice upon irradiation, which along with the desorption of large cations gives a new insight into X-ray induced photochemistry in CO ice.

We investigate the generation of Primordial Black Holes (PBHs) with the aid of gravitationally increased friction mechanism originated from the NonMinimal field Derivative Coupling (NMDC) to gravity framework, with the quartic potential. Applying the coupling parameter as a two-parted function of inflaton field and fine-tuning of four parameter assortments we can acquire ultra slow-roll phase to slow down the inflaton field due to high friction. This enables us to achieve enough enhancement in the amplitude of curvature perturbations power spectra to generate PBHs with different masses. The reheating stage is considered to obtain criteria for PBHs generation during radiation dominated era. We demonstrate that two cases of asteroid mass PBHs ($10^{-13}M_{\odot}$ and $10^{-15}M_{\odot})$ can be very interesting candidates for comprising $98.3\%$ and $99.1\%$ of the total Dark Matter (DM) content of the universe. Moreover, we analyse the production of induced Gravitational Waves (GWs), and illustrate that their spectra of current density parameter $(\Omega_{\rm GW_0})$ for all parameter cases foretold by our model have climaxes which cut the sensitivity curves of GWs detectors, ergo the veracity of our outcomes can be tested in light of these detectors. At last, our numerical results exhibit that the spectra of $\Omega_{\rm GW_0}$ behave as a power-law function with respect to frequency, $\Omega_{\rm GW_0} (f) \sim (f/f_c)^{n} $, in the vicinity of climaxes. Also, in the infrared regime $f\ll f_{c}$, the power index satisfies the relation $n=3-2/\ln(f_c/f)$.

Solar wind electrons play an important role in the energy balance of the solar wind acceleration by carrying energy into interplanetary space in the form of electron heat flux. The heat flux is stored in the complex electron velocity distribution functions (VDFs) shaped by expansion, Coulomb collisions, and field-particle interactions. We investigate how the suprathermal electron deficit in the anti-strahl direction, which was recently discovered in the near-Sun solar wind, drives a kinetic instability and creates whistler waves with wave vectors that are quasi-parallel to the direction of the background magnetic field. We combine high-cadence measurements of electron pitch-angle distribution functions and electromagnetic waves provided by Solar Orbiter during its first orbit. Our case study is based on a burst-mode data interval from the Electrostatic Analyser System (SWA-EAS) at a distance of 112 $R_S$ (0.52 au) from the Sun, during which several whistler wave packets were detected by Solar Orbiter's Radio and Plasma Waves (RPW) instrument. The sunward deficit creates kinetic conditions under which the quasi-parallel whistler wave is driven unstable. We directly test our predictions for the existence of these waves through solar wind observations. We find whistler waves that are quasi-parallel and almost circularly polarised, propagating away from the Sun, coinciding with a pronounced sunward deficit in the electron VDF. The cyclotron-resonance condition is fulfilled for electrons moving in the direction opposite to the direction of wave propagation, with energies corresponding to those associated with the sunward deficit.

We construct a Bayesian inference deep learning machine for parameter estimation of gravitational wave events of binaries of black hole coalescence. The structure of our deep Bayseian machine adopts the conditional variational autoencoder scheme by conditioning both the gravitational wave strains and the variations of amplitude spectral density of the detector noise. We show that our deep Bayesian machine is capable of yielding the posteriors compatible with the ones from the nest sampling method, and of fighting against the noise outliers. We also apply our deep Bayesian machine to the LIGO/Virgo O3 events, and find that conditioning detector noise to fight against its drifting is relevant for the events with medium signal-to-noise ratios.

Gravitational waves provide us with a new window into our Universe, and have already been used to place strong constrains on the existence of light scalar fields, which are a common feature in many alternative theories of gravity. However, spin effects are still relatively unexplored in this context. In this work, we construct an effective point-particle action for a generic spinning body that can couple both conformally and disformally to a real scalar field, and we show that requiring the existence of a self-consistent solution automatically implies that if a scalar couples to the mass of a body, then it must also couple to its spin. We then use well-established effective field theory techniques to conduct a comprehensive study of spin-orbit effects in binary systems to leading order in the post-Newtonian (PN) expansion. Focusing on quasicircular nonprecessing binaries for simplicity, we systematically compute all key quantities, including the conservative potential, the orbital binding energy, the radiated power, and the gravitational-wave phase. We show that depending on how strongly each member of the binary couples to the scalar, the spin-orbit effects that are due to a conformal coupling first enter into the phase at either 0.5PN or 1.5PN order, while those that arise from a disformal coupling start at either 3.5PN or 4.5PN order. This suppression by additional PN orders notwithstanding, we find that the disformal spin-orbit terms can actually dominate over their conformal counterparts due to an enhancement by a large prefactor. Accordingly, our results suggest that upcoming gravitational-wave detectors could be sensitive to disformal spin-orbit effects in double neutron star binaries if at least one of the two bodies is sufficiently scalarised.