Papers for Friday, Feb 11 2022

An exact turbulence universal scaling law-8/3 for the density fluctuations of cosmic ray (CR) is obtained on the basis of a new analytical compressible turbulence theory and known two-fluid model of the CR dynamics. It is shown that the origin of this scaling law may be due to the breaking of the nonlinear simple waves in CR medium near the scale of their Larmor radii as for the space plasma of solar wind and magnetosheath.

Observations that require large physical instrument dimensions and/or a considerable amount of cryogens, as it is for example the case for high spatial resolution far infrared astronomy, currently still face technological limits for their execution from space. The high cost and finality of space missions furthermore call for a very low risk approach and entail long development times. For certain spectral regions, prominently including the mid- to far-infrared as well as parts of the ultraviolet, stratospheric balloons offer a flexible and affordable complement to space telescopes, with short development times and comparatively good observing conditions. Yet, the entry burden to use balloon-borne telescopes is high, with research groups typically having to shoulder part of the infrastructure development as well. Aiming to ease access to balloon-based observations, we present the efforts towards a community-accessible balloon-based observatory, the European Stratospheric Balloon Observatory (ESBO). ESBO aims at complementing space-based and airborne capabilities over the next 10-15 years and at adding to the current landscape of scientific ballooning activities by providing a service-centered infrastructure for broader astronomical use, performing regular flights and offering an operations concept that provides researchers with a similar proposal-based access to observation time as practiced on ground-based observatories. We present details on the activities planned towards the goal of ESBO, the current status of the STUDIO (Stratospheric UV Demonstrator of an Imaging Observatory) prototype platform and mission, as well as selected technology developments with extensibility potential to space missions undertaken for STUDIO.

Asteroid-mass primordial black holes (PBH) can explain the observed dark matter abundance while being consistent with the current indirect detection constraints. These PBH can produce gamma-ray signals from Hawking radiation that are within the sensitivity of future measurements by the AMEGO and e-ASTROGAM experiments. PBH which give rise to such observable gamma-ray signals have a cosmic origin from large primordial curvature fluctuations. There must then be a companion, stochastic gravitational wave (GW) background produced by the same curvature fluctuations. We demonstrate that the resulting GW signals will be well within the sensitivity of future detectors such as LISA, DECIGO, BBO, and the Einstein Telescope. The multi-messenger signal from the observed gamma-rays and GW will allow a precise measurement of the primordial curvature perturbation that produces the PBH. Indeed, we argue that the resulting correlation between the two types of observations can provide a smoking-gun signal of PBH.

Despite the impressive success of the standard cosmological model, several anomalies defy its triumph. Among them is the so-called lensing anomaly: the Planck satellite observes stronger CMB gravitational lensing than expected. The role of neutrinos in this anomaly has been mostly overlooked, despite their key role in CMB lensing, because in the standard scenario they tend to increase the tension. Here, we show that this strongly depends on the assumed neutrino equation of state. We demonstrate that if neutrinos have yet undiscovered long-range interactions, the lensing pattern is significantly affected, rendering the lensing anomaly as a pure statistical fluctuation. Our results thus open up a window to link anomalous CMB lensing with present and future cosmological, astrophysical, and laboratory measurements of neutrino properties.

Triangulum, M33, is a low mass, relatively undisturbed spiral galaxy that offers a new regime in which to test models of dynamical heating. In spite of its proximity, the dynamical heating history of M33 has not yet been well constrained. In this work, we present the TREX Survey, the largest stellar spectroscopic survey across the disk of M33. We present the stellar disk kinematics as a function of age to study the past and ongoing dynamical heating of M33. We measure line of sight velocities for ~4,500 disk stars. Using a subset, we divide the stars into broad age bins using Hubble Space Telescope and Canada-France-Hawaii-Telescope photometric catalogs: massive main sequence stars and helium burning stars (~80 Myr), intermediate mass asymptotic branch stars (~1 Gyr), and low mass red giant branch stars (~4 Gyr). We compare the stellar disk dynamics to that of the gas using existing HI, CO, and Halpha kinematics. We find that the disk of M33 has relatively low velocity dispersion (~16 km/s), and unlike in the Milky Way and Andromeda galaxies, there is no strong trend in velocity dispersion as a function of stellar age. The youngest disk stars are as dynamically hot as the oldest disk stars and are dynamically hotter than predicted by most M33 like low mass simulated analogs in Illustris. The velocity dispersion of the young stars is highly structured, with the large velocity dispersion fairly localized. The cause of this high velocity dispersion is not evident from the observations and simulated analogs presented here.

We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multi-zone galactic chemical evolution model. We adopt a metallicity-independent ("primary") N yield from massive stars and a metallicity-dependent ("secondary") N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]-[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By $t\approx5$ Gyr, the model approaches an equilibrium [N/O]-[O/H] relation, which traces the radial oxygen gradient. We find good agreement with the [N/O]-[O/H] trend observed in extra-galactic systems if we adopt an IMF-averaged massive star yield $y_\text{N}^\text{CC}=3.6\times10^{-4}$, consistent with predictions for rapidly rotating progenitors, and a fractional AGB yield that is linear in mass and metallicity $y_\text{N}^\text{AGB}=(9\times10^{-4})(M_*/M_\odot)(Z_*/Z_\odot)$. This model reproduces the [N/O]-[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]-[O/H] relation, but the AGB time-delay is required to match the APOGEE stellar age and [O/Fe] trends. If we add $\sim$40\% oscillations to the star formation rate, the model reproduces the scatter in gas-phase [N/O] vs. [O/H] observed in external galaxies by MaNGA. We also construct models using published AGB yields and examine their empirical successes and shortcomings. For all AGB yields we consider, simple stellar populations release half their N after only $\sim$250 Myr.

The massive stars that survive the phase of red supergiants (RSGs) spend the rest of their life in extremity. Their unstable atmospheres facilitate the formation and episodic ejection of shells that alter the stellar appearance and surroundings. In the present study, we revise the evolutionary state of eight hypergiants in the Magellanic Clouds, four of early-A type and four of FG type, and complement the short list of the eruptive post-RSGs termed as yellow hypergiants (YHGs). We refine the outdated temperatures and luminosities of the stars by means of high-resolution spectroscopy with FEROS. The A-type stars are suggested to be in their early, post-main sequence phase, showing spectrophotometric characteristics of redward evolving supergiants. On the other hand, the FG-type stars manifest themselves through the enhanced atmospheric activity that is traced by emission filling in H$\alpha$ and the dynamical modulation of the low-excitation BaII line. Of these stars, the dusty HD269723 is suggested to have recently departed from a cool phase. We identify double-peaked emission in the FEROS data of HD269953 that emerges from an orbiting disk-hosting companion. The highlight of the study is an episode of enhanced mass loss of HD271182 that manifests as a dimming event in the lightcurve and renders the star "modest" analogue to $\rho$ Cas. The luminosity $\log(L/L_{\odot})= 5.6$ of HD271182 can serve as an updated threshold for the luminosity of stars exhibiting a post-RSG evolution in the Large Magellanic Cloud.

We present the discovery of neutral gas detected in both damped Ly$\alpha$ absorption (DLA) and HI 21-cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an HI bridge connecting two interacting dwarf galaxies (log$(M_{\text{star}}/\text{M}_{\odot}) = 8.5 \pm 0.2$) that host a $z = 0.026$ DLA with log[$N$(HI)/cm$^{-2}$]$ = 20.60 \pm 0.05$ toward the QSO J2339-5523 ($z_{\text{QSO}} = 1.35$). At impact parameters of $d = 6$ and $33$ kpc, the dwarf galaxies have no companions more luminous than $\approx 0.05L_{*}$ within at least $\Delta v = \pm 300$ km s$^{-1}$ and $d \approx 350$ kpc. HI 21-cm emission is spatially coincident with the DLA at the 2-3$\sigma$ level per spectral channel over several adjacent beams. However, HI 21-cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with spin temperature to covering factor ratio $T_{s}/f_{c} > 1880$ K). VLT-MUSE observations demonstrate that the $\alpha$-element abundance of the ionized ISM is consistent with the DLA ($\approx 10$% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf-dwarf interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and HI 21-cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments.

Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood. To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe.

We study star formation driven outflows in a $z\sim0.02$ starbursting disk galaxy, IRAS08339+6517, using spatially resolved measurements from the Keck Cosmic Web Imager (KCWI). We develop a new method incorporating a multi-step process to determine whether an outflow should be fit in each spaxel, and then subsequently decompose the emission line into multiple components. We detect outflows ranging in velocity, $v_{\rm out}$, from $100-600$ km s$^{-1}$ across a range of star formation rate surface densities, $\Sigma_{\rm SFR}$, from $\sim$0.01-10 M$_\odot$ yr$^{-1}$ kpc$^{-2}$ in resolution elements of a few hundred parsec. Outflows are detected in $\sim100\%$ of all spaxels within the half-light radius, and $\sim70\%$ within $r_{90}$, suggestive of a high covering fraction for this starbursting disk galaxy. Around $2/3$ of the total outflowing mass originates from the star forming ring, which corresponds to $<10\%$ of the total area of the galaxy. We find that the relationship between $v_{\rm out}$ and the $\Sigma_{\rm SFR}$, as well as between the mass loading factor, $\eta$, and the $\Sigma_{\rm SFR}$, are consistent with trends expected from energy-driven feedback models. We study the resolution effects on this relationship and find stronger correlations above a re-binned size-scale of $\sim500$ pc. Conversely, we do not find statistically significant consistency with the prediction from momentum-driven winds.

Molecular outflows are often detected originating from both protostellar and extragalactic sources. Studies of low-mass, isolated high-mass, and extragalactic sources reveal scaling relations connecting the force carried by an outflow and the properties of the source that drives it. The aim of this work is to examine the effects of clustered star formation on the protostellar outflows and their scaling relations and to explore the possibility that outflows varying greatly scale and energetics are consistent with being launched by the same physical processes. High-angular resolution CO J = 3-2 observations were used of ten high-mass protostars in the Cygnus-X molecular cloud, obtained at the SMA as part of the PILS-Cygnus survey. From these data, the outflow force was measured. In addition, an extended sample of protostellar and extragalactic outflow-force measurements was assembled from existing literature to allow for a direct comparison of the scaling relations of the two types of outflows. Molecular outflows were detected originating from all ten sources of the PILS-Cygnus survey, and their outflow forces are found in close agreement with measurements from the literature. The comparison of the protostellar and extragalactic sources reveals, with 95\% confidence, that Class 0 protostars and extragalactic sources follow the same outflow force--bolometric luminosity correlation. The close agreement between the Cygnus-X sources and sources of similar envelope mass and bolometric luminosity suggests that clustered star formation has no significant effect on protostellar outflows. We find a strong indication that protostellar and extragalactic outflows are consistent with having a similar launch mechanism. The existence of such a mechanism would enable the development of a single universal outflow launch model, although more observations are required in order to verify this connection.

In this work, we test the robustness of the constancy of the Supernova absolute magnitude $M_B$ using Non-parametric Reconstruction Techniques (NRT). We isolate the luminosity distance parameter $d_L(z)$ from the Baryon Acoustic Oscillations (BAO) data set and cancel the expansion part from the observed distance modulus $\mu(z)$. Consequently, the degeneracy between the absolute magnitude with the Hubble constant $H_0$, is replaced by a degeneracy with the sound horizon at drag epoch $r_d$. When imposing the $r_d$ value, this yields the $M_B(z) = M_B + \delta M_B(z)$ value from NRT. We perform the respective reconstructions using the model independent Artificial Neural Network (ANN) technique and Gaussian processes (GP) regression. For the ANN we infer $M_B = -19.22\pm0.20$, and for the GP we get $M_B = -19.25\pm0.39$ when using the sound horizon from late time measurements. These estimations provide a $1\,\sigma$ possibility of a nuisance parameter presence $\delta M_B(z)$ at higher redshifts. We also tested different known nuisance models with the Markov Chain Monte Carlo (MCMC) technique which showed a strong preference for the constant model, but it was not possible not single out a best fit nuisance model.

A weak continuous line has been recently discovered onboard Voyager 1 in the interstellar medium, whose origin raised two major questions. First, how can this line be produced by plasma quasi-thermal noise on the Voyager short antenna? Second, why does this line emerge at some distance from the heliopause? We provide a simple answer to these questions, which elucidates the origin of this line. First, a minute quantity of supra-thermal electrons, as generally present in plasmas, whence the qualifier quasi-thermal, can produce a small plasma frequency peak on a short antenna, of amplitude independent of the concentration of these electrons; furthermore, the detection required long spectral averages, alleviating the smallness of the peak compared to the background. We therefore attribute the observed line to a minute proportion of fast electrons that contribute negligibly to the pressure. Second, we suggest that, up to some distance from the heliopause, the large compressive fluctuations ubiquitous in this region prevent the line to emerge from the statistical fluctuations of the receiver noise because it is blurred out by the averaging required for detection,especially in the presence of short-wavelength density fluctuations. These results open up novel perspectives for interstellar missions, by showing that a minute proportion of fast electrons may be sufficient to measure the density even with a relatively short antenna, because the quietness of the medium enables a large number of spectra to be averaged.

This paper describes a method by which the metrology system of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray space observatory, which uses two lasers to characterize the relative motion of the optics and focal plane benches, can be approximated should one laser fail. The two benches are separated by a ten-meter-long rigid mast that undergoes small amounts of thermal flexing which need to be compensated for in order to produce a non-blurred image. We analyze the trends of mast motion by archival observation parameters in order to discover whether the mast motion in future observations can be predicted. We find that, by using the solar aspect angle (SAA), observation date, and orbital phase, we can simulate the motion of one laser by translating the track produced by the other and applying modifications to the resulting mast aspect solution, allowing the reconstruction of a minimally distorted point spread function in most cases. We will implement the generation of simulated mast files alongside the usual NuSTAR data reduction pipeline for contingency purposes. This work has implications for reducing the risk of implementing laser metrology systems on future missions that use deployable masts to achieve the long focal lengths required in high-energy astronomy by mitigating the impact of a metrology laser failure in the extended phase of a mission.

In traditional lucky imaging (TLI), many consecutive images of the same scene are taken with a high frame-rate camera, and all but the sharpest images are discarded before constructing the final shift-and-add image. Here we present an alternative image analysis pipeline -- The Thresher -- for these kinds of data, based on online multi-frame blind deconvolution. It makes use of all available data to obtain a best estimate of the astronomical scene in the context of reasonable computational limits; it does not require prior estimates of the point-spread functions in the images, or knowledge of point sources in the scene that could provide such estimates. Most importantly, the scene it aims to return is the optimum of a justified scalar objective based on the likelihood function. Because it uses the full set of images in the stack, The Thresher outperforms TLI in signal-to-noise; as it accounts for the individual-frame PSFs, it does this without loss of angular resolution. We demonstrate the effectiveness of our algorithm on both simulated data and real Electron-Multiplying CCD images obtained at the Danish 1.54m telescope (hosted by ESO, La Silla). We also explore the current limitations of the algorithm, and find that for the choice of image model presented here, non-linearities in flux are introduced into the returned scene. Ongoing development of the software can be viewed at https://github.com/jah1994/TheThresher.

We present the first results of a multi-year program to map the orbits of M dwarf multiples within 25 parsecs. The observations were conducted primarily during 2019 - 2020 using speckle interferometry at the Southern Astrophysical Research (SOAR) Telescope in Chile, using the High-Resolution Camera mounted on the adaptive optics module (HRCam+SAM). The sample of nearby M dwarfs is drawn from three sources: multiples from the RECONS long-term astrometric monitoring program at the SMARTS 0.9m, known multiples for which these new observations will enable or improve orbit fits, and candidate multiples flagged by their astrometric fits in Gaia Data Release 2 (DR2). We surveyed 333 of our 338 M dwarfs via 830 speckle observations, detecting companions for 63% of the stars. Most notably, this includes new companions for 76% in the subset selected from Gaia DR2. In all, we report the first direct detections of 97 new stellar companions to the observed M dwarfs. Here we present the properties of those detections, the limits of each non-detection, and five orbits with periods 0.67 - 29 yr already observed as part of this program. Companions detected have projected separations of 0.024 - 2.0 arcsec (0.25 - 66 AU) from their primaries and have $\Delta I \lesssim 5.0$ mag. This multi-year campaign will ultimately map complete orbits for nearby M dwarfs with periods up to 3 yr, and provide key epochs to stretch orbital determinations for binaries to 30 yr.

The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [CII] 158 $\mu$m observations revealed an expanding bubble, the Veil shell, being powered by stellar winds and ionization feedback. We have identified a protrusion-like substructure in the Northwest portion of the Orion Veil Shell that may indicate additional feedback mechanisms that are highly directional. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. We use the [CII] 158 $\mu$m map of the Orion Nebula obtained with the upGREAT instrument onboard SOFIA. The spectral and spatial resolution of the observations are 0.3 km/s and 16 arcsec, respectively. We consider three possible origins for this protrusion: Fossil outflow cavities created by jets/outflows during the protostellar accretion phase, pre-existing clumpiness in the OMC-1 core, and the stellar wind during the main sequence phase. Based on the energetics and the morphology, we conclude that the northwestern part of the pre-existing cloud was locally perturbed by outflows ejected from massive protostars in the Trapezium cluster. This suggests that the protrusion of the Veil is the result of mechanical rather than radiative feedback. Furthermore, we argue that the location of the protrusion is a suitable place to break the Orion Veil owing to the photo-ablation from the walls of the protrusion. We conclude that the outflows of massive protostars can influence the morphology of the future \hii\,region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core pre-processed by the protostellar jet is important.

Impact crater counts on the Saturnian satellites are a key element for estimating their surface ages and placing constraints on their impactor population. The Cassini mission radar observations allowed crater counts to be made on the surface of Titan, revealing an unexpected scarcity of impact craters that show high levels of degradation. Following previous studies on impact cratering rates on the Saturnian satellites, we modeled the cratering process on Titan to constrain its surface chronology and to assess the role of centaur objects as its main impactors. A theoretical model previously developed was used to calculate the crater production on Titan, considering the centaur objects as the main impactors and including two different slopes for the size-frequency distribution (SFD) of the smaller members of their source population. A simple model for the atmospheric shielding effects is considered within the cratering process and our results are then compared with other synthetic crater distributions and updated observational crater counts. This comparison is then used to compute Titan's crater retention age for each crater diameter. The cumulative crater distribution produced by the SFD with a differential index of $s_2 = 3.5$ is found to consistently predict large craters (D > 50 km) on the surface of Titan, while it overestimates the number of smaller craters. As both the modeled and observed distributions flatten for craters $D \lesssim 25 $ km due to atmospheric shielding, the difference between them can be considered as a proxy for the scale to which erosion processes have acted on the surface of Titan throughout the Solar System age. Our results for the surface chronology of Titan indicate that craters with D > 50 km can prevail over the Solar System age, whereas smaller craters may be completely obliterated due to erosion processes acting globally.

Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.71\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs.

We present an analysis of XMM-Newton observations of four stars in the young (670 Myr) open cluster Praesepe. The planets hosted by these stars all lie close in radius-period space to the radius-period valley and/or the Neptunian desert, two features that photoevaporation by X-ray and extreme ultraviolet (EUV) photons could be driving. Although the stars are no longer in the saturated regime, strong X-ray and extreme ultraviolet irradiation is still ongoing. Based on EUV time evolution slopes we derived in a previous paper, in all four cases, two-thirds of their EUV irradiation is still to come. We compare the XMM-Newton light curves to those simultaneously measured with K2 at optical wavelengths, allowing us to search for correlated variability between the X-ray and optical light curves. We find that the X-ray flux decreases and flattens off while the optical flux rises throughout for K2-100, something that could result from active regions disappearing from view as the star spins. Finally, we also investigate possible futures for the four planets in our sample with simulations of their atmosphere evolution still to come, finding that complete photoevaporative stripping of the envelope of three of the four planets is possible, depending on the current planet masses.

Lava worlds belong to a class of short orbital period planets reaching dayside temperatures high enough to melt their silicate crust. Theory predicts that the resulting lava oceans outgas their volatile components, attaining equilibrium with the overlying vapour. This creates a tenuous, silicate-rich atmosphere that may be confined to the permanent dayside of the planet. The launch of JWST will provide the much needed sensitivity and spectral coverage to characterise these worlds. In this paper, we assess the observability of characterisable spectral features by self-consistently modelling silicate atmospheres for all the currently confirmed targets having sufficient substellar temperatures (> 1500 K). To achieve this we use outgassed equilibrium chemistry and radiative transfer methods to compute temperature-pressure profiles, atmospheric chemical compositions and emission spectra. We explore varying melt compositions, free of highly-volatile elements, accounting for possible atmospheric evolution. Our models include a large number of neutral and ionic species, as well as all up-to-date opacities. The results indicate that SiO and SiO2 infrared features are the best, unique identifiers of silicate atmospheres, detectable using MIRI instrument of JWST. Detection of these two species in emission would allow for strong constraints on atmospheric thermal structure and possibly the composition of the melt. We also propose that certain species, e.g., TiO, may be directly tied to different classes of melts, possibly revealing surface and interior dynamics. Currently, there are nearly a dozen confirmed lava planets ideal for characterisation of silicate atmospheres using JWST, with two of these already accepted for the initial General Observers program.

Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and non-linear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented due to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field due to a stronger Alfv\'{e}n speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfv\'{e}n speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D, but suggest that as one goes to weak reconnecting field (or strong guide field), the system becomes more 2.5D like and the 2.5D and 3D rates converge. These results have relevance to situations like nano-flare heating and flare current sheets in the corona.

The assembly of massive black holes in the early universe remains a poorly constrained open question in astrophysics. The merger and accretion of light seeds (remnants of Population III stars with mass below $\sim 1000M_{\odot}$) or heavy seeds (in the mass range $10^4-10^6 M_{\odot}$) could both explain the formation of massive black holes, but the abundance of seeds and their merging mechanism are highly uncertain. In the next decades, the gravitational-wave observatories coming online are expected to observe very high-redshift mergers, shedding light on the seeding of the first black holes. In this Letter we explore the potential and limitations for LISA, Cosmic Explorer and Einstein Telescope to constrain the mixture ratio of light and heavy seeds as well as the probability that central black holes in merging galaxies merge as well. Since the third generation ground-based gravitational-wave detectors will only observe light seed mergers, we demonstrate two scenarios in which the inference of the seed mixture ratio and merging probability can be limited. The synergy of multi-band gravitational-wave observations and electromagnetic observations will likely be necessary in order to fully characterize the process of high-redshift black hole formation.

The stellar cataclysms producing astronomical transients have long been modeled as either a point-like explosion or jet-like engine ignited at the center of a spherically symmetric star. However, many stars are observed, or are expected on theoretical grounds, not to be precisely spherically symmetric, but rather to have a slightly flattened geometry similar to that of an oblate spheroid. Here we present axisymmetric two-dimensional hydrodynamical simulations of the dynamics of point-like explosions initiated at the center of an aspherical massive star with a range of oblateness. We refer to these exploding aspherical stars as "ellipsars" in reference to the elliptical shape of the iso-density contours of their progenitors in the two-dimensional axisymmetric case. We find that ellipsars are capable of accelerating expanding rings of relativistic ejecta which may lead to the production of astronomical transients including low-luminosity GRBs, relativistic supernovae, and Fast Blue Optical Transients (FBOTs.)

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system.

In the manuscript, we discuss properties of the SDSS J1042-0018 which is a broad line AGN but mis-classified as a HII galaxy in the BPT diagram (SDSS J1042-0018 called as a mis-classified broad line AGN). The emission lines around H$\alpha$ and around H$\beta$ are well described by different model functions, considering broad Balmer lines to be described by Gaussian or Lorentz functions. Different model functions lead to different determined narrow emission line fluxes, but the different narrow emission line flux ratios lead the SDSS J1042-0018 as a HII galaxy in the BPT diagram. In order to explain the unique properties of the mis-classified broad line AGN SDSS J1042-0018, two methods are proposed, the starforming contributions and the compressed NLRs with high electron densities near to critical densities of forbidden emission lines. Fortunately, the strong starforming contributions can be preferred in the SDSS J1042-0018. The mis-classified broad line AGN SDSS J1042-0018, well explained by starforming contributions, could provide further clues on the applications of BPT diagrams to the normal broad line AGN.

The Nancy Grace Roman Space Telescope Coronagraph Instrument is a high-contrast imager, polarimeter, and spectrometer that will enable the study of exoplanets and circumstellar disks at visible wavelengths ($\sim$550--850~nm) at contrasts 2--3 orders of magnitude better than can currently be achieved by ground or space-based direct imaging facilities. To capitalize on this sensitivity, precise flux calibration will be required. The Roman Coronagraph, like other space-based missions, will use on-orbit flatfields to measure and correct for phenomena that impact the measured total effective throughput. However, the Coronagraph does not have internal lamp sources, therefore we have developed a method to perform flatfield calibrations using observations of extended sources, such as Uranus and Neptune, using a combination of rastering the Coronagraph's Fast Steering Mirror, tiling the planet across the field of view, and matched-filter image processing. Here we outline the process and present the results of simulations using images of Uranus and Neptune from the Hubble Space Telescopes Wide Field Camera 3, in filters approximate to the Coronagraph's Band 1 and Band 4. The simulations are performed over the Coronagraph's direct imaging and polarimetric modes. We model throughput effects in 3 different spatial frequency regimes including 1) high spatial frequency detector pixel-to-pixel quantum efficiency variations, 2) medium spatial frequency "measles" caused by particle deposition on the detector or other focal-plane optics post-launch, and 3) low spatial frequency detector fringing caused by self-interference due to internal reflections in the detector substrate as well as low spatial frequency vignetting at the edges of the Coronagraph's field of view. We show that Uranus and Neptune can be used as astrophysical flat sources with high precision ($\sim$0.5% relative error)

The ratio of nitrogen isotopes in the Martian atmosphere is a key constraint on the planet's atmospheric evolution. However, enrichment of the heavy isotope expected due to atmospheric loss from sputtering and photochemical processes is greater than measurements. A massive, multi-bar early CO2-dominated atmosphere and recent volcanic outgassing have been proposed to explain this discrepancy, and many previous models have assumed atmospheric nitrogen rapidly reached a steady state where loss to space balanced volcanic outgassing. Here we show using time-dependent models that the abundance and isotopic composition of nitrogen in the Martian atmosphere can be explained by a family of evolutionary scenarios in which the initial partial pressure of nitrogen is sufficiently high that a steady state is not reached and nitrogen levels gradually decline to present-day values over 4 billion years. Our solutions do not require a multi-bar early CO2 atmosphere and are consistent with volcanic outgassing indicated by both geologic mapping and the atmospheric 36Ar/38Ar ratio. Monte Carlo simulations that include these scenarios estimate that the partial pressure of N2 was 60 - 740 mbar (90% confidence, with a median value of 310 mbar) at 3.8 billion years ago when the valley networks formed. We suggest that such a high nitrogen partial pressure could have contributed substantially to warming on early Mars.

Changing-look Active Galactic Nuclei (CLAGN) are AGN that change type as their broad emission lines appear or disappear, which is usually accompanied by strong flux changes inn their blue featureless continuum. We search for Turn-On CLAGN by selecting type-2 AGN from the spectroscopic 6dF Galaxy Survey, whose photometry, as observed over a decade later by the SkyMapper Southern Survey, is consistent with type-1 AGN. Starting from a random sample of 235 known type-2 AGN we select 18 candidates and confirm 13 AGN to have changed into type-1 spectra; observations of an incomplete sample reveal nine further Turn-On CLAGN. While our search was not intended to reliably discover Turn-Off CLAGN, we discover four such cases as well. This result suggests a Turn-On CLAGN rate of 12% over ~15 years and imply a total CLAGN rate of ~25% over this period. Finally, we present observations of AGN that are atypical for the CLAGN phenomenology, including J1109146 - a CLAGN that did not appear as an AGN in 6dFGS; J1406507 - the second reported Changinglook NLS1; and J1340153 - a CLAGN with a change timescale of three months.

Thirteen yr observation data of 4FGL J0824.9+3915 with the Large Area Telescope on board the Fermi Gamma Ray Space Telescope (Fermi/LAT) are analyzed for revisiting whether 4C +39.23B, a compact steep-spectrum (CSS) source closed to a flat-spectrum radio quasar (FSRQ) 4C +39.23A in the $\gamma$-ray emitting region of 4FGL J0824.9+3915, is a $\gamma$-ray emitter. We find that the time-integrated $\gamma$-ray emission of 4FGL J0824.9+3915 is overwhelmingly dominated by 4C +39.23A. It shows significant variability at a 6.7$\sigma$ confidence level and the average $\gamma$-ray flux in the 0.1--300 GeV energy band is $(1.60\pm 0.15)\times10^{-8}$ ph cm$^{-2}$ s$^{-1}$ with a power-law photon spectral index of $2.48\pm0.05$. During MJD 57500--58500, 4FGL J0824.9+3915 is in a low state with a steady $\gamma$-ray flux. Analyzed the Fermi/LAT observation data in this time interval, it is found that the TS values of the $\gamma$-ray emission from 4C +39.23A and 4C +39.23B are $\sim5$ and $\sim 31$, respectively, indicating that the $\gamma$-ray emission in this time interval is dominated by the CSS 4C +39.23B. The derived average flux in this time interval for 4C +39.23B is $(9.40\pm4.10)\times 10^{-9}$ ph cm $^{-2}$ s$^{-1}$ with $\Gamma_{\gamma}=2.45\pm0.17$. Attributing the spectral energy distribution (SED) of 4C +39.23B to the radiations from its core and extended region, we show that the SED can be represented with a two-zone leptonic model. Its $\gamma$-ray emission is contributed by the core region. The derived magnetic field strength and Doppler boosting factor of the core are 0.13 G and 6.5. Comparing 4C +39.23B with other $\gamma$-emitting CSSs and compact symmetric objects (CSOs) in the $\Gamma_{\gamma}-L_{\gamma}$ plane, it resembles CSSs.

The last decade has seen a revolution in our knowledge of the Galaxy thanks to the Gaia and asteroseismic space missions and the ground-based spectroscopic surveys. To complete this picture, it is necessary to map the ages of its stellar populations. During recent years, the dependence on time of abundance ratios involving slow (s) neutron-capture and $\alpha$ elements (called chemical-clocks) has been used to provide estimates of stellar ages, usually in a limited volume close to the Sun. We aim to analyse the relations of chemical clocks in the Galactic disc extending the range to R$_{\rm GC}\sim$6-20~kpc. Using the sixth internal data release of the Gaia-ESO survey, we calibrated several relations between stellar ages and abundance ratios [s/$\alpha$] using a sample of open clusters, the largest one so far used with this aim. Thanks to their wide galactocentric coverage, we investigated the radial variations of the shape of these relations, confirming their non-universality. We estimated our accuracy and precision in recovering the global ages of open clusters, and the ages of their individual members. We applied the multi-variate relations with the highest correlation coefficients to the field star population. We confirm that there is no single age-chemical clock relationship valid for the whole disc, but that there is a dependence on the galactocentric position, which is related to the radial variation of the star formation history combined with the non-monotonic dependence on metallicity of the yields of the s-process elements from low- and intermediate-mass stars. Finally, the abundance ratios [Ba/$\alpha$] are more sensitive to age than those with [Y/$\alpha$] for young disc stars, and their slopes vary less with galactocentric distance.

We argue that double inflation may occur when a spectator field is non-minimally coupled to gravity. As a concrete example, we study a two-field inflationary model where the initial spectator field is non-minimally coupled to gravity while the initial inflaton field is minimally coupled. The non-minimal coupling results in the growth of the spectator field which, in turn, drives the second stage of inflation in a significant region of parameter space. The isocurvature fluctuations originating from the spectator field source adiabatic ones, and hence the spectator non-minimal coupling can modify the inflationary predictions for the spectral index and the tensor-to-scalar ratio even though the initial inflaton field is minimally coupled to gravity. We explicitly show that quadratic chaotic inflation can become viable by the introduction of the spectator non-minimal coupling.

We report the results of AstroSat observations of Cygnus X-2 during February 2016. The source's power density spectrum generated using LAXPC data revealed the presence of a prominent Quasi-periodic Oscillation (QPO) at $\sim42$ Hz with broadband continuum noise at lower frequencies at $\sim10$ Hz. The large effective area of LAXPC at $\gtrsim$30 keV allowed for an unprecedented study of the energy dependence of the QPO and the broad noise continuum. The fractional r.m.s increases with energy, and its shape is similar for both the QPO and the continuum noise, suggesting a common radiative origin. However, while the QPO exhibits hard time lags, with the high energy photons lagging the low ones by a few milliseconds, the continuum noise shows the opposite behavior. The photon spectrum from SXT and LAXPC in $0.7-30$ keV band comprises the soft component from a disc and a hard Comptonized component from a hot corona. While the energy dependence of the r.m.s shows that the QPO and the continuum noise variability are dominated by the Comptonized component, the change in sign of the time-lag suggests that the dynamic origin of the QPO may be in the disk while the noise continuum may originate from the corona.

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) which were both on board the Hubble Space Telescope (HST). The FUV-UV colour-magnitude diagram (CMD) shows a main sequence (MS) turnoff at FUV $\approx$ 22 mag and FUV-UV $\approx$ 3 mag. The MS extends 4 mag below the turnoff, and a prominent horizontal branch (HB) and blue straggler (BS) sequence can be seen. A total of 1218 MS stars, 185 red giant branch stars, 47 BS stars and 41 HB stars are identified, along with 78 sources blueward of the MS which consist of white dwarfs (WDs) and objects in the gap between the WDs and the MS that include potential cataclysmic variable (CV) candidates. The radial distribution of the BS population is concentrated towards the cluster centre, indicating that mass segregation has occurred. The blue and red sub-populations of the double BS sequence appear mixed in the ultraviolet CMD, and no significant central concentration of CV candidates is seen in this cluster.

A complete description is given of two-mirror telescopes with a flat medial focal surface, on which the images of stars are circles of least confusion. Particular attention is paid to aplanats, since their field of view is noticeably larger than that of classical systems. Two sets of appropriate solutions correspond to Schwarzschild and Gregorian telescopes. As a result, it becomes possible to use flat light detectors with wide-field two-mirror telescopes. New designs are of particular interest when as few reflective surfaces as possible are required, which is typical for space exploration and non-optical observations.

We select 456 galaxies with kinematically misaligned gas and stellar components from 9546 parent galaxies in MaNGA, and classify them into 72 star-forming galaxies, 142 green-valley galaxies and 242 quiescent galaxies. Comparing the spatial resolved properties of the misaligned galaxies with control samples closely match in the D$_n$4000 and stellar velocity dispersion, we find that: (1) the misaligned galaxies have lower values in $V_{\rm gas}/{\sigma}_{\rm gas}$ and $V_{\rm star}/{\sigma}_{\rm star}$ (the ratio between ordered to random motion of gas and stellar components) across the entire galaxies than their control samples; (2) the star-forming and green-valley misaligned galaxies have enhanced central concentrated star formation than their control galaxies. The difference in stellar population between quiescent misaligned galaxies and control samples is small; (3) gas-phase metallicity of the green valley and quiescent misaligned galaxies are lower than the control samples. For the star forming misaligned galaxies, the difference in metallicity between the misaligned galaxies and their control samples strongly depends on how we select the control samples. All these observational results suggest external gas accretion influences the evolution of star forming and green valley galaxies, not only in kinematics/morphologies, but also in stellar populations. However, the quiescent misaligned galaxies have survived from different formation mechanisms.

We suggest that the recently discovered, enigmatic pulsar with a period of 18.18 minutes, GLEAM-X J162759.5-523504.3, is most likely a hot subdwarf (proto white dwarf). A magnetic dipole model explains the observed period and period-derivative for a highly magnetized ($\sim 10^8$G), hot subdwarf of typical mass $\sim 0.5M_\odot$ and radius $\sim 0.3R_\odot$, and an age of $\sim 3\times 10^4$yr. The subdwarf spin is close to its breakup speed and its spindown luminosity is near its Eddington limit, likely as a result of accretion from a companion.

Molecules in space are synthesized via a large variety of gas-phase reactions, and reactions on dust-grain surfaces, where the surface acts as a catalyst. Especially, saturated, hydrogen-rich molecules are formed through surface chemistry. Astrochemical models have developed over the decades to understand the molecular processes in the interstellar medium, taking into account grain surface chemistry. However, essential input information for gas-grain models, such as binding energies of molecules to the surface, have been derived experimentally only for a handful of species, leaving hundreds of species with highly uncertain estimates. Moreover, some fundamental processes are not well enough constrained to implement these into the models. The proceedings gives three examples how computational chemistry techniques can help answer fundamental questions regarding grain surface chemistry.

In the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068, it has been found that the molecular gas in the Circum-nuclear Disk (CND) is outflowing, which is a manifestation of ongoing AGN feedback. The outflowing gas has a large spread of velocities, which likely drive different shock chemistry signatures at different locations in the CND. We perform a multi-line molecular study using two shock tracers, SiO and HNCO, with the aim to determine the gas properties traced by these two species, and explore the possibility of reconstructing the shock history in the CND. Five SiO transitions and three HNCO transitions were imaged at high resolution $0''.5-0''.8$ with the Atacama Large Millimeter/submillimeter Array (ALMA). We performed both LTE and non-LTE radiative transfer analysis coupled with Bayesian inference process in order to characterize the gas properties, such as molecular gas density and gas temperature. We found clear evidence of chemical differentiation between SiO and HNCO, with the SiO/HNCO ratio ranging from greater than one on the east of CND to lower than one on the west side. The non-LTE radiative transfer analysis coupled with Bayesian inference confirms that the gas traced by SiO has different densities - and possibly temperatures - than that traced by HNCO. We find that SiO traces gas affected by fast shocks while the gas traced by HNCO is either just affected by slow shocks or not shocked at all. A distinct differentiation between SiO and HNCO has been revealed in our observations and the further analysis of the gas properties traced by both species, which confirms the results from previous chemical modelings.

The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies.

Our aim is to test the accretion scenario of lambda Boo stars. This model predicts that a binary system with two early-type stars passing through a diffuse cloud should both display the same superficial peculiarity. We carried out a detailed abundance determination of three multiple systems hosting a candidate lambda Boo star: the remarkable triple system HD 15164/65/65C and the binary systems HD 193256/281 and HD 198160/161. The abundance analysis of HD 15164/65/65C shows a clear lambda Boo object (HD 15165) and two objects with near solar composition (HD 15164 and 15165C). Notably, the presence of a lambda Boo star (HD 15165) together with a near solar early-type object (HD 15164) is difficult to explain under the accretion scenario. Also, the solar-like composition derived for the late-type star of the system (HD 15165C) could be used, for the first time, as a proxy for the initial composition of the lambda Boo stars. Then, by reviewing abundance analysis of all known binary systems with candidate lambda Boo stars from literature and including the systems analyzed here, we find no binary/multiple system having two clear "bonafide" lambda Boo stars, as expected from the accretion scenario. The closer candidates to show two lambda Boo-like stars are HD 84948, HD 171948 and HD 198160; however, in our opinion they show mild rather than clear lambda Boo patterns. Our results brings little support to the accretion scenario. Then, there is an urgent need of additional binary and multiple systemsto be analyzed through a detailed abundance analysis.[abridged]

In this paper, we present results from long-term $V(RI)_{c}$ photometric observations of the pre-main-sequence stars V2764 Ori and LkH$\alpha$ 301, located in the field of the McNeil's Nebula within the Orion star-forming complex. Our observations were performed in the period from August 2004 to November 2021 with three telescopes and seven different types of CCD cameras. Photometric observations, especially concerning the long-term behavior of the stars, are missing in the literature. We present the first photometric monitoring for them, which cover 17 years. Our data indicate that the variability of both stars is typical for classical T Tauri stars.

Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was initiated. In its first iteration in June 2021, ten rare event search collaborations contributed to this initiative via talks and discussions. The contributing collaborations were CONNIE, CRESST, DAMIC, EDELWEISS, MINER, NEWS-G, NUCLEUS, RICOCHET, SENSEI and SuperCDMS. They presented data about their observed energy spectra and known backgrounds together with details about the respective measurements. In this paper, we summarize the presented information and give a comprehensive overview of the similarities and differences between the distinct measurements. The provided data is furthermore publicly available on the workshop's data repository together with a plotting tool for visualization.

Dust plays a pivotal role in determining the observed spectral energy distribution (SED) of galaxies. Yet our understanding of dust attenuation is limited and our observations suffer from the dust-metallicity-age degeneracy in SED fitting (single galaxies), large individual variances (ensemble measurements), and the difficulty in properly dealing with uncertainties (statistical considerations). In this study, we create a population Bayesian model to rigorously account for correlated variables and non-Gaussian error distributions and demonstrate the improvement over a simple Bayesian model. We employ a flexible 5-D linear interpolation model for the parameters that control dust attenuation curves as a function of stellar mass, star formation rate (SFR), metallicity, redshift, and inclination. Our setup allows us to determine the complex relationships between dust attenuation and these galaxy properties simultaneously. Using Prospector fits of nearly 30,000 3D-HST galaxies, we find that the attenuation slope ($n$) flattens with increasing optical depth ($\tau$), though less so than in previous studies. $\tau$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $\tau$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $\tau$ than face-on galaxies, but only for $\log~M_*\gtrsim 10$, reflecting the lack of triaxiality for low-mass galaxies. Redshift evolution of dust attenuation is strongest for low-mass, low-SFR galaxies, with higher optical depths but flatter curves at high redshift. Finally, $n$ has a complex relationship with stellar mass, highlighting the intricacies of the star-dust geometry. We have publicly released software (https://github.com/Astropianist/DustE) for users to access our population model.

Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system's inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555-2908. By modelling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is $397\pm4$ km s$^{-1}$ (errors quoted at 95\% confidence interval), as well as a binary inclination of $>75^{\rm o}$. Combined with $\gamma$-ray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{-0.09}$ M$_\odot$ and a companion mass of 0.060$^{+0.005}_{-0.003}$ M$_\odot$, placing PSR J1555-2908 at the observed upper limit of what is considered a black widow system.

We revise the dynamical properties of a class of cosmological models where the dark sector interacts through an interacting term that changes sign during evolution. In particular, we obtain the critical points and we investigate the existence and stability conditions for cosmological solutions, describing radiation, matter and dark energy dominated eras. We find that all the studied models admit a stable critical point corresponding to an accelerated phase. We use background data to find the best fit parameters for one of the studied models, resulting an interacting parameter with a definite sign within $1\sigma$ confidence level, consistent with the results of the dynamical system analysis. We also compute the statefinder parameters and plot the $r-q$ and $r-s$ planes, where we observe different trajectories when we vary the interaction parameter for a specific model and when we vary the interacting scenario. We can in this sense distinguish among models, including $\Lambda$CDM.

We provide an overview of the observational properties of the four major Taurid showers, namely the Northern and Southern Taurids (#017 NTA and #002 STA), the Beta Taurids (#173 BTA) and the Zeta Perseids (#172 ZPE). Analysing more than two decades of meteor observations from visual, optical and radar measurements we present the Taurids average activity, annual variations in strength, radiant drift and orbital variations as a function of solar longitude and particle size. The Taurid showers are detected over several weeks in the spring and autumn, but their annual activity level is generally low (less than 15 visual meteors per hour). We find the STA to be predominant in autumn, while its twin the ZPE dominates over the BTA in spring. Due to their long duration, the position of each shower's radiant and orbital elements are variable with time. Optical measurements have previously recorded enhanced STA activity and increased fireball rates caused by the return of a swarm of meteoroids trapped in the 7:2 mean motion resonance with Jupiter. However, we find no presence of the swarm in radar data, suggesting that small meteoroids are removed from the resonance faster than fireball-producing meteoroids. We also find the STA to be enriched in smaller particles early in their activity period. The differences we identify in our analysis between the showers at different particle sizes provide strong observational constraints to future dynamical modelling of the Taurid Meteoroid Complex.

We propose a precessing transient magnetar model for the recently discovered radio source GLEAM-X J162759.5-523504.3. We identify the observed period of $\sim 1$ ks as the precession period of the magnetar deformed due to its strong ($B_{\phi} \sim 10^{16}$ G) toroidal field. The resulting deformation of order $10^{-4}$ implies a spin period of $P_{\rm s} = 0.1$ s. Assuming a strong dipole field of $B_{\rm d} \sim 10^{14}$ G we predict a period derivative of $\dot{P}_{\rm s}\sim 10^{-11}$ s/s. We also predict that the precession period of the magnetar can be observed in the hard X-ray band just as the other three galactic magnetars exhibit precession.

We study the tidal response of a superfluid neutron star in a binary system, focussing on Newtonian models with superfluid neutrons present throughout the star's core and the inner crust. Within the two-fluid formalism, we consider the main aspects that arise from the presence of different regions inside the star, with particular focus on the various interfaces. Having established the relevant theory, we determine the tidal excitation of the most relevant oscillation modes during binary inspiral. Our results suggest that superfluid physics has a negligible impact on the static tidal deformation. The overwhelming contribution to the Love number is given by, as for normal matter stars, the ordinary fundamental mode (f-mode). Strong entrainment, here described by a phenomenological expression which mimics the large effective neutron mass expected at the bottom of the crust, is shown to have significant impact on the superfluid modes, but our results for the dynamical tide are nevertheless similar to the static limit: the fundamental modes are the ones most significantly excited by the tidal interaction, with the ordinary f-mode dominating the superfluid one. We also discuss the strain built up in the star's crust during binary inspiral, showing that the superfluid f-mode may (depending on entrainment) reach the limit where the crust breaks, although it does so after the ordinary f-mode. Overall, our results suggest that the presence of superfluidity may be difficult to establish from binary neutron star gravitational-wave signals.

Magnetars are X-ray emitting neutron stars that are deemed too luminous to be spin-powered. The prevailing view that magnetars' X-ray luminosities exceed their spin-down luminosities is however based on the assumption that the decay with distance of the flux of the X-rays received from magnetars obeys the inverse-square law. The results presented here, of testing the hypothesis of independence of luminosities and distances of magnetars by means of the Efron-Petrosian statistic, imply that the observational data in the McGill Magnetar Catalog are consistent with the dependence $S\propto D^{-3/2}$ of the flux densities $S$ of these objects on their distances $D$ at substantially higher levels of significance than they are with the dependence $S\propto D^{-2}$. This is not incompatible with the requirements of the conservation of energy because the radiation process described in Ardavan [ Mon. Not. R. Astron. Soc., 507, 4530-4563 (2021)], by which the superluminally moving current sheet in the magnetosphere of a neutron star is shown to generate the observed X-ray pulses, is intrinsically transient: the difference in the fluxes of power across any two spheres centred on the star is balanced by the change with time of the energy contained inside the shell bounded by those spheres. Once their over-estimation is rectified, the ratios of X-ray to spin-down luminosities of known magnetars turn out to be invariably lower than $1$. A magnetar differs from a spin-powered pulsar only in that it is observed along a privileged latitudinal direction: the closer is the line of sight to a direction in which the radiation from the current sheet is focused, the higher the frequency and the lower the decay rate with distance of the observed radiation.

Gravitational-wave (GW) detections of merging neutron star-black hole (NSBH) systems probe astrophysical neutron star (NS) and black hole (BH) mass distributions, especially at the transition between NS and BH masses. Of particular interest are the maximum NS mass, minimum BH mass, and potential mass gap between them. While previous GW population analyses assumed all NSs obey the same maximum mass, if rapidly spinning NSs exist, they can extend to larger maximum masses than nonspinning NSs. In fact, several authors have proposed that the $\sim2.6\,M_\odot$ object in the event GW190814 -- either the most massive NS or least massive BH observed to date -- is a rapidly spinning NS. We therefore infer the NSBH mass distribution jointly with the NS spin distribution, modeling the NS maximum mass as a function of spin. Using 4 LIGO-Virgo NSBH events including GW190814, if we assume that the NS spin distribution is uniformly distributed up to the maximum (breakup) spin, we infer the maximum non-spinning NS mass is $2.7^{+0.5}_{-0.4}\,M_\odot$ (90\% credibility), while assuming only nonspinning NSs, the NS maximum mass must be $>2.53 M_\odot$ (90\% credibility). The data support the mass gap's existence, with a minimum BH mass at $5.4^{+0.7}_{-1.0} M_\odot$. With future observations, under simplified assumptions, 150 NSBH events may constrain the maximum nonspinning NS mass to $\pm0.02\,M_\odot$, and we may even measure the relation between the NS spin and maximum mass entirely from GW data. If rapidly rotating NSs exist, their spins and masses must be modeled simultaneously to avoid biasing the NS maximum mass.

Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough evaluation of the presence of a third low-mass planetary companion, which started emerging during a previous campaign. Radial velocities (RVs) were calculated using both a cross-correlation function (CCF) and a template matching approach. The RV analysis includes a component to model Proxima's activity using a Gaussian process (GP). We use the CCF's full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit.

The energy emitted by active galactic nuclei (AGN) may provide a self-regulating process (AGN feedback) that shapes the evolution of galaxies. This is believed to operate along two modes: on galactic scales by clearing the interstellar medium via outflows, and on circumgalactic scales by preventing the cooling and accretion of gas onto the host galaxy. Radio jets associated with radiatively-inefficient AGN are known to contribute to the latter mode of feedback. However, such jets could also play a role on circum-nuclear and galactic scales, blurring the distinction between the two modes. We have discovered a spatially-resolved, massive molecular outflow, carrying $\sim75\%$ of the gas in the central region of the host galaxy of a radiatively-inefficient AGN. The outflow coincides with the radio jet 540 pc offset from the core, unambiguously pointing to the jet as the driver of this phenomenon. The modest luminosity of the radio source ($L\rm_{1.4 GHz}=2.1 \times 10\rm^{23}~\rm W~\rm Hz^{-1}$) confirms predictions of simulations that jets of low-luminosity radio sources carry enough power to drive such outflows. Including kpc-scale feedback from such sources -- comprising of the majority of the radio AGN population -- in cosmological simulations may assist in resolving some of their limitations.

We compute the evolution of the grain size distribution (GSD) in a suite of numerical simulations of an isolated Milky-Way-like galaxy using the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4}. The full GSD is sampled on a logarithmically spaced grid with 30 bins, and its evolution is calculated self-consistently with the hydrodynamical and chemical evolution of the galaxy using a state-of-the-art star formation and feedback model. In previous versions of this model, the GSD tended to be slightly biased towards larger grains and the extinction curve had a tendency to be flatter than the observations. This work addresses these issues by considering the diffusion of dust and metals through turbulence on subgrid scales and introducing a multi-phase subgrid model that enables a smoother transition from diffuse to dense gas. We show that diffusion can significantly enhance the production of small grains and improve the agreement with the observed dust extinction curve in the Milky Way.

Accretion disks surrounding compact objects, and other environmental factors, deviate satellites from geodesic motion. Unfortunately, setting up the equations of motion for such relativistic trajectories is not as simple as in Newtonian mechanics. Here, we propose a simple method aimed at generating physically accurate covariant 4-forces. We apply this method to several conservative and dissipative forces. In particular, we compute the drag due to gravitational and hard-sphere collisions in dust, gas and radiation media. We recover and covariantly extend known forces such as Epstein drag, Chandrasekhar's dynamical friction and Poynting-Robertson drag. Variable-mass effects are also considered, namely Hoyle-Lyttleton accretion and the variable-mass rocket. We conclude with two applications: 1. The free-falling spring. We find that Hooke's law amounts to an effective Anti-de Sitter tidal force; 2. Black hole infall with drag. We numerically compute some trajectories on a Schwarzschild background supporting a dust-like accretion disk.

We study the phenomenology of the Higgs-Dilaton model in the context of Einstein-Cartan gravity, focusing on the separate impact of the Holst and Nieh-Yan terms on the inflationary observables. Using analytical and numerical techniques, we show the predictions of these scenarios to display an attractor-like behaviour intrinsically related to the curvature of the field-space manifold in the metric formulation of the theory. Beyond that, the analysis of the Nieh-Yan case reveals the existence of an additional attractor solution induced by a cubic pole in the inflaton kinetic term that becomes relevant at large dilaton couplings. This constitutes a unique feature of the Einstein-Cartan formulation as compared to the metric and Palatini counterparts.

The nature of neutrinos, whether Dirac or Majorana, is hitherto not known. Assuming that the neutrinos are Dirac, which needs $B-L$ to be an exact symmetry, we make an attempt to explain the observed proportionality between the relic densities of dark matter (DM) and baryonic matter in the present Universe ${\it i.e.,}\,\, \Omega_{\rm DM} \approx 5\, \Omega_{\rm B}$. Assuming the existence of heavy $SU(2)_L$ scalar doublet $(X= (X^0, X^-)^T)$ in the early Universe, an equal and opposite $B-L$ asymmetry can be generated in left and right-handed sectors by the CP-violating out-of-equilibrium decay $X^0 \to \nu_L \nu_R$ since $B-L$ is an exact symmetry. We ensure that $\nu_L-\nu_R$ equilibration does not occur until below the electroweak (EW) phase transition during which a part of the lepton asymmetry gets converted to dark matter asymmetry through a dimension eight operator, which conserves $B-L$ symmetry and is in thermal equilibrium. The remaining $B-L$ asymmetry then gets converted to a net B-asymmetry through EW-sphalerons which are active at a temperature above 100 GeV. To alleviate the small-scale anomalies of $\Lambda$CDM, we assume the DM to be self-interacting via a light mediator, which not only depletes the symmetric component of the DM, but also paves a way to detect the DM at terrestrial laboratories through scalar portal mixing.

In this contribution, we summarize our recent studies on the chiral invariant mass and the chiral condensates in neutron star matter. We construct a unified equations of state assuming the crossover phase transition from hadronic matter described by a parity doublet model to quark matter by an Nambu--Jona-Lasinio type quark model. We first show that the chiral invariant mass is constrained to be 600 MeV $\lesssim m_0 \lesssim$ 900 MeV from recent observations of neutron stars. We then determine the density dependence of the chiral condensate in the crossover description, and show that the chiral condensates are actually smoothly connected from the hadronic matter where the change is driven by the positive chiral scalar charge in a nucleon, to the quark matter where the change is by the modification of the quark Dirac sea, reflecting the hadron-quark crossover.

We discuss cosmology based on a cuscuta-galileon gravity theory, which preserves just two degrees of freedom. Although there exists no additional degrees of freedom, introduction of a potential of a scalar field changes the dynamics. The scalar field is completely determined by matter fields. Giving an exponential potential as an example, we discuss the cosmological dynamics. The gravitational "constant" $G_{\rm F}$ appeared in the effective Friedmann equation becomes time dependent. We also present how to construct a potential when we know the evolution of the Hubble parameter. When we assume the $\Lambda$CDM cosmology for the background evolution, we find the potential form. We then analyze the density perturbations, which equation is characterized only by a change of the gravitational "constant"$G_{\rm eff}$, which also becomes time dependent. From the observational constraints such as the constraint from the big-bang nucleosynthesis and the constraint on time-variation of gravitational constant, we restrict the parameters in our models. The time dependence of the gravitational constant in the effective Friedmann equation, we may have a chance to explain the Hubble tension problem.

In recent years, deep learning techniques have been introduced into the field of trajectory optimization to improve convergence and speed. Training such models requires large trajectory datasets. However, the convergence of low thrust (LT) optimizations is unpredictable before the optimization process ends. For randomly initialized low thrust transfer data generation, most of the computation power will be wasted on optimizing infeasible low thrust transfers, which leads to an inefficient data generation process. This work proposes a deep neural network (DNN) classifier to accurately identify feasible LT transfer prior to the optimization process. The DNN-classifier achieves an overall accuracy of 97.9%, which has the best performance among the tested algorithms. The accurate low-thrust trajectory feasibility identification can avoid optimization on undesired samples, so that the majority of the optimized samples are LT trajectories that converge. This technique enables efficient dataset generation for different mission scenarios with different spacecraft configurations.

We study 1-loop corrections to the primordial stochastic background of gravitational waves produced during inflation. While in single-clock, at the leading order in slow-roll, quantum corrections keep the amplitude scale-free this is not the case when the pattern of symmetry breaking is different. In particular, when spatial diffeomorphisms are also broken during inflation as for solid inflation, a log-running in the external momentum is generated. We relate the appearance of a log-running to the spontaneous breaking of dilatation invariance. The running could be instrumental to distinguish single-clock from alternative models of inflation in future high sensitivity CMB polarisation and GWs experiments.

Deep learning methods have gained popularity in high energy physics for fast modeling of particle showers in detectors. Detailed simulation frameworks such as the gold standard Geant4 are computationally intensive, and current deep generative architectures work on discretized, lower resolution versions of the detailed simulation. The development of models that work at higher spatial resolutions is currently hindered by the complexity of the full simulation data, and by the lack of simpler, more interpretable benchmarks. Our contribution is SUPA, the SUrrogate PArticle propagation simulator, an algorithm and software package for generating data by simulating simplified particle propagation, scattering and shower development in matter. The generation is extremely fast and easy to use compared to Geant4, but still exhibits the key characteristics and challenges of the detailed simulation. We support this claim experimentally by showing that performance of generative models on data from our simulator reflects the performance on a dataset generated with Geant4. The proposed simulator generates thousands of particle showers per second on a desktop machine, a speed up of up to 6 orders of magnitudes over Geant4, and stores detailed geometric information about the shower propagation. SUPA provides much greater flexibility for setting initial conditions and defining multiple benchmarks for the development of models. Moreover, interpreting particle showers as point clouds creates a connection to geometric machine learning and provides challenging and fundamentally new datasets for the field. The code for SUPA is available at https://github.com/itsdaniele/SUPA.

In this contribution, we present for the first time a scenario according to which early quark deconfinement in compact stars is triggered by the Bose-Einstein condensation (BEC) of a light sexaquark (S) with a mass $m_S<2054$ MeV, that has been suggested as a candidate particle to explain the baryonic dark matter in the Universe. The onset of S BEC marks the maximum mass of hadronic neutron stars and it occurs when the condition for the baryon chemical potential $\mu=m_S/2$ is fulfilled in the center of the star, corresponding to $M_{\rm onset}\lesssim 0.7~M_\odot$. In the gravitational field of the star the density of the BEC of the S increases until a new state of the matter is attained, where each of the S-states got dissociated into a triplet of color-flavor-locked (CFL) diquark states. These diquarks are the Cooper pairs in the color superconducting CFL phase of quark matter, so that the developed scenario corresponds to a Bose-Einstein condensation - Bardeen-Cooper-Schrieffer (BEC-BCS) transition in strongly interacting matter. For the description of the CFL phase, we develop here for the first time the three-flavor extension of the density-functional formulation of a chirally symmetric Lagrangian model of quark matter where confining properties are encoded in a divergence of the scalar self-energy at low densities and temperatures.

The advent of gravitational-wave (GW) astronomy has presented us with a completely new means for observing the Universe, allowing us to probe its structure and evolution like never before. In this thesis, we explore three distinct but complementary avenues for using GW observations to gain new insights into cosmology and fundamental physics. In chapter 1, we study the astrophysical GW background (AGWB): the cumulative GW signal arising from a large number of compact binary coalescences (CBCs) throughout the Universe. Since these compact binaries reside in galaxies, the AGWB contains anisotropies that trace out the large-scale structure of the cosmic matter distribution. We investigate the angular power spectrum of the AGWB, with the goal of developing predictions that can be confronted with directional AGWB searches. In chapter 2, we calculate the nonlinear GW memory emitted by cusps and kinks on cosmic string loops, which are among the most promising cosmological sources of GWs. We show that, surprisingly, the cusp memory signal diverges for sufficiently large loops, indicating a breakdown in the validity of the weak-field description of the cusp. We then present one tentative possible solution to this divergence, in which the portion of the string surrounding the cusp collapses to form a primordial black hole (PBH). Finally, in chapter 3 we develop a powerful new method for GW detection based on precision measurements of the orbits of binary systems. In the presence of a stochastic GW background (GWB) the trajectories of the binary's components are perturbed, giving rise to a random walk in the system's orbital parameters over time. We calculate the sensitivity of binary pulsars and lunar laser ranging to the GWB through this effect, and show that present data are already sensitive enough to place the strongest constraints to date in the $\mu$Hz frequency band.

Nuclear matter properties based on a relativistic approach suitable for the description of multi-component systems are calculated. We use a set of nuclear relativistic mean-field models that satisfy acceptable nuclear matter properties and neutron star observations. The effects of the density dependence of the symmetry energy and of the Landau quantization due to the presence of a strong external magnetic field are discussed. Properties such as the proton fraction, the Landau mass, Landau parameters and entrainment matrix, the adiabatic index and speed of sound are calculated for cold $\beta$-equilibrium matter. A large dispersion on the calculated properties is obtained at two to three times saturation density $\rho_0 $. The proton Landau mass can be as low as one third of the vacuum nucleon mass at 2-3$~\rho_0 $. Similar effects are obtained for the Landau parameters, in particular, the ones involving protons, where the relative dispersion of $F^0_{pp}$ and $F^1_{pp}$ is as high as 30\% to 50\% at 2-3$~\rho_0 $. These parameters are particularly sensitive to the symmetry energy. The effect of the magnetic field on the nuclear properties is small for fields as high as 10$^{18}$G except for a small range of densities just above the crust-core transition. Tables with the EoS, and the parameters, are provided in the Supplementary Material section.

In this letter we discuss infinite statistics and motivate its role in quantum gravity. Then, we connect infinite statistics to a dynamical form of dark energy, and we obtain an expression for the evolution of the Hubble parameter that we compare to observation. The equation of state parameter $w_{eff} < -1$ in this framework.