Papers for Thursday, Sep 01 2022

We present photometric and spectroscopic data of SN 2018lab, a low luminosity type II-P supernova (LLSN) with a V-band peak luminosity of $-15.1\pm0.1$ mag. SN 2018lab was discovered by the Distance Less Than 40 Mpc (DLT40) SNe survey only 0.73 days post-explosion, as determined by observations from the Transiting Exoplanet Survey Satellite (TESS). TESS observations of SN 2018lab yield a densely sampled, fast-rising, early time light curve likely powered by circumstellar medium (CSM) interaction. The blue-shifted, broadened flash ionization lines in the earliest spectra ($<$2 days) of SN 2018lab provide further evidence for ejecta-CSM interaction. The early emission features in the spectra of SN 2018lab are well described by models of a red supergiant progenitor with an extended envelope and close-in CSM. As one of the few LLSNe with observed flash ionization features, SN 2018lab highlights the need for more early spectra to explain the diversity of flash feature morphology in type II SNe.

We use Gaia DR3 data to study the Collinder 132-Gulliver 21 region via the machine learning algorithm StarGO, and find eight subgroups of stars (ASCC 32, Collinder 132 gp 1--6, Gulliver 21) located in close proximity. Three co-moving populations were identified among these eight subgroups: (i) a coeval 25 Myr-old moving group (Collinder 132); (ii) an intermediate-age (50--100 Myr) group; and (iii) the 275 Myr-old dissolving cluster Gulliver 21. These three populations form parallel diagonal stripe-shape over-densities in the U--V distribution, which differ from open clusters and stellar groups in the solar neighborhood. We name this kinematic structure the Collinder 132-Gulliver 21 stream, as it extends over 270 pc in the 3D space. The oldest population Gulliver21 is spatially surrounded by the Collinder 132 moving group and the intermediate-age group. Stars in the Collinder 132-Gulliver 21 stream have an age difference up to 250 Myr. Metallicity information shows a variation of 0.3 dex between the youngest and oldest populations. The formation of the Collinder132-Gulliver 21 stream involves both star formation and dynamical heating. The youngest population (Collinder 132 moving group) with homogeneous metallicity is probably formed through filamentary star formation. The intermediate-age and the oldest population were then scatted by the Galactic bar or spiral structure resonance to intercept Collinder 132's orbit. Without mutual interaction between each population, the three populations are flying by each other currently and will become distinct three groups again in approximately ~50Myr.

Many stripped envelope supernovae (SNe) present a signature of high-velocity material responsible for broad absorption lines in the observed spectrum. These include SNe that are associated with long gamma-ray bursts (LGRBs) and low-luminosity GRBs (llGRBs), and SNe that are not associated with GRBs. Recently it was suggested that this high velocity material originates from a cocoon that is driven by a relativistic jet. In LGRBs this jet breaks out successfully from the stellar envelope, while in llGRBs and SNe that are not associated with GRBs the jet is choked. Here we use numerical simulations to explore the velocity distribution of an outflow that is driven by a choked jet and its dependence on the jet and progenitor properties. We find that in all cases where the jet is not choked too deep within the star, the outflow carries a roughly constant amount of energy per logarithmic scale of proper velocity over a wide range of velocities, which depends mostly on the cocoon volume at the time of its breakout. This is a universal property of jets driven outflows, which does not exist in outflows of spherically symmetric explosions or when the jets are choked very deep within the star. We therefore conclude that jets that are choked (not too deep) provide a natural explanation to the fast material seen in the early spectra of stripped envelope SNe that are not associated with LGRBs and that properties of this material could reveal information on the otherwise hidden jets.

Supermassive black holes have been known to disrupt passing stars producing outbursts called Tidal Disruption Events offering a unique view on the early stages of accretion disk and jet formation. The advent of large scale optical time-domain surveys has significantly increased the number of known events and challenged our understanding of their dynamics and emission processes. Here, we present the linear polarization curve of the most polarized tidal disruption without any indication of contribution from a jet to the emission. Our observations demonstrate that optical TDE emission can be powered by tidal stream shocks.

We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV 191b using backlighting by the superimposed elliptical system VV 191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using the James Webb Space Telescope (JWST) and the Hubble Space Telescope (HST) spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Distinct dust lanes are found to continue well beyond the bright spiral arms, and trace a complex web, not always following the spiral arm pitch, with a very sharp radial cutoff near 1.7 Petrosian radii. We present attenuation profiles and coverage statistics in each band at radii 14-21 kpc. We derive the attenuation law with wavelength; largely thanks to the leverage provided by the 1.5-micron measurements, the data both within and between the dust lanes clearly favor a steeper reddening behavior (R_V ~ 3.3 between 0.6 and 1.5 microns) than found for starbursts and star-forming regions of galaxies. The value of R appears to shift to larger values at increasing wavelengths (R~3.5 between 0.9 and 1.5 microns). For Milky Way-like grain populations, these low values suggests that most of the dust structure is spatially resolved at the 60-pc level. Mixing regions of different column density would flatten the wavelength behavior(larger R_V). The NIRCam images reveal a lens arc from a distant background galaxy at z~3.2$, spanning 90 degrees azimuthally at 2.8" from the foreground elliptical galaxy nucleus, plus a counterimage and additional weakly-lensed background galaxy. The lens model and imaging data give a mass/light ratio M/L_B=6.6 in solar units within the Einstein radius 1.8 kpc.

We present the Living Swift-XRT Point Source Catalogue (LSXPS) and real-time transient detector. Uniquely among X-ray catalogues, LSXPS is updated in near real-time, making this the first up-to-date record of the point sources detected by a sensitive X-ray telescope: the Swift-X-ray Telescope (XRT). The associated upper limit calculator likewise makes use of all available data allowing contemporary upper limits to be rapidly produced. The combination of these two products, along with the rapid timescale on which XRT data are available, enables us to carry out for the first time low-latency searches for new transient X-ray events fainter than those available to the current generation of wide-field imagers. In this paper we describe the system and the classification and dissemination of transient alerts it generates.

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state.

The origin of cosmic rays is a pivotal open issue of high-energy astrophysics. Supernova remnants are strong candidates to be the Galactic factory of cosmic rays, their blast waves being powerful particle accelerators. However, supernova remnants can power the observed flux of cosmic rays only if they transfer a significant fraction of their kinetic energy to the accelerated particles, but conclusive evidence for such efficient acceleration is still lacking. In this scenario, the shock energy channeled to cosmic rays should induce a higher post-shock density than that predicted by standard shock conditions. Here we show this effect, and probe its dependence on the orientation of the ambient magnetic field, by analyzing deep X-ray observations of the Galactic remnant of SN 1006. By comparing our results with state-of-the-art models, we conclude that SN 1006 is an efficient source of cosmic rays and obtain an observational support for the quasi-parallel acceleration mechanism.

One of the largest uncertainties in stellar structure and evolution theory is the transport of angular momentum in the stellar interiors. Asteroseismology offers a powerful tool for measuring the internal rotation frequencies of pulsating stars, but the number of such measurements has remained few for $\gtrsim 3\,{\rm M}_\odot$ main-sequence stars. In this work, we compile a list of 52 slowly pulsating B stars for which the interior rotation has been measured asteroseismically. The measurements of the spin parameters, which describe the relative importance of rotation, for the gravito-inertial mode oscillations show that for 40 of the stars the oscillations fall within the sub-inertial regime. We find that the core rotation frequencies of the stars decrease as a function of age, and show evidence of angular momentum transport occurring on the main-sequence. Finally, we derive the inclination angles of the stars, showing that they are generally consistent with the expectations from surface cancellation effects for the given oscillation modes.

James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) images of the luminous infrared (IR) galaxy VV 114 are presented. This redshift ~ 0.020 merger has a western component (VV 114W) rich in optical star clusters and an eastern component (VV 114E) hosting a luminous mid-IR nucleus hidden at UV and optical wavelengths by dust lanes. With MIRI, the VV 114E nucleus resolves primarily into bright NE and SW cores separated by 630 pc. This nucleus comprises 45% of the 15um light of VV 114, with the NE and SW cores having IR luminosities, L_ IR (8-1000um) ~ 8+/-0.8x10^10 L_sun and ~ 5+/-0.5x10^10 L_sun, respectively, and IR densities, Sigma_IR >~ 2+/-0.2x10^13 L_sun / kpc^2 and >~ 7+/-0.7x10^12 L_sun / kpc^2, respectively -- in the range of Sigma_IR for the Orion star-forming core and the nuclei of Arp 220. The NE core, previously speculated to have an Active Galactic Nucleus (AGN), has starburst-like mid-IR colors. In contrast, the VV 114E SW has AGN-like colors. Approximately 40 star-forming knots with L_IR ~ 0.02-5x10^10 L_sun are identified, 25% of which have no optical counterpart. Finally, diffuse emission accounts for 40-60% of the mid-IR emission. Mostly notably, filamentary Poly-cyclic Aromatic Hydrocarbon (PAH) emission stochastically excited by UV and optical photons accounts for half of the 7.7um light of VV 114. This study illustrates the ability of JWST to detect obscured compact activity and distributed PAH emission in the most extreme starburst galaxies in the local Universe.

The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument operating with an f/32 beam can be deployed using GNAO. The GNAO project includes the development of a new laser guide star facility which will consist of four side-launched laser beams supporting the two primary AO modes of GNAO: a wide-field mode providing an improved image quality over natural seeing for a 2-arcminute circular field-of-view and a narrow-field mode providing near diffraction-limited performance over a 20x20 arcsecond square field-of-view. The GNAO wide field mode will enable GIRMOS's multi-IFU configuration in which the science beam to each individual IFU will be additionally corrected using multi-object AO within GIRMOS. The GNAO narrow field mode will feed the GIRMOS tiled IFU configuration in which all IFUs are combined into a "super"-IFU in the center of the field. GNAO also includes the development of a new Real Time Controller, a new GNAO Facility System Controller and finally the development of a new AO Bench. We present in this paper an overview of the GNAO facility and provide a status update of each product.

The transiting planet HD80606b undergoes a 1000-fold increase in insolation during its 111-day orbit due to it being highly eccentric (e=0.93). The planet's effective temperature increases from 400K to over 1400K in a few hours as it makes a rapid passage to within 0.03AU of its host star during periapsis. Spectroscopic observations during the eclipse (which is conveniently oriented a few hours before periapsis) of HD80606b with the James Webb Space Telescope (JWST) are poised to exploit this highly variable environment to study a wide variety of atmospheric properties, including composition, chemical and dynamical timescales, and large scale atmospheric motions. Critical to planning and interpreting these observations is an accurate knowledge of the planet's orbit. We report on observations of two full-transit events: 7 February 2020 as observed by the TESS spacecraft and 7--8 December 2021 as observed with a worldwide network of small telescopes. We also report new radial velocity observations which when analyzed with a coupled model to the transits greatly improve the planet's orbital ephemeris. Our new orbit solution reduces the uncertainty in the transit and eclipse timing of the JWST era from tens of minutes to a few minutes. When combined with the planned JWST observations, this new precision may be adequate to look for non-Keplerian effects in the orbit of HD80606b.

We present the characterization of the low-gravity M6 dwarf 2MASS J0619-2903 previously identified as an unusual field object based on its strong IR excess and variable near-IR spectrum. Multiple epochs of low-resolution (R~150) near-IR spectra show large-amplitude (~0.1-0.5 mag) continuum variations on timescales of days to 12 years, unlike the small-amplitude variability typical for field ultracool dwarfs. The variations between epochs are well-modeled as changes in the relative extinction ($\Delta{A_V}\approx2$ mag). Likewise, Pan-STARRS optical photometry varies on timescales as long as 11 years (and possibly as short as an hour) and implies similar amplitude $A_V$ changes. NEOWISE mid-IR light curves also suggest changes on 6-month timescales, with amplitudes consistent with the optical/near-IR extinction variations. However, near-IR spectra, near-IR photometry, and optical photometry obtained in the past year indicate the source can also be stable on hourly and monthly timescales. From comparison to objects of similar spectral type, the total extinction of 2MASS J0619-2903 seems to be $A_V\approx4-6$ mag, with perhaps epochs of lower extinction. Gaia EDR3 finds that 2MASS J0619-2903 has a wide-separation (1.2' = 10450 AU) stellar companion, with an isochronal age of $31^{+22}_{-10}$ Myr and a mass of $0.30^{+0.04}_{-0.03}$ Msun. Adopting this companion's age and EDR3 distance (145.2$\pm$0.6 pc), we estimate a mass of 0.11-0.17 Msun for 2MASS J0619-2903. Altogether, 2MASS J0619-2903 appears to possess an unusually long-lived primordial circumstellar disk, perhaps making it a more obscured analog to the "Peter Pan" disks found around a few M dwarfs in nearby young moving groups.

Because of the diversity of stellar masses and orbital sizes of binary systems, and the complex interaction between star-star, star-planet and planet-planet, it has been difficult to fully characterize the planetary systems associated with binary systems. Here, we report high-precision astrometric observations of the low-mass binary system GJ 896AB, revealing the presence of a Jupiter-like planetary companion (GJ 896Ab). The planetary companion is associated to the main star GJ 896A, with an estimated mass of 2.3 Jupiter masses and an orbit period of 284.4 days. A simultaneous analysis of the relative astrometric data obtained in the optical and infrared with several telescopes, and the absolute astrometric data obtained at radio wavelengths with the Very Long Baseline Array (VLBA), reveals, for the first time, the fully characterized three-dimensional (3D) orbital plane orientation of the binary system and the planetary companion. The planetary and binary orbits are found to be in a retrograde configuration and with a large mutual inclination angle ($\Phi$ = 148 deg) between both orbital planes. Characterizing the 3D orbital architecture of binary systems with planets is important in the context of planet formation, as it could reveal whether the systems were formed by disk fragmentation or turbulence fragmentation, as well as the origin of spin-orbit misalignment. Furthermore, since most stars are in binary or multiple systems, our understanding of systems such as this one will help to further understand the phenomenon of planetary formation in general.

The Athena X-ray Integral Unit (X-IFU) is the high resolution X-ray spectrometer, studied since 2015 for flying in the mid-30s on the Athena space X-ray Observatory, a versatile observatory designed to address the Hot and Energetic Universe science theme, selected in November 2013 by the Survey Science Committee. Based on a large format array of Transition Edge Sensors (TES), it aims to provide spatially resolved X-ray spectroscopy, with a spectral resolution of 2.5 eV (up to 7 keV) over an hexagonal field of view of 5 arc minutes (equivalent diameter). The X-IFU entered its System Requirement Review (SRR) in June 2022, at about the same time when ESA called for an overall X-IFU redesign (including the X-IFU cryostat and the cooling chain), due to an unanticipated cost overrun of Athena. In this paper, after illustrating the breakthrough capabilities of the X-IFU, we describe the instrument as presented at its SRR, browsing through all the subsystems and associated requirements. We then show the instrument budgets, with a particular emphasis on the anticipated budgets of some of its key performance parameters. Finally we briefly discuss on the ongoing key technology demonstration activities, the calibration and the activities foreseen in the X-IFU Instrument Science Center, and touch on communication and outreach activities, the consortium organisation, and finally on the life cycle assessment of X-IFU aiming at minimising the environmental footprint, associated with the development of the instrument. It is expected that thanks to the studies conducted so far on X-IFU, along the design-to-cost exercise requested by ESA, the X-IFU will maintain flagship capabilities in spatially resolved high resolution X-ray spectroscopy, enabling most of the original X-IFU related scientific objectives of the Athena mission to be retained (abridged).

We obtained 838 Sloan r-band images (~28 hrs) of the GOODS-North field with the Large Binocular Camera (LBC) on the Large Binocular Telescope in order to study the presence of extended, low surface brightness features in galaxies and investigate the trade-off between image depth and resolution. The individual images were sorted by effective seeing, which allowed for optimal resolution and optimal depth mosaics to be created with all images with seeing FWHM < 0.9" and FWHM < 2.0", respectively. Examining bright galaxies and their substructure as well as accurately deblending overlapping objects requires the optimal resolution mosaic, while detecting the faintest objects possible (to a limiting magnitude of $m_{AB}$ ~ 29.2 mag) requires the optimal depth mosaic. The better surface brightness sensitivity resulting from the larger LBC pixels, compared to those of extant WFC3/UVIS and ACS/WFC cameras aboard the Hubble Space Telescope (HST) allows for unambiguous detection of both diffuse flux and very faint tidal tails. We created azimuthally-averaged radial surface brightness profiles for the 360 brightest galaxies in the mosaics. We find little difference in the majority of the light profiles from the optimal resolution and optimal depth mosaics. However, $\lesssim$ 15% of the profiles show excess flux in the galaxy outskirts down to surface brightness levels of $\mu^{AB}_{r} $ $\simeq$ 31 mag arcsec $^{-2}$. This is relevant to Extragalactic Background Light (EBL) studies as diffuse light in the outer regions of galaxies are thought to be a major contribution to the EBL. While some additional diffuse light exists in the optimal depth profiles compared to the shallower, optimal resolution profiles, we find that diffuse light in galaxy outskirts is a minor contribution to the EBL overall in the r-band.

Localized magnetic reconnection at the dayside magnetopause leads to the production of Flux Transfer Events (FTEs). The magnetic field within the FTEs exhibit complex helical flux-rope topologies. Leveraging the Adaptive Mesh Refinement (AMR) strategy, we perform a 3-dimensional magnetohydrodynamic simulation of the magnetosphere of an Earth-like planet and study the evolution of these FTEs. For the first time, we detect and track the FTE structures in 3D and present a complete volumetric picture of FTE evolution. The temporal evolution of thermodynamic quantities within the FTE volumes confirm that continuous reconnection is indeed the dominant cause of active FTE growth as indicated by the deviation of the P-V curves from an adiabatic profile. An investigation into the magnetic properties of the FTEs show a rapid decrease in the perpendicular currents within the FTE volume exhibiting the tendency of internal currents toward being field aligned. An assessment on the validity of the linear force-free flux rope model for such FTEs show that the structures drift towards a constant-$\alpha$ state but continuous reconnection inhibits the attainment of a purely linear force-free configuration. Additionally, the flux enclosed by the selected FTEs are computed to range between 0.3-1.5 MWb. The FTE with the highest flux content constitutes $\sim$ 1% of the net dayside open flux. These flux values are further compared against the estimates provided by the linear force-free flux-rope model. For the selected FTEs, the linear force-free model underestimated the flux content by up to 40% owing to the continuous reconnected flux injection.

Temporal variations of apparent magnitude, called light curves, are observational statistics of interest captured by telescopes over long periods of time. Light curves afford the exploration of Space Domain Awareness (SDA) objectives such as object identification or pose estimation as latent variable inference problems. Ground-based observations from commercial off the shelf (COTS) cameras remain inexpensive compared to higher precision instruments, however, limited sensor availability combined with noisier observations can produce gappy time-series data that can be difficult to model. These external factors confound the automated exploitation of light curves, which makes light curve prediction and extrapolation a crucial problem for applications. Traditionally, image or time-series completion problems have been approached with diffusion-based or exemplar-based methods. More recently, Deep Neural Networks (DNNs) have become the tool of choice due to their empirical success at learning complex nonlinear embeddings. However, DNNs often require large training data that are not necessarily available when looking at unique features of a light curve of a single satellite. In this paper, we present a novel approach to predicting missing and future data points of light curves using Gaussian Processes (GPs). GPs are non-linear probabilistic models that infer posterior distributions over functions and naturally quantify uncertainty. However, the cubic scaling of GP inference and training is a major barrier to their adoption in applications. In particular, a single light curve can feature hundreds of thousands of observations, which is well beyond the practical realization limits of a conventional GP on a single machine. Consequently, we employ MuyGPs, a scalable framework for hyperparameter estimation of GP models that uses nearest neighbors sparsification and local cross-validation. MuyGPs...

We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be suppressed by nucleon superfluidity. In the present study, focusing on the DU process and nucleon superfluidity, we investigate the effects of NS cooling on the X-ray bursts using a general-relativistic stellar-evolution code. We find that the DU process leads find the longer recurrence time and the higher peak luminosity, which could be obstructed by the neutrons superfluidity. We also apply our burst models to the comparison with {\it Clocked burster} GS 1826$-$24, and to the recurrence time of superburst triggered by carbon ignition. These effects are significant within a certain range of binary parameters and uncertainty of the NS equation of state.

The Primordial Disk of small icy planetesimals, once located at 15-30 AU from the Sun, was disrupted by giant planet migration in the early Solar System. The Primordial Disk thereby became the source region of objects in the current-day Kuiper Belt, Scattered Disk, and Oort Cloud. I present the thermophysics code "Numerical Icy Minor Body evolUtion Simulator", or NIMBUS, and use it to study the thermophysical evolution of planetesimals in the Primordial Disk prior to its disruption. Such modelling is mandatory in order to understand the behaviour of dynamically new comets from the Oort Cloud, as well as the activity of Centaurs and short-period comets from the Scattered Disk, that return pre-processed to the vicinity of the Sun. I find that bodies in the midst of the Primordial Disk with diameters ranging 4-200 km lost all their CO ice on time-scales of order 0.1-10 Myr depending on size, through a combination of protosolar and long-lived radionuclide heating. CO and other hypervolatiles therefore require a less volatile host for their storage. I consider two possible hosts: amorphous water ice and CO2 ice. Because of the high luminosity of the protosun, some Primordial Disk bodies may have sustained significant crystallisation, CO:CO2 segregation, and CO2 sublimation in the uppermost few tens of meters. I discuss how this may affect coma abundance ratios and distant activity in dynamically new comets.

We report on the galaxy density environment around a high-z radio galaxy (HzRG) at z=4.72, HSC J083913.17+011308.1 (HSC J0839+0113), probed using an r-dropout Lyman break galaxy (LBG) sample from the Hyper Suprime-Cam Subaru Strategic Program data. We find that HSC J0839+0113 resides in the outskirt of an overdense region identified by the r-dropout galaxies at a 4.7 sigma significance level. The projected distance between HSC J0839+0113 and the peak position of the overdense region is 0.4 physical Mpc which is shorter than the typical protocluster radius in this epoch. According to the extended Press Schechter and the light cone models, the HSC J0839+0113-hosted overdense region is expected to evolve into a halo > 10^14 Msun at z=0 with a high probability of >80 %. These findings suggest that HSC J0839+0113 is associated with a protocluster. The HSC J0839+0113 rich-system is the most overdense region of LBGs among the known protoclusters with LBGs in the same cosmic epoch.

How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on to the Milky Way and Andromeda galaxies, is related to the cosmic web. The analysis considers the three high-resolution simulations available in the HESTIA simulation suite, as well as the derived velocity shear and tidal tensors. We notice two eras in the Local Group accretion history, delimited by an epoch around $z \approx 0.7$. We also find that satellites can travel up to $\sim 4$ Mpc, relative to their parent halo before crossing its viral radius $R_{200}$. Finally, we observe a strong alignment of the infall direction with the axis of slowest collapse $\vec{e_3}$ of both tidal and shear tensors, implying satellites of the Local Group originated from one particular region of the cosmic web and were channeled towards us via the process of accretion.This alignment is dominated by the satellites that enter during the early infall era, i.e $z>0.7$.

The ALPACA experiment is a project aiming to observe sub-PeV gamma rays for the first time in the southern hemisphere. The main goal of ALPACA is to identify PeVatrons, the accelerators of Galactic PeV cosmic rays, by observing sub-PeV pion-decay gamma rays generated in interactions between PeV cosmic rays and the interstellar medium. This new air shower experiment is located at an altitude of 4,740 m above sea level in the middle of Mt. Chakartaya in Bolivia. The air shower array consists of 401 scintillation counters covering an 83,000 m$^2$ surface area. In addition, a water-Cherenkov-type muon detector array with an area of 3,700 m$^2$ is installed to discriminate gamma rays from background cosmic rays. The prototype array ALPAQUITA will start data taking in 2022 and will extend to ALPACA in 2024. We report on a general introduction to ALPACA, the progress of the project, and the sensitivity to sub-PeV gamma rays.

The LOPES experiment was a radio interferometer built at the existing air shower array KASCADE-Grande in Karlsruhe, Germany. The last configuration of LOPES was called LOPES 3D and consisted of ten tripole antennas. Each of these antennas consisted of three crossed dipoles east-west, north-south, and vertically aligned. With this, LOPES 3D had the unique possibility to study the benefits of measurements with vertically aligned antennas in the environment of the well understood and calibrated particle detector array KASCADE-Grande. The measurements with three spatially coincident antennas allows a redundant reconstruction of the electric field vector. Several methods to exploit the redundancy were developed and tested. Furthermore, for the first time in LOPES, the background noise could be studied polarization- and direction dependent. With LOPES 3D it could be demonstrated that radio detection reaches a higher efficiency for inclined showers when including measurements with vertically aligned antennas and that the vertical component gets more important for the measurement of inclined showers. In this contribution we discuss a weighting scheme for the best combination of three redundant reconstructed electric field vectors. Furthermore, we discuss the influence of these weighting schemes on the ability to reconstruct air showers using the radio method. We show an estimate of the radio efficiency for inclined showers with focus on the benefits of measurements with vertically aligned antennas and we present the direction dependent noise in the different polarizations.

The Principal Component Analysis (PCA) method and the Singular Value Decomposition (SVD) method are widely used for foreground subtraction in 21 cm intensity mapping experiments. We show their equivalence, and point out that the condition for completely clean separation of foregrounds and cosmic 21 cm signal using the PCA/SVD is unrealistic. We propose a PCA-based foreground subtraction method, dubbed "Singular Vector Projection (SVP)" method, which exploits a priori information of the left and/or right singular vectors of the foregrounds. We demonstrate with simulation tests that this new, semi-blind method can reduce the error of the recovered 21 cm signal by orders of magnitude, even if only the left and/or right singular vectors in the largest few modes are exploited. The SVP estimators provide a new, effective approach for 21 cm observations to remove foregrounds and uncover the physics in the cosmic 21 cm signal.

Our photometric and spectroscopic monitoring shows that starting with 2020 June 4, day +217 from optical maximum and well into its advanced nebular stage, Nova Sct 2019 began displaying a series of nine large amplitude flares (up to Delta(m)~1.7 mag), characterized by a rapid rise to peak (=<10 hours) and a fast exponential decline (e-folding time =<50 hours). The time interval Delta(t) between flares follows an ordered sequence, declining from 8.43 to 4.90 days, that safely allows to exclude that any other flare occured without being recorded by the observations. When the sequence of flares was over by 2020 July 28 (day +271), Nova Sct 2019 slowed its overall decline rate from Delta(m)=0.0067 mag/day to 0.0027 mag/day. The flares were caused by material expelled at high velocity (~1000 km/s) from the still burning WD. The cooler pseudo-photosphere forming at each flare in the expelled material, resulted in a recombination wave to spread through the original nova ejecta (at ~170 AU from the WD), quenching emission from [FeX] and [FeVII] and boosting that from lower ionization species. After each flare, once the small amount of expelled material had turned optically thin, the original nova ejecta resumed displaying [FeX] and [FeVII] emission lines, a fact that clearly proves the direct photo-ionization action exerted on the ejecta by the burning WD. While the other known flaring novae (V458 Vul, V4745 Sgr, and V5588 Sgr) presented the flares close to maximum brightness and with increasing Delta(t), Nova Sct 2019 is unique in having displayed them during the advanced nebular stage and with decreasing Delta(t).

Supra-arcade downflows (SADs) are dark structures descending towards post-reconnection flare loops observed in extreme ultraviolet or X-ray observations and are closely related to magnetic reconnection during solar flares. Due to the lack of statistical study on SADs in a single flare, evolutions of kinematic and thermal properties of SADs during the flare process still remain obscure. In this work, we identified 81 SADs in a flare that occurred on 2013 May 22 using observations of the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO). The kinematic properties of each SAD, including the appearance time, height, projective velocity, and acceleration were recorded. We found that the appearance heights of SADs become larger during the flare, which is likely due to the lift of the bottom of the plasma sheet. In the flare decay phase, the region where SADs mainly appear moves from the north part to the south side possibly related to a secondary eruption in the south side. The trajectories of most SADs can be fitted by one or two deceleration processes, while some special ones have positive accelerations during the descent. For the thermal properties, we selected 54 SADs, whose front and body could be clearly distinguished from the surrounding during the entire descent, to perform Differential Emission Measure analysis. It is revealed that the temperatures of the SAD front and body tend to increase during their downward courses, and the relationship between the density and temperature indicates that the heating is mainly caused by adiabatic compression.

We put constraints on the peak absolute magnitude of type Ia supernova using the Pantheon sample for type Ia supernova observations and the cosmic chronometers data for the Hubble parameter by a model independent and non-parametric approach. Our analysis is based on the Gaussian process regression. We find percent level bounds on the peak absolute magnitude. For completeness and to check the consistency of the results, we also include the Baryon acoustic oscillation data and the prior of the comoving sound horizon from Planck 2018 cosmic microwave background observations. The inclusion of these two data gives tighter constraints on it at the sub-percent level. The mean values of peak absolute magnitude from all these data are consistent with each other and the values are approximately equal to -19.4.

The low-luminosity Active Galactic Nuclei M87, archetype of Fanaroff-Riley I radio-galaxies, was observed in a historically quiet state in 2017. While one-zone leptonic jet models alone cannot explain the core radio-to-gamma-ray spectrum, we explore a hybrid jet-disc scenario. In this work, we model the overall spectral energy distribution of M87's core with a dominating one-zone lepto-hadronic jet component, coupled with the contribution from the accretion flow. We find close-to-equipartition parameter sets for which the jet component fits the radio-to-optical data as well as the gamma-ray band, while the accretion flow mainly contributes to the X-ray band. The effects of gamma-ray absorption by the Extragalactic Background Light during the propagation towards Earth are probed and are found to be negligible for this model. The neutrino flux produced by such scenarios is also calculated, but remains below the current instruments' sensitivity.

There are no strong constraints placed thus far on the amplitude of internal gravity waves (IGWs) that are stochastically excited in the radiative interiors of solar-type stars. Late F-type stars have relatively thin convective envelopes with fast convective flows and tend to be fast rotators compared to solar-type stars of later spectral types. These two elements are expected to directly impact the IGW excitation rates and properties. We want to estimate the amplitude of stochastically excited gravity modes (g-modes) in F-type stars for different rotational regimes. We used the ASH code to perform 3D simulations of deep-shell models of 1.3 $M_\odot$ F-type solar-type stars, including the radiative interior and the shallow convective envelope. The IGWs are excited by interface interactions between convective plumes and the top of the radiative interior. We were able to characterise the IGWs and g-mode properties in the radiative interior, and we compared these properties using the computation from the 1D oscillation code GYRE. The amplitude of low-frequency modes is significantly higher in fast-rotating models and the evolution of the period spacing of consecutive modes exhibits evidence of a behaviour that is modified by the influence of the Coriolis force. For our fastest rotating model, we were able to detect the intermediate degree g-mode signature near the top of the simulation domain. Nevertheless, the predicted luminosity perturbations from individual modes still remain at small amplitudes. We obtained mode amplitudes that are several orders of magnitude higher than those of prior 3D simulations of solar models. Our simulations suggest that g-mode signatures could be detectable in late F-type stars. [abridged]

Emission of dust up to a few micrometer in size by impacts of sand grains during saltation is thought to be one source of dust within the Martian atmosphere. To study this dust fraction, we carried out laboratory impact experiments. Small numbers of particles of about 200\textmu{}m in diameter impacted a simulated Martian soil (bimodal \textit{Mars Global Simulant}). Impacts occurred at angles of $\sim 18^\circ$ in vacuum with an impact speed of $\sim 1 \rm m/s$. Ejected dust was captured on adjacent microscope slides and the emitted particle size distribution (PSD) was found to be related to the soil PSD. We find that the ejection of clay sized dust gets increasingly harder the smaller these grains are. However, in spite of strong cohesive forces, individual impacts emit dust of 1\textmu{}m and less, i.e. dust in the size range that can be suspended in the Martian atmosphere. More generally, the probability of ejecting dust of a given size can be characterized by a power law in the size range between 0.5\textmu{}m and 5\textmu{}m (diameter).

We present a detailed study of the barium star at the heart of the planetary nebula Abell 70. Time-series photometry obtained over a period of more than ten years demonstrates that the barium-contaminated companion is a rapid rotator with temporal variability due to spots. The amplitude and phasing of the photometric variability changes abruptly, however there is no evidence for a change in the rotation period (P = 2.06~d) over the course of the observations. The co-addition of 17 high-resolution spectra obtained with VLT-UVES allow us to measure the physical and chemical properties of the companion, confirming it to be a chromospherically-active, late G-type sub-giant with more than +1~dex of barium enhancement. We find no evidence of radial velocity variability in the spectra, obtained over the course of approximately 130~d with a single additional point some 8 years later, with the radial velocities of all epochs approximately $-$10 \kms{} from the previously measured systemic velocity of the nebula. This is perhaps indicative that the binary has a relatively long period (P $\gtrsim$ 2~yr) and high eccentricity ($e\gtrsim$ 0.3), and that all the observations were taken around radial velocity minimum. However, unless the binary orbital plane is not aligned with the waist of the nebula or the systemic velocity of the binary is not equal to the literature value for the nebula, this would imply an unfeasibly large mass for the nebular progenitor.

We present a new method based on information theory to find the optimal number of bands required to measure the physical properties of galaxies with a desired accuracy. As a proof of concept, using the recently updated COSMOS catalog (COSMOS2020), we identify the most relevant wavebands for measuring the physical properties of galaxies in a Hawaii Two-0 (H20)- and UVISTA-like survey for a sample of $i<25$ AB mag galaxies. We find that with available $i$-band fluxes, $r$, $u$, IRAC/$ch2$ and $z$ bands provide most of the information regarding the redshift with importance decreasing from $r$-band to $z$-band. We also find that for the same sample, IRAC/$ch2$, $Y$, $r$ and $u$ bands are the most relevant bands in stellar mass measurements with decreasing order of importance. Investigating the inter-correlation between the bands, we train a model to predict UVISTA observations in near-IR from H20-like observations. We find that magnitudes in $YJH$ bands can be simulated/predicted with an accuracy of $1\sigma$ mag scatter $\lesssim 0.2$ for galaxies brighter than 24 AB mag in near-IR bands. One should note that these conclusions depend on the selection criteria of the sample. For any new sample of galaxies with a different selection, these results should be remeasured. Our results suggest that in the presence of a limited number of bands, a machine learning model trained over the population of observed galaxies with extensive spectral coverage outperforms template-fitting. Such a machine learning model maximally comprises the information acquired over available extensive surveys and breaks degeneracies in the parameter space of template-fitting inevitable in the presence of a few bands.

We present detailed timing and spectral analyses of the transient X-ray pulsar RX J0209.6$-$7427 in the Small Magellanic Cloud during its 2019 giant outburst. With a better known distance than most galactic X-ray pulsars, its peak luminosity is determined to be $(1.11\pm0.06)\times 10^{39}\, \rm erg\ s^{-1}$; it is thus a {\it bonda fide} pulsating ultraluminous X-ray source (PULX). Owing to the broad energy band of \textit{Insight}-HXMT, its pulsed X-ray emission was detected from 1 keV up to the 130$-$180 keV band, which is the highest energy emission detected from any PULXs outside the Milky Way. This allows us to conclude that its main pulsed X-ray emission is from the "fan beam" of the accretion column, and its luminosity is thus intrinsic. We also estimate its magnetic field of (4.8$-$8.6)$\times10^{12}$ G or (1.7$-$2.2)$\times10^{13}$ G, from its spin evolution or transition in the accretion column structure during the outburst; we suggest that the two values of the magnetic field strength correspond to the dipole and multipole magnetic fields of the neutron star, similar to the recent discovery in the Galactic PULX Swift J0243.6+6124. Therefore, the nature of the neutron star and its ULX emission can be understood within the current theoretical frame of accreting neutron stars. This may have implications for understanding the nature of those farther away extragalactic PULXs.

Context. Filamentary structures appear to be ubiquitous in the interstellar medium. Being able to detect and characterize them is the first step toward understanding their origin, their evolution, and their role in the Galactic cycle of matter. Aims. We present a new method, called FilDReaMS, to detect and analyze filaments in a given image. This method is meant to be fast, user-friendly, multi-scale, and suited for statistical studies. Methods. The input image is scanned with a rectangular model bar, which makes it possible to uncover structures that can be locally approximated by this bar and to derive their orientations. The bar width can be varied over a broad range of values to probe filaments of different widths. Results. We performed several series of tests to validate the method and to assess its sensitivity to the level of noise, the filament aspect ratios, and the dynamic range of filament intensities. We found that the method exhibits very good performance at recovering the orientation of the filamentary structures, with an accuracy of 0.5{\deg} in nominal conditions, up to 3{\deg} in the worst case scenario with high level of noise. The filaments width is recovered with uncertainties better than 0.5 px (pixels) in most of the cases, which could extend up to 3 px in case of low signal-to-noise ratios. Some attempt to build a correspondence between Plummer-type filament profiles and outcomes of the method is proposed, but remains sensitive to the local environment. Conclusions. Our method is found to be robust and adapted to the identification and the reconstruction of filamentary structures in various environments, from diffuse to dense medium. It allows us to explore the hierarchical scales of these filamentary structures with a high reliability, especially when dealing with their orientation.

RISTRETTO is a visible high-resolution spectrograph fed by an extreme adaptive optics (XAO) system, to be proposed as a visitor instrument on ESO VLT. The main science goal of RISTRETTO is the detection and atmospheric characterization of exoplanets in reflected light, in particular the temperate rocky planet Proxima b. RISTRETTO will be able to measure albedos and detect atmospheric features in a number of exoplanets orbiting nearby stars for the first time. It will do so by combining a high-contrast AO system working at the diffraction limit of the telescope to a high-resolution spectrograph, via a 7-spaxel integral-field unit (IFU) feeding single-mode fibers. Further science cases for RISTRETTO include the study of accreting protoplanets such as PDS 70 b & c through spectrally-resolved H-alpha emission; and spatially-resolved studies of Solar System objects such as icy moons and the ice giants Uranus and Neptune. The project is in an advanced design phase for the spectrograph and IFU/fiber-link sub-systems, and a preliminary design phase for the AO front-end. Construction of the spectrograph and IFU/fiber-link will start at the end of 2022. RISTRETTO is a pathfinder instrument in view of similar developments at ESO ELT, in particular the SCAO-IFU mode of ELT-ANDES and the future ELT-PCS instrument.

Both simulations and observations of the interstellar medium show that the study of the relative orientations between filamentary structures and the magnetic field can bring new insight into the role played by magnetic fields in the formation and evolution of filaments and in the process of star formation. We provide a first application of FilDReaMS, the new method to detect and analyze filaments. Our goal is to investigate the relative orientations between the detected filaments and the magnetic field. We apply FilDReaMS to a sample of four Herschel fields (G210, G300, G82, G202) characterized by different Galactic environments and evolutionary stages. First, we examine the networks formed by filaments of different bar widths as well as their hierarchical organization. Second, we compare the filament orientations to the magnetic field orientation inferred from Planck polarization data and, for the first time, we study the statistics of the relative orientation angle as functions of both spatial scale and ${\rm H_2}$ column density. We find preferential relative orientations in the four Herschel fields: small filaments with low column densities are slightly more parallel than perpendicular to the magnetic field, while large filaments (higher column densities) are oriented nearly perpendicular (or, in the case of G202, nearly parallel). In the two nearby fields (G210 and G300), we observe a transition from mostly parallel to mostly perpendicular at an $N_{\rm H_2} \simeq 1.1$ and $1.4\times10^{21}\,$cm$^{-2}$, respectively, consistent with the results of previous studies. Our results confirm the existence of a coupling between magnetic fields at cloud scales and filaments at smaller scale. They also illustrate the potential of combining Herschel and Planck observations, and they call for further statistical analyses with our dedicated method.

Observations of comet C/2016 R2 (PanSTARRS) have revealed exceptionally bright emission bands of N$_2^+$, the strongest ever observed in a comet spectrum. Alternatively, it appears to be poor in CN compared to other comets, and remarkably depleted in H$_2$O. Here we quantify the N$_2$ production rate from N$_2^+$ emission lines using the Haser model. We derived effective parent and daughter scalelengths for N2 producing N2+. This is the first direct measurement of such parameters. Using a revised fluorescence efficiency for N2+, the resulting production rate of molecular nitrogen is inferred to be Q(N$_2$) ~ 1 $\times 10^{28}$ molecules.s-1 on average for 11, 12, and 13 Feb. 2018, the highest for any known comet. Based on a CO production rate of Q(CO) ~ 1.1 $\times 10^{29}$ molecules.s-1, we find Q(N$-2$)/Q(CO)~0.09, which is consistent with the N$_2^+$/CO$^+$ ratio derived from the observed intensities of N$_2^+$ and CO$^+$ emission lines. We also measure significant variations in this production rate between our three observing nights, with Q(N$_2$) varying by plus or minus 20% according to the average value

Polars are highly magnetic cataclysmic variables which have been long observed to have both high and low brightness states. The duration of these states has been previously seen to vary from a number of days up to years. Despite this; these states and their physical origin has not been explained in a consistent manner. We present observations of the shortest duration states of a number of Polars observed by ZTF and TESS. This has allowed us to determine that short duration states are a relatively common feature across the population of Polars. Furthermore we have been able to generalise the model of star spot migration to explain both short lived high and low states in Polars by incorporating the interaction between the magnetic field of the white dwarf and that of the star spots.

The RISTRETTO project is aiming to build an instrument that will detect the reflected light from close-by exoplanet. It is a two stage instrument: An extreme AO system in the visible, followed by a seven spaxel single mode High resolution Spectrograph. In this paper we present the design of this spectrograph: a classical echelle spectrograph fed with single mode fibers. Standard single mode fibers have been chosen and are forming a long tilted slit in order to have the right order spacing on the detector. The instrument will be under vacuum and thermally controlled in order to make it stable.

In the framework of the GIARPS@TNG High-resolution Observations of T Tauri stars (GHOsT) project, we study the accretion properties of 37 Classical T Tauri Stars of the Taurus-Auriga star forming region (SFR) with the aim of characterizing their relation with the properties of the central star, of jets and disk winds, and of the global disk structure, in synergy with complementary ALMA millimiter observations. We derive stellar parameters, optical veiling, accretion luminosity ($\rm L_{acc}$) and mass accretion rate ($\rm \dot M_{acc}$) in a homogeneous and self-consistent way using high-resolution spectra acquired at the Telescopio Nazionale Galileo with the HARPS-N and GIANO spectrographs, and flux-calibrated based on contemporaneous low-resolution spectroscopic and photometric ancillary observations. The $\rm L_{acc}$-$\rm L_{\star}$, $\rm \dot{M}_{acc}$-$\rm M_{\star}$ and $\rm \dot{M}_{acc}$-$\rm M_{disk}$ relationships of the Taurus sample are provided and compared with those of the coeval SFRs of Lupus and Chamaeleon I. Our results demonstrate the potential of contemporaneous optical and near-infrared high-resolution spectroscopy to simultaneously provide precise measurements of stellar and accretion/wind properties of young stars.

The Gravitational-wave Optical Transient Observer (GOTO) is a wide-field telescope project focused on detecting optical counterparts to gravitational-wave sources. Each GOTO robotic mount holds eight 40 cm telescopes, giving an overall field of view of 40 square degrees. As of 2022 the first two GOTO mounts have been commissioned at the Roque de los Muchachos Observatory on La Palma, Canary Islands, and construction of the second node with two additional 8-telescope mounts has begin at Siding Spring Observatory in New South Wales, Australia. Once fully operational each GOTO mount will be networked to form a robotic, multi-site observatory, which will survey the entire visible sky every two nights and enable rapid follow-up detections of transient sources.

Strongly-magnetized ($B\geq10^{12}$ G) accreting neutron stars (NSs) are prime targets for studying the launching of jets by objects with a solid surface; while classical jet-launching models predict that such NSs cannot launch jets, recent observations and models argue otherwise. Transient Be/X-ray binaries (BeXRBs) are critical laboratories for probing this poorly-explored parameter space for jet formation. Here, we present the coordinated monitoring campaigns of three BeXRBs across four outbursts: giant outbursts of SAX 2103.5+4545, 1A 0535+262, and GRO J1008-57, as well as a Type-I outburst of the latter. We obtain radio detections of 1A 0535+262 during ten out of twenty observations, while the other targets remained undetected at typical limits of $20$-$50$ $\mu$Jy. The radio luminosity of 1A 0535+262 positively correlates with its evolving X-ray luminosity, and inhabits a region of the $L_X$-$L_R$ plane continuing the correlation observed previously for the BeXRB Swift J0243.6+6124. We measure a BeXRB $L_X$-$L_R$ coupling index of $\beta = 0.86 \pm 0.06$ ($L_R \propto L_X^\beta$), similar to the indices measured in NS and black hole low-mass X-ray binaries. Strikingly, the coupling's $L_R$ normalisation is $\sim 275$ and $\sim 6.2\times10^3$ times lower than in those two comparison samples, respectively. We conclude that jet emission likely dominates during the main peak of giant outbursts, but is only detectable for close-by or super-Eddington systems at current radio sensitivities. We discuss these results in the broader context of X-ray binary radio studies, concluding that our results suggest how supergiant X-ray binaries may host a currently unidentified additional radio emission mechanism.

The ESA/NASA Solar Orbiter space mission has been successfully launched in February 2020. Onboard is the Polarimetric and Helioseismic Imager (SO/PHI), which has two telescopes, a High Resolution Telescope (HRT) and the Full Disc Telescope (FDT). The instrument is designed to infer the photospheric magnetic field and line-of-sight velocity through differential imaging of the polarised light emitted by the Sun. It calculates the full Stokes vector at 6 wavelength positions at the Fe I 617.3 nm absorption line. Due to telemetry constraints, the instrument nominally processes these Stokes profiles onboard, however when telemetry is available, the raw images are downlinked and reduced on ground. Here the architecture of the on-ground pipeline for HRT is presented, which also offers additional corrections not currently available on board the instrument. The pipeline can reduce raw images to the full Stokes vector with a polarimetric sensitivity of $10^{-3}\cdot I_{c}$ or better.

Every moment, countless meteoroids enter our atmosphere unseen. The detection and measurement of meteors offer the unique opportunity to gain insights into the composition of our solar systems' celestial bodies. Researchers, therefore, carry out a wide-area-sky-monitoring to secure 360-degree video material, saving every single entry of a meteor. Existing machine intelligence cannot accurately recognize events of meteors intersecting the earth's atmosphere due to a lack of high-quality training data publicly available. This work presents four reusable open source solutions for researchers trained on data we collected due to the lack of available labeled high-quality training data. We refer to the proposed dataset as the NightSkyUCP dataset, consisting of a balanced set of 10,000 meteor- and 10,000 non-meteor-events. Our solutions apply various machine learning techniques, namely classification, feature learning, anomaly detection, and extrapolation. For the classification task, a mean accuracy of 99.1\% is achieved. The code and data are made public at figshare with DOI: 10.6084/m9.figshare.16451625

Observations of CH$^+$ are used to trace the physical properties of diffuse clouds, but this requires an accurate understanding of the underlying CH$^+$ chemistry. Until this work, the most uncertain reaction in that chemistry was dissociative recombination (DR) of CH$^+$. Using an electron-ion merged-beams experiment at the Cryogenic Storage Ring, we have determined the DR rate coefficient of the CH$^+$ electronic, vibrational, and rotational ground state applicable for different diffuse cloud conditions. Our results reduce the previously unrecognized order-of-magnitude uncertainty in the CH$^+$ DR rate coefficient to $\sim \pm 20\%$ and are applicable at all temperatures relevant to diffuse clouds, ranging from quiescent gas to gas locally heated by processes such as shocks and turbulence. Based on a simple chemical network, we find that DR can be an important destruction mechanism at temperatures relevant to quiescent gas. As the temperature increases locally, DR can continue to be important up to temperatures of $ \sim 600\,\mathrm{K} $ if there is also a corresponding increase in the electron fraction of the gas. Our new CH$^+$ DR rate coefficient data will increase the reliability of future studies of diffuse cloud physical properties via CH$^+$ abundance observations.

We present simulations of the formation and evolution of clusters in spiral arms. The simulations follow two different spiral arm regions, and the total gas mass is varied to produce a range of different mass clusters. We find that including photoionizing feedback produces the observed cluster mass radius relation, increasing the radii of clusters compared to without feedback. Supernovae have little impact on cluster properties. We find that in our high density, high gas mass simulations, star formation is less affected by feedback, as star formation occurs rapidly before feedback has much impact. In our lowest gas density simulation, the resulting clusters are completely different (e.g. the number of clusters and their masses) to the case with no feedback. The star formation rate is also significantly suppressed. The fraction of stars in clusters in this model decreases with time flattening at about 20\%. In our lowest gas simulation model, we see the formation of a star forming group with properties similar to an OB association, in particular similar to Orion Ia. We suggest that low densities, and stronger initial dynamics are conducive to forming associations rather than clusters. In all models cluster formation is complex with clusters merging and splitting. The most massive clusters which form have tended to undergo more mergers.

The theory of the formation of the first stars in the Universe, the so-called Population III (Pop III), has until now largely neglected the impact of magnetic fields. Complementing a series of recent studies of the magneto-hydrodynamic (MHD) aspects of Pop III star formation, we here carry out a suite of idealized numerical experiments where we ascertain how the fragmentation properties of primordial protostellar discs are modified if MHD effects are present. Specifically, starting from cosmological initial conditions, we focus on the central region in a select minihalo at redshift $z\sim$ 25, inserting a magnetic field at an intermediate evolutionary stage, normalized to a fraction of the equipartition value. To explore parameter space, we consider different field geometries, including uniform, radial, toroidal, and poloidal field configurations, with the toroidal configuration being the most realistic. The collapse of the gas is followed for $\sim$8 orders of magnitude in density after the field was inserted, until a maximum of $10^{15}{\rm \,cm}^{-3}$ is reached. We find that the magnetic field leads to a delay in the collapse of the gas. Moreover, the toroidal field has the strongest effect on the collapse as it inhibits the fragmentation of the emerging disc surrounding the central core and leads to the formation of a more massive core. The full understanding of the formation of Pop~III stars and their mass distribution thus needs to take into account the effect of magnetic fields. We further conclude that ideal MHD is only a first step in this endeavor, to be followed-up with a comprehensive treatment of dissipative effects, such as ambipolar diffusion and Ohmic dissipation.

Shocks and torques produced by non-axisymmetric structures such as spiral arms and bars may transport gas to galaxy central regions. We test this hypothesis by studying the dependence of concentration of CO luminosity ($C_{CO}$), molecular gas ($C_{mol}$), and star formation rate ($C_{SFR}$) in central $\sim$ 2 kpc on the $strength$ of non-axisymmetric disk structure using a sample of 57 disk galaxies selected from the EDGE-CALIFA survey. $C_{mol}$ is calculated using a CO-to-H$_2$ conversion factor that decreases with higher metallicity and higher stellar surface density. We find that $C_{mol}$ is systematically 0.22 dex lower than $C_{CO}$. We confirm that high $C_{mol}$ and strong non-axisymmetric disk structure are more common in barred galaxies than in unbarred galaxies. However, we find that spiral arms also increase $C_{mol}$. We show that there is a good correlation between $C_{mol}$ and the $strength$ of non-axisymmetric structure (which can be due to a bar, spiral arms, or both). This suggests that the stronger the bars and spirals, the more efficient the galaxy is at transporting cold gas to its center. Despite the small subsample size, $C_{mol}$ of the four Seyferts are not significantly reduced compared to inactive galaxies of similar disk structure, implying that the AGN feedback in Seyferts may not notably affect the molecular gas distribution in the central $\sim$2kpc. We find that $C_{SFR}$ tightly correlates with $C_{mol}$ in both unbarred and barred galaxies. Likewise, elevated $C_{SFR}$ is found in galaxies with strong disk structure. Our results suggest that the disk structure, either spirals or bars, can transport gas to the central regions, with higher inflow rates corresponding to stronger structure, and consequently boost central star formation. Both spirals and bars play, therefore, an essential role in the secular evolution of disk galaxies.

Theoretical works have looked into the various topologies and amplitudes, as well as the stability of the magnetic field that is expected to be present in the radiative interior of stars evolving after the main sequence. Such internal magnetic fields have never been observed in evolved stars. As a result, there is a major piece missing from our global picture of stars as dynamical bodies. Asteroseismology opened a window onto stellar internal dynamics, as oscillation frequencies, amplitudes, and lifetimes are affected by processes that are taking place inside the star. In this scope, magnetic signatures on mixed-mode frequencies have recently been characterized, but the task of detection remains challenging as the mixed-mode frequency pattern is highly complex and affected by rotational effects, while modes of different radial orders are often intertwined. In this work, we aim to build a bridge between theoretical prescriptions and complex asteroseismic data analysis to facilitate a future search and characterization of internal magnetism with asteroseismology. We investigate the effect of magnetic fields inside evolved stars with solar-like oscillations on the estimation of the period spacing of gravity-mode components of simulated mixed gravito-acoustic modes. We derived a new corrected stretching function of the power spectrum density to account for the presence of magnetic signatures on their frequencies. We demonstrate that the strong dependency of the amplitude of the magnetic signature with mixed-mode frequencies leads to biased estimates of period spacings towards lower values. We also show that a careful analysis of the oscillation frequency pattern through various period spacing estimates and across a broad frequency range might lead to the first detection of magnetic fields inside red giants and at the same time, we adjust the measured value of gravity-mode period spacing.

The migration history of Jupiter in the sun's natal disk remains poorly constrained. Here we consider how Jupiter's migration affects small-body reservoirs and how this constrains its original orbital distance from the Sun. We study the implications of large-scale and inward radial migration of Jupiter for the inner solar system while considering the effects of collisional evolution of planetesimals. We use analytical prescriptions to simulate the growth and migration of Jupiter in the gas disk. We assume the existence of a planetesimal disk inside Jupiter's initial orbit. This planetesimal disk received an initial total mass and size-frequency distribution (SFD). Planetesimals feel the effects of aerodynamic gas drag and collide with one another, mostly while shepherded by the migrating Jupiter. Our main goal is to measure the amount of mass in planetesimals implanted into the main asteroid belt (MAB) and the SFD of the implanted population. We also monitor the amount of dust produced during planetesimal collisions. We find that the SFD of the planetesimal population implanted into the MAB tends to resemble that of the original planetesimal population interior to Jupiter. We also find that unless very little or no mass existed between 5 au and Jupiter's original orbit, it would be difficult to reconcile the current low mass of the MAB with the possibility that Jupiter migrated from distances beyond 15 au. This is because the fraction of the original disk mass that gets implanted into the MAB is very large. Finally, we discuss the implications of our results in terms of dust production to the so-called NC-CC isotopic dichotomy.

I start by describing some updates to the pPXF method, which has been used to measure stellar and gas kinematics as well as the formation history (SFH) and chemical composition of galaxies. I outline the novel linearly-constrained least-squares optimization algorithm used by pPXF and I illustrate the changes I made to be able to include photometric measurements together with full-spectrum fitting in pPXF. Then I present an application of the revised pPXF method to the study of the non-parametric SFH and metallicity $[M/H]$ of a sample of 3200 galaxies at redshift $0.6<z<1$ (median $z=0.76$), complete above a stellar mass $M_\ast>3\times10^{10}$ M$_\odot$, with spectroscopy from the LEGA-C survey and 28-bands photometry from two alternative catalogues. I extract and compare the stellar population using three independent stellar population synthesis (SPS) methods and both photometric catalogues. I find robust trends in the global light-weighted ages and $[M/H]$ consistent and of similar quality as those from nearby galaxy surveys, with the well-known main dependence on the galaxies' stellar velocity dispersion $\sigma_\ast$ (or alternative measures of central density). The recovered SFH indicate a sharp and strikingly clear boundary from star formation to quenching at $\lg(\sigma_\ast/km s^{-1})\approx2.3$, similar to what is invoked by some models. Equally clear quenching boundaries are seen at $[M/H]\approx-0.1$ and for a Sersic exponent $n_{\rm Ser}\approx0.5$. Results are consistent with two SPS methods and both photometric catalogues, but the third SPS method displays significant differences, highlighting the importance of comparing model assumptions. The pPXF software is available from https://pypi.org/project/ppxf/.

We study the generation of a localized peak in the primordial spectrum of curvature perturbations from a transient dissipative phase during inflation, leading to a large population of primordial black holes. The enhancement of the power spectrum occurs due to stochastic thermal noise sourcing curvature fluctuations. We solve the stochastic system of Einstein equations for many realizations of the noise and obtain the distribution for the curvature power spectrum. We then propose a method to find its expectation value using a deterministic system of differential equations. In addition, we find a single stochastic equation whose analytic solution helps to understand the main features of the spectrum. Finally, we derive a complete expression and a numerical estimate for the energy density of the stochastic background of gravitational waves induced at second order in perturbation theory. This includes the gravitational waves induced during inflation, during the subsequent radiation epoch and their mixing. Our scenario provides a novel way of generating primordial black hole dark matter with a peaked mass distribution and a detectable stochastic background of gravitational waves from inflation.

Reconstruction of gravitational lensing effects in the CMB from current and upcoming surveys is still dominated by temperature anisotropies. Extragalactic foregrounds in temperature maps can induce significant biases in the lensing power spectrum obtained with the standard quadratic estimators. Techniques such as masking cannot remove these foregrounds fully, and the residuals can still lead to large biases if unaccounted for. In this paper, we study the "shear-only" estimator, an example of a class of geometric methods that suppress extragalactic foreground contamination while making only minimal assumptions about foreground properties. The shear-only estimator has only been formulated in the flat-sky limit and so is not easily applied to wide surveys. Here, we derive the full-sky version of the shear-only estimator and its generalisation to an $m=2$ multipole estimator that has improved performance for lensing reconstruction on smaller scales. The multipole estimator is generally not separable, and so is expensive to compute. We explore separable approximations based on a singular-value decomposition, which allow efficient evaluation of the estimator with real-space methods. Finally, we apply these estimators to simulations that include extragalactic foregrounds and verify their efficacy in suppressing foreground biases.

Spitzer/IRAC imaging has revealed that the brightest $z\sim7-8$ galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency among early galaxies. However, IRAC's limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here we use CEERS JWST/NIRCam data to characterize 118 UV-faint (median M$_{UV}=-19.6$) $z\sim6.5-8$ galaxies. We find that the SEDs are typically dominated by young ($\sim$10-50 Myr), low-mass ($M_\ast\sim10^8\ M_\odot$) stellar populations, with no need for extreme masses ($\sim10^{11} M_\odot$) among our sample in contrast to recent findings in CEERS. Considering previous studies of UV-bright (M$_{UV}\sim-22$) $z\sim7-8$ galaxies, we find evidence for a strong (5-10$\times$) increase in specific star formation rate toward lower luminosities (median sSFR=82 Gyr$^{-1}$ in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median $\xi_{ion}=10^{25.7}$ erg$^{-1}$ Hz). In spite of their much larger sSFRs, we find no significant increase in [OIII]$+$H$\beta$ EWs towards fainter M$_{UV}$ (median $\approx$780 $\mathring{A}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities ($\lesssim$3% $Z_\odot$) where [OIII] emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our CEERS galaxies. While the majority of our objects are very blue (median $\beta=-2.0$), we identify a significant tail of very dusty galaxies ($\beta\sim-1$) at $\approx$0.5$L_{UV}^\ast$ which may contribute significantly to the $z\sim7-8$ star formation rate density.

We investigate shock structures driven by merger events in high-resolution simulations that result in a galaxy with a virial mass M ~ 1e12 Msol. We find that the sizes and morphologies of the internal shocks resemble remarkably well those of the newly-detected class of odd radio circles (ORCs). This would highlight a so-far overlooked mechanism to form radio rings, shells and even more complex structures around elliptical galaxies. Mach numbers of M = 2-3 for such internal shocks are in agreement with the spectral indices of the observed ORCs. We estimate that ~5 percent of galaxies could undergo merger events which occasionally lead to such prominent structures within the galactic halo during their lifetime, explaining the low number of observed ORCs. At the time when the shock structures are matching the physical sizes of the observed ORCs, the central galaxies are typically classified as early-type galaxies, with no ongoing star formation, in agreement with observational findings. Although the energy released by such mergers could potentially power the observed radio luminosity already in Milky-Way-like halos, our predicted luminosity from a simple, direct shock acceleration model is much smaller than the observed one. Considering the estimated number of candidates from our cosmological simulations and the higher observed energies, we suggest that the proposed scenario is more likely for halo masses around 1e13 Msol in agreement with the observed stellar masses of the galaxies at the center of ORCs.

It has recently been demonstrated that Reissner-Nordstr\"om black holes in composed Einstein-Maxwell-scalar field theories can support static scalar field configurations with a non-minimal negative coupling to the Maxwell electromagnetic invariant of the charged spacetime. We here reveal the physically interesting fact that scalar field configurations with a non-minimal {\it positive} coupling to the spatially-dependent Maxwell electromagnetic invariant ${\cal F}\equiv F_{\mu\nu}F^{\mu\nu}$ can also be supported in black-hole spacetimes. Intriguingly, it is explicitly proved that the positive-coupling black-hole spontaneous scalarization phenomenon is induced by a non-zero combination $a\cdot Q\neq0$ of {\it both} the spin $a\equiv J/M$ and the electric charge $Q$ of the central supporting black hole. Using analytical techniques we prove that the regime of existence of the positive-coupling spontaneous scalarization phenomenon of Kerr-Newman black holes with horizon radius $r_+(M,a,Q)$ and a non-zero electric charge $Q$ (which, in principle, may be arbitrarily small) is determined by the {\it critical onset line} $(a/r_+)_{\text{critical}}=\sqrt{2}-1$. In particular, spinning and charged Kerr-Newman black holes in the composed Einstein-Maxwell-scalar field theory are spontaneously scalarized by the positively coupled fields in the dimensionless charge regime $0<{{Q}\over{M}}\leq\sqrt{2\sqrt{2}-2}$ if their dimensionless spin parameters lie above the critical onset line ${{a(Q)}\over{M}}\geq \big[{{a(Q)}\over{M}}\big]_{\text{critical}}={{1+\sqrt{1-2(2-\sqrt{2}){(Q/M)}^2}}\over{2\sqrt{2}}}$.

In this work, we investigate the existence of string theory solutions with a $d$-dimensional (quasi-) de Sitter spacetime, for $3 \leq d \leq 10$. Considering classical compactifications, we derive no-go theorems valid for general $d$. We use them to exclude (quasi-) de Sitter solutions for $d \geq 7$. In addition, such solutions are found unlikely to exist in $d=6,5$. For each no-go theorem, we further compute the $d$-dependent parameter $c$ of the swampland de Sitter conjecture, $M_p \frac{|\nabla V|}{V} \geq c$. Remarkably, the TCC bound $c \geq \frac{2}{\sqrt{(d-1)(d-2)}}$ is then perfectly satisfied for $d \geq 4$, with several saturation cases. However, we observe a violation of this bound in $d=3$. We finally comment on related proposals in the literature, on the swampland distance conjecture and its decay rate, and on the so-called accelerated expansion bound.

The centers of the core-collapse supernovae are one of the densest environments in the universe. Under such conditions, it is conceivable that a first-order phase transition from ordinary nuclear matter to the quark-gluon plasma occurs. This transition releases a large amount of latent heat that can drive a supernova explosion and may imprint a sharp signature in the neutrino signal. We show how this snap feature, if observed at large-scale neutrino detectors, can set competitive limits on the neutrino masses and assist the localization of the supernova via triangulation. The 95%C.L. limit on the neutrino mass can reach 0.16 eV in Ice-Cube, 0.22 eV in Hyper-Kamiokande, and 0.58 eV in DUNE, for a supernova at a distance of 10 kpc. For the same distance and in the most optimistic neutrino conversion case, the triangulation method can constrain the $1\sigma$ angular uncertainty of the supernova localization within $\sim 0.3-9.0$ deg in the considered pairs of the detectors, leading to an improvement up to an order of magnitude with respect to the often considered in the literature rise time of the neutronization burst.

We explore the possibility of using superfluid helium for direct detection of sub-GeV dark matter (DM). We discuss the relevant phenomenology resulting from the scattering of an incident dark matter particle on a Helium nucleus. Rather than directly exciting quasi-particles, DM in this mass range will interact with a single He atom, triggering an atomic cascade which eventually also includes emission and thermalization of quasi-particles. We present in detail the analytical framework needed for modeling these processes and determining the resulting flux of quasi-particles. We propose a novel method for detecting this flux with modern force-sensitive devices, such as nanoelectro-mechanical system (NEMS) oscillators, and derive the sensitivity projections for a generic sub-GeV DM detection experiment using such sensors.

We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit.

We propose a scenario where superheavy dark matter (DM) can be produced via symmetry restoration first-order phase transition during inflation triggered by the evolution of the inflaton field. The phase transition happens in a spectator sector coupled to the inflaton field. During the phase transition, the spectator field tunnels from a symmetry-broken vacuum to a symmetry-restored vacuum. The massive particles produced after bubble collisions are protected against decaying by the restored symmetry and may serve as a DM candidate in the later evolution of the Universe. We show that the latent heat released during the phase transition can be sufficient to produce the DM relic abundance observed today. In addition, accompanied with the super heavy DM, this first-order phase transition also produces gravitational waves detectable via future gravitational wave detectors.

Tilt-to-length (TTL) noise from angular jitter in LISA is projected to be the dominant noise source in the milli-Hertz band unless corrected in post-processing. The correction is only possible after removing the overwhelming laser phase noise using time-delay interferometry (TDI). We present here a frequency domain model that describes the effect of angular motion of all three spacecraft on the interferometric signals after propagating through TDI. We then apply a Fisher information matrix analysis to this model to calculate the minimum uncertainty with which TTL coupling coefficients may be estimated. Furthermore, we show the impact of these uncertainties on the residual TTL noise in the gravitational wave readout channel, and compare it to the impact of the angular witness sensors' readout noise. We show that the residual TTL noise post-subtraction in the TDI variables for a case using the LISA angular jitter requirement and integration time of one day is limited to the 8\,pm/$\sqrt{\rm Hz}$ level by angular sensing noise. However, using a more realistic model for the angular jitter we find that the TTL coupling uncertainties are 70 times larger, and the noise subtraction is limited by these uncertainties to the 14\,pm/$\sqrt{\rm Hz}$ level.

We present a comprehensive dynamical systems analysis of homogeneous and isotropic Friedmann-La\^{i}matre-Robertson-Walker cosmologies in the Hu-Sawicki $f(R)$ dark energy model for the parameter choice $\{n,C_1\}=\{1,1\}$. For a generic $f(R)$ theory, we outline the procedures of compactification of the phase space, which in general is 4-dimensional. We also outline how, given an $f(R)$ model, one can determine the coordinate of the phase space point that corresponds to the present day universe and the equation of a surface in the phase space that represents the $\Lambda$CDM evolution history. Next, we apply these procedures to the Hu-Sawicki model under consideration. We identify some novel features of the phase space of the model such as the existence of invariant submanifolds and 2-dimensional sheets of fixed points. We determine the physically viable region of the phase space, the fixed point corresponding to possible matter dominated epochs and discuss the possibility of a non-singular bounce, re-collapse and cyclic evolution. We also provide a numerical analysis comparing the $\Lambda$CDM evolution and the Hu-Sawicki evolution.