Papers for Monday, May 09 2022

Several recent studies have shown that velocity differences of very wide binary stars, measured to high precision with GAIA, can potentially provide an interesting test for modified-gravity theories which attempt to emulate dark matter. These systems should be entirely Newtonian according to standard dark-matter theories, while the predictions for MOND-like theories are distinctly different, if the various observational issues can be overcome. Here we provide an updated version of our 2019 study using the recent GAIA EDR3 data release: we select a large sample of 73 159 candidate wide binary stars with distance <300 parsec and magnitudes G<17 from GAIA EDR3, and estimate component masses using a main-sequence mass-luminosity relation. We then examine the frequency distribution of pairwise relative projected velocity (relative to circular-orbit value) as a function of projected separation, compared to simulations; as before, these distributions show a clear peak at a value close to Newtonian expectations, along with a long 'tail' which extends to much larger velocity ratios and may well be caused by hierarchical triple systems with an unresolved or unseen third star. We then fit these observed distributions with a simulated mixture of binary, triple and flyby populations, for GR or MOND orbits, and find that standard gravity is somewhat preferred over one specific implementation of MOND; though we have not yet explored the full parameter space of triple population models and MOND versions. Improved data from future GAIA releases, and followup of a subset of systems to better characterise the triple population, should allow wide binaries to become a decisive test of GR vs MOND in the future.

C/2017 K2 (PANSTARRS) is an Oort cloud comet previously observed to be active at heliocentric distances r>20 au on what is likely its first passage through the inner solar system. We observed the comet on 2021 March 19-20 at r=6.82 au pre-perihelion and 8.35 deg phase angle with the Hubble Space Telescope (HST), and obtained high spatial resolution photometry and polarimetry mapping the properties of dust over the coma prior to the onset of water ice sublimation activity on the nucleus. We found clear radial gradients in the color and polarization of the coma: the F475W-F775W (g'-i') reflectance slope increased from ~4.5% per 100 nm within ~10,000 km of the nucleus up to ~7% per 100 nm by ~50,000 km, while the negative polarization in F775W (i') strengthened from about -2% to -3.5% over the same range. The radial intensity profiles moreover strongly deviate from profiles simulated for stable dust grains. Near-infrared imagery obtained with the Palomar Hale Telescope on 2021 May 18 at r=6.34 au revealed a continued absence of micron-sized grains in the tail, but showed no clear spatial gradient in JHKs colors. The observed patterns collectively appear consistent with the inner coma being optically dominated by sublimating, micron-sized water ice grains, unlike the tail of more stable, millimeter-sized grains. Finally, we evaluated these results alongside other Oort cloud comets, and found in a reanalysis of HST observations of C/2012 S1 (ISON) that the near-nucleus polarimetric halo reported for that comet is likely an observational artifact.

The interstellar medium (ISM) of galaxies very often contains a gas component that reaches the temperature of several million degrees, whose physical and chemical properties can be investigated through imaging and spectroscopy in the X-rays. We review the current knowledge on the origin and retention of the hot ISM in star-forming and early-type galaxies, from a combined theoretical and observational standpoint. As a complex interplay between gravitational processes, environmental effects, and feedback mechanisms contributes to its physical conditions, the hot ISM represents a key diagnostic of the evolution of galaxies.

We aim to develop a simple prescription for migration and accretion in 1D disc models, calibrated with results of 3D hydrodynamic simulations. Our focus lies on non-self-gravitating discs, but we also discuss to what degree our prescription could be applied when the discs are self-gravitating. We study migration using torque densities. Our model for the torque density is based on existing fitting formulas, which we subsequently modify to prevent premature gap-opening. At higher planetary masses, we also apply torque densities from hydrodynamic simulations directly to our 1D model. These torque densities allow modelling the orbital evolution of an initially low-mass planet that undergoes runaway-accretion to become a massive planet. The two-way exchange of angular momentum between disc and planet is included. This leads to a self-consistent treatment of gap formation that only relies on directly accessible disc parameters. We present a formula for Bondi- and Hill- gas accretion in the disc-limited regime. This formula is self-consistent in the sense that mass is removed from the disc in the location from where it is accreted. We find that the resulting evolution in mass and semi-major axis in the 1D framework is in good agreement with those from 3D hydrodynamical simulations for a range of parameters. Our prescription is valuable for simultaneously modelling migration and accretion in 1D-models. We conclude that it is appropriate and beneficial to apply torque densities from hydrodynamic simulations in 1D models, at least in the parameter space we study here. More work is needed to in order to determine whether our approach is also applicable in an even wider parameter space and in situations with more complex disc thermodynamics, or when the disc is self-gravitating.

The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations.

There have been relatively few published long-duration, uninterrupted light curves of magnetic cataclysmic variable stars in which the accreting white dwarf's rotational frequency is slightly desynchronized from the binary orbital frequency (asynchronous polars). We report Kepler K2 and TESS observations of two such systems. The first, SDSS J084617.11+245344.1, was observed by the Kepler spacecraft for 80 days during Campaign 16 of the K2 mission, and we identify it as a new asynchronous polar with a likely 4.64 h orbital period. This is significantly longer than any other asynchronous polar, as well as all but several synchronous polars. Its spin and orbital periods beat against each other to produce a conspicuous 6.77 d beat period, across which the system's accretion geometry gradually changes. The second system in this study, Paloma, was observed by TESS for one sector and was already known to be asynchronous. Until now, there had been an ambiguity in its spin period, but the TESS power spectrum pinpoints a spin period of 2.27 h. During the resulting 0.7 d spin-orbit beat period, the light curve phased on the spin modulation alternates between being single- and double-humped. We explore two possible explanations for this behavior: the accretion flow being diverted from one of the poles for part of the beat cycle, or an eclipse of the emitting region responsible for the second hump.

Astronomers have discovered a handful of exoplanets with rocky bulk compositions but orbiting so close to their host star that the surface of the planet must be at least partially molten. It is expected that the dayside of such "lava planets" harbors a rock vapor atmosphere that flows quickly towards the airless nightside -- this partial atmosphere is critical to the interpretation of lava planet observations, but transports negligible heat towards the nightside. As a result, the surface temperature of the magma ocean may range from 3000~K near the sub-stellar point down to 1500~K near the day-night terminator. We use simple models incorporating the thermodynamics and geochemistry of partial melt to predict the physical and chemical properties of the magma ocean as a function of the distance from the sub-stellar point. Our two principal findings are that 1) the dayside magma ocean is much deeper than previously thought, probably extending down to the core-mantle boundary in some locations, and 2) much of the dayside is only partially molten, leading to gradients in the surface chemistry of the magma ocean. These findings have important implications for the dynamics of the magma ocean as well as the composition and dynamics of the atmosphere.

Measuring a pulsar's rotational evolution is crucial to understanding the nature of the pulsar. Here we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER data. For the newly-discovered fast energetic young pulsar, PSR J0058-7218, we increase the baseline of its timing model from 1.4 days to 8 months and not only measure more precisely its spin-down rate nudot = (-6.2324+/-0.0001)x10^-11 Hz s^-1 but also for the first time the second time derivative of spin rate nuddot = (4.2+/-0.2)x10^-21 Hz s^-2. For the fastest and most energetic young pulsar, PSR J0537-6910, we detect 4 more glitches, for a total of 15 glitches over 4.5 years of NICER monitoring, and show that its spin-down behavior continues to set this pulsar apart from all others, including a long-term braking index n = -1.234+/-0.009 and interglitch braking indices that asymptote to <~ 7 for long times after a glitch. For PSR J1101-6101, we measure a much more accurate spin-down rate that agrees with a previous value measured without phase-connection. For PSR J1412+7922 (also known as Calvera), we extend the baseline of its timing model from our previous 1-year model to 4.4 years, and for PSR J1849-0001, we extend the baseline from 1.5 years to 4.7 years. We also present a long-term timing model of the energetic pulsar, PSR J1813-1749, by fitting previous radio and X-ray spin frequencies from 2009-2019 and new ones measured here using 2018 NuSTAR and 2021 Chandra data.

(Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible.

Planetary obliquity is a first order control on planetary climate and seasonal contrast, which has a number of cascading consequences for life. How moderately high obliquity (obliquities greater than Earth's current obliquity up to 45$^{\circ}$) affects a planet's surface physically has been studied previously, but we lack an understanding of how marine life will respond to these conditions. We couple the ROCKE-3D general circulation model to the cGENIE 3D biogeochemical model to simulate the ocean biosphere's response to various planetary obliquities, bioessential nutrient inventories, and biospheric structure. We find that the net rate of photosynthesis increased by 35$\%$ and sea-to-air flux of biogenic oxygen doubled between the 0$^{\circ}$ and 45$^{\circ}$ obliquity scenarios, which is an equivalent response to doubling bioessential nutrients. Our results suggest that moderately high-obliquity planets have higher potential for biospheric oxygenation than their low-obliquity counterparts and that life on moderately high-obliquity habitable planets may be easier to detect with next generation telescopes. These moderately high-obliquity habitable planets may also be more conducive to the evolution of complex life.

We present new HST imaging photometry for the site of the Type IIn supernova (SN) 2009ip taken almost a decade after explosion. The optical source has continued to fade steadily since the SN-like event in 2012. In the F606W filter, which was also used to detect its luminous blue variable (LBV) progenitor 13~yr before the SN, the source at the position of SN2009ip is now 1.2mag fainter than that quiescent progenitor. It is 6-7mag fainter than the pre-SN outbursts in 2009--2011. This definitively rules out a prediction that the source would return to its previous state after surviving the 2012 event. Instead, the late-time fading matches expectations for a terminal explosion. The source fades at a similar rate in all visual-wavelength filters without significant color changes, therefore also ruling out the hypothesis of a luminous dust-obscured survivor or transition to a hotter post-LBV survivor. The late-time continuum with steady color and strong H$\alpha$ emission detected in a narrow F657N filter are, however, entirely expected for ongoing shock interaction with circumstellar material in a decade-old core-collapse SN. Interestingly, the ultraviolet flux has stayed nearly constant since 2015, supporting previous conjectures that the F275W light traces main-sequence OB stars in an underlying young star cluster. We expect that the visual-wavelength continuum will eventually level off, tracing this cluster light. Without any additional outbursts, it seems prudent to consider the 2012 event as a terminal SN explosion, and we discuss plausible scenarios.

The field of machine learning has drawn increasing interest from various other fields due to the success of its methods at solving a plethora of different problems. An application of these has been to train artificial neural networks to solve differential equations without the need of a numerical solver. This particular application offers an alternative to conventional numerical methods, with advantages such as lower memory required to store the solutions, parallelization, and in some cases less overall computational cost than its numerical counterparts. In this work, we train artificial neural networks to represent a bundle of solutions of the differential equations that govern the background dynamics of the Universe for four different models. The models we have chosen are $\Lambda \mathrm{CDM}$, the Chevallier-Polarski-Linder parametric dark energy model, a quintessence model with an exponential potential, and the Hu-Sawicki $f\left(R\right)$ model. We used the solutions that the networks provide to perform statistical analyses to estimate the values of each model's parameters with observational data; namely, estimates of the Hubble parameter from Cosmic Chronometers, the Supernovae type Ia data from the Pantheon compilation, and measurements from Baryon Acoustic Oscillations. The results we obtain for all models match similar estimations done in the literature using numerical solvers. In addition, we estimated the error of the solutions by comparing them to the analytical solution when there is one, or to a high-precision numerical solution when there is not. Through those estimations we found that the error of the solutions was at most $\sim1\%$ in the region of the parameter space that concerns the $95\%$ confidence regions that we found using the data, for all models and all statistical analyses performed in this work.

We describe a semi-analytic model to predict the triaxial shapes of dark matter halos utilizing the sequences of random merging events captured in merger trees to follow the evolution of each halo's energy tensor. When coupled with a simple model for relaxation toward a spherical shape, we find that this model predicts distributions of halo axis length ratios which approximately agree with those measured from cosmological N-body simulations in their medians and inter-quartile ranges. We demonstrate the predictive and explanatory power of this model by considering conditioned distributions of axis length, and the mass-dependence of halo shapes, finding these to be in excellent agreement with N-body results. This model provides both insight into the physics driving the evolution of halo triaxial shapes, and rapid quantitative predictions for the statistics of triaxiality connected directly to the formation history of the halo.

The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the constant-cross-section gravothermal solution to models with arbitrary velocity dependence in the cross section. We show that the gravothermal solutions after mapping are in good agreement with Arepo N-Body simulation results. We illustrate the power of this approach by applying this fast mapping method to a halo hosting a low surface brightness galaxy UGC 128. We show that this fast mapping method can be used to constrain free parameters in a physically motivated cross-section model and illustrate parameter space favored by the rotation curve measurement.

Ultralight bosonic particles are fascinating candidates of dark matter (DM). It behaves as classical waves in our Galaxy due to its large number density. There have been various methods proposed to search for the wave-like DM, such as methods utilizing interferometric gravitational-wave detectors. Understanding the characteristics of DM signals is crucial to extract the properties of DM from data. While the DM signal is nearly monochromatic with the angular frequency of its mass, the amplitude and phase are gradually changing due to the velocity dispersion of DMs in our Galaxy halo. The stochastic amplitude and phase should be properly taken into account to accurately constrain the coupling constant of DM from data. Previous works formulated a method to obtain the upper bound on the coupling constant incorporating the stochastic effects. One of these works compared the upper bound with and without the stochastic effect in a measurement time that is much shorter than the variation time scale of the amplitude and phase. In this paper, we extend their formulation to arbitrary measurement time and evaluate the stochastic effects. Moreover, we investigate the velocity-dependent signal for dark photon DM including an uncertainly of the velocity. We demonstrate that our method accurately estimates the upper bound on the coupling constant with numerical simulations. We also estimate the expected upper bound of the coupling constant of axion DM and dark photon DM from future experiments in a semi-analytic way. The stochasticity especially affects constraints on a small mass region. Our formulation offers a generic treatment of the ultralight bosonic DM signal with the stochastic effect.

The search for extraterrestrial intelligence at radio frequencies has largely been focused on continuous-wave narrowband signals. We demonstrate that broadband pulsed beacons are energetically efficient compared to narrowband beacons over longer operational timescales. Here, we report the first extensive survey searching for such broadband pulsed beacons towards 1883 stars as a part of the Breakthrough Listen's search for advanced intelligent life. We conducted 233 hours of deep observations across 4 to 8 GHz using the Robert C. Byrd Green Bank Telescope and searched for three different classes of signals with artificial (or negative) dispersion. We report a detailed search -- leveraging a convolutional neural network classifier on high-performance GPUs -- deployed for the very first time in a large-scale search for signals from extraterrestrial intelligence. Due to the absence of any signal-of-interest from our survey, we place a constraint on the existence of broadband pulsed beacons in our solar neighborhood: $\lesssim$1 in 1000 stars have transmitter power-densities $\gtrsim$10$^5$ W/Hz repeating $\leq$500 seconds at these frequencies.

Gamma-ray bursts are divided into short gamma-ray bursts and long gamma-ray bursts based on the bimodal distribution of their durations. Long bursts and short bursts are typically characterized by different statistical characteristics. Nevertheless, there are some samples that challenge such a framework, such as GRB 060614, a long-duration burst with short burst characteristics. Furthermore, gamma-ray bursts are generally considered to be an event with no periodic or repetitive behavior, since the progenitors usually undergo destructive events, such as massive explosions or binary compact star mergers. In this work, we investigated Fermi data for possible quasi-periodic oscillations and repetitive behaviors of gamma-ray bursts using timing analysis methods and report a special event GRB 201104A, which is a long-duration burst with characteristics of a short burst, and it exhibits a "repetitive" behavior. We propose that such a situation may arise from lensed short gamma-ray bursts and attempt to verify it by Bayesian analysis. In addition, we extend the spectrum analysis to Bayesian inference. In spite of the existence of at least two distinct time periods with similar spectral distributions, there is no strong evidence that they result from a lensing gamma-ray burst. Taking the gravitational-lensing scenario out of consideration, a long burst would resemble a short burst in its repetitive behavior, which presents a challenge for the current classification scheme.

In this article, we measure the mean magnetic shear from the morphological evolution of flare ribbons, and examine the evolution of flare thermal and non-thermal X-ray emissions during the progress of flare reconnection. We analyze three eruptive flares and three confined flares ranging from GOES class C8.0 to M7.0. They exhibit well-defined two ribbons along the magnetic polarity inversion line (PIL), and have been observed by the Atmospheric Imaging Assembly and the Ramaty High Energy Solar Spectroscopic Imager from the onset of the flare throughout the impulsive phase. The analysis confirms the strong-to-weak shear evolution in the core region of the flare, and the flare hard X-ray emission rises as the shear decreases. In eruptive flares in this sample, significant non-thermal hard X-ray emission lags the ultraviolet emission from flare ribbons, and rises rapidly when the shear is modest. In all flares, we observe that the plasma temperature rises in the early phase when the flare ribbons rapidly spread along the PIL and the shear is high. We compare these results with prior studies, and discuss their implications, as well as complications, related to physical mechanisms governing energy partition during flare reconnection.

With the availability of large-scale surveys like Kepler and TESS, there is a pressing need for automated methods to classify light curves according to known classes of variable stars. We introduce a new algorithm for classifying light curves that compares 7000 time-series features to find those which most effectively classify a given set of light curves. We apply our method to Kepler light curves for stars with effective temperatures in the range 6500--10,000K. We show that the sample can be meaningfully represented in an interpretable five-dimensional feature space that separates seven major classes of light curves (delta Scuti stars, gamma Doradus stars, RR Lyrae stars, rotational variables, contact eclipsing binaries, detached eclipsing binaries, and non-variables). We achieve a balanced classification accuracy of 82% on an independent test set of Kepler stars using a Gaussian mixture model classifier. We use our method to classify 12,000 Kepler light curves from Quarter 9 and provide a catalogue of the results. We further outline a confidence heuristic based on probability density with which to search our catalogue, and extract candidate lists of correctly-classified variable stars.

NGC 5548 is a very well-studied Seyfert 1 galaxy in broad wavelengths. Previous multiwavelength observation campaigns have indicated that its multiple absorbers are highly variable and complex. A previous study applied a two-zone partial covering model with different covering fractions to explain the complex X-ray spectral variation and reported a correlation between one of the covering fractions and the photon index of the power-law continuum. However, it is not straightforward to physically understand such a correlation. In this paper, we propose a model to avoid this unphysical situation; the central X-ray emission region is partially covered by clumpy absorbers composed of double layers. These "double partial coverings" have precisely the same covering fraction. Based on our model, we have conducted an extensive spectral study using the data taken by XMM-Newton, Suzaku, and NuSTAR in the range of 0.3-78 keV for 16 years. Consequently, we have found that the X-ray spectral variations are mainly explained by independent changes of the following three components; (1) the soft excess spectral component below ~1 keV, (2) the cut-off power-law normalization, and (3) the partial covering fraction of the clumpy absorbers. In particular, spectral variations above ~1 keV are mostly explained only by the changes of the partial covering fraction and the power-law normalization. In contrast, the photon index and all the other spectral parameters are not significantly variable.

Wind-fed supergiant X-ray binaries are precious laboratories not only to study accretion under extreme gravity and magnetic field conditions, but also to probe still highly debated properties of massive star winds. These includes the so-called clumps, originated from the inherent instability of line driven winds, and larger structures. In this paper, we report on the results of the last (and not yet published) monitoring campaigns that our group has been carrying out since 2007 with both XMM-Newton and the Swift Neil Gehrels observatory. Data collected with the EPIC cameras on-board XMM-Newton allow us to carry out a detailed hardness ratio-resolved spectral analysis that can be used as an efficient way to detect spectral variations associated to the presence of clumps. Long-term observations with the XRT on-board Swift, evenly sampling the X-ray emission of supergiant X-ray binaries over many different orbital cycles, are exploited to look for the presence of large scale structures in the medium surrounding the compact objects. The results reported in this paper represent the outcomes of the concluded observational campaigns we carried out on the supergiant X-ray binaries 4U 1907+09, IGR J16393-4643, IGR J19140+0951, and XTE J1855-026, as well as the supergiant fast X-ray transients IGR J17503-2636, IGR J18410-0535, and IGR J11215-5952. All results are discussed in the context of wind-fed supergiant X-ray binaries and shall ideally serve to optimally shape the next observational campaigns aimed at sources in the same classes. We show in one of the paper appendices that IGR J17315-3221, preliminary classified in the literature as a possible supergiant X-ray binary discovered by INTEGRAL, is the product of a data analysis artifact and should thus be disregarded for future studies.

The optical spectra of the B-supergiant LS III+52 24 (IRAS 22023+5249) obtained at the 6-meter telescope BTA with a resolution R$\ge$60000 in 2010-2021 revealed signs of wind variability and velocity stratification in the extended atmosphere. The H$\alpha$ and H$\beta$ lines have a P Cyg type profile; their wind absorption changes position in the range from $-270$ to $-290$ km/s. The intensity of the H$\alpha$ emission reaches record values with respect to the local continuum: I/Icont $\approx$70. The stationary radial velocity according to the positions of symmetric forbidden emissions and permitted metal emissions was taken as the systemic velocity Vsys=$-149.6\pm$0.6 km/s. Based on the positions of absorptions of NII and OII ions, a time variability of the radial velocity in the range from $-127.2$ to $-178.3$ km/s was found for the first time for this star. This variability indicates the possible presence of a companion and/or atmospheric pulsations. The change of the oxygen triplet profile OI 7775 A due to the occurrence of unstable emission was registered. The set of interstellar absorptions of the NaI D-lines profile in the range from $-10.0$ to $-167.2$ km/s is formed in the Local Arm and subsequent arms of the Galaxy. The distance to the star, d$\ge$5.3 kpc, combined with the high systemic velocity, indicates that the star is located in the Galaxy beyond the Scutum-Crux arm.

We describe a new software package for simulating channelised, high-time resolution data streams from radio telescopes. The software simulates data from the telescope and observing system taking into account the observation strategy, receiver system and digitisation. The signatures of pulsars, fast radio bursts and flare stars are modelled, including frequency-dependent effects such as scattering and scintillation. We also simulate more generic signals using spline curves and images. Models of radio frequency interference include signals from satellites, terrestrial transmitters and impulsive, broadband signals. The simulated signals can also be injected into real data sets. Uses of this software include the production of machine learning training data sets, development and testing of new algorithms to search for anomalous patterns and to characterise processing pipelines.

Coronal mass ejections (CMEs) are large clouds of magnetized plasma ejected from the Sun, and are often associated with acceleration of electrons that can result in radio emission via various mechanisms. However, the underlying mechanism relating the CMEs and particle acceleration still remains a subject of heated debate. Here, we report multi-instrument radio and extreme ultraviolet (EUV) imaging of a solar eruption event on 24 September 2011. We determine the emission mechanism of a moving radio burst, identify its three-dimensional (3D) location with respect to a rapidly expanding EUV wave, and find evidence for CME shocks that produce quasiperiodic acceleration of electron beams.

The determination of the velocities, accelerations and the gravitational field intensity at a given location in a galaxy could potentially be achieved in an unexpected manner with the environment of the observer, for instance, the local mean mass density in the galaxy. This idea, mathematically supported by the asymmetric distance concept, is illustrated here by a study regarding the rotation of spiral galaxies. This suggestion is new in the astrophysics field (in the following, it is called the \k{appa}-model) and could help to mimic the main effects seen in modified Newtonian dynamics (MOND) theory, modified gravity (MOG) models, or other related models built with the aim of eliminating dark matter that are already well-established theories. Thus, starting from two selected examples of galaxies, in section 5, we show that there is an equivalence between MOND and the \k{appa}-model. In particular, on the opposite side, we have the speculative nature of the dominant paradigm, the elusive dark matter, a matter whose properties always remain undefined despite intense theoretical, experimental and observational efforts for over 50 years.

We develop a novel semi-analytic spectral fitting approach to quantify the star-formation histories (SFHs) and chemical enrichment histories (ChEHs) of individual galaxies. We construct simple yet general chemical evolution models that account for gas inflow and outflow processes as well as star formation, to investigate the evolution of merger-free star-forming systems. These models are fitted directly to galaxies' absorption-line spectra, while their emission lines are used to constrain current gas phase metallicity and star formation rate. We apply this method to spiral galaxies selected from the SDSS-IV MaNGA survey. By fitting the co-added absorption-line spectra for each galaxy, and using the emission-line constraints on present-day metallicity and star formation, we reconstruct both the SFHs and the ChEHs for all objects in the sample. We can use these reconstructions to obtain archaeological measures of derived correlations such as the mass--metallicity relation at any redshift, which compare favourably with direct observations. We find that both the SFHs and ChEHs have strong mass dependence: massive galaxies accumulate their stellar masses and become enriched earlier. This mass dependence causes the observed flattening of the mass--metallicity relation at lower redshifts. The model also reproduces the observed gas-to-stellar mass ratio and its mass dependence. Moreover, we are able to determine that more massive galaxies have earlier gas infall times and shorter infall time-scales, and that the early chemical enrichment of low-mass galaxies is suppressed by strong outflows, while outflows are not very significant in massive galaxies.

We explore indirect methods to detect Polycyclic Aromatic Hydrocarbons (PAHs) in gas-rich, absorption-selected galaxies at high redshift. We look at the optical VLT/X-shooter observations of an intervening, extremely strong damped Lyman-alpha absorber (or ESDLA, with log(N(HI))>~21.7)) towards QSO SDSS J1143+1420 at redshift, z(ESDLA)=2.323. Literature studies have shown that this ESDLA contains signatures of dust and diffuse molecular hydrogen and it was specifically chosen for our study due to its close spatial proximity (impact parameter, rho=0.6+/-0.3 kpc) with its associated galaxy. There is no direct detection of PAHs emission in the limited observations of infrared(IR)-spectra along this sight-line. Hence, we use CLOUDY numerical simulation modelling to indirectly probe the presence of PAH in the ESDLA. We note that PAHs need to be included in the models to reproduce the observed column densities of warm H2 and CI. Thus, we infer the presence of PAHs indirectly in our ESDLA, with an abundance of PAH/H = 10^(-7.046). We also measure a low 2175 A bump strength (E(bump)~0.03-0.19 mag) relative to star-forming galaxies by modelling extinction of QSO spectra by dust at the absorber rest-frame. This is consistent with the low PAH abundance obtained indirectly using CLOUDY modelling. Our study highlights the usage of CLOUDY modelling to indirectly detect PAH in high-redshift gas-rich absorption-selected galaxies.

Oscillatory reconnection is a process that has been suggested to underlie several solar and stellar phenomena, and is likely to play an important role in transient events such as flares. Quasi-periodic pulsations (QPPs) in flare emissions may be a manifestation of oscillatory reconnection, but the underlying mechanisms remain uncertain. In this paper, we present 2D magnetohydrodynamic (MHD) simulations of two current-carrying magnetic flux ropes with an out-of-plane magnetic field undergoing oscillatory reconnection in which the two flux ropes merge into a single flux rope. We find that oscillatory reconnection can occur intrinsically without an external oscillatory driver during flux rope coalescence, which may occur both during large-scale coronal loop interactions and the merging of plasmoids in fragmented current sheets. Furthermore, we demonstrate that radially propagating non-linear waves are produced in the aftermath of flux rope coalescence, due to the post-reconnection oscillations of the merged flux rope. The behaviour of these waves is found to be almost independent of the initial out-of-plane magnetic field. It is estimated that the waves emitted through merging coronal loops and merging plasmoids in loop-top current sheets would have a typical phase speed of 90 km/s and 900 km/s respectively. It is possible that the properties of the waves emitted during flux rope coalescence could be used as a diagnostic tool to determine physical parameters within a coalescing region.

Spectropolarimetric data allow for simultaneous monitoring of stellar chromospheric $\log{R^{\prime}_{\rm{HK}}}$ activity and the surface-averaged longitudinal magnetic field, $B_l$, giving the opportunity to probe the relationship between large-scale stellar magnetic fields and chromospheric manifestations of magnetism. We present $\log{R^{\prime}_{\rm{HK}}}$ and/or $B_l$ measurements for 954 mid-F to mid-M stars derived from spectropolarimetric observations contained within the PolarBase database. Our magnetically active sample complements previous stellar activity surveys that focus on inactive planet-search targets. We find a positive correlation between mean $\log{R^{\prime}_{\rm{HK}}}$ and mean $\log|B_l|$, but for G stars the relationship may undergo a change between $\log{R'_{\rm{HK}}}\sim-4.4$ and $-4.8$. The mean $\log{R^{\prime}_{\rm{HK}}}$ shows a similar change with respect to the $\log{R^{\prime}_{\rm{HK}}}$ variability amplitude for intermediately-active G stars. We also combine our results with archival chromospheric activity data and published observations of large-scale magnetic field geometries derived using Zeeman Doppler Imaging. The chromospheric activity data indicate a slight under-density of late-F to early-K stars with $-4.75\leq\log{R'_{\rm HK}}\leq-4.5$. This is not as prominent as the original Vaughan-Preston gap, and we do not detect similar under-populated regions in the distributions of the mean $|B_l|$, or the $B_l$ and $\log{R'_{\rm HK}}$ variability amplitudes. Chromospheric activity, activity variability and toroidal field strength decrease on the main sequence as rotation slows. For G stars, the disappearance of dominant toroidal fields occurs at a similar chromospheric activity level as the change in the relationships between chromospheric activity, activity variability and mean field strength.

From the Sun we know that coronal mass ejections (CMEs) are a transient phenomenon, often correlated with flares. They have an impact on solar mass- and angular momentum loss, and therefore solar evolution, and make a significant part of space weather. The same is true for stars, but stellar CMEs are still not well constrained, although new methodologies have been established, and new detections presented in the recent past. So far, probable detections of stellar CMEs have been presented, but their physical parameters which are not directly accessible from observations, such as electron density, optical thickness, temperature, etc., have been so far not determined for the majority of known events. We apply cloud modeling, as commonly used on the Sun, to a known event from the literature, detected on the young dMe star V374 Peg. This event manifests itself in extra emission on the blue side of the Balmer lines. By determining the line source function from 1D NLTE modeling together with the cloud model formulation we present distributions of physical parameters of this event. We find that except for temperature and area all parameters are at the upper range of typical solar prominence parameters. The temperature and the area of the event were found to be higher than for typical solar prominences observed in Balmer lines. We find more solutions for the filament than for the prominence geometry. Moreover we show that filaments can appear in emission on dMe stars contrary to the solar case.

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

We present an improved background model for the Large Area X-ray Proportional Counter (LAXPC) detectors on-board AstroSat. Because of the large collecting area and high pressure, the LAXPC instrument has a large background count rate, which varies during the orbit. Apart from the variation with latitude and longitude during the orbit there is a prominent quasi-diurnal variation which has not been modelled earlier. Using over 5 years of background observations, we determined the period of the quasi-diurnal variation to be 84495 s and using this period, it is possible to account for the variation and also identify time intervals where the fit is not good. These lead to a significant improvement in the background model. The quasi-diurnal variation can be ascribed to the changes in charged particle flux in the near Earth orbit.

The eclipsing time variations in post-common-envelope binaries were proposed to be due to the time-varying component of the stellar quadrupole moment. This is suggested to be produced by changes in the stellar structure due to an internal redistribution of the angular momentum and the effect of the centrifugal force. We examine this hypothesis presenting 3D simulations of compressible magneto-hydrodynamics (MHD) performed with the Pencil Code, modeling the stellar dynamo for a solar mass star with angular velocities of 20 and 30 times solar. We include and vary the strength of the centrifugal force, comparing with reference simulations without the centrifugal force and including a simulation where its effect is enhanced. The centrifugal force is causing perturbations in the evolution of the stars, so that the outcome in the details becomes different as a result of non-linear evolution. While the average density profile is unaffected by the centrifugal force, a relative change in density difference between high altitudes and the equator of order $\sim 10^{-4}$ is found. The power spectrum of the convective velocity is found to be more sensitive to the angular velocity than the strength of the centrifugal force. The quadrupole moment of the stars includes a fluctuating and a time-independent component which varies with the rotation rate. As very similar behavior is produced in the absence of the centrifugal force, we conclude that it is not the main ingredient for producing the time-averaged quadrupole moment of the star. In a real physical system, we thus expect contributions from both components, that is, due to the time-dependent gravitational force from the variation in the quadrupole term and due to spin-orbit coupling due to the persistent part of the quadrupole.

High-resolution observations of horizontal motions in the penumbra are needed to complement the concept of penumbrae obtained from spectropolarimetry. Time series of intensity images of a large sunspot in AR 10634 acquired with the Swedish Solar Telescope in the G band and red continuum are analysed. The two simultaneous time series last six hours and five minutes. Horizontal motions of penumbral grains (PGs), structures in dark bodies of filaments, the outer penumbral border, and G-band bright points are measured in time slices that cover the whole width of the penumbra and the neighbouring granulation. The spatial and temporal resolutions are 90 km and 20.1 s, respectively. In the inner penumbra, PGs move toward the umbra (inwards) with a mean speed of -0.7 km s-1. The direction of motion changes from inwards to outwards at approximately 60% of the penumbral width and the mean speed increases gradually in the outer penumbra, approaching 0.5 km/s. This speed is also typical of an expansion of the penumbra-granulation border during periods that typically last one hour and are followed by a fast contraction. The majority of the G-band bright points moves away from the sunspot, with a typical speed of 0.6 km/s. High outward speeds, 3.6 km/s on average, are observed in dark bodies of penumbral filaments. According to the model of penumbral filaments, it is suggested that the speeds detected in the dark bodies of filaments are associated with the Evershed flow and that the opposite directions of PG motions in the inner and outer penumbrae may be explained by the interaction of rising plasma in filament heads with a surrounding, differently inclined magnetic field.

The evolution of molecular interstellar clouds is a complex, multi-scale process. The power-law density exponent describes the steepness of density profiles, and it has been used to characterize the density structures of the clouds yet its usage is usually limited to spherically symmetric systems. Importing the Level-Set Method, we develop a new formalism that generates robust maps of a generalized density exponent $k_{\rho}$ at every location for complex density distributions. By applying it to high fidelity, high dynamical range map of the Perseus molecular cloud constructed using data from the Herschel and Planck satellites, we find that the density exponent exhibits a surprisingly wide range of variation ($-3.5 \lesssim k_{\rho} \lesssim -0.5$). Regions at later stages of gravitational collapse are associated with steeper density profiles. Inside a region, gas located in the vicinities of dense structures has very steep density profiles with $k_{\rho} \approx -3$, which forms because of depletion. This density exponent analysis reveals diverse density structures, forming a coherent picture that gravitational collapse leads to a continued steepening of the density profile. We expect our method to be effective in studying other power-law-like density structures, including granular materials and the Large-Scale Structure of the Universe.

Recently studies discovered that part of the Gould Belt belongs to a 2.7 kpc-long {coherent, thin} wave consisting of a chain of clouds, where a damped undulation pattern has been identified from the spatial arrangement of the clouds. We use the proper motions of Young Stellar Objects (YSOs) anchored inside the clouds to study the kinematic structure of the Radcliffe Wave in terms of $v_z$, and identify a damped, wave-like pattern from the $v_z$ space, which we call "velocity undulation". We propose a new formalism based on the Ensemble Empirical Mode Decomposition (EEMD) to determine the amplitude, period, and phase of the undulation pattern, and find that the spatial and the velocity undulation share an almost identical spatial frequency of about 1.5 kpc, and both are damped when measured from one side to the other. Measured for the first cycle, they exhibit a phase difference of around $2\pi/3$. The structure is oscillating around the midplane of the Milky Way disk with an amplitude of $\sim\,130\,\pm\,20\,\rm pc$. The vertical extent of the Radcliffe Wave exceeds the thickness of the molecular disk, suggesting that the undulation of the undulation signature might originate from a perturbation, e.g. the passage of a dwarf galaxy.

On 2021 August 8, the recurrent nova RS Ophiuchi erupted again, after an interval of 15.5 yr. Regular monitoring by the Neil Gehrels Swift Observatory began promptly, on August 9.9 (0.37 day after the optical peak), and continued until the source passed behind the Sun at the start of November, 86 days later. Observations then restarted on day 197, once RS Oph emerged from the Sun constraint. This makes RS Oph the first Galactic recurrent nova to have been monitored by Swift throughout two eruptions. Here we investigate the extensive X-ray datasets from 2006 and 2021, as well as the more limited data collected by EXOSAT in 1985. The hard X-rays arising from shock interactions between the nova ejecta and red giant wind are similar following the last two eruptions. In contrast, the early super-soft source (SSS) in 2021 was both less variable and significantly fainter than in 2006. However, 0.3-1 keV light-curves from 2021 reveal a 35 s quasi-periodic oscillation consistent in frequency with the 2006 data. The Swift X-ray spectra from 2021 are featureless, with the soft emission typically being well parametrized by a simple blackbody, while the 2006 spectra showed much stronger evidence for superimposed ionized absorption edges. Considering the data after day 60 following each eruption, during the supersoft phase the 2021 spectra are hotter, with smaller effective radii and lower wind absorption, leading to an apparently reduced bolometric luminosity. We explore possible explanations for the gross differences in observed SSS behaviour between the 2006 and 2021 outbursts.

Optical clusters identified from red-sequence galaxies suffer from projection effects, where interloper galaxies along the line-of-sight to a cluster are mistaken as genuine members of the cluster. In the previous study (Sunayama et al. 2020), we found that the projection effects cause the boost on the amplitudes of clustering and lensing on large scale compared to the expected amplitudes in the absence of any projection effects. These boosts are caused by preferential selections of filamentary structure aligned to the line-of-sight due to distance uncertainties in photometric surveys. We model the projection effects with two simple assumptions and develop a novel method to quantify the size of the boost using cluster-galaxy cross-correlation functions. We validate our method using mock cluster catalogs built from cosmological N-body simulations and find that we can obtain unbiased constraints on the boost parameter with our model. We then apply our analysis on the SDSS redMaPPer clusters and find that the size of the boost is roughly 20% for all the richness bins except the cluster sample with the richness bin $\lambda \in [30,40]$. This is the first study to constrain the boost parameter independent from cluster cosmology studies and provides a self-consistency test for the projection effects.

Star cluster formation in giant molecular clouds involves the local collapse of the cloud into small gas-rich subclusters, which can then subsequently collide and merge to build up the final star cluster(s). In this paper, we simulate collisions between these subclusters, using coupled smooth particle hydrodynamics for the gas and N-body dynamics for the stars. We are guided by previous radiation hydrodynamics simulations of molecular cloud collapse which provide the global properties of the colliding clusters, such as their stellar and gas masses, and their initial positions and velocities. The subclusters in the original simulation were treated as sink particles which immediately merged into a single entity after the collision. We show that the more detailed treatment provides a more complex picture. At collisional velocities above ~ 10 km/s, the stellar components of the cluster do not form a monolithic cluster within 3 Myr, although the gas may do so. At lower velocities, the clusters do eventually merge but over timescales that may be longer than the time for a subsequent collision. The structure of the resultant cluster is not well-fit by any standard density distribution, and the clusters are not in equilibrium but continue to expand over our simulation time. We conclude that the simple sink particle treatment of subcluster mergers in large-scale giant molecular cloud simulations provides an upper limit on the final cluster properties.

We analyze the physical conditions and chemical composition of the photoionized Herbig-Haro object HH~514, which emerges from the proplyd 170-337 in the core of the Orion Nebula. We use high-spectral resolution spectroscopy from UVES at the Very Large Telescope and IFU-spectra from MEGARA at the Gran Telescopio de Canarias. We observe two components of HH~514, the jet base and a knot, with $n_{\rm e}= (2.3 \pm 0.1) \times 10^5 \text{cm}^{-3}$ and $n_{\rm e}= (7 \pm 1) \times 10^4 \text{cm}^{-3}$, respectively, both with $T_{\rm e}\approx 9000 \text{ K}$. We show that the chemical composition of HH~514 is consistent with that of the Orion Nebula, except for Fe, Ni and S, which show higher abundances. The enhanced abundances of Fe and Ni observed in HH objects compared with the general interstellar medium is usually interpreted as destruction of dust grains. The observed sulphur overabundance (more than two times solar) is challenging to explain since the proplyd photoevaporation flow from the same disk shows normal sulphur abundance. If the aforementioned S-overabundance is due to dust destruction, the formation of sulfides and/or other S-bearing dust reservoirs may be linked to planet formation processes in protoplanetary disks, which filter large sulfide dust grains during the accretion of matter from the disk to the central star. We also show that published kinematics of molecular emission close to the central star are not consistent with either a disk perpendicular to the optical jet, nor with an outflow that is aligned with it.

The clustering of dark matter halos depends on the assembly history of halos at fixed halo mass; a phenomenon referred to as \textit{halo assembly bias}. Halo assembly bias is readily observed in cosmological simulations of dark matter. However, it is difficult to detect it in observations. The identification of galaxy or cluster properties correlated with the formation time of the halo at fixed halo mass and the ability to select galaxy clusters free from projection effects are the two most significant hurdles in the observational detection of halo assembly bias. The latter, in particular, can cause a misleading detection of halo assembly bias by boosting the amplitude of lensing and clustering on large scales. This study uses twelve different properties of central galaxies of SDSS redMaPPer clusters derived from spectroscopy to divide the clusters into sub-samples. We test the dependence of the clustering amplitude on these properties at fixed richness. We first infer halo mass and bias using weak lensing signals around the clusters using shapes of galaxies from the SDSS survey. We validate the bias difference between the two subsamples using cluster-galaxy cross-correlations. This methodology allows us to decouple the contamination due to the projection effects from the halo assembly bias signals. We do not find any significant evidence of a difference in the clustering amplitudes correlated with any of our explored properties. Our results indicate that central galaxy properties may not correlate significantly with the halo assembly histories at fixed richness.

Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($\lambda$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions.

The discovery of the chaotic motion of the planets in the Solar System dates back more than 30 years. Still, no analytical theory has satisfactorily addressed the origin of chaos so far. Implementing canonical perturbation theory in the computer algebra system TRIP, we systematically retrieve the secular resonances at work along the orbital solution of a forced long-term dynamics of the inner planets. We compare the time statistic of their half-widths to the ensemble distribution of the maximum Lyapunov exponent and establish dynamical sources of chaos in an unbiased way. New resonances are predicted by the theory and checked against direct integrations of the Solar System. The image of an entangled dynamics of the inner planets emerges.

We present a charge distribution based model that computes the infrared spectrum of polycyclic aromatic hydrocarbon (PAH) molecules using recent measurements or quantum chemical calculations of specific PAHs. The model is applied to a sample of well-studied photodissociation regions (PDRs) with well-determined physical conditions (the radiation field strength, $G_{0}$, electron density $n_{e}$, and the gas temperature, $T_{\rm gas}$). Specifically, we modelled the emission of five PAHs ranging in size from 18 to 96 carbon atoms, over a range of physical conditions characterized by the ionization parameter $\gamma = G_{0}\times T_{\rm gas}^{1/2}/n_{e}$. The anions emerge as the dominant charge carriers in low $\gamma $ ($< 2\times 10^{2}$) environments, neutrals in the intermediate $\gamma$ ($10^{3} - 10^{4}$) environments, and cations in the high $\gamma$ ($ > 10^{5}$) environments. Furthermore, the PAH anions and cations exhibit similar spectral characteristics. The similarity in the cationic and anionic spectra translates into the interpretation of the 6.2/(11.0+11.2) band ratio, with high values of this ratio associated with large contributions from either cations or anions. The model's predicted values of 6.2/(11.0+11.2) and 3.3/6.2 compared well to the observations in the PDRs NGC 7023, NGC 2023, the horsehead nebula, the Orion bar, and the diffuse ISM, demonstrating that changes in the charge state can account for the variations in the observed PAH emission. We also reassess the diagnostic potential of the 6.2/(11.0+11.2) vs 3.3/(11.0+11.2) ratios and show that without any prior knowledge about $\gamma$, the 3.3/(11.0+11.2) can predict the PAH size, but the 6.2/(11.0+11.2) cannot predict the $\gamma$ of the astrophysical environment.

Understanding the formation and evolution of high mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are reliable and representative of clusters real parameters and the impact of line-of-sight orientation. In this paper we derive parameters for a simulated cluster formed from a cloud-cloud collision with the full 6D phase space, and compare them with those derived from three different 2D line-of-sight orientations for the cluster. We show the same qualitative conclusions can be reached when viewing clusters in 2D versus 3D, but that drawing quantitative conclusions when viewing in 2D is likely to be inaccurate. The greatest divergence occurs in the perceived kinematics of the cluster, which in some orientations appears to be expanding when the cluster is actually contracting. Increases in the cluster density compounds pre-existing perspective issues, reducing the relative accuracy and consistency of properties derived from different orientations. This is particularly problematic for determination of the number, and membership, of subclusters present in the cluster. We find the fraction of subclusters correctly identified in 2D decreases as the cluster evolves, reaching less than 3.4% at the evolutionary end point for our cluster.

In a substantial number of core-collapse supernovae, early-time interaction indicates a dense circumstellar medium (CSM) that may be produced by outbursts from the progenitor star. Wave-driven mass loss is a possible mechanism to produce these signatures, with previous work suggesting this mechanism is most effective for low mass ($ \sim \! 11 \, M_\odot$) SN progenitors. Using one-dimensional hydrodynamic simulations with MESA, we study the effects of this wave heating in SN progenitors of masses $M_{\rm ZAMS} = 10-13\, M_{\odot}$. This range encompasses stars that experience semi-degenerate, central neon burning and more degenerate, off-center neon ignition. We find that central Ne ignition at $M_{\rm ZAMS} = 11\, M_{\odot}$ produces a burst of intense wave heating that transmits $\sim10^{47}$ erg of energy at $10$ years before core collapse, whereas other masses experience smaller levels of wave heating. Wave heating does not hydrodynamically drive mass loss in any of our models and is unlikely to produce a very massive CSM on its own. However, wave heating can cause large radial expansion (by more than an order of magnitude), photospheric cooling, and luminosity brightening by up to $\sim 10^6\, L_{\odot}$ in hydrogen-poor stripped star models. Some type Ib/c progenitors could drastically change their appearance in the final years of their lives, with brightness in visual bands increasing by nearly 3 mags. Moreover, interaction with a close binary companion could drive intense mass loss, with implications for type Ibn and other interaction-powered SNe.

Stars in massive star clusters exhibit intrinsic variations in some light elements (MPs) that are difficult to explain in a coherent formation scenario. In recent years, HST photometry has led to the characterisation of the global properties of th MPs in an unparalleled level of detail. In particular, the colour-colour diagrams known as chromosome maps have been proven to be very efficient at separating cluster stars with field-like metal abundance distribution (1P) from object with distinctive light-element abundance anti-correlations (2P). The unexpected wide colour ranges covered by the 1P group in the chromosome maps of the majority of the investigated Galactic GCs have been recently attributed to intrinsic metallicity variations up to ~0.30 dex, from the study of SGB stars in two metal rich Galactic GCs by employing appropriate HST filter combinations. On the other hand, high-resolution spectroscopy of small samples of 1P stars in NGC 3201 and NGC 2808 have so far provided conflicting results, with a spread of metallicity detected in NGC 3201 but not in NGC 2808. We present a new method that employs HST near-UV and optical photometry of RGB stars, to independently confirm these results. Our approach has been firstly validated using observational data for M2, a GC hosting a small group of 1P stars with enhanced (by ~0.5 dex) metallicity with respect to the main component. We have then applied our method to three clusters that cover a much larger metallicity range, and have well populated, extended first population sequences in their chromosome maps, namely M92, NGC2808, and NGC6362. We confirm that metallicity spreads are present among 1P stars in these clusters, thus solidifying the case for the existence of unexpected variations up to a factor of two of metal abundances in most GCs. We also confirm the complex behaviour of the mean metallicity differences between 1P and 2P stars.

We present a novel global 3-D coronal MHD model called COCONUT, polytropic in its first stage and based on a time-implicit backward Euler scheme. Our model boosts run-time performance in comparison with contemporary MHD-solvers based on explicit schemes, which is particularly important when later employed in an operational setting for space weather forecasting. It is data-driven in the sense that we use synoptic maps as inner boundary input for our potential field initialization as well as an inner boundary condition in the further MHD time evolution. The coronal model is developed as part of the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) and will replace the currently employed, more simplistic, empirical Wang-Sheeley-Arge (WSA) model. At 21.5 Rs where the solar wind is already supersonic, it is coupled to EUHFORIA's heliospheric model. We validate and benchmark our coronal simulation results with the explicit-scheme Wind-Predict model and find good agreement for idealized limit cases as well as real magnetograms, while obtaining a computational time reduction of up to a factor 3 for simple idealized cases, and up to 35 for realistic configurations, and we demonstrate that the time gained increases with the spatial resolution of the input synoptic map. We also use observations to constrain the model and show that it recovers relevant features such as the position and shape of the streamers (by comparison with eclipse white-light images), the coronal holes (by comparison with EUV images) and the current sheet (by comparison with WSA model at 0.1 AU).

Supernovae (SNe) are among the most energetic events in the universe still far from being fully understood. An early and prompt detection of neutrinos is a one-time opportunity for the realization of the first multi-messenger observation of these events. In this work, we present the prospects of detecting neutrinos produced before (pre-SN) and during a SN while running an advanced cryogenic detector. The recent advancements of the cryogenic detector technique and the discovery of coherent elastic neutrino-nucleus scattering offer a wealth of opportunities in neutrino detection. The combination of the excellent energy resolution of this experimental technique, with the high cross section of this detection channel and its equal sensitivity to all neutrino flavors enables the realization of highly sensitive cm-scale neutrino telescopes, as the newly proposed RES-NOVA experiment. We present a detailed study on the detection promptness of pre-SN and SN neutrino signals, with direct comparisons among different classes of test statistics. While the well-established Poisson test offers in general best performance under optimal conditions, the non-parametric Recursive Product of Spacing statistical test (RPS) is more robust and ideal for triggering astrophysical neutrino signals with no specific prior knowledge. Based on our statistical tests the RES-NOVA experiment is able to identify SN neutrino signals at a 15 kpc distance with 95% of success rate, and pre-SN signal as far as 480 pc with a pre-warn time of the order of 10 s. These results demonstrate the potential of RPS for the identification of neutrino signals and the physics reach of the RES-NOVA experiment.

The search for a particle dark matter signal in terms of radiation produced by dark matter annihilation or decay has to cope with the extreme faintness of the predicted signal and the presence of masking astrophysical backgrounds. It has been shown that using the correlated information between the dark matter distribution in the Universe with the fluctuations of the cosmic radiation fields has the potential to allow setting apart a pure dark matter signal from astrophysical emissions, since spatial fluctuations in the radiation field due to astrophysical sources and dark matter emission have different features. The cross-correlation technique has been proposed and adopted for dark matter studies by looking at dark matter halos (over-densities). In this paper we extend the technique by focusing on the information on dark matter distribution offered by cosmic voids, and by looking specifically at the gamma-ray dark matter emission: we show that, while being under-dense and therefore producing a reduced emission as compared to halos, nevertheless in voids the relative size of the cross-correlation signal due to decaying dark matter vs. astrophysical sources is significantly more favourable, producing signal-to-background ratios $S/B$ (even significantly) larger than 1 for decay lifetimes up to $2 \times 10^{30}$ s. This is at variance with the case of halos, where $S/B$ is typically (even much) smaller than 1. We show that forthcoming galaxy surveys such as Euclid combined with future generation gamma-ray detectors with improved specifications have the ability to provide a hint of such a signal with a predicted significance up to $4.2\sigma$ for galaxies and $2.7\sigma$ for the cosmic shear. The bound on the dark matter lifetime attainable exploiting voids is predicted to improve on current bounds in a mass range for the WIMP of $20\div200$ GeV.

In recent decades our understanding of solar active regions (ARs) has improved substantially due to observations made with better angular resolution and wider spectral coverage. While prior AR observations have shown that these structures were always brighter than the quiet Sun at centimeter wavelengths, recent observations at millimeter and submillimeter wavelengths have shown ARs with well defined dark umbrae. Given this new information, it is now necessary to update our understanding and models of the solar atmosphere in active regions. In this work, we present a data-constrained model of the AR solar atmosphere, in which we use brightness temperature measurements of NOAA 12470 at three radio frequencies: 17 (NoRH), 100 and 230 GHz (ALMA). Based on our model, which assumes that the radio emission originates from thermal free-free and gyroresonance processes, we calculate radio brightness temperature maps that can be compared with the observations. The magnetic field at distinct atmospheric heights was determined in our modelling process by force-free field extrapolation using photospheric magnetograms taken by HMI/SDO. In order to determine the best plasma temperature and density height profiles necessary to match the observations, the model uses a genetic algorithm that modifies a standard quiet Sun atmospheric model. Our results show that the height of the transition region (TR) of the modelled atmosphere varies with the type of region being modelled: for umbrae the TR is located at 1080 +/- 20 km above the solar surface; for penumbrae, the TR is located at 1800 +/- 50 km; and for bright regions outside sunspots, the TR is located at 2000 +/- 100 km. With these results, we find good agreement with the observed AR brightness temperature maps. Our modelled AR can be used to estimate the emission at frequencies without observational coverage.

We study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and \Moller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides a good approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides an approximation to differential self-scattering for most of the halo evolution. However, it can significantly underestimate the growth rate of the central density at late stages of the collapse phase. This indicates that constant and velocity-dependent dark matter self-interactions are fundamentally different, as for the latter the cross section evolves with the halo dynamically, boosting the collapse. This feature may help test different self-interacting dark matter models.

In the era of precision cosmology it has became crucial to find new and competitive probes to estimate cosmological parameters, in an effort of finding answers to the current cosmological tensions/discrepancies. In this work, we show the possibility of using observations of Super Massive Black Hole (SMBH) shadows as an anchor for the distance ladder, substituting the sources usually exploited for such purpose, such as Cepheid variable stars. Compared to the standard approaches, the use of SMBH has the advantage of not needing to be anchored with distance calibrators outside the Hubble flow since the shadows physical size can be estimated knowing the mass of the SMBH. Furthermore, SMBH are supposed to inhabit the center of all galaxies which, in principle, means that we can measure the size of the shadows in any Supernova type Ia host galaxy. Under the assumption that the mass of the SMBH can be accurately and reliably estimated, we find that the Hubble constant can be constrained with a $\approx10\%$ precision even considering current experimental design of ground-based interferometers. By constructing a SMBH catalogue based on a specific choice of the SMBH Mass Function (BHMF), we forecast the constraints on the Hubble constant, finding that a precision of $\approx4\%$ may be within reach of future interferometers.

In recent observations of extremely metal-poor low-mass starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. We investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as Pair-Instability Supernovae (PISN). We run chemical evolution models with yields that account for the contribution by PISN. We use both the recent non-rotating stellar yields from Goswami et al. 2021, and new yields from rotating very massive stars calculated on purpose in this work. We also search for the best initial mass function (IMF) that is able to reproduce the observations. We can reproduce the observations by adopting a bi-modal IMF and by including an initial burst of rotating very massive stars. Only with a burst of very massive stars can we reproduce the almost solar Fe/O ratios at the estimated young ages. We also confirm that rotation is absolutely needed to concomitantly reproduce the observed N/O ratios. These results stress the importance of very massive stars in galactic chemical evolution studies and strongly support a top-heavy initial mass function in the very early evolutionary stages of metal poor starburst galaxies.

We study the scenario of dark photon Dark Matter where the mass is generated through the Higgs mechanism rather than the constant Stueckelberg mass. In this construction the dark sector contains necessarily extra degrees of freedom and interactions that lead to non-trivial dynamics including thermalization, phase transitions, cosmic string production. As a consequence the predictions of Stueckelberg theories are vastly modified, strongly depending on the couplings to curvature and on the scale of inflation $H_I$ compared to the scale $f$ of spontaneous symmetry breaking. We find in particular that only in extreme regions of parameter space the phenomenology of Stueckelberg dark photon is reproduced. These scenarios are strongly constrained by isocurvature perturbations unless the dark sector is approximately Weyl invariant.

It is proposed that the rapid observed homogeneous nucleation of ice dust in a cold, weakly-ionized plasma depends on the formation of negative hydroxyl ions by fast electrons impacting water molecules. These OH$^{-}$ ions attract neutral water molecules because of the high dipole moment of the water molecules and so hydrates of the form (OH)$^{-}$(H$_{2}$O)$_{n}$ are formed. The hydrates continuously grow in the cold environment to become macroscopic ice grains. These ice grains are negatively charged as a result of electron impact and so continue to attract water molecules. Because the hydroxyl ions are negative, unlike positive ions they do not suffer recombination loss from collision with plasma electrons. Recombination with positive ions is minimal because positive ions are few in number (weak ionization) and slow-moving as result of being in thermal equilibrium with the cold background gas.

We discuss a minimal extension of the Standard Model (SM) where a single real scalar field serves as both inflaton and dark matter. The corresponding Lagrangian contains the renormalizable interactions of the inflaton field. Quantum effects generally induce a non-minimal coupling to gravity which facilitates inflation consistent with the PLANCK constraints. A large fraction of the inflaton quanta produced after inflation must be converted into the SM radiation reheating the Universe and the rest remains dark matter today. We consider thermal and non-thermal production of inflaton dark matter. In the non-thermal case, we take into account collective effects with the help of lattice simulations. Combining analytic and numerical results with the unitarity consideration, we find that the inflaton dark matter model is viable only in the thermal case where the inflaton mass is near half the Higgs mass.

Gravitational entropy is an elusive concept. Various theoretical proposals have been presented, initially based on Penrose's Weyl Curvature Hypothesis, and variations of it. A more recent proposal by Clifton, Ellis, and Tavakol (CET) considered a novel approach by defining such entropy from a Gibbs equation constructed from an effective stress-energy tensor that emerges from the 'square root' algebraic decomposition of the Bel-Robinson tensor, the simplest divergence-less tensor related to the Weyl tensor. Since, so far all gravitational entropy proposals have been applied to highly restrictive and symmetric spacetimes, we probe in this paper the CET proposal for a class of much less idealized spactimes (the Szekeres class I models) capable of describing the joint evolution of arrays of arbitrary number of structures: overdensities and voids, all placed on selected spatial locations in an asymptotic $\Lambda$CDM backgound. By using suitable covariant variables and their fluctuations, we find the necessary and sufficient conditions for a positive CET entropy production to be a negative sign of the product of the density and Hubble expansion fluctuations. To examine the viability of this theoretical result we examine numerically the CET entropy production for two elongated over dense regions surrounding a central spheroidal void, all evolving jointly from initial linear perturbations at the last scattering era into present day Mpc-size CDM structures. We show that CET entropy production is positive for all times after last scattering at the precise spatial locations where structure growth occurs and where the exact density growing mode is dominant. The present paper provides the least idealized (and most physically robust) probe of a gravitational entropy proposal in the context of structure formation.

This study investigates the properties of waves in relativistic extended magnetohydrodynamics (RXMHD), which includes Hall and electron thermal inertia effects. We focus on the case when the electron temperature is ultrarelativistic, and thus, the electron thermal inertia becomes finite at near the proton inertial scale. We derive the linear dispersion relation of RXMHD and find that the Hall and electron thermal inertia effects couple with the displacement current, giving rise to three superluminous waves in addition to the slow, fast, and Alfv\'en waves. We also show that the phase- and group-velocity surfaces of fast and Alfv\'en waves are distorted by the Hall and electron thermal inertia effects. There is a range of scales where the group velocity of fast wave is smaller than that of the Alfv\'en and slow waves. These findings are applicable to a region near the funnel base of low-luminosity accretion flows where electrons can be ultrarelativistic.

We study disc galaxies in the framework of General Relativity to focus on the possibility that, even in the low energy limit, there are relevant corrections with respect to the purely Newtonian approach. Our analysis encompasses the model both considering a low energy expansion and exact solutions, making it clear the connection between these different approaches. In particular, we focus on two different limits: the well known gravitomagnetic analogy and a new limit, called strong gravitomagnetism, which has corrections in c of the same order as the Newtonian terms. We show that these two limits of the general class of solution can account for the observed flat velocity profile, contrary to what happens using Newtonian models, where a dark matter contribution is required. Hence, we suggest a geometrical origin for a certain amount of dark matter effects.

Explosive astrophysical systems - such as supernovae or compact star binary mergers - provide conditions where exotic degrees of freedom can be populated. Within the covariant density functional theory of nuclear matter we build several general purpose equations of state which, in addition to the baryonic octet, account for $\Delta(1232)$ resonance states. The thermodynamic stability of $\Delta$-admixed nuclear matter is investigated in the limiting case of vanishing temperature for charge fractions $Y_Q=0.01$ and $Y_Q=0.5$ and wide ranges of the coupling constants to the scalar and vector mesonic fields. General purpose equation of state models with exotica presently available on the \textsc{CompOSE} database are further reviewed; for a selection of them we then investigate thermal properties for thermodynamic conditions relevant for core-collapse supernovae and binary neutron star mergers. Modifications induced by hyperons, $\Delta(1232)$, $K^-$, pions and quarks are discussed.

Extended theories of gravity have been extensively investigated during the last thirty years, aiming at fixing infrared and ultraviolet shortcomings of General Relativity and of the associated $\Lambda$CDM cosmological model. Recently, non-local theories of gravity have drawn increasing attention due to their potential to ameliorate both the ultraviolet and infrared behavior of gravitational interaction. In particular, Integral Kernel theories of Gravity provide a viable mechanism to explain the late time cosmic acceleration so as to avoid the introduction of any form of unknown dark energy. On the other hand, these models represent a natural link towards quantum gravity. Here, we study a scalar-tensor equivalent model of General Relativity corrected with non-local terms, where corrections are selected by the existence of Noether symmetries. After performing the weak field limit and generalizing the results to extended mass distributions, we analyse the non-local model at galaxy cluster scales, by comparing the theoretical predictions with gravitational lensing observations from the CLASH program. We obtain agreement with data at the same level of statistical significance as General Relativity. We also provide constraints for the Navarro--Frenk--White parameters and lower bounds for the non-local length scales. The results are finally compared with those from the literature.