Papers for Thursday, Nov 10 2022

Uncertainty in the initial-final mass relation (IFMR) has long been a problem in understanding the final stages of massive star evolution. One of the major challenges of constraining the IFMR is the difficulty of measuring the mass of non-luminous remnant objects (i.e. neutron stars and black holes). Gravitational wave detectors have opened the possibility of finding large numbers of compact objects in other galaxies, but all in merging binary systems. Gravitational lensing experiments using astrometry and photometry are capable of finding compact objects, both isolated and in binaries, in the Milky Way. In this work we improve the PopSyCLE microlensing simulation code in order to explore the possibility of constraining the IFMR using the Milky Way microlensing population. We predict that the Roman Space Telescope's microlensing survey will likely be able to distinguish different IFMRs based on the differences at the long end of the Einstein crossing time distribution and the small end of the microlensing parallax distribution, assuming the small ($\pi_E \lesssim 0.02$) microlensing parallaxes characteristic of black hole lenses are able to be measured accurately. We emphasize that future microlensing surveys need to be capable of characterizing events with small microlensing parallaxes in order to place the most meaningful constraints on the IFMR.

The sudden release of magnetic energy on the Sun drives powerful solar flares and coronal mass ejections. The key issue is the difficulty in predicting the occurrence time and location of strong solar eruptions, i.e., those leading to the high impact space weather disturbances at the near-Earth environment. Solar radio imaging helps identify the magnetic field characteristics of active regions susceptible to intense flares and energetic coronal mass ejections. Mapping of the Sun at X-band (8.1 -- 9.3 GHz) with the 12-m radio telescope at the Arecibo Observatory allows monitoring of the evolution of the brightness temperature of active regions in association with the development of magnetic complexity, which can lead to strong eruptions. For a better forecasting strategy in the future, such ground-based radio observations of high-spatial and temporal resolution, along with a full polarization capability, would have tremendous potential not only to understand the magnetic activity of solar eruptions, but also for revealing the particle acceleration mechanism and additional exciting science.

This chapter of the Planetary Exploration Horizon 2061 Report reviews the way the six key questions about planetary systems, from their origins to the way they work and their habitability, identified in chapter 1, can be addressed by means of solar system exploration, and how one can find partial answers to these six questions by flying to the different provinces to the solar system: terrestrial planets, giant planets, small bodies, and up to its interface with the local interstellar medium. It derives from this analysis a synthetic description of the most important space observations to be performed at the different solar system objects by future planetary exploration missions. These observation requirements illustrate the diversity of measurement techniques to be used as well as the diversity of destinations where these observations must be made. They constitute the base for the identification of the future planetary missions we need to fly by 2061, which are described in chapter 4. Q1- How well do we understand the diversity of planetary systems objects? Q2- How well do we understand the diversity of planetary system architectures? Q3- What are the origins and formation scenarios for planetary systems? Q4- How do planetary systems work? Q5- Do planetary systems host potential habitats? Q6- Where and how to search for life?

Changing-look (CL) AGN are unique probes of accretion onto supermassive black holes (SMBHs), especially when simultaneous observations in complementary wavebands allow investigations into the properties of their accretion flows. We present the results of a search for CL behaviour in 412 Swift-BAT detected AGN with multiple epochs of optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS). 125 of these AGN also have 14-195 keV ultra-hard X-ray light-curves from Swift-BAT which are contemporaneous with the epochs of optical spectroscopy. Eight CL events are presented for the first time, where the appearance or disappearance of broad Balmer line emission leads to a change in the observed Seyfert type classification. Combining with known events from the literature, 21 AGN from BASS are now known to display CL behaviour. Nine CL events have 14-195 keV data available, and five of these CL events can be associated with significant changes in their 14-195 keV flux from BAT. The ultra-hard X-ray flux is less affected by obscuration and so these changes in the 14-195 keV band suggest that the majority of our CL events are not due to changes in line-of-sight obscuration. We derive a CL rate of 0.7-6.2 per cent on 10-25 year time-scales, and show that many transitions happen within at most a few years. Our results motivate further multi-wavelength observations with higher cadence to better understand the variability physics of accretion onto SMBHs.

The Vera Rubin Observatory, set to begin observations in mid-2024, will increase our discovery rate of supernovae to well over one million annually. There has been a significant push to develop new methodologies to identify, classify and ultimately understand the millions of supernovae discovered with the Rubin Observatory. Here, we present the first simulation-based inference method using normalizing flows, trained to rapidly infer the parameters of toy supernovae model in multivariate, Rubin-like datastreams. We find that our method is well-calibrated compared to traditional inference methodologies (specifically MCMC), requiring only one-ten-thousandth of the CPU hours during test time.

We present a multifrequency analysis of the radio galaxy 3CR 196.1 ($z = 0.198$), associated with the brightest galaxy of the cool core cluster CIZAJ0815.4-0303. This nearby radio galaxy shows a hybrid radio morphology and an X-ray cavity, all signatures of a turbulent past activity, potentially due to merger events and AGN outbursts. We present results of the comparison between $Chandra$ and VLT/MUSE data for the inner region of the galaxy cluster, on a scale of tens of kpc. We discovered H$\alpha$ + [N II]$\lambda6584$ emission spatially associated with the X-ray cavity (at $\sim$10 kpc from the galaxy nucleus) instead of with its rim. This result differs from previous discoveries of ionized gas surrounding X-ray cavities in other radio galaxies harbored in galaxy clusters and could represent the first reported case of ionized gas filling an X-ray cavity, either due to different AGN outbursts or to the cooling of warm ($10^4<T\leq10^7$ K) AGN outflows. We also found that the H$\alpha$, [N II]$\lambda\lambda6548,6584$ and [S II]$\lambda\lambda6718,6733$ emission lines show an additional redward component, at $\sim$1000 km$\,$s$^{-1}$ from rest frame, with no detection in H$\beta$ or [O III]$\lambda\lambda4960,5008$. We believe the most likely explanation for this redward component is the presence of a background gas cloud since there appears to be a discrete difference in velocities between this component and the rest frame.

With the advent of high cadence, all-sky automated surveys, supernovae (SNe) are now discovered closer than ever to their dates of explosion. However, young pre-maximum light follow-up spectra of Type Ic supernovae (SNe Ic), probably arising from the most stripped massive stars, remain rare despite their importance. In this paper we present a set of 49 optical spectra observed with the Las Cumbres Observatory through the Global Supernova Project for 6 SNe Ic, including a total of 17 pre-maximum spectra, of which 8 are observed more than a week before V-band maximum light. This dataset increases the total number of publicly available pre-maximum light SN Ic spectra by 25% and we provide publicly available SNID templates that will significantly aid in the fast identification of young SNe Ic in the future. We present detailed analysis of these spectra, including Fe II 5169 velocity measurements, O I 7774 line strengths, and continuum shapes. We compare our results to published samples of stripped supernovae in the literature and find one SN in our sample that stands out. SN 2019ewu has a unique combination of features for a SN Ic: an extremely blue continuum, high absorption velocities, a P-cygni shaped feature almost 2 weeks before maximum light that TARDIS radiative transfer modeling attributes to C II rather than H$\alpha$, and weak or non-existent O I 7774 absorption feature until maximum light.

We use high resolution numerical simulations in order to analyze the stellar bar evolution in spinning dark matter (DM) halos. Previous works have shown that the halo spin has a substantial effect on the bar evolution and can lead to bar dissolution following the vertical buckling instability. Here, we invoke the DM spin sequence, $\lambda=0-0.09$, and study the effect of DM density along this $\lambda$-sequence by varying the compactness of DM halo. We find that (1) varying the DM density has a profound effect on the stellar bar evolution along the $\lambda$-sequence, namely, on its amplitude, pattern speed, buckling time, etc.; (2) For $\lambda\gtrsim 0.04$, the buckling instability has been delayed progressively, and does not occur when the bar has reached its maximal strength; (3) Instead, stellar bars remain near maximal strength, and their amplitude plateau stage extends over $\sim 1-7$ Gyr, terminating with the buckling instability; (4) Although stellar bars remain strong during the plateau, their pattern speed stays nearly constant. The reason for this unusual behavior of stellar bars follows from the highly reduced gravitational torques which they experience due to the DM bar being aligned with the stellar bar. The performed orbital analysis shows that the delayed buckling results from a slow evolution of stellar oscillations along the bar major and vertical axes -- thus postponing the action of the vertical 2:1 resonance which pumps the rotational energy into vertical motions; (5) Peanut/boxy shaped bulges form at the beginning of the plateau and grow with time; (6) Strong stellar bars in spinning halos can avoid fast braking, resolving the long standing discrepancy between observations and $N$-body simulations. This behavior of stellar bars along the $\lambda$- and DM density-sequences, reveals a wealth of stellar bar properties which require additional study.

In this work, we combine the semi-analytic model of galaxy formation and evolution SAG with the $102$ relaxed simulated galaxy clusters from The Three Hundred project, and we study the link between the quenching of star formation (SF) and the physical processes that galaxies experience through their dynamical history in and around clusters. We classify galaxies in four populations based on their orbital history: recent and ancient infallers, and backsplash and neighbouring galaxies. We find that $\sim 85$ per cent of the current population of quenched galaxies located inside the clusters are ancient infallers with low or null content of hot and cold gas. The fraction of quenched ancient infallers increases strongly between the first and second pericentric passage, due to the removal of hot gas by the action of ram-pressure stripping (RPS). The majority of them quenches after the first pericentric passage, but a non-negligible fraction needs a second passage, specially galaxies with $M_\star \leq 10^{10.5} \, {\rm M_\odot}$. Recent infallers represent $\sim 15$ per cent of the quenched galaxies located inside the cluster and, on average, they contain a high proportion of hot and cold gas; moreover, pre-processing effects are the responsible for quenching the recent infallers prior to infall onto the main cluster progenitor. The $\sim 65$ per cent of quenched galaxies located around clusters are backsplash galaxies, for which the combination of RPS acting during a pre-processing stage and inside the cluster is necessary for the suppression of SF in this population.

The dynamical excitation of asteroids due to mean motion resonant interactions with planets is enhanced when their parent star leaves the main sequence. However, numerical investigation of resonant outcomes within post-main-sequence simulations is computationally expensive, limiting the extent to which detailed resonant analyses have been performed. Here, we combine the use of a high-performance computer cluster and the general semianalytical libration width formulation of Gallardo et al. (2021) in order to quantify resonant stability, strength and variation instigated by stellar evolution for a single-planet system containing asteroids on both crossing and non-crossing orbits. We find that resonant instability can be accurately bound with only main-sequence values by computing a maximum libration width as a function of asteroid longitude of pericentre. We also quantify the relative efficiency of mean motion resonances of different orders to stabilize versus destabilize asteroid orbits during both the giant branch and white dwarf phases. The 4:1, 3:1 and 2:1 resonances represent efficient polluters of white dwarfs, and even when in the orbit-crossing regime, both the 4:3 and 3:2 resonances can retain small reservoirs of asteroids in stable orbits throughout giant branch and white dwarf evolution. This investigation represents a preliminary step in characterising how simplified extrasolar Kirkwood gap structures evolve beyond the main-sequence.

We use two zoom-in $\Lambda$CDM hydrodynamical simulations of massive disk galaxies to study the possible existence of fixed satellite groups showing a kinematically-coherent behaviour across evolution (angular momentum conservation and clustering). We identify three such groups in the two simulations, defining kinematically-coherent, time-persistent planes (KPPs) that last at least from virialization to $z=0$ (more than 7 Gyrs). This proves that orbital pole clustering is not necessarily set in at low redshift, representing a long-lived property of galaxy systems. KPPs are thin and oblate, represent $\sim25-40\%$ of the total number of satellites in the system, and are roughly perpendicular to their corresponding central disk galaxies during certain periods, consistently with Milky Way $z=0$ data. KPP satellite members are statistically distinguishable from satellites outside KPPs: they show higher specific orbital angular momenta, orbit more perpendicularly to the central disk galaxy, and have larger pericentric distances, than the latter. We numerically prove, for the first time, that KPPs and the best-quality positional planes share the same space configuration across time, such that KPPs act as `skeletons' preventing the latter of being washed out in short timescales. In one of the satellite-host systems, we witness the late capture of a massive dwarf galaxy endowed with its own satellite system, also organized into a KPP configuration prior to its capture. We briefly explore the consequences this event has on the host's KPP, and on the possible enhancement of the asymmetry in the number of satellites rotating in one sense or the opposite within the KPP.

Over the past decade, the disparity between the value of the cosmic expansion rate directly determined from measurements of distance and redshift or instead from the standard $\Lambda$CDM cosmological model calibrated by measurements from the early Universe, has grown to a level of significance requiring a solution. Proposed systematic errors are not supported by the breadth of available data (and "unknown errors" untestable by lack of definition). Simple theoretical explanations for this "Hubble tension" that are consistent with the majority of the data have been surprisingly hard to come by, but in recent years, attention has focused increasingly on models that alter the early or pre-recombination physics of $\Lambda$CDM as the most feasible. Here, we describe the nature of this tension, emphasizing recent developments on the observational side. We then explain why early-Universe solutions are currently favored and the constraints that any such model must satisfy. We discuss one workable example, early dark energy, and describe how it can be tested with future measurements. Given an assortment of more extended recent reviews on specific aspects of the problem, the discussion is intended to be fairly general and understandable to a broad audience.

Spicules are plasma jets, observed in the dynamic interface region between the visible solar surface and the hot corona. At any given time, it is estimated that about 3 million spicules are present on the Sun. We find an intriguing parallel between the simulated spicular forest in a solar-like atmosphere and the numerous jets of polymeric fluids when both are subjected to harmonic forcing. In a radiative magnetohydrodynamic numerical simulation with sub-surface convection, solar global surface oscillations are excited similarly to those harmonic vibrations. The jets thus produced match remarkably well with the forests of spicules detected in observations of the Sun. Taken together, the numerical simulations of the Sun and the laboratory fluid dynamics experiments provide insights into the mechanism underlying the ubiquity of jets: the nonlinear focusing of quasi-periodic waves in anisotropic media of magnetized plasma as well as polymeric fluids under gravity is sufficient to generate a forest of spicules on the Sun.

We present a 6-D map of the Orphan-Chenab (OC) stream by combining the data from 5 years of Southern Stellar Stream Spectroscopic Survey $S^5$ observations with Gaia EDR3 data. We reconstruct the proper motion, radial velocity, distance and on-sky track of stream stars with spline models and extract the stellar density along the stream. The stream has a total luminosity of $M_V=-8.2$ and an average metallicity of $[Fe/H]=-1.9$, similar to classical MW satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1 kpc, a constant line of sight velocity dispersion of 5 km/s, but an increase in the velocity dispersion along the stream near pericenter to $\sim$ 10 km/s. Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6-D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended LMC potential. This allows us to constrain the mass profile of the MW within the distance range 15.6 < r < 55.5 kpc, with the best measured enclosed mass of $(2.85\pm 0.1)\times10^{11}\,M_\odot$ within 32.4 kpc. With the OC stream's closest approach distance to the LMC of $\sim 21$ kpc, our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of the LMC's enclosed mass being at 32.8 kpc with $M=(7.02\pm 0.9)\times10^{10}\, {M}_\odot$. We confidently detect that the LMC DM halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars on the OC stream show large energy and angular momentum spreads caused by the LMC perturbation and revealing the limitations of orbital invariants for substructure identification in the MW halo.

Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular structures observed in protoplanetary disks, starting with sub-micron-sized dust grains. We connect the models of dust coagulation and drift, planetesimal formation in the streaming instability, gravitational interactions between planetesimals, pebble accretion, and planet migration, into one uniform framework. We find that planetesimals forming early at the massive end of the size distribution grow quickly dominantly by pebble accretion. These few massive bodies grow on the timescales of ~100 000 years and stir the planetesimals formed later preventing the emergence of further planetary cores. Additionally, a migration trap occurs allowing for retention of the growing cores. Pressure bumps are favourable locations for the emergence and rapid growth of planetary cores by pebble accretion as the dust density and grain size are increased and the pebble accretion onset mass is reduced compared to a smooth-disk model.

Binary asteroids probe thermal-radiation effects on the main-belt asteroids' evolution. We discuss the possibility of detecting binary minor planet systems by the astrometric wobble of the center-of-light around the center-of-mass. This method enables the exploration of the phase-space of binary asteroids, which is difficult to explore using common detection techniques. We describe a forward model that projects the center-of-light position with respect to the center-of-mass, as it is seen by the observer. We study the performance of this method using simulated Gaia-like data. We apply the astrometric method to a subset of the Gaia DR2 Solar System catalog and find no significant evidence of binary asteroids. This is likely because the Gaia DR2 removed astrometric outliers, which in our case may be due to astrophysical signals. Applying this method to binary asteroid (4337) Arecibo, for which Gaia DR3 reported a possible astrometric signal with a period of P = 32.85+/-0.38 hr, reveals a possible 2.2-sigma solution with a period of 16.26 hr (about half the reported period). We find a small, marginally significant, excess of astrometric noise in the known binary asteroid population from Pravec et al. relative to the entire asteroid population in the Gaia DR2 Solar System catalog. We also discuss some caveats like precession and asteroid rotation.

Many energetic supernovae are thought to be powered by the rotational-energy of a highly-magnetized, rapidly-rotating neutron star. The emission from the associated luminous pulsar wind nebula (PWN) can photoionize the supernova ejecta, leading to a nebular spectrum of the ejecta with signatures possibly revealing the PWN. SN 2012au is hypothesized to be one such supernova. We investigate the impact of different ejecta and PWN parameters on the supernova nebular spectrum, and test if any photoionization models are consistent with SN 2012au. We study how constraints from the nebular phase can be linked into modelling of the diffusion phase and the radio emission of the magnetar. We present a suite of late-time (1-6y) spectral simulations of SN ejecta powered by an inner PWN. Over a large grid of 1-zone models, we study the behaviour of the SN physical state and line emission as PWN luminosity ($L_{\rm PWN}$), injection SED temperature ($T_{\rm PWN}$), ejecta mass ($M_{\rm ej}$), and composition (pure O or realistic) vary. We discuss the resulting emission in the context of the observed behaviour of SN 2012au, a strong candidate for a PWN-powered SN. The supernova nebular spectrum varies as $T_{\rm PWN}$ varies, as the ejecta become less ionized as $T_{\rm PWN}$ increases. Low ejecta mass models at high PWN power obtain runaway ionization for O I and, in extreme cases, also O II, causing a sharp decrease in their ion fraction over a small change in the parameter space. Certain models can reproduce the oxygen lines luminosities of SN 2012au reasonably well at individual epochs, but we find no model that fits over the whole time evolution; this is likely due to the simple model setup. Using our derived constraints from the nebular phase, we predict that the magnetar powering SN 2012au had an initial rotation period $\sim$ 15 ms, and should be a strong radio source (F > 100 mJy) for decades.

In galactic disks, the Parker instability results when non-thermal pressure support exceeds a certain threshold. The non-thermal pressures considered in the Parker instability are cosmic ray pressure and magnetic pressure. This instability takes a long time to saturate $(>500 \, \mathrm{Myr})$ and assumes a background with fixed cosmic ray pressure to gas pressure ratio. In reality, galactic cosmic rays are injected into localized regions $(< 100 \,\mathrm{pc})$ by events like supernovae, increasing the cosmic ray pressure to gas pressure ratio. In this work, we examine the effect of such cosmic ray injection on large scales $ (\sim 1\,\mathrm{kpc})$ in cosmic ray magnetohydrodynamic simulations using the \texttt{Athena++} code. We vary the background properties, dominant cosmic ray transport mechanism, and injection characteristics between our simulation runs. We find the injection will disrupt the interstellar medium on shorter timescales than the Parker instability. If cosmic ray transport by advection is dominant, cosmic ray injection disrupts the disk on short time scales $(<100\,\mathrm{Myr})$. If cosmic ray transport by the streaming instability is dominant, the injection creates a buoyant flux tube long after the initial injection $(>150\,\mathrm{Myr})$. Finally, when cosmic ray transport by diffusion dominates, the injected cosmic rays make an entire flux tube over pressured in a short time $(\sim 10 \, \mathrm{Myr})$. This over pressure pushes gas off the tube and drives buoyant rise on time scales similar to the advection dominated case.

We report the results of our continuous UBVRI-band photometry of Betelgeuse from 1999 to 2022 using the same photometric system. There are two advantages in our observation: (1) we used a photodiode as a detector to avoid saturation, and (2) our data set includes U-band light curve, which is not widely observed in recent CCD photometries. Using our light curves, we conducted the periodicity analysis, and found ~405- and ~2160-day periods. We also discuss the tentative detection of a long-period variation over 20 years or longer. Finally, we discuss the peculiar variation of the U-B color index during the "Great Dimming" event between late 2019 and early 2020.

We report on optical observations and modeling of HD96670, a single-line spectroscopic binary in the Carina OB2 association. We collected 10 epochs of optical spectroscopy, and optical photometry on 17 non-consecutive nights on the source. We construct a radial velocity curve from the spectra, and update the orbital period of the binary to be $P = 5.28388 \pm 0.00046$ days. The spectra show oxygen and helium absorption, consistent with an O-type primary. We see no evidence for spectral lines from the secondary star in the binary. We model the optical light curve and radial velocity curve simultaneously using the Wilson-Devinney code and find a best fit mass of $M_1 = 22.7^{+5.2}_{-3.6} M_\odot$ for the primary, and $M_2 = 6.2^{+0.9}_{-0.7} M_\odot$ for the secondary. An object of this mass is consistent with either a B-type star, or a black hole. Given that we see no absorption lines from the secondary, in combination with an observed hard power-law X-ray spectrum with $\Gamma = 2.6$ detected past 10 keV, maybe produced by wind accretion onto the secondary, we conclude that the secondary is most likely a black hole. We see asymmetrical helium lines with a shape consistent with the presence of a third star. If the secondary is indeed a black hole, this system would add to the small sample of only four possible black hole high mass X-ray binaries in the galaxy.

It is shown that the relativistic jets associated with the growth and past activity of the supermassive black hole in the Andromeda galaxy could be the main source of cosmic rays with energies above $10^{15}$ eV. Most of the cosmic ray energy is related to a bow shock of the jet that produces multi-PeV cosmic rays with light composition. The highest energy cosmic rays with heavy composition are produced in the jet itself. The spectra of energetic particles produced in Andromeda galaxy and propagated to the Earth are calculated and compared with observations.

Line intensity mapping (LIM) experiments probing the nearby universe can expect a considerable amount of cosmic infrared background (CIB) contiuum emission coming from near and far-infrared galaxies. For the purpose of using the LIM data to constrain the star formation rate (SFR), we argue that the CIB continuum - traditionally treated as contamination - can be combined with the LIM signal to enhance the SFR constraints achievable. We first present a power spectrum model that is capable of joining continuum and line emissions that assume the same prior SFR model. We subsequently analyze the effectiveness of the joint model in the context of the EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM), which utilizes the [CII] molecular line to study the SFR. We numerically compute the theoretical power spectra according to our model and the EXCLAIM survey specifics, and perform Fisher analysis to obtain SFR parameter constraints. We find that although the joint model has no considerable advantage over LIM alone assuming the current survey level of EXCLAIM, its effects become significant when we consider more optimistic values of survey resolution and angular span that are expected of future LIM experiments. By manipulating the Fisher formalism, we show that the CIB is not only an additional SFR sensitive signal, but also serves to break the SFR parameter degeneracy that naturally emerges from the [CII] Fisher matrix. For this reason, addition of the CIB will allow improvements in the survey parameters to be better reflected in the SFR constraints, and can be effectively utilized by future LIM experiments.

Close encounters between two initially unbound objects can result in a binary system if enough energy is released as gravitational waves (GWs). We address the scenario in which such encounters occur in merging elliptical galaxies. There is evidence that elliptical galaxies can harbor intermediate-mass black holes. Therefore, these systems are potentially the breeding grounds of sources of gravitational waves corresponding to inspiraling compact objects onto a massive black hole due to the dynamics, the large densities, and the number of compact remnants they contain. We show that this process is efficient for intermediate-mass black holes (IMBHs) with masses ranging from M $\in (10^3,10^5)$ M$_{\odot}$ and results in the formation of intermediate mass-ratio inspirals (IMRIs). We consider a set of IMBHs and smaller black holes with masses $m_2 \in (10,10^3)$ M$_{\odot}$ to estimate the IMRI formation rate. We find rates ranging between 10$^{-8}$ yr$^{-1}$, and 10$^{-5}$ yr$^{-1}$, and the IMRI formation rate per comoving volume in merging galaxies as a function of the redshift. The peak frequencies of the gravitational radiation emitted when these IMRIs are formed are within the detection band of space-borne detectors such as LISA and TianQin; taking into account the observable volume of these detectors, the total amount of IMRI detections per year is significant.

The Milky Way's stellar halo, which extends to $>100$ kpc, encodes the evolutionary history of our Galaxy. However, most studies of the halo to date have been limited to within a few kpc of the Sun. Here, we characterize differences between this local halo and the stellar halo in its entirety. We construct a composite stellar halo model by combining observationally motivated N-body simulations of the Milky Way's nine most massive disrupted dwarf galaxies that account for almost all of the mass in the halo. We find that (1) the representation by mass of different dwarf galaxies in the local halo compared to the whole halo can be significantly overestimated (e.g., the Helmi Streams) or underestimated (e.g., Cetus) and (2) properties of the overall halo (e.g., net rotation) inferred via orbit integration of local halo stars are significantly biased, because e.g., highly retrograde debris from Gaia-Sausage-Enceladus is missing from the local halo. Therefore, extrapolations from the local to the global halo should be treated with caution. From analysis of a sample of 11 MW-like simulated halos, we identify a population of recently accreted ($\lesssim5$ Gyrs) and disrupted galaxies on high angular momenta orbits that are entirely missing from local samples, and awaiting discovery in the outer halo. Our results motivate the need for surveys of halo stars extending to the Galaxy's virial radius.

Evidence of chemical disequilibria and other anomalous observations in the Venusian atmosphere motivate the search for life within the planet's temperate clouds. To find signs of a Venusian aerial biosphere, a dedicated astrobiological space mission is required. Venus Life Finder (VLF) missions encompass unique mission concepts with specialized instruments to search for habitability indicators, biosignatures and even life itself. A key in the search for life is direct capture, concentration and visualization of particles of biological potential. Here, we present a short overview of Fluid-Screen (FS) technology, a recent advancement in the dielectrophoretic (DEP) microbial particle capture, concentration and separation. FS is capable of capturing and separating biochemically diverse particles, including multicellular molds, eukaryotic cells, different species of bacteria and even viruses, based on particle dielectric properties. In this short communication, we discuss the possible implementation of Fluid-Screen in the context of the VLF missions, emphasizing the unique science output of the Fluid-Screen instrument. FS can be coupled with other highly sophisticated instruments such as an autofluorescence microscope or a laser desorption mass spectrometer. We discuss possible configurations of Fluid-Screen that upon modification and testing, could be adapted for Venus. We discuss the unique science output of the FS technology that can capture biological particles in their native state and hold them in the focal plane of the microscope for the direct imaging of the captured material. We discuss the challenges for the proposed method posed by the concentrated sulfuric acid environment of Venus' clouds. While Venus' clouds are a particularly challenging environment, other bodies of the solar system, e.g., with liquid water present, might be especially suitable for Fluid-Screen application.

For the upcoming PLAnetary Transits and Oscillation of stars (PLATO) satellite mission, a large number of target stars are required to yield a statistically significant number of planet transits. Locating the centres of the long duration observational phase (LOP) fields closer to the Galactic plane will increase the target star numbers but also the astrophysical false positives (FPs) from blended eclipsing binary systems. We utilise the Binary Stellar Evolution and Population Synthesis (BiSEPS) code, to create a complete synthetic stellar and planetary population for the proposed southern LOP field (LOPS0), as well as for a representative portion of the northern LOP field (LOPNsub). For LOPS0 we find an overall low FP rate for planets smaller than Neptunes. The FP rate generally shows little variation with Galactic longitude (l), and a modest increase with decreasing Galactic latitude (|b|). The location of the LOPS field centre within the current allowed region is not strongly constrained by FPs. Analysis of LOPNsub suggests a markedly increased number of FPs across the full range of planet radii at low |b| resulting in approximately twice the percent FP rate in the LOPNsub compared to the corresponding southern field segment in the planet radius range -0.2 < log(R/Rsun) <= 0.4. However, only a few percent of fully eclipsing FPs in LOPS0 in this radius range have periods

KELT-9b is an ultra-hot Jupiter observed to be undergoing extreme mass loss. Its A0-type host star has a radiative envelope, which makes its surface layers prone to retaining recently accreted material. To search for potential signs of planetary material polluting the stellar surface, we carry out the most comprehensive chemical characterisation of KELT-9 to-date. New element detections include Na and Y, which had previously been detected in the ultra-hot Jupiter but not studied in the star; these detections complete the set of nine elements measured in both star and planet. In comparing KELT-9 with similar open cluster stars we find no strong anomalies. This finding is consistent with calculations of photospheric pollution accounting for stellar mixing and using observationally estimated KELT-9b mass loss rates. We also rule out recent, short-lived intensive mass transfer such as the stellar ingestion of an Earth-mass exomoon.

Discoveries of planet- and stellar remnant-hosting pulsars challenge our understanding as the violent supernova explosion that forms the pulsar presumably destabilizes the system. Type II supernova explosions lead to the formation of eccentric bound systems, free-floating planets, neutron stars, pulsars, and white dwarfs. Analytical and numerical studies of high mass-loss rate systems based on perturbation theory so far have focused mainly on planet-star systems. In this paper, we extend our understanding of the fate of planet-star and binary systems by assuming a homologous envelope expansion model using a plausible ejection velocity ($1000-10000\,\mathrm{km/s}$), envelope- and neutron star masses. The investigation covers secondary masses of 1-10MJup for planetary, and 1-20MSun for stellar companions. We conduct and analyze over 2.5 million simulations assuming different semi-major axes (2.23 - 100au), eccentricities (0-0.8), and true-anomalies (0-2pi) for the companion. In a homologous expansion scenario, we confirm that the most probable outcome of the explosion is the destabilization of the system, while the retention of a bound system requires a highly eccentric primordial orbit. In general, a higher ejecta velocity results in a lower eccentricity orbit independent of secondary mass. The explanation of close-in pulsar planets requires exotic formation scenarios, rather than survival through the type II supernova explosion model. Post-explosion bound star systems gain a peculiar velocity (<100\,km/s), even though the explosion model is symmetric. The applied numerical model allows us to derive velocity components for dissociating systems. The peculiar velocities of free-floating planets and stellar corpses are in the range of 10^-6-275km/s.

Dual AGN are important for understanding galaxy-merger-triggered fueling of black holes and hierarchical growth of structures. The least explored type of dual AGN are those associated with mergers of two dwarf galaxies. According to observations and cosmological simulations, dwarf galaxies are the most abundant type of galaxies in the early Universe and the galaxy merger rate is dominated by dwarfs. However, these mergers are generally too distant to be directly observed, and low-redshift dwarf-dwarf merger-related dual AGN are notoriously hard to find. In this paper, we present the first results of our large-scale search for this elusive type of object and the first two candidates for dual AGN in dwarf-dwarf mergers. Both objects exhibit tidal features (tails and bridges) characteristic of galaxy mergers/interactions. One object is apparently in a late-stage merger with an AGN separation of < 5kpc, while the second is in an early-stage merger with interacting galaxies having established a tidal bridge. Both objects have dual, luminous X-ray sources that are most likely due to actively accreting massive black holes. Also, both objects have infrared counterparts, with colors consistent with being AGN. Follow-up observations will provide us a glimpse into key processes that govern the earliest phases of growth of galaxies, their central black holes, and merger-induced star formation.

We present a subgrid model for supernova feedback designed for simulations of galaxy formation. The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the OWLS and EAGLE simulations, respectively. In the thermal channel, the energy is distributed statistically isotropically and injected stochastically in large amounts per event, which minimizes spurious radiative energy losses. In the kinetic channel, we inject the energy in small portions by kicking gas particles in pairs in opposite directions. The implementation of kinetic feedback is designed to conserve energy, linear momentum and angular momentum, and is statistically isotropic. To test and validate the model, we run simulations of isolated Milky Way-mass and dwarf galaxies, in which the gas is allowed to cool down to 10 K. Using the thermal and kinetic channels together, we obtain smooth star formation histories and powerful galactic winds with realistic mass loading factors. Furthermore, the model produces spatially resolved star formation rates and velocity dispersions that are in agreement with observations. We vary the numerical resolution by several orders of magnitude and find excellent convergence of the global star formation rates and the mass loading of galactic winds. We show that large thermal-energy injections generate a hot phase of the interstellar medium (ISM) and modulate the star formation by ejecting gas from the disc, while the low-energy kicks increase the turbulent velocity dispersion in the neutral ISM, which in turn helps suppress star formation.

We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e. $\rho<10^5\ {\rm cm^{-3}}$), including photoheating, photoelectric heating, H$_2$ heating/cooling, cosmic-ray heating, H/He cooling, metal-line cooling, CO cooling, and dust cooling (recombination and gas-grain collisions). We validate the model by comparing 1D equilibrium simulations across six dex in metallicity to existing 1D ISM models in the literature. We apply PRISM to high-resolution (4.5 pc) isolated dwarf galaxy simulations that include state-of-the-art models for star formation and stellar feedback to take an inventory of which cooling and heating processes dominate each different gas phase of a galaxy and to understand the importance of non-equilibrium effects. We show that most of the ISM gas is either close to thermal equilibrium or exhibits a slight cooling instability, while from a chemical perspective, the non-equilibrium electron fraction is often more than three times higher or lower than the equilibrium value, which impacts cooling, heating, and observable emission lines. Electron enhancements are attributed to recombination lags while deficits are shown to be due to rapid cosmic-ray heating. The PRISM model and its coupling to RAMSES-RTZ is applicable to a wide variety of astrophysical scenarios, from cosmological simulations to isolated giant molecular clouds, and is particularly useful for understanding how changes to ISM physics impact observable quantities such as metallic emission lines.

Tidal disruption events (TDEs) show a correlation between the UV to X-ray spectral index and the Eddington ratio, with non-thermal X-ray emission at the low Eddington ratio. We consider the corona surrounding the accretion disc as a non-thermal X-ray source. We construct a time-dependent and non-relativistic advective accretion disc-corona model for TDEs. The infalling debris is assumed to form a seed disc in time $t_c$, that evolves due to the mass gain from the infalling debris at the constant outer radius with a mass fallback rate $\dot{M}_{\rm fb}$ and the mass loss through accretion onto the black hole. The viscous stress in our model depends on gas ($P_g$) and total ($P_t$) pressures as $\tau_{r\phi} \propto P_g^{1-\mu} P_t^{\mu}$, where $\mu$ is a constant. We find that the mass accretion rate $\dot{M}_a$ evolves from Eddington to sub-Eddington accretion with a late-time evolution close to $t^{-5/3}$, where $t$ is the time. We find that the bolometric disc luminosity follows a late-time evolution close to $t^{-5/3}$. The ratio of total X-ray luminosity from corona to bolometric disc luminosity increases with $\mu$ and increases at late times for $\mu \neq 1$. We obtain the X-ray blackbody temperature of the disc that agrees with the temperature from X-ray observations ($\sim~10^5~{\rm K}$). We find the radiative efficiency of the disc increases with time and decreases for a disc when the corona is included. We have neglected the outflow, and our model is more applicable for near-to-sub-Eddington accretion and when $\dot{M}_{\rm fb}$ is sub-Eddington.

The ionospheric path delay impacts single-band very long baseline interferometry (VLBI) group delays, which limits their applicability for absolute astrometry. I consider two important cases: when observations are made simultaneously at two bands, but delays at only one band are available for a subset of observations and when observations are made at one band only by design. I developed optimal procedures of data analysis for both cases using Global Navigation Satellite System (GNSS) ionosphere maps, provided a stochastic model that describes ionospheric errors, and evaluated their impact on source position estimates. I demonstrate that the stochastic model is accurate at a level of 15%. I found that using GNSS ionospheric maps as is introduces serious biases in estimates of declinations and I developed the procedure that almost eliminates them. I found serendipitously that GNSS ionospheric maps have multiplicative errors and have to be scaled by 0.85 in order to mitigate the declination bias. A similar scale factor was found in comparison of the vertical total electron contents from satellite altimetry against GNSS ionospheric maps. I favor interpretation of this scaling factor as a manifestation of the inadequacy of the thin shell model. I showed in this study that we are able to model the ionospheric path delay to the extent that no systematic errors emerge and we are able to adequately assess the contribution of the ionosphere-driven random errors on source positions. This makes single-band absolute astrometry a viable option that can be used for source position determination.

The temperature predicted by photoionization models for the Narrow Line Region of Seyfert 2 galaxies is lower than the value inferred from the observed [O III] {\lambda}4363A/{\lambda}5007A line ratio. We explore the possibility of considering a harder ionizing continuum than typically assumed. The spectral ionizing energy distribution, which can generate the observed {\lambda}4363A/{\lambda}5007A ratio, is characterized by a secondary continuum peak at 200 eV.

Active galactic nuclei (AGNs) make up about 35 per cent of the more than 250 sources detected in very-high-energy (VHE) gamma rays to date with Imaging Atmospheric Cherenkov Telescopes. Apart from four nearby radio galaxies and two AGNs of unknown type, all known VHE AGNs are blazars. Knowledge of the cosmological redshift of gamma-ray blazars is key to enabling the study of their intrinsic emission properties, as the interaction between gamma rays and the extragalactic background light (EBL) results in a spectral softening. Therefore, the redshift determination exercise is crucial to indirectly placing tight constraints on the EBL density and to studying blazar population evolution across cosmic time. Due to the powerful relativistic jets in blazars, most of their host galaxies' spectral features are outshined, and dedicated high signal-to-noise spectroscopic observations are required. Deep medium- to high-resolution spectroscopy of 33 gamma-ray blazar optical counterparts was performed with the European Southern Observatory New Technology Telescope, Keck II telescope, Shane 3-meter telescope and the Southern African Large Telescope. From the sample, spectra from 25 objects display spectral features or are featureless and have high signal-to-noise. The other eight objects have low quality featureless spectra. We systematically searched for absorption and emission features and estimated, when possible, the fractional host galaxy flux in the measured total flux. Our measurements yielded 14 firm spectroscopic redshifts, ranging from 0.0838 to 0.8125, one tentative redshift, and two lower limits: one at z > 0.382 and the other at z > 0.629.

Empirical correlations between various key parameters have been extensively explored ever since the discovery of gamma-ray bursts (GRBs) and have been widely used as standard candles to probe the Universe. The Amati relation and the Yonetoku relation are two good examples, which have been paid special attention to. The former reflects the connection between the peak photon energy (Ep) and the isotropic $\gamma$-ray energy release (Eiso), while the latter links Ep with the isotropic peak luminosity (Lp), both in the form of a power law function. Most GRBs are found to well follow these correlations, but a theoretical interpretation is still lacking. Meanwhile, there are also some obvious outliers, which may be off-axis GRBs and may follow different correlations as compared with the on-axis ones. Here we present a simple analytical derivation for the Amati relation and the Yonetoku relation in the framework of the standard fireball model, the correctness of which are then confirmed by numerical simulations. The off-axis Amati relation and Yonetoku relation are also derived, which differ from the corresponding on-axis relation markedly. Our results reveal the intrinsic physics lying behind the radiation processes of GRBs, and highlight the importance of viewing angle in the empirical correlations of GRBs.

One of the puzzles to have emerged from the Kepler and TESS missions is the existence of unexplained dips in the lightcurves of a small fraction of rapidly-rotating M dwarfs in young open clusters and star-forming regions. We present a theoretical investigation of one possible explanation - that these are caused by dust clouds trapped in the stellar magnetic fields. The depth and duration of the observed dips allow us to estimate directly the linear extent of the dust clouds and their distances from the rotation axis. The dips are found to be between 0.4-4.8%. We find that their distance is close to the co-rotation radius: the typical location for stable points where charged particles can be trapped in a stellar magnetosphere. We estimate the charge acquired by a dust particle due to collisions with the coronal gas and hence determine the maximum grain size that can be magnetically supported, the stopping distance due to gas drag and the timescale on which dust particles can diffuse out of a stable point. Using the observationally-derived magnetic field of the active M dwarf V374 Peg, we model the distribution of these dust clouds and produce synthetic light curves. We find that for 1 micron dust grains, the light curves have dips of 1% - 3% and can support masses of order of $10^{12}$kg. We conclude that magnetically-trapped dust clouds (potentially from residual disc accretion or tidally-disrupted planetesimal or cometary bodies) are capable of explaining the periodic dips in the Kepler and TESS data.

The study on the dynamic evolution of young supernova remnants (SNRs) is an important way to understand the density structure of the progenitor's circumstellar medium. We have reported the acceleration or deceleration, proper motion and brightness changes of 260 compact radio features in the second youngest known SNR Cas A at 5\,GHz based on the VLA data of five epochs from 1987 to 2004. The radio expansion center locates at $\alpha(1950)=23^{\rm h}21^{\rm m}9^{\rm s}_{\cdot}7 \pm 0^{\rm s}_{\cdot}29, \delta(1950)=+58^{\circ}32^{\prime}25^{\prime\prime}_{\cdot}2 \pm 2^{\prime\prime}_{\cdot}2$. Three-quarters of the compact knots are decelerating, this suggests that there are significant density fluctuations in the stellar winds of the remnant's progenitor. We have verified that the acceleration or deceleration of compact knots are not related with the distribution of brightness. The brightening, fading, disappearing or new appearing of compact radio features in Cas A suggests that the magnetic field in the remnant is changing rapidly.

SN 2018ivc is an unusual type II supernova (SN II). It is a variant of SNe IIL, which might represent a transitional case between SNe IIP with a massive H-rich envelope, and IIb with only a small amount of the H-rich envelope. However, SN 2018ivc shows an optical light curve evolution more complicated than canonical SNe IIL. In this paper, we present the results of prompt follow-up observations of SN 2018ivc with the Atacama Large Millimeter/submillimeter Array (ALMA). Its synchrotron emission is similar to that of SN IIb 1993J, suggesting that it is intrinsically an SN IIb-like explosion of a He star with a modest (~0.5 - 1 Msun) extended H-rich envelope. Its radio, optical, and X-ray light curves are explained primarily by the interaction between the SN ejecta and the circumstellar material (CSM); we thus suggest that it is a rare example (and the first involving the `canonical' SN IIb ejecta) for which the multi-wavelength emission is powered mainly by the SN-CSM interaction. The inner CSM density, reflecting the progenitor activity in the final decade, is comparable to that of SN IIb 2013cu that showed a flash spectral feature. The outer CSM density, and therefore the mass-loss rate in the final ~200 years, is larger than that of SN 1993J by a factor of ~5. We suggest that SN 2018ivc represents a missing link between SNe IIP and IIb/Ib/Ic in the binary evolution scenario.

Millimetre-wave astronomy is an important tool for cosmology. In this context, the Line Intensity Mapping (LIM) technique has been proposed to map in three dimensions the specific intensity due to line (e.g. [CII], CO) emission. This mapping is particularly interesting to study the primordial galaxies as a function of redshift. LIM observations are typically carried out on single dish telescopes. Hyper-spectral integrated devices have the potential to replace the current Fourier transform, or the planned Fabry-Perot based instruments to achieve efficient, i.e. large field-of-view, line intensity mapping at millimetre and sub-millimetre wavelengths. The aim is to perform hyper-spectral mapping, with a spectral resolution 100-1000, over large, i.e. thousands of beams, instantaneous patches of the Sky. The device that we have developed allows avoiding moving parts, complicated or dispersive optics or tunable filters to be operated at cryogenics temperatures. We designed, fabricated and tested an innovative integrated hyperspectral device sensitive in the 80-90GHz range. It contains nineteen horns and sixteen spectral-imaging channels, each selecting a frequency band of about 0.1GHz. A conical horn antenna, coupled to a planar superconducting filter, collects the radiation. A capacitively coupled Lumped Element Kinetic Inductance Detector (LEKID) is then in charge of dissipating and sensing the super-current established in the filter. The prototype is fabricated with only two photo-lithography steps over a commercial mono-crystalline sapphire substrate. It exhibits a spectral resolution R=800. The optical noise equivalent power is in the observational relevant 1e-17 W per sqrt(Hz) range. The device is polarisation-sensitive, paving the way to spectro-polarimetry measurements over very large field-of-views.

A large-scale, horse-shoe-like filament was investigated and the magnetic field around it was reconstructed. This is an intermediate filament (IF) that appeared on the solar disk for the first time at 02:00 UT on 2015 November 7, and took 8 days to move to the central median on the solar disk. The active region AR 12452 around which the filament occurred was diffuse so that the magnetic field nearby was weak, the average field strength is 106 G. Therefore, the existing approaches to extrapolating the coronal magnetic field and to constructing the filament configuration in the region with strong background field do not work well here. On the basis of the regularized Biot-Savart laws method, we successfully constructed a data-constrained, non-linear force-free field configuration for this IF observed on 2015 November 14. The overall IF configuration obtained in this way matches well the morphology suggested by a 304~\AA \ image taken by the Atmospheric Imaging Assembly on board Solar Dynamics Observatory. Magnetic dips in the configuration were coincident in space with the H$\alpha$ features of the filament, which is lower in altitude than the features seen in 304~\AA. This suggests that the cold plasma fills the lower part of the filament, and hot plasma is situated in the higher region. A quasi-separatrix layer wraps the filament, and both the magnetic field and the electric current are stronger near the inner edge of the filament.

We study the role the the p-mode-like vertical oscillation on the photosphere in driving solar winds in the framework of Alfven-wave-driven winds. By performing one-dimensional magnetohydrodynamical numerical simulations from the photosphere to the interplanetary space, we discover that the mass-loss rate is raised up to 4 times as the amplitude of longitudinal perturbations at the photosphere increases. When the longitudinal fluctuation is added, transverse waves are generated by the mode conversion from longitudinal waves in the chromosphere, which increases Alfvenic Poynting flux in the corona. As a result, the coronal heating is enhanced to yield higher coronal density by the chromospheric evaporation, leading to the increase of the mass-loss rate. Our findings clearly show the importance of the p-mode oscillation in the photosphere and the mode conversion in the chromosphere in determining the basic properties of the wind from the sun and solar-type stars.

Once a cometary plasma cloud has been created through ionisation of the cometary neutrals, it presents an obstacle to the solar wind and the magnetic field within it. The acceleration and incorporation of the cometary plasma by the solar wind is a complex process that shapes the cometary plasma environment and is responsible for the creation of boundaries such as a bow shock and diamagnetic cavity boundary. It also gives rise to waves and electric fields which in turn contribute to the acceleration of the plasma. This chapter aims to provide an overview of how the solar wind is modified by the presence of the cometary plasma, and how the cometary plasma is incorporated into the solar wind. We will also discuss models and techniques widely used in the investigation of the plasma environment in the context of recent findings by Rosetta. In particular, this chapter highlights the richness of the processes and regions within this environment and how processes on small scales can shape boundaries on large scales. It has been fifteen years since the last book on Comets was published and since then we have made great advances in the field of cometary research. But many open questions remain which are listed and discussed with particular emphasis on how to advance the field of cometary plasma science through future space missions.

Theoretical models predict that the obscuration of stellar radiation by irregularities on the surface of a protoplanetary disk can cause self-generating waves traveling towards the star. However, this process is traditionally simulated using the 1+1D approach, the key approximations of which - vertical hydrostatic equilibrium of the disk and vertical diffusion of IR radiation - can distort the picture. This article presents a two-dimensional radiative hydrodynamic model of the evolution of an axially symmetric gas and dust disk. Within this model, but using simplified assumptions from 1+1D models, we have reproduced the spontaneous generation and propagation of thermal surface waves. The key conclusion of our work is that taking into account two-dimensional hydrodynamics and diffusion of IR radiation suppresses the spontaneous generation and development of thermal waves observed in the 1+1D approximation. The search for the possibility of the existence of surface thermal waves should be continued by studying the problem for various parameters of protoplanetary disks.

The Hubble diagram of type-Ia supernovae (SNe-Ia) provides cosmological constraints on the nature of dark energy with an accuracy limited by the flux calibration of currently available spectrophotometric standards. The StarDICE experiment aims at establishing a 5-stage metrology chain from NIST photodiodes to stars, with a targeted accuracy of \SI{1}{mmag} in $griz$ colors. We present the first two stages, resulting in the calibration transfer from NIST photodiodes to a demonstration \SI{150}{Mpixel} CMOS sensor (Sony IMX411ALR as implemented in the QHY411M camera by QHYCCD). As a side-product, we provide full characterization of this camera. A fully automated spectrophotometric bench is built to perform the calibration transfer. The sensor readout electronics is studied using thousands of flat-field images from which we derive stability, high resolution photon transfer curves and estimates of the individual pixel gain. The sensor quantum efficiency is then measured relative to a NIST-calibrated photodiode. Flat-field scans at 16 different wavelengths are used to build maps of the sensor response. We demonstrate statistical uncertainty on quantum efficiency below \SI{0.001}{e^-/\gamma} between \SI{387}{nm} and \SI{950}{nm}. Systematic uncertainties in the bench optics are controlled at the level of \SI{1e-3}{e^-/\gamma}. Uncertainty in the overall normalization of the QE curve is 1\%. Regarding the camera we demonstrate stability in steady state conditions at the level of \SI{32.5}{ppm}. Homogeneity in the response is below \SI{1}{\percent} RMS across the entire sensor area. Quantum efficiency stays above \SI{50}{\percent} in most of the visible range, peaking well above \SI{80}{\percent} between \SI{440}{nm} and \SI{570}{nm}. Differential non-linearities at the level of \SI{1}{\percent} are detected. A simple 2-parameter model is proposed to mitigate the effect.

The fields of cosmic ray astrophysics, gamma-ray astrophysics, and neutrino astrophysics have diverged somewhat. But for the air showers in the GeV and TeV energy ranges, the ground-based detector techniques have considerable overlaps. VHE gamma-ray astronomy is the observational study measuring the directions, flux, energy spectra, and time variability of the sources of these gamma rays. With the low flux of gamma rays, and the background of charged particle cosmic rays, the distinguishing characteristic of gamma-ray air shower detectors is large size and significant photon to charge particle discrimination. Air shower telescopes for gamma-ray astronomy consist of an array of detectors capable of measuring the passage of particles through the array elements. To maximize signal at energies of a TeV or so, the array needs to be built at high altitude as the maximum number of shower particles is high in the atmosphere. These detectors have included sparse arrays of shower counters, dense arrays of scintillators or resistive plate counters (RPC), buried muon detectors in concert with surface detectors, or many-interaction-deep Water Cherenkov Detectors (WCD). In general these detectors are sensitive over a large field of view, the whole of the sky is a typical sensitivity and perhaps 2/3 of the sky selected for clean analysis, but with only moderate resolution in energy, typically due to shower-to-shower fluctuations and the intrinsic sampling of the detector. These telescopes though, operate continuously, despite weather, moonlight, day or night, and without needing to be pointed to a specific target for essentially a 100\% duty cycle. In this chapter we will examine the performance and characteristics of such detectors. These are contrasted with the Imaging Air Cherenkov Telescopes which also operate in this energy range, and both current and future proposed experiments are described.

Context: The orientation and rotation of Mars, which can be described by a set of Euler angles, is estimated from radioscience data and is then used to infer Mars internal properties. The data are analyzed using a modeling expressed within the Barycentric Celestial Reference System (BCRS). Aims: We provide new and more accurate (to the $0.1$ mas level) estimations of the relativistic corrections to be included in the BCRS model of the orientation and rotation of Mars to avoid a misinterpretation of the data. Methods: There are two types of relativistic contributions in Mars rotation and orientation: (i) those that directly impact the Euler angles and (ii) those resulting from the time transformation between a local Mars reference frame and BCRS. The former correspond essentially to the geodetic effect. We compute them assuming that Mars evolves on a Keplerian orbit. As for the latter, we compute the effect of the time transformation and compare the rotation angle corrections obtained using realistic orbits as described by ephemerides. Results: The relativistic correction in longitude comes mainly from the geodetic effect and results in the geodetic precession (6.754mas/yr) and the geodetic annual nutation (0.565 mas amplitude). For the rotation angle, the correction is dominated by the effect of the time transformation. The main annual, semi-annual, and ter-annual terms have amplitudes of 166.954 mas, 7.783 mas, and 0.544mas, respectively. The amplitude of the annual term differs by about 9 mas from the estimate usually considered by the community. We identify new terms at the Mars-Jupiter and Mars-Saturn synodic periods (0.567 mas and 0.102 mas amplitude) that are relevant considering the current level of uncertainty of the measurements, as well as a contribution to the rotation rate (7.3088 mas/day). There is no significant correction that applies to the obliquity.

Giant planets around young stars serve as a clue to unveiling their formation history and orbital evolution. CI Tau is a 2\,Myr-old classical T-Tauri star hosting an eccentric hot Jupiter, CI Tau\,b. The standard formation scenario of a hot Jupiter predicts that planets formed further out and migrated inward. A high eccentricity of CI Tau b may be suggestive of high-$e$ migration due to secular gravitational perturbations by an outer companion. Also, ALMA 1.3\,mm-continuum observations show that CI Tau has at least three annular gaps in which unseen planets may exist. We present high-contrast imaging around CI Tau taken from Keck/NIRC2 $L^{\prime}$-band filter and vortex coronagraph that allows us to search for an outer companion. We did not detect any outer companion around CI Tau from angular differential imaging (ADI) using two deep imaging data sets. The detection limits from ADI-reduced images rule out the existence of an outer companion beyond $\sim30$\,au that can cause the Kozai-Lidov migration of CI Tau\,b. Our results suggest that CI Tau\,b may have experienced Type II migration from $\lesssim 2$\,au in Myrs. We also confirm that no planets with $\geq 2-4\,M_\mathrm{Jup}$ are hidden in two outer gaps.

Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if such mergers have very short delay times and/or beneficial masses of the compact objects at low metallicities. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models.

The so-called great dimming event of alpha Ori in late 2019 and early 2020 sparked our interest in the behaviour of chromospheric activity during this period. To study the timeline of chromospheric activity, we derive a S_MWO time series of TIGRE and Mount Wilson values, and we compare this long time series with photometric data from the AAVSO database. In addition, we determine the absolute and normalised excess flux of the Ca II H&K lines. To do so, we estimate the changing effective temperature from TIGRE spectra and find a clear drop of about 80 K between November 2019 and February 2020, which coincides with the minimum of visual brightness. During the same period, the S-index increased significantly, yet this is a mere contrast effect, because the normalised excess flux of the Ca II H&K lines did not change significantly. However, the latter dropped immediately after this episode. Comparing the combined S_MWO values and visual magnitude time series, we find a similar increase in the S-index during another noticeable decrease in the visual magnitude of alpha Ori, which took place in 1984 and 1985. To also probe the dynamics of the upper photosphere, we analysed the lines in 6251-6263 A and found core distance varies which shows a relation with the great dimming event. This type of variation could be caused by rising and sinking cool plumes as a temporary spill-over of convection on alpha Ori. Based on our study, we conclude that the cause for the great dimming is located in the photosphere. Furthermore, the long-term spectroscopic and photometric time series suggests that this great dimming does not appear to be a unique phenomenon, but rather that such dimmings do occur more frequently, which motivates further monitoring of alpha Ori with facilities such as TIGRE.

We identify and investigate a possible correlation between the $\rm{[CII]} 158{\mu}m$ luminosity and linewidth in the $\rm{[CII]}$-detected galaxies. Observationally, the strength of the $\rm{[CII]} 158{\mu}m$ emission line is usually stronger than that of the CO emission line and this $\rm{[CII]}$ line has been used as another tracer of the galactic characteristics. Moreover, many $\rm{[CII]}$-detected galaxies are identified in $z > 4$. Motivated by previous studies of the CO luminosity - FWHM correlation relation (LFR) and the available new $\rm{[CII]}$ measurements, we compile samples of the $\rm{[CII]}$-detected galaxies in the literature and perform the linear regression analysis. The $\rm{[CII]}$ LFR is confirmed at a robust level. We also demonstrate the possible application of the $\rm{[CII]}$ LFR by utilizing it on the distance measurement of the high-$z$ galaxy. As a result, we extend the cosmic spatial scale beyond the redshift $z$ of $7$. With the outcome of the distance measurement, we constrain the cosmology parameters in the Chevallier-Polarski-Linder model, which considers the evolution of dark energy. Consequently, the uncertainties of the $\textit{w}_{0}$ and $\textit{w}_{a}$ are reduced significantly when the measured distance data of the $\rm{[CII]}$-detected galaxies are included in the cosmological parameter constraint, exemplifying the potential of using the $\rm{[CII]}$-detected galaxies as a tracer to constrain the cosmological parameters.

Phosphine could be a key molecule in the understanding of exotic chemistry happening in (exo)planetary atmospheres. While it has been detected in the Solar System's giant planets, it has not been observed in exoplanets yet. In the exoplanetary context however it has been theorized as a potential biosignature molecule. The goal of our study is to identify which illustrative science cases for PH3 chemistry are observable with a space-based mid-infrared nulling interferometric observatory like the LIFE (Large Interferometer For Exoplanets) concept. We identified a representative set of scenarios for PH3 detections in exoplanetary atmospheres varying over the whole dynamic range of the LIFE mission. We used chemical kinetics and radiative transfer calculations to produce forward models of these informative, prototypical observational cases for LIFEsim, our observation simulator software for LIFE. In a detailed, yet first order approximation it takes a mission like LIFE: (i) about 1h to find phosphine in a warm giant around a G star at 10 pc, (ii) about 10 h in H2 or CO2 dominated temperate super-Earths around M star hosts at 5 pc, (iii) and even in 100h it seems very unlikely that phosphine would be detectable in a Venus-Twin with extreme PH3 concentrations at 5 pc. Phosphine in concentrations previously discussed in the literature is detectable in 2 out of the 3 cases and about an order of magnitude faster than comparable cases with JWST. We show that there is a significant number of objects accessible for these classes of observations. These results will be used to prioritize the parameter range for the next steps with more detailed retrieval simulations. They will also inform timely questions in the early design phase of a mission like LIFE and guide the community by providing easy-to-scale first estimates for a large part of detection space of such a mission.

This document describes the Netherlands' decadal strategic planning process for the current decade. We give the scientific rationale for our prioritization of research areas and the facility choices that follow from our scientific priorities. We also describe actions needed for the sustainability of our community and our work, and the budgets needed to fulfil our stated ambitions. The names listed as authors are in fact the editors of this paper, which results from the work of the entire Netherlands astronomy community.

We explore the possibility of using deep learning to generate multifield images from state-of-the-art hydrodynamic simulations of the CAMELS project. We use a generative adversarial network to generate images with three different channels that represent gas density (Mgas), neutral hydrogen density (HI), and magnetic field amplitudes (B). The quality of each map in each example generated by the model looks very promising. The GAN considered in this study is able to generate maps whose mean and standard deviation of the probability density distribution of the pixels are consistent with those of the maps from the training data. The mean and standard deviation of the auto power spectra of the generated maps of each field agree well with those computed from the maps of IllustrisTNG. Moreover, the cross-correlations between fields in all instances produced by the emulator are in good agreement with those of the dataset. This implies that all three maps in each output of the generator encode the same underlying cosmology and astrophysics.

Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region.

The G230LB grating used with STIS's CCD detector scatters red light. In red objects, the scattered light mingles with the ultraviolet signal, causing spurious short-wavelength flux and weakening absorption features. Recent calibration observations characterize the scattered light using duplicate observations with the MAMA detector and similar grating G230L. The full two-dimensional spectrum contains little helpful information to mitigate the scattered light problem. For one-dimensional, extracted spectra, the scattered light can be approximately modeled as a ramped pedestal whose amplitude is proportional to the object's V-band flux. We present formulae for scattered light corrections. For stars warmer than G0 spectral type, correction is superfluous. Off-slit-center positioning appears not to affect the properties of the scattered light. Therefore, we are able to extrapolate correction formulae for extended objects from the point source formulae. Polynomials for flux corrections due to off-center slit positioning in the 0.2 arcsec slit are also tabulated.

Short-term variability and multiple velocity components in the powerful highly ionized wind of the archetypal UFO PG 1211+143 are indicative of inner disc instabilities or short-lived accretion events. The detection of a high velocity (~ 0.3c) inflow of highly ionized matter, located at 20 R_g, offered the first direct observational support for the latter scenario, where matter approaching at a high inclination to the black hole spin plane may result in warping and tearing of the inner accretion disc, with subsequent inter-ring collisions producing shocks, loss of rotational support and rapid mass infall. Simultaneous soft x-ray spectra reveal a lower velocity (~ 0.1c) inflow of less ionized matter, identified as 'upstream' at 200 R_g, with a line of sight through matter converging on the supermassive black hole. We discuss here why ultrafast ionized winds are relatively common in luminous Seyfert galaxies, while detection of the 0.3c inflow in PG 1211+143 remains a rare example.

In the procedure of constraining the cosmological parameters with the observational Hubble data and the type Ia supernova data, the combination of Masked Autoregressive Flow and Denoising Autoencoder can perform a good result. The above combination extracts the features from OHD with DAE, and estimates the posterior distribution of cosmological parameters with MAF. We ask whether we can find a better tool to compress large data in order to gain better results while constraining the cosmological parameters. Information maximising neural networks, a kind of simulation-based machine learning technique, was proposed at an earlier time. In a series of numerical examples, the results show that IMNN can find optimal, non-linear summaries robustly. In this work, we mainly compare the dimensionality reduction capabilities of IMNN and DAE. We use IMNN and DAE to compress the data into different dimensions and set different learning rates for MAF to calculate the posterior. Meanwhile, the training data and mock OHD are generated with a simple Gaussian likelihood, the spatially flat {\Lambda}CDM model and the real OHD data. To avoid the complex calculation in comparing the posterior directly, we set different criteria to compare IMNN and DAE.

Integrated-light models that incorporate common types of chemically peculiar (CP) stars are assembled using synthetic spectra. Selected spectral features encode significant age information for populations with ages $\sim$50 Myr $<$ age $< \sim$2 Gyr. Due to the alleviation of template mismatch, the inclusion of CP star features in model spectra improves the accuracy of recovered stellar population parameters, but we are not able to show that new or unique age information can be extracted from the weak CP features compared to continuum fitting and strong-feature strengths, at least at the present state of the art. An age-extraction routine that recovers 2- and 3-burst age structures is employed to analyze the spectra of local group galaxies. NGC 224 (M31) has a stellar population too old for the types of CP stars we examine. NGC 221 (M32) also shows no CP spectral features. It appears to contain a component at age $\sim$1 Gyr at 1% by mass in addition to its dominant 4.7 Gyr population. Unlike SDSS galaxy spectrum averages, NGC 205 (M110) contains no features due to HgMn stars. This excludes the age range associated with HgMn production, and its near-nuclear spectrum is best fit by a 68$\pm$2 Myr population superimposed on an older population with a 1.85$\pm$0.1 Gyr component. Both NGC 205 and NGC 221 have an ancient component whose mass is not easy to constrain given the overwhelming light-dominance of the younger populations.

Radio interferometers designed to probe the 21cm signal from Cosmic Dawn and the Epoch of Reionisation must contend with systematic effects that make it difficult to achieve sufficient dynamic range to separate the 21cm signal from foreground emission and other effects. For instance, the instrument's chromatic response modulates the otherwise spectrally smooth foregrounds, making them difficult to model, while a significant fraction of the data must be excised due to the presence of radio frequency interference (RFI), leaving gaps in the data. Errors in modelling the (modulated and gappy) foregrounds can easily generate spurious contamination of what should otherwise be 21cm signal-dominated modes. Various approaches have been developed to mitigate these issues by (e.g.) using non-parametric reconstruction of the foregrounds, in-painting the gaps, and weighting the data to reduce the level of contamination. We present a Bayesian statistical method that combines these approaches, using the coupled techniques of Gaussian constrained realisations (GCR) and Gibbs sampling. This provides a way of drawing samples from the joint posterior distribution of the 21cm signal modes and their power spectrum in the presence of gappy data and an uncertain foreground model in a computationally scalable manner. The data are weighted by an inverse covariance matrix that is estimated as part of the inference, along with a foreground model that can then be marginalised over. We demonstrate the application of this technique on a simulated HERA-like delay spectrum analysis, comparing three different approaches for accounting for the foreground components.

We present CoSANG (Coupling Semi-Analytic and N-body Galaxies), a new hybrid model for cosmological dark matter and stellar halo simulations. In this approach a collisionless model (Gadget3) is used for gravitational interactions while a coupled semi-analytic model (SAGE) calculates baryonic effects at each time-step. This live self-consistent interaction at each time-step is the key difference between CoSANG and traditional semi-analytic models that are mainly used for post-processing. By accounting for the gravitational effect of the baryons, CoSANG can overcome some of the deficiencies of pure N-body simulations, while being less computationally expensive than hydrodynamic simulations. Moreover CoSANG can produce stellar halo populations via tagging tracer dark matter particles. We demonstrate the performance and dynamical accuracy of this approach using both controlled test simulations and a set of three cosmological zoom-in simulations of Milky Way (MW) mass halos. We simulate each target halo both without the coupling (dark matter only, hereafter DMO) and with the coupling (CoSANG). We compare the internal structure (subhalo distribution, shape and orientation) of the halos. The following changes are observed in the CoSANG model compared to the DMO model: 1) the total number of subhaloes close to the center of the halo is reduced, 2) the $V_{\mathrm{max}}$ distribution peaks at a lower value and lies below the DMO model, 3) the axis ratio is smaller. The difference between DMO and CoSANG simulations is more significant in early forming halos.

Volatiles, notably water, are key to the habitability of rocky planets. The presence of water in planetary material can be inferred from the atmospheric oxygen abundances of polluted white dwarfs, but this interpretation is often complex. We study the accretion process, and find that ices may sublimate and accrete before more refractory minerals reach the star. As a result, a white dwarf's relative photospheric abundances may vary with time during a single accretion event, and do not necessarily reflect the bulk composition of a pollutant. We offer two testable predictions for this hypothesis: 1. cooler stars will more often be inferred to have accreted wet pollutants, and 2. there will be rare occurrences of accretion events with inferred volatile levels far exceeding those of pristine comets. To observationally test these predictions, we statistically constrain the water content of white dwarf pollutants. We find that in the current sample, only three stars show statistically significant evidence of water at the 2$\sigma$ level, due to large typical uncertainties in atmospheric abundances and accretion states. In the future, an expanded sample of polluted white dwarfs with hydrogen-dominated atmospheres will allow for the corroboration of our theoretical predictions. Our work also shows the importance of interpreting pollutant compositions statistically, and emphasizes the requirement to reduce uncertainties on measured abundances to allow for statistically significant constraints on their water content.

Polluted white dwarfs serve as astrophysical mass spectrometers - their photospheric abundances are used to infer the composition of planetary objects that accrete onto them. We show that due to asymmetries in the accretion process, the composition of the material falling onto a star may vary with time during the accretion of a single planetary body. Consequently, the instantaneous photospheric abundances of white dwarfs do not necessarily reflect the bulk composition of their pollutants, especially when their diffusion timescales are short. In particular, we predict that when an asteroid with an iron core tidally disrupts around a white dwarf, a larger share of its mantle is ejected, and that the core/mantle fraction of the accreting material varies with time during the event. Crucially, this implies that the core fraction of differentiated pollutants cannot be determined for white dwarfs with short diffusion timescales, which sample only brief episodes of longer accretion processes. The observed population of polluted white dwarfs backs up the proposed theory. More white dwarfs have accreted material with high Fe/Ca than low Fe/Ca relative to stellar abundance ratios, indicating the ejection of mantle material. Additionally, we find tentative evidence that the accretion rate of iron decreases more rapidly than that of magnesium or calcium, hinting at variability of the accreted composition. Further corroboration of the proposed theory will come from the up-coming analysis of large samples of young white dwarfs.

Galactic outflows driven by supernovae (SNe) are thought to be a powerful regulator of a galaxy's star-forming efficiency. Mass, energy, and metal outflows ($\eta_M$, $\eta_E$, and $\eta_Z$, here normalized by the star formation rate, the SNe energy and metal production rates, respectively) shape galaxy properties by both ejecting gas and metals out of the galaxy and by heating the circumgalactic medium (CGM), preventing future accretion. Traditionally, models have assumed that galaxies self-regulate by ejecting a large fraction of the gas which enters the interstellar medium (ISM), even though such high mass-loadings are in growing tension with observations. To better understand how the relative importance of ejective (i.e. high mass-loading) vs preventative (i.e. high energy-loading) feedback affects the present-day properties of galaxies, we develop a simple gas-regulator model of galaxy evolution, where the stellar mass, ISM, and CGM are modeled as distinct reservoirs which exchange mass, metals, and energy at different rates within a growing halo. Focusing on the halo mass range from $10^{10}$ to $10^{12} M_{\odot}$, we demonstrate that, with reasonable parameter choices, we can reproduce the stellar-to-halo mass relation and the ISM-to-stellar mass relation with low mass-loaded ($\eta_M \sim 0.1-10$) but high energy-loaded ($\eta_E \sim 0.1-1$) winds, with self-regulation occurring primarily through heating and cooling of the CGM. We show that the model predictions are robust against changes to the mass-loading of outflows but are quite sensitive to our choice of the energy-loading, preferring $\eta_E \sim 1$ for the lowest mass halos and $\sim 0.1$ for Milky Way-like halos.

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion-neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at $\sim 261$ au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion-neutral decoupling condition, and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfv\'{e}nic, and its largest possible injection scale is $\sim 194$ au.

Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron star mergers, or if they appear at even higher densities that are not realized in the Cosmos. In this article, we perform numerical-relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiralling binary component. Our simulations reveal that shortly before the merger it is possible that such stars undergo a "reverse phase transition", i.e., densities decrease and the quark core inside the star disappears leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to a rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the postmerger phase transition, e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition.

In the early universe, Dirac neutrino magnetic moments due to their chirality-flipping nature could lead to thermal production of right-handed neutrinos, which would make a significant contribution to the effective neutrino number, $N_{\rm eff}$. We present in this paper a dedicated computation of the neutrino chirality-flipping rate in the thermal plasma. With a careful and consistent treatment of soft scattering and the plasmon effect in finite temperature field theories, we find that neutrino magnetic moments above $3.3\times 10^{-12}\mu_B$ have been excluded by current CMB and BBN measurements of $N_{\rm eff}$, assuming three right-handed neutrinos were thermalized. This limit is stronger than the latest bounds from XENONnT and LUX-ZEPLIN experiments and comparable with those from stellar cooling considerations.

The direct discovery of gravitational waves (GWs) from the coalescence of compact binary components by the LIGO/Virgo/KAGRA Collaboration provides an unprecedented opportunity for exploring the underlying theory of gravity that drives the coalescence process in the strong and highly-dynamical field regime of gravity. In this paper, we consider the observational effects of spatial covariant gravities on the propagation of GWs in the cosmological background and obtain the observational constraints on coupling coefficients in the action of spatial covariant gravities from GW observations. We first decompose the GWs into the left-hand and right-hand circular polarization modes and derive the effects of the spatial covariant gravities on the propagation equation of GWs. We find that these effects can be divided into three classes: (1) frequency-independent effects on GW speed and friction, (2) parity-violating amplitude and velocity birefringences, and (3) Lorentz-violating damping rate and dispersion of GWs. With these effects, we calculate the corresponding modified waveform of GWs generated by the coalescence of compact binaries. By comparing these new effects with the publicly available posterior samples or results from various tests of gravities with LIGO/Virgo/KAGRA data in the literature, we derive the observational constraints on coupling coefficients of the spatial covariant gravities. These results represent the most comprehensive constraints on the spatial covariant gravities in the literature.

We compute the Lema{\^i}tre-Hubble diagram for axial Bianchi IX universes with comoving dust. We motivate our choice by defining a {\it minimal} symmetry breaking of the cosmological principle. This criterium admits only two possibilities: the axial Bianchi I and IX universes. The latter have positive curvatures and reduce to the former in the zero curvature limit. Remarkably, negative curvatures are excluded by this minimal symmetry breaking in presence of comoving dust.

We present the novel finite-temperature FSU2H$^*$ equation-of-state model that covers a wide range of temperatures and lepton fractions for the conditions in proto-neutron stars, neutron star mergers and supernovae. The temperature effects on the thermodynamical observables and the composition of the neutron star core are stronger when the hyperonic degrees of freedom are considered. We pay a special attention to the temperature and density dependence of the thermal index in the presence of hyperons and conclude that the true thermal effects cannot be reproduced with the use of a constant $\Gamma$ law.

We present here the new results obtained with the INPOP planetary ephemerides and BepiColombo radio-science simulations. We give new constraints for the classic General Relativity tests in terms of violation of the PPN parameters $\beta$ and $\gamma$ and the time variation of the gravitational constant G. We also present new limits for the mass of the graviton and finally we obtain new acceptable intervals for the dilaton parameters $\alpha_{0}$, $\alpha_{T}$ and $\alpha_{G}$. Besides these tests of gravitation, we also study the possibility of detecting the Sun core rotation.

We demonstrate the consistency of the quark deconfinement phase transition parameters in the beta-stable neutron star matter and in the nearly symmetric nuclear matter formed in heavy-ion collisions (HICs). We investigate the proton and $\Lambda$ flow in Au+Au collisions at 3 and 4.5 GeV/nucleon incident beam energies with the pure hadron cascade version of a multi-phase transport model. The phase transition in HICs and neutron stars is described based on a class of hybrid equations of state from the quark mean-field model for the hadronic phase and a constant-speed-of-sound parametrization for the high-density quark phase. The measurements of the anisotropic proton flow at 3 GeV/nucleon by the STAR collaboration favor a relatively low phase transition density lower than $\sim 2.5$ times saturation density indicated by the gravitational wave and electromagnetic observations of neutron stars. And the proton flow data at the higher energy of 4.5 GeV/nucleon can be used to effectively constrain the softness of high-density quark matter equations of state. Finally, compared to the proton flow, the $\Lambda$ flow is found to be less sensitive and not constraining to the equations of state.

AEROS aims to develop a nanosatellite as a precursor of a future system of systems, which will include assets and capabilities of both new and existing platforms operating in the Ocean and Space, equipped with state-of-the-art sensors and technologies, all connected through a communication network linked to a data gathering, processing and dissemination system. This constellation leverages scientific and economic synergies emerging from New Space and the opportunities in prospecting, monitoring, and valuing the Ocean in a sustainable manner, addressing the demand for improved spatial, temporal, and spectral coverage in areas such as coastal ecosystems management and climate change assessment and mitigation. Currently, novel sensors and systems, including a miniaturized hyperspectral imager and a flexible software-defined communication system, are being developed and integrated into a new versatile satellite structure, supported by an innovative on-board software. Additional sensors, like the LoRaWAN protocol and a wider field of view RGB camera, are under study. To cope with data needs, a Data Analysis Centre, including a cloud-based data and telemetry dashboard and a back-end layer, to receive and process acquired and ingested data, is being implemented to provide tailored-to-use remote sensing products for a wide range of applications for private and institutional stakeholders.

In this work, we explored warm inflation in the background of $f(R,T)$ gravity in the strong dissipation regime. Considering scalar field for FLRW universe, we derived modified field equations. We then deduced slow-roll parameters under slow-roll approximations followed by power spectrum for scalar and tensor perturbations and their corresponding spectral indices. We have considered Chaotic and Natural potentials and estimated scalar spectral index and tensor-to-scalar ratio for constant as well as variable dissipation factor $\Gamma$. We found that both the rejected potentials can be revived under the context of $f(R,T)$ gravity with suitable choice of the model parameters. Further, it is seen that within the warm inflationary scenario both the potentials are consistent with Planck 2018 bounds at the Planckian and sub Planckian energy scales.

We compute the axion-pion scattering $a \pi \to \pi \pi$, relevant for the axion thermalization rate in the early universe, within unitarized NLO chiral perturbation theory. The latter extends the range of validity of the chiral expansion of axion-pion scattering and thus provides a crucial ingredient for the reliable determination of the relic density of thermal axions, whenever the axion decoupling temperature is below that of the QCD phase transition. Implications for cosmological observables are briefly discussed.

In our previous studies we have examined solar wind and magnetospheric plasmas turbulence, including Markovian character on large inertial magneto-hydrodynamic scales. Here we present the results of statistical analysis of magnetic field fluctuations in the Earth's magnetosheath based on Magnetospheric Multiscale mission at much smaller kinetic scales. Following our results on spectral analysis with very large slopes of about -16/3, we apply Markov processes approach to turbulence in this kinetic regime. It is shown that the Chapman-Kolmogorov equation is satisfied and the lowest-order Kramers-Moyal coefficients describing drift and diffusion with a power-law dependence are consistent with a generalized Ornstein-Uhlenbeck process. The solutions of the Fokker-Planck equation agree with experimental probability density functions, which exhibit a universal global scale invariance through the kinetic domain. In particular, for moderate scales we have the kappa distribution described by various peaked shapes with heavy tails, which with large values of kappa parameter are reduced to the Gaussian distribution for large inertial scales. This shows that the turbulence cascade can be described by the Markov processes also on very small scales. The obtained results on kinetic scales may be useful for better understanding of the physical mechanisms governing turbulence