Papers for Wednesday, Nov 30 2022

This chapter reviews for each province and destination of the Solar System the representative space missions that will have to be designed and implemented by 2061 to address the six key science questions about the diversity, origins, workings and habitability of planetary systems (described in chapter 1) and to perform the critical observations that have been described in chapters 3 and partly 2. It derives from this set of future representative missions, some of which will have to be flown during the 2041-2061 period, the critical technologies and supporting infrastructures that will be needed to fly these challenging missions, thus laying the foundation for the description of technologies and infrastructures for the future of planetary exploration that is given in chapters 5 and 6, respectively.

The final mass distribution of primordial black holes is sensitive to the equation of state of the Universe at the scales accessible by the power spectrum. Motivated by the presence of phase transitions in several beyond the Standard Model theories, some of which are strongly coupled, we analyse the production of primordial black holes during such phase transitions, which we model using the gauge/gravity duality. We focus in the (often regarded as physically uninteresting) case for which the phase transition is just a smooth crossover. We find an enhancement of primordial black hole production in the range $M_{\rm{PBH}}\in[10^{-16},10^{-6}]M_{\odot}$.

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations from Kov\'acs et al. (2019) used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS.

High-contrast imaging has afforded astronomers the opportunity to study light directly emitted by adolescent (tens of Myr) and ``proto" ($<$10Myr) planets still undergoing formation. Direct detection of these planets is enabled by empirical Point Spread Function (PSF) modeling and removal algorithms. The computational intensity of such algorithms, and their multiplicity of tunable input parameters, has led to the prevalence of ad-hoc optimization approaches to high-contrast imaging results. In this work, we present a new, systematic approach to optimization vetted using data of the high-contrast stellar companion HD 142527 B from the Magellan Adaptive Optics (MagAO) Giant Accreting Protoplanet Survey (GAPlanetS). More specifically, we present a grid search technique designed to explore three influential parameters of the PSF-subtraction algorithm pyKLIP -- annuli, movement, and KL modes. We consider multiple metrics for post-processed image quality in order to optimally recover at H$\alpha$ (656nm) synthetic planets injected into contemporaneous continuum (643nm) images. These metrics include: peak (single-pixel) SNR, average (multi-pixel average) SNR, 5$\sigma$ contrast, and false-positive fraction. We apply continuum-optimized KLIP reduction parameters to six H$\alpha$ direct detections of the low-mass stellar companion HD142527 B, and recover the companion at a range of separations. Relative to a single-informed, non-optimized set of KLIP parameters applied to all datasets uniformly, our multi-metric grid search optimization led to improvements in companion SNR of up to 1.2$\sigma$, with an average improvement of 0.6$\sigma$. Since many direct imaging detections lie close to the canonical 5$\sigma$ threshold, even such modest improvements may result in higher yields in future imaging surveys.

We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a $v_\mathrm{max}$-$r_\mathrm{max}$ relation. We then estimate the dark matter subhalo sensitivity in $16{,}000$ simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that, with a $3\sigma$ detection threshold, $2.35$ per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass $M_\mathrm{max}\leq 10^{10}M_\odot$, $0.03$ per cent are sensitive to $M_\mathrm{max}\leq 10^{9}M_\odot$, and, the limit of sensitivity is found to be $M_\mathrm{max}=10^{8.8\pm0.2}M_\odot$. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield $1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$ detectable subhaloes per lens in the entire sample, but this increases to $35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$ per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction $f_\mathrm{sub}^{-1}$. Assuming $f_\mathrm{sub}=0.01$, one in every $70$ lenses in general should yield a detection, or one in every $\sim$ three lenses in the most sensitive sample. From $170,000$ new strong lenses detected by Euclid, we expect $\sim 2500$ new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e., those with half-mode masses $M_\mathrm{hm}>10^8M_\odot$.

The dependence of star formation efficiency (SFE) on galactic structures, especially whether the SFE in the bar region is lower than those in the other regions, has recently been debated. We report the SFEs of 18 nearby gas-rich massive star-forming barred galaxies with a large apparent bar major axis ($\geqq 75^{\prime\prime}$). We statistically measure the SFE by distinguishing the center, bar-end, and bar regions for the first time. The molecular gas surface density is derived from archival CO(1-0) and/or CO(2-1) data by assuming a constant CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$), and the star formation rate surface density is derived from a linear combination of far-ultraviolet and mid-infrared intensities. The angular resolution is $15^{\prime\prime}$, which corresponds to $0.3 - 1.8~\rm kpc$. We find that the ratio of the SFE in the bar to that in the disk was systematically lower than unity (typically $0.6-0.8$), which means that the star formation in the bar is systematically suppressed. Our results are inconsistent with similar recent statistical studies that reported that SFE tends to be independent of galactic structures. This inconsistency can be attributed to the differences in the definition of the bar region, spatial resolution, $\alpha_{\rm CO}$, and sample galaxies. Furthermore, we find a negative correlation between SFE and velocity width of the CO spectrum, which is consistent with the idea that the large dynamical effects, such as strong shocks, large shear, and fast cloud-cloud collisions caused by the noncircular motion of the bar, result in a low SFE.

Hickson Compact Groups (HCGs) are dense configurations of 4 to 10 galaxies, whose HI (neutral gas) morphology appears to follow an evolutionary sequence of three phases, with gas initially confined to galaxies, then significant amounts spread throughout the intra-group medium, and finally with almost no gas remaining in the galaxies themselves. The HI deficiency of HCGs is expected to increase as the HI morphological phase progresses along this sequence, potentially making it a useful proxy for evolutionary phase. We test this hypothesis for the first time with a large sample of 38 HCGs with VLA HI observations that are uniformly reduced and analysed with a purpose-built pipeline. However, we find little evidence that HI deficiency can be used as a proxy for the evolutionary phase of a HCG in either of the first two phases, with the distribution of HI deficiency being consistent in both, although it does greatly increase in the third phase. This appears to be the result to three factors: a) there is already a broad range of HI deficiencies in Phase 1 HCGs, possibly due to their differing locations relative to large scale structures; b) the timescale for major interactions and morphological changes is, in general, considerably shorter than the timescale for the destruction or consumption of HI gas; and c) some groups have their HI content rejuvenated by the late addition of a new gas-rich member (for which we added a new sub-phase, 3c, to the established evolutionary sequence). Finally, across all HCGs studied, we identify only a few cases where there is strong evidence for the existence of a previously proposed diffuse HI component in the intra-group medium, which might be detectable with improved observations. This work was completed with considerable attention paid to scientific reproducibility, and all reduction and analysis has been made public via Github and Zenodo. (Abridged)

Non-equilibrium chemistry is a key process in the study of the InterStellar Medium (ISM), in particular the formation of molecular clouds and thus stars. However, computationally it is among the most difficult tasks to include in astrophysical simulations, because of the typically high (>40) number of reactions, the short evolutionary timescales (about $10^4$ times less than the ISM dynamical time) and the characteristic non-linearity and stiffness of the associated Ordinary Differential Equations system (ODEs). In this proof of concept work, we show that Physics Informed Neural Networks (PINN) are a viable alternative to traditional ODE time integrators for stiff thermo-chemical systems, i.e. up to molecular hydrogen formation (9 species and 46 reactions). Testing different chemical networks in a wide range of densities ($-2< \log n/{\rm cm}^{-3}< 3$) and temperatures ($1 < \log T/{\rm K}< 5$), we find that a basic architecture can give a comfortable convergence only for simplified chemical systems: to properly capture the sudden chemical and thermal variations a Deep Galerkin Method is needed. Once trained ($\sim 10^3$ GPUhr), the PINN well reproduces the strong non-linear nature of the solutions (errors $\lesssim 10\%$) and can give speed-ups up to a factor of $\sim 200$ with respect to traditional ODE solvers. Further, the latter have completion times that vary by about $\sim 30\%$ for different initial $n$ and $T$, while the PINN method gives negligible variations. Both the speed-up and the potential improvement in load balancing imply that PINN-powered simulations are a very palatable way to solve complex chemical calculation in astrophysical and cosmological problems.

We present a new suite of cosmological zoom-in hydrodynamical ($\approx 20\, \mathrm{pc}$) simulations of Milky-Way mass galaxies to study how a varying mass ratio for a Gaia-Sausage-Enceladus (GSE) progenitor impacts the $z=0$ chemodynamics of halo stars. Using the genetic modification approach, we create five cosmological histories for a Milky-Way-mass dark matter halo ($M_{200} \approx 10^{12} \, M_\mathrm{\odot}$), incrementally increasing the stellar mass ratio of a $z\approx2$ merger from 1:25 to 1:2, while fixing the galaxy's final dynamical, stellar mass and large-scale environment. We find markedly different morphologies at $z=0$ following this change in early history, with a growing merger resulting in increasingly compact and bulge-dominated galaxies. Despite this structural diversity, all galaxies show a radially-biased population of inner halo stars like the Milky-Way's GSE which, surprisingly, has a similar magnitude, age, $\rm [Fe/H]$ and $\rm [\alpha/Fe]$ distribution whether the $z\approx2$ merger is more minor or major. This arises because a smaller ex-situ population at $z\approx2$ is compensated by a larger population formed in an earlier merger-driven starburst, with both populations strongly overlapping in the $\rm [Fe/H]-\rm [\alpha/Fe]$ plane. Our study demonstrates that multiple high-redshift histories can lead to similar $z=0$ chemodynamical features in the halo, highlighting the need for additional constraints to distinguish them, and the importance of considering the full spectrum of progenitors when interpreting $z=0$ data to reconstruct our Galaxy's past.

Interactions between photons and electrons are ubiquitous in astrophysics. Photons can be down scattered (Compton scattering) or up scattered (inverse Compton scattering) by moving electrons. Inverse Compton scattering, in particular, is an essential process for the production of astrophysical gamma rays. Computations of inverse Compton emission typically adopts an isotropic or an ultrarelativistic assumption to simplify the calculation, which makes them unable to broadcast the formula to the whole phase space of source particles. In view of this, we develop a numerical scheme to compute the interactions between anisotropic photons and electrons without taking ultrarelativistic approximations. Compared to the ultrarelativistic limit, our exact results show major deviations when target photons are down scattered or when they possess energy comparable to source electrons. We also consider two test cases of high-energy inverse Compton emission to validate our results in the ultrarelativistic limit. In general, our formalism can be applied to cases of anisotropic electron-photon scattering in various energy regimes, and for computing the polarizations of the scattered photons.

The formation histories of compact binary mergers, especially stellar-mass binary-black hole mergers, have recently come under increased scrutiny and revision. In this paper we revisit the question of the dominant formation channel and efficiency of forming binary neutron-star mergers. We use the stellar and binary evolution code MESA and implement an up-to-date and detailed method for common envelope and mass transfer. We preform simulations for donor masses between 8-20 solar masses with a neutron star companion of 1.4 and 2.0 solar masses, at two metallicities, using varying common envelope efficiencies, and two prescriptions for electron-capture supernovae. In contrast to the case of binary-black hole mergers, for a neutron star companion of 1.4 solar masses, all binary neutron star mergers are formed following a common envelope phase, while for a neutron star mass of 2.0 solar masses we identify a small subset of mergers following only stable mass transfer if the neutron star receives a large natal kick. Regardless of neutron star companion mass, we find that large supernova natal kicks are favored in the formation of binary neutron star mergers, and find more mergers at subsolar metallicity compared to solar.

Close-in planets orbiting around low-mass stars are exposed to intense energetic photon and particle radiation and harsh space weather. We have modeled such conditions for Proxima Centauri b (Garraffo et al. 2016b), a rocky planet orbiting in the habitable-zone of our closest neighboring star, finding a stellar wind pressure three orders of magnitude higher than the solar wind pressure on Earth. At that time, no Zeeman-Doppler observations of the surface magnetic field distribution of Proxima Cen were available and a proxy from a star with similar Rossby number to Proxima was used to drive the MHD model. Recently, the first ZDI observation of Proxima Cen became available (Klein et al. 2021). We have modeled Proxima b's space weather using this map and compared it with the results from the proxy magnetogram. We also computed models for a high-resolution synthetic magnetogram for Proxima b generated by a state-of-the-art dynamo model. The resulting space weather conditions for these three scenarios are similar with only small differences found between the models based on the ZDI observed magnetogram and the proxy. We conclude that our proxy magnetogram prescription based on Rossby number is valid, and provides a simple way to estimate stellar magnetic flux distributions when no direct observations are available. Comparison with models based on the synthetic magnetogram show that the exact magnetogram details are not important for predicting global space weather conditions of planets, reinforcing earlier conclusions that the large-scale (low-order) field dominates, and that the small-scale field does not have much influence on the ambient stellar wind.

The earliest stages of star formation, when young stars are still deeply embedded in their natal clouds, represent a critical phase in the matter cycle between gas clouds and young stellar regions. Until now, the high-resolution infrared observations required for characterizing this heavily obscured phase (during which massive stars have formed, but optical emission is not detected) could only be obtained for a handful of the most nearby galaxies. One of the main hurdles has been the limited angular resolution of the Spitzer Space Telescope. With the revolutionary capabilities of the JWST, it is now possible to investigate the matter cycle during the earliest phases of star formation as a function of the galactic environment. In this Letter, we demonstrate this by measuring the duration of the embedded phase of star formation and the implied time over which molecular clouds remain inert in the galaxy NGC628 at a distance of 9.8Mpc, demonstrating that the cosmic volume where this measurement can be made has increased by a factor of $>100$ compared to Spitzer. We show that young massive stars remain embedded for $5.1_{-1.4}^{+2.7}$Myr ($2.3_{-1.4}^{+2.7}$Myr of which being heavily obscured), representing $\sim20\%$ of the total cloud lifetime. These values are in broad agreement with previous measurements in five nearby ($D < 3.5$Mpc) galaxies and constitute a proof of concept for the systematic characterization of the early phase of star formation across the nearby galaxy population with the PHANGS-JWST survey.

There are excesses of sub-Neptunes just wide of period commensurabilities like the 3:2 and 2:1, and corresponding deficits narrow of them. Any theory that explains this period ratio structure must also explain the strong transit timing variations (TTVs) observed near resonance. Besides an amplitude and a period, a sinusoidal TTV has a phase. Often overlooked, TTV phases are effectively integration constants, encoding information about initial conditions or the environment. Many TTVs near resonance exhibit non-zero phases. This observation is surprising because dissipative processes that capture planets into resonance also damp TTV phases to zero. We show how both the period ratio structure and the non-zero TTV phases can be reproduced if pairs of sub-Neptunes capture into resonance in a gas disc while accompanied by a third non-resonant body. Convergent migration and eccentricity damping drives pairs to orbital period ratios wide of commensurability and, after the disk clears, secular forcing by the third-body perturber phase-shifts the TTVs. The scenario predicts that resonant planets are apsidally aligned and possess eccentricities up to an order of magnitude larger than previously thought.

We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, 56Ni mass synthesized, host extinction, and the time of explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (t less than 25 days, t less than 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of explosion, steepness parameter beta, and host extinction are reasonably constrained with incomplete light curve data, whereas mass-loss rate, kinetic energy and 56Ni mass estimates generally require complete light curves spanning greater than 100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow up observations. Our findings are relevant to identify, prioritize, and coordinate efficient follow up of transients discovered by Vera C. Rubin Observatory.

We present new GRIFFIN project hydrodynamical simulations that model the formation of galactic star cluster populations in low-metallicity ($Z=0.00021$) dwarf galaxies, including radiation, supernova and stellar wind feedback of individual massive stars. In the simulations, stars are sampled from the stellar initial mass function (IMF) down to the hydrogen burning limit of 0.08 M$_\odot$. Mass conservation is enforced within a radius of 1 pc for the formation massive stars. We find that massive stars are preferentially found in star clusters and follow a correlation set at birth between the highest initial stellar mass and the star cluster mass, in agreement with observations. With a fully sampled IMF, star clusters loose mass in the galactic tidal field according to mass-loss rates observed in nearby galaxies. Of the released stellar feedback, 60\% of the supernova material and up to 35\% of the wind material reside either in the hot interstellar medium (ISM) or in gaseous, metal enriched outflows. While stellar winds (instantaneously) and supernovae (delayed) start enriching the ISM right after the first massive stars form, the formation of supernova-enriched stars is significantly delayed (by $>50$ Myr) compared to the formation of stars enriched by stellar winds. The first forming star clusters are therefore solely enriched by stellar wind material. Overall, supernova ejecta dominate the enrichment by mass, while the number of enriched stars is determined by continuous stellar winds. These results present a concept for the formation of chemically distinct populations of stars in bound star clusters, reminiscent of multiple populations in globular clusters.

AT 2000ch is a highly variable massive star and supernova imposter in NGC 3432 first detected in 2000. It is similar and often compared to SN 2009ip, and it is therefore expected to undergo a core-collapse supernova -- a SN imposter of similar brightness -- in the near future. We characterize the long-term variability of AT 2000ch in the radio and optical regimes with archival data reaching back to the year 1984. We use these newly reduced observations in addition to observations in the literature to restrict the mass-loss rates of AT 2000ch at multiple epochs based on different approaches, and to infer the general properties of its circumstellar nebula with respect to the detected radio brightness. We extend the known optical light curve of AT 2000ch up to the beginning of 2022 by performing point spread function photometry on archival data from the Palomar Transient Factory and the Zwicky Transient Facility. We reduced archival radio continuum observations obtained with the Very Large Array using standard calibration and imaging methods and complemented these with pre-reduced \changes observations as well as observations obtained with the Westerbork Synthesis Radio Telescope and LOw Frequency ARray. For the analysis of AT 2000ch, we consider the optical light curve and color evolution, its radio continuum brightness at different frequencies and times, and the corresponding spectral indices. We estimated mass-loss rates and optical depths based on radio continuum brightnesses and Ha fluxes. We report two newly detected outbursts of AT 2000ch similar to those found in the 2000s and 13 re-brightening events, of which at least four are not conclusively detected because of insufficient sampling of the light curve. The dates of all outbursts and significant, well-sampled re-brightening events are consistent with a period of $\sim 201 \pm 12\,$days over a total time-span of two decades. Such a behavior has never been found for any SN imposter, especially not for candidate SN~2009ip analogs. During 2010 to 2012 and 2014 to 2018, we only have a few detections, which is insufficient to come to any conclusion as to a possible less eruptive phase of the transient. We find steady dimming after the most recent re-brightening events and possible evidence of porosity in the circumstellar envelope, suggesting AT~2000ch may currently be in transition to a state of relative calm. We identified a second, unrelated source at a projected distance of $\sim 23\,$pc ($\sim0.5^{\prime\prime}$) that has contaminated the optical measurements of AT~2000ch at its minimum luminosity over the last two decades probably on a $5\%-10\,\%$ level, but this does not affect our overall findings and is negligible during re-brightening. We are able to restrict the mass-loss rate of AT~2000ch to range between several $10^{-6}\,\textrm{M}_{\odot}/\textrm{yr}$ and several $10^{-5}\,\textrm{M}_{\odot}/\textrm{yr}$. The fresh ejecta appear to be optically thick to radio continuum emission at least within the first $\sim 25\,$days after significant re-brightening.

We present the first high spectral resolution mid-infrared survey in the Orion BN/KL region, covering 7.2 to 28.3 micron. With SOFIA/EXES we target the enigmatic source Orion IRc2. While this is in the most prolifically studied massive star-forming region, longer wavelengths and molecular emission lines dominated previous spectral surveys. The mid-infrared observations in this work access different components and molecular species in unprecedented detail. We unambiguously identify two new kinematic components, both chemically rich with multiple molecular absorption lines. The "blue clump" has vLSR = -7.1 \pm 0.7 km/s and the "red clump" 1.4 \pm 0.5 km/s. While the blue and red clumps have similar temperatures and line widths, molecular species in the blue clump have higher column densities. They are both likely linked to pure rotational H2 emission also covered by this survey. This work provides evidence for the scenario that the blue and red clumps are distinct components unrelated to the classic components in the Orion BN/KL region. Comparison to spectroscopic surveys towards other infrared targets in the region show that the blue clump is clearly extended. We analyze, compare, and present in depth findings on the physical conditions of C2H2, 13CCH2, CH4, CS, H2O, HCN, H13CN, HNC, NH3, and SO2 absorption lines and an H2 emission line associated with the blue and red clumps. We also provide limited analysis of H2O and SiO molecular emission lines towards Orion IRc2 and the atomic forbidden transitions [FeII], [SI], [SIII], and [NeII].

The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the respective electron and positron spectra. For the latter, those would coincide in energy and would contrast predictions of smooth spectra as from particle dark matter. In this work, using a library of pulsar population models for the local part of the Milky Way, we perform a power-spectrum analysis on the cosmic-ray positron fraction. We also develop a technique to cross-correlate the electron and positron fluxes. We show that both such analyses, can be used to search statistically for the presence of spectral wiggles in the cosmic-ray data. For a significant fraction of our pulsar simulations, those techniques are already sensitive enough to give a signal for the presence of those features above the regular noise, with forthcoming observations making them even more sensitive. Finally, by cross-correlating the AMS-02 electron and positron spectra, we find an intriguing first hint for a positive correlation between them, of the kind expected by a population of local pulsars.

We constructed an oxygen abundance map and N/O ratio map of the unusually low excitation dwarf irregular galaxy NGC 4163 based on publicly available spectroscopy obtained by the MaNGA survey. We detected auroral emission line [OII]$\lambda\lambda$7320,7330 which allows us to measure chemical abundance by direct T$_e$ method. We found that the scatter of the oxygen abundance derived by the strong line method is large. The oxygen abundances 12 + log(O/H) derived by strong line method vary from ~7.3 to ~7.8 with a mean value of ~7.55. The oxygen abundances derived in two apertures of 2 arcseconds by the direct T$_e$ method using our measurements of the O$^+$ auroral line is about 7.8 dex. The nitrogen-to-oxygen ratio log(N/O) of about -1.5 is typical value for a low metallicity galaxy, maybe slightly shifted towards higher N/O ratios with respect to the N/O values in the HII regions in nearby galaxies. An unusual negative trend between log(N/O) and oxygen abundance is detected. NGC 4163 is a gas-poor galaxy with a neutral atomic gas mass fraction of around 0.25. The oxygen abundance in the galaxy is only around 0.1 of the oxygen abundance potentially attainable in a galaxy with such a gas mass fraction. The low metallicity coupled with the low gas mass fraction implies that either the metallicity of the interstellar medium of the galaxy was reduced by pristine gas infall in the recent epoch or the evolution of this galaxy was accompanied by strong galactic winds.

Neutron stars are rich laboratories of multiple branches of modern physics. These include gravitational physics, nuclear and particle physics, (quantum) electrodynamics, and plasma astrophysics. In this chapter, we present the pioneering theoretical studies and the pivotal historical observations on which our understanding of neutron stars is based on. Then, we discuss the usage of neutron stars as probes of fundamental theories of physics.

We report the statistical properties of stars and brown dwarfs obtained from three radiation hydrodynamical simulations of star cluster formation with metallicities of 1, 1/10 and 1/100 of the solar value. The star-forming clouds are subjected to cosmic microwave background radiation that is appropriate for star formation at a redshift z=5. The results from the three calculations are compared to each other, and to similar previously published calculations that had levels of background radiation appropriate for present-day (z=0) star formation. Each of the calculations treat dust and gas temperatures separately and include a thermochemical model of the diffuse interstellar medium. We find that whereas the stellar mass distribution is insensitive to the metallicity for present-day star formation, at z=5 the characteristic stellar mass increases with increasing metallicity and the mass distribution has a deficit of brown dwarfs and low-mass stars at solar metallicity compared to the Galactic initial mass function. We also find that the multiplicity of M-dwarfs decreases with increasing metallicity at z=5. These effects are a result of metal-rich gas being unable to cool to as low temperatures at z=5 compared to at z=0 due to the hotter cosmic microwave background radiation, which inhibits fragmentation at high densities.

Some Hydrogen-poor superluminous supernovae are likely powered by a magnetar central engine, making their luminosity larger than common supernovae. Although a significant amount of X-ray flux is expected from the spin down of the magnetar, direct observational evidence is still to be found, giving rise to the "missing energy" problem. Here we present NuSTAR observations of nearby SN 2018hti 2.4y (rest frame) after its optical peak. We expect that, by this time, the ejecta have become optically thin for photons more energetic than about 15keV. No flux is detected at the position of the supernova down to $F_{\rm{10-30keV}} = 9.0\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, or an upper limit of $7.9 \times 10^{41}$ erg s$^{-1}$ at a distance of 271Mpc. This constrains the fraction of bolometric luminosity from the putative spinning down magnetar to be $f_{\rm X} \lesssim 36$% in the 10-30keV range in a conservative case, $f_{\rm X} \lesssim 11$% in an optimistic case.

The 3UCubed project is a 3U CubeSat being jointly developed by the University of New Hampshire, Sonoma State University, and Howard University as a part of the NASA Interstellar Mapping and Acceleration Probe, IMAP, student collaboration. This project comprises of a multidisciplinary team of undergraduate students from all three universities. The mission goal of the 3UCubed is to understand how Earths polar upper atmosphere the thermosphere in Earths auroral regions, responds to particle precipitation and solar wind forcing, and internal magnetospheric processes. 3UCubed includes two instruments with rocket heritage to achieve the science mission: an ultraviolet photomultiplier tube, UVPMT, and an electron retarding potential analyzer ERPA. The spacecraft bus consists of the following subsystems: Attitude Determination and Control, Command and Data Handling, Power, Communication, Structural, and Thermal. Currently, the project is in the post-PDR stage, starting to build and test engineering models to develop a FlatSat prior to critical design review in 2023. The goal is to launch at least one 3U CubeSat to collect science data close to the anticipated peak of Solar Cycle 25 around July 2025. Our mother mission, IMAP, is also projected to launch in 2025, which will let us jointly analyze the science data of the main mission, providing the solar wind measurements and inputs to the magnetosphere with that of 3UCubed, providing the response of Earths cusp to these inputs.

Aims. Lightning mapper sensors on board weather satellites can be successfully used to observe fireballs. These sensors use a very narrow spectral band at 777nm, which is only a small fraction of the total fireball radiation. In this spectral band, the oxygen O I-1 triplet is dominant for fast meteors and the Planck continuum can prevail in slow meteors. It is possible to estimate the meteor brightness in the visible spectral range from this narrowband radiation, but it is vital to first study the dependence of this radiation on the meteor velocity. Methods. We used observations from the well-established European Fireball Network with newly developed digital spectral cameras that allowed us to study the oxygen triplet in meteor spectra and its relation to the meteor velocity and altitude. In addition, we studied strong magnesium and sodium lines. Results. We developed a method for calibration of fireball observation reported by Geostationary Lightning Mapper (GLM) sensors on board the Geostationary Operational Environmental Satellite (GOES) weather satellites. We confirm that in slow meteors, the radiation of the Planck continuum dominates, but for faster meteors, a correction on velocity is needed. We observe that the altitude where the oxygen line was recorded can also affect the radiation at 777 nm. In addition, determining whether or not the meteor showed a bright flare could also lead to a similar effect. Thus, the meteor brightness estimate may be impacted by these characteristics. We derived simple corrections on the altitude and on the meteor brightness that helped to improve the overall precision of the magnitude estimate of our sample. This allowed us to estimate the magnitude of meteors observed by GLM with an accuracy of ~ 1 in magnitude. The Na/Mg line intensity ratio was found to be constant for velocities above 25 km/s and increasing toward lower velocities.

Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy. Recent papers suggest that extreme rotation is required to retrieve the photometric split. Our new grids of MESA models and the prevalent SYCLIST models show, however, that initial slow (0-35% of the linear Keplerian rotation velocities) and intermediate (50-65% of the Keplerian rotation velocities) rotation are adequate to explain the photometric split. These values are consistent with the recent spectroscopic measurements of cluster and field stars, and are likely to reflect the birth spin distributions of upper main-sequence stars. A fraction of the initially faster-rotating stars may be able to reach near-critical rotation at the end of their main-sequence evolution and produce Be stars in the turn-off region of young star clusters. However, we find that the presence of Be stars up to two magnitudes below the cluster turnoff advocates for a crucial role of binary interaction in creating Be stars. We argue that surface chemical composition measurements may help distinguish these two Be star formation channels. While only the most rapidly rotating, and therefore nitrogen-enriched, single stars can evolve into Be stars, slow pre-mass-transfer rotation and inefficient accretion allows for mild or no enrichment even in critically rotating accretion-induced Be stars. Our results shed new light on the origin of the spin distribution of young and evolved B-type main sequence stars.

Massive stars are the progenitors of black holes and neutron stars, the mergers of which can be detected with gravitational waves (GW). The expansion of massive stars is one of the key factors affecting their evolution in close binary systems, but it remains subject to large uncertainties in stellar astrophysics. For population studies and predictions of GW sources, the stellar expansion is often simulated with the analytic formulae from Hurley et al. (2000). These formulae need to be extrapolated for stars beyond 50 solar masses and are often considered outdated. In this work we present five different prescriptions developed from 1D stellar models to constrain the maximum expansion of massive stars. We adopt these prescriptions to investigate how stellar expansion affects mass transfer interactions and in turn the formation of GW sources. We show that limiting radial expansion with updated 1D stellar models, when compared to the use of Hurley et al. (2000) radial expansion formulae, does not significantly affect GW source properties (rates and masses). This is because most mass transfer events leading to GW sources are initialised before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (total mass > 50 solar masses) formed from stars that may evolve beyond the Humphreys-Davidson limit, whose radial expansion is the most uncertain. We conclude that understanding the expansion of massive stars and the origin of the Humphrey-Davidson limit is a key factor for the study of GW sources.

Generative adversarial networks are a promising tool for image generation in the astronomy domain. Of particular interest are conditional generative adversarial networks (cGANs), which allow you to divide images into several classes according to the value of some property of the image, and then specify the required class when generating new images. In the case of images from Imaging Atmospheric Cherenkov Telescopes (IACTs), an important property is the total brightness of all image pixels (image size), which is in direct correlation with the energy of primary particles. We used a cGAN technique to generate images similar to whose obtained in the TAIGA-IACT experiment. As a training set, we used a set of two-dimensional images generated using the TAIGA Monte Carlo simulation software. We artificiallly divided the training set into 10 classes, sorting images by size and defining the boundaries of the classes so that the same number of images fall into each class. These classes were used while training our network. The paper shows that for each class, the size distribution of the generated images is close to normal with the mean value located approximately in the middle of the corresponding class. We also show that for the generated images, the total image size distribution obtained by summing the distributions over all classes is close to the original distribution of the training set. The results obtained will be useful for more accurate generation of realistic synthetic images similar to the ones taken by IACTs.

Intervening metal absorption lines in the spectra of z > 6 quasars are fundamental probes of the ionization state and chemical composition of circumgalactic and intergalactic gas near the end of the reionization epoch. Large absorber samples are required to robustly measure typical absorber properties and to refine models of the synthesis, transport, and ionization of metals in the early Universe. The "Ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" (XQR-30) has obtained high signal-to-noise spectra of 30 luminous quasars, nearly quadrupling the existing sample of 12 high quality z~6 quasar spectra. We use this unprecedented sample to construct a catalog of 778 systems showing absorption in one or more of MgII (360 systems), FeII (184), CII (46), CIV (479), SiIV (127), and NV (13) which span 2 < z < 6.5. This catalog significantly expands on existing samples of z > 5 absorbers, especially for CIV and SiIV which are important probes of the ionizing photon background at high redshift. The sample is 50% (90%) complete for rest-frame equivalent widths W > 0.03AA (0.09AA). We publicly release the absorber catalog along with completeness statistics and a Python script to compute the absorption search path for different ions and redshift ranges. This dataset is a key legacy resource for studies of enriched gas from the era of galaxy assembly to cosmic noon, and paves the way for even higher redshift studies with the James Webb Space Telescope and 30m-class telescopes.

We propose a semi-analytic model that is developed to understand the cosmological evolution of the mean metallicity in the Universe. In particular, we study the contributions of Population III (Pop III) and Population II (Pop II) stars to the production of $\mathrm{Fe,~Si,~Zn, ~Ni,~P, ~Mg, ~Al, ~S, ~C, ~N}$, and $\mathrm{~O}$. We aim to quantify the roles of two different models in the chemical enrichment of the Universe. The first model (A) considers both stars with Pop III and Pop II yields. For the second model (B), the yields involved are only for Pop II stars. We start by describing the cosmic star formation rate (CSFR) through an adaptation of a scenario developed within the hierarchical scenario of structure formation with a Press-Schechter-like formalism. We adapt the formalism to implement the CSFR to the standard chemical evolution scenario to investigate the course of chemical evolution on a cosmological basis. Calculations start at redshift $z\sim 20$, and we compare the results of our two models with data from damped Lyman-$\alpha$ systems (DLAs), and globular clusters (GCs). Our main results find that metal production in the Universe occurred very early, quickly increasing with the formation of the first stars. When comparing results for [Fe/H] with observations from GCs, yields of Pop II stars are not enough to explain the observed chemical abundances, requiring stars with physical properties similar those expected from Pop III stars. Our semi-analytic model can deliver consistent results for the evolution of cosmic metallicities. Our results show that the chemical enrichment in the early Universe is rapid, and at redshift $\sim 12.5$, the metallicity reaches $10^{-4}\, Z_{\odot}$ for the model that includes Pop III stars. In addition, we explore values for the initial mass function (IMF) within the range $[0.85, 1.85]$.

We calculate the evolution of a star-forming cloud core using a three-dimensional resistive magnetohydrodynamics simulation, treating dust grains as Lagrangian particles, to investigate the dust motion in the early star formation stage. We prepare six different-sized set of dust particles in the range $a_{\rm d}=0.01$--$1000\,\mu$m, where $a_{\rm d}$ is the dust grain size. In a gravitationally collapsing cloud, a circumstellar disk forms around a protostar and drives a protostellar outflow. Almost all the small dust grains ($a_{\rm d} \lesssim 10$--$100\,\mu$m) initially distributed in the region $\theta_0 \lesssim 45^\circ$ are ejected from the center by the outflow, where $\theta_0$ is the initial zenith angle relative to the rotation axis, whereas only a small number of the large dust grains ($a_{\rm d} \gtrsim 100\,\mu$m) distributed in the region are ejected. All other grains fall onto either the protostar or disk without being ejected by the outflow. Regardless of the dust grain size, the behavior of the dust motion is divided into two trends after dust particles settle into the circumstellar disk. The dust grains reaching the inner disk region from the upper envelope preferentially fall onto the protostar, while those reaching the outer disk region or disk outer edge from the envelope can survive without an inward radial drift. These surviving grains can induce dust growth. Thus, we expect that the outer disk regions could be a favored place of planet formation.

Excess noise from scattered light poses a persistent challenge in the analysis of data from gravitational wave detectors such as LIGO. We integrate a physically motivated model for the behavior of these "glitches" into a standard Bayesian analysis pipeline used in gravitational wave science. This allows for the inference of the free parameters in this model, and subtraction of these models to produce glitch-free versions of the data. We show that this inference is an effective discriminator of the presence of the features of these glitches, even when those features may not be discernible in standard visualizations of the data.

Using various latest cosmological datasets including Type-Ia supernovae, cosmic microwave background radiation, baryon acoustic oscillations, and estimations of the Hubble parameter, we test some dark energy models with parameterized equations of state and try to distinguish or select observation-preferred models. We obtain the best fitting results of the six models and calculate their values of the Akaike Information Criteria and Bayes Information Criterion. And we can distinguish these dark energy models from each other by using these two information criterions. However, the $\Lambda $CDM model remains the best fit model. Furthermore, we perform geometric diagnostics including statefinder and Om diagnostics to understand the geometric behaviour of the dark energy models. We find that the six DE models can be distinguished from each other and from $\Lambda$CDM, Chaplygin gas, quintessence models after the statefinder and Om diagnostics were performed. Finally, we consider the growth factor of the dark energy models with comparison to $\Lambda $CDM model. Still, we find the models can be distinguished from each other and from $\Lambda $CDM model through the growth factor approximation.

Delta Scuti ($\delta$ Sct) variables are intermediate mass stars that lie at the intersection of the main sequence and the instability strip on the Hertzsprung-Russel diagram. Various lines of evidence indicate that nonlinear mode interactions shape their oscillation spectra, including the particularly compelling detection of resonantly interacting mode triplets in the $\delta$ Sct star KIC 8054146. Motivated by these observations, we use the theory of three-mode coupling to study the strength and prevalence of nonlinear mode interactions in fourteen $\delta$ Sct models that span the instability strip. For each model, we calculate the frequency detunings and nonlinear coupling strengths of $\sim 10^4$ unique combinations of mode triplets. We find that all the models contain at least $\sim 100$ well-coupled triplets whose detunings and coupling strengths are consistent with the triplets identified in KIC 8054146. Our results suggest that resonant mode interactions can be significant in $\delta$ Sct stars and may explain why many exhibit rapid changes in amplitude and oscillation period.

Tidal disruption events (TDEs) around super massive black holes (SMBHs) are a potential laboratory to study super-Eddington accretion disks and sometimes result in powerful jets or outflows which may shine in the radio and sub millimeter bands. In this work, we model the thermal synchrotron emission of jets from general relativistic radiation magneto-hydrodynamics (GRRMHD) simulations of a BH accretion disk/jet system which assumes the TDE resulted in a magnetized accretion disk around a BH accreting at $\sim 12-25$ times the Eddington accretion rate. Through synthetic observations with the Next Generation Event Horizon Telescope (ngEHT) and an image reconstruction analysis, we demonstrate that TDE jets may provide compelling targets, within the context of the models explored in this work. In particular, we find that jets launched by a SANE super-Eddington disk around a spin $a_*=0.9$ reach the ngEHT detection threshold at large distances (up to 100 Mpc in this work). A two-temperature plasma in the jet or weaker jets, such as a spin $a_*=0$ model, requires a much closer distance as we demonstrate detection at 10 Mpc for limiting cases of $a_*=0,\,\mathcal{R}=1$ or $a_*=0.9,\, \mathcal{R}=20$. We also demonstrate that TDE jets may appear as superluminal sources if the BH is rapidly rotating and the jet is viewed nearly face on.

We investigate a novel satellite avoidance strategy to mitigate the impact of large commercial satellite constellations in low-Earth orbit on the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We simulate the orbits of currently planned Starlink and OneWeb constellations ($\sim$40,000 satellites) to test how effectively an upgraded Rubin scheduler algorithm can avoid them, and assess how the overall survey is affected. Given a reasonably accurate satellite orbit forecast, we find it is possible to adjust the scheduler algorithm to effectively avoid some satellites. Overall, sacrificing 10% of LSST observing time to avoid satellites reduces the fraction of LSST visits with streaks by a factor of two. Whether such a mitigation will be required depends on the overall impact of streaks on science, which is not yet well quantified. This is due to a lack of adequate information about satellite brightness distributions as well as the impact of glints and low surface brightness residuals on alert purity and systematic errors in cosmological parameter estimation. A significant increase in the number of satellites or their brightness during Rubin Operations may make implementing this satellite avoidance strategy worthwhile.

It is well established that solar eruptions are powered by free magnetic energy stored in current-carrying magnetic field in the corona. It has also been generally accepted that magnetic flux ropes (MFRs) are a critical component of many coronal mass ejections (CMEs). What remains controversial is whether MFRs are present well before the eruption. Our aim is to identify progenitors of MFRs, and investigate pre-eruptive magnetic properties associated with these progenitors. Here we analyze 28 MFRs erupting within 45 deg from the disk center from 2010 to 2015. All MFRs'feet are well identified by conjugate coronal dimmings. We then calculate magnetic properties at the feet of the MFRs, prior to their eruptions, using Helioseismic and Magnetic Imager (HMI) vector magnetograms. Our results show that only 8 erupting MFRs are associated with significant non-neutralized electric currents, 4 of which also exhibit pre-eruptive dimmings at the foot-prints. Twist and current distributions are asymmetric at the two feet of these MFRs. The presence of pre-eruption dimmings associated with non-neutralized currents suggests the pre-existing MFRs. Furthermore, evolution of conjugate dimmings and electric currents within the foot-prints can provide clues about the internal structure of MFRs and their formation mechanism.

We present new mass estimates for Andromeda (M31) using the orbital angular momenta of four satellite galaxies (M33, NGC 185, NGC 147, IC 10) derived from existing proper motions, distances, and line-of-sight velocities. We infer two masses for M31: $M_{\rm vir}= 2.81^{+1.48}_{-0.75}\times10^{12}\, M_{\odot}$ using satellite galaxy phase space information derived with HST-based M31 proper motions and $M_{\rm vir}=2.99^{+1.32}_{-0.68}\times10^{12}\, M_{\odot}$ using phase space information derived with the weighted average of HST+Gaia-based M31 proper motions. The precision of our new M31 mass estimates (23-53%) improves by a factor of two compared to previous mass estimates using a similar methodology with just one satellite galaxy and places our results amongst the highest precision M31 estimates in recent literature. Furthermore, our results are consistent with recently revised estimates for the total mass of the Local Group (LG), with the stellar mass--halo mass relation, and with observed kinematic data for both M31 and its entire population of satellites. An M31 mass $> 2.5 \times 10^{12}\, M_{\odot}$ could have major implications for our understanding of LG dynamics, M31's merger and accretion history, and our understanding of LG galaxies in a cosmological context.

In this paper, we perform a comprehensive study for the physical properties of SN 2018gk which is a luminous type IIb supernova (SN). We find that the early-time photospheric velocity vary from a larger value to a smaller value before the photosphere reach a temperature floor. We generalize the photosphere modulus and fit the multiband light curves (LCs) of SN 2018gk. We find that the $^{56}$Ni mass model require $\sim$1.1 M$_\odot$ of $^{56}$Ni which is larger than the derived ejecta mass ($\sim$0.58 M$_\odot$) as well as the derived $^{56}$Ni mass ($\sim$0.4 M$_\odot$) in the literature. Alternatively, we use the magnetar plus $^{56}$Ni and the fallback plus $^{56}$Ni models to fit the LCs of SN 2018gk, finding that the two models can can fit the LCs. We favor the magnetar plus $^{56}$Ni since the parameters are well constrained and rather reasonable ($M_{\rm ej} =$ 2.80 M$_\odot$, $M_{\rm Ni}=0.16 $ M$_\odot$ which is smaller than the upper limit of the value of the $^{56}$Ni mass can by synthesized by the neutrino-powered core collapse SNe, $B=8.47\times10^{14}$ G which is comparable to those of luminous and superluminous SNe studied in the literature, and $P_0=$ 10.83 ms which is comparable to those of luminous SNe). Therefore, we suggest that SN 2018gk might be a SNe IIb mainly powered by a central engine. Finally, we confirm the NIR excesses of the SEDs of SN 2018gk at some epochs and constrain the physical properties of the putative dust using the blackbody plus dust emission model.

The star formation rate in galaxies is well known to correlate with stellar mass (the `star-forming main sequence'). Here we extend this further to explore any additional dependence on galaxy surface brightness, a proxy for stellar mass surface density. We use a large sample of low redshift ($z \leq 0.08$) galaxies from the GAMA survey which have both SED derived star formation rates and photometric bulge-disc decompositions, the latter providing measures of disc surface brightness and disc masses. Using two samples, one of galaxies fitted by a single component with S\'{e}rsic index below 2 and one of the discs from two-component fits, we find that once the overall mass dependence of star formation rate is accounted for, there is no evidence in either sample for a further dependence on stellar surface density.

The Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope system is located on the Canary Island of La Palma and inspects the very high-energy (VHE, few tens of GeV and above) gamma-ray sky. MAGIC consists of two imaging atmospheric Cherenkov telescopes (IACTs), which capture images of the air showers originating from the absorption of gamma rays and cosmic rays by the atmosphere, through the detection of Cherenkov photons emitted in the shower. The sensitivity of IACTs to gamma-ray sources is mainly determined by the ability to reconstruct the properties (type, energy, and arrival direction) of the primary particle generating the air shower. The state-of-the-art IACT pipeline for shower reconstruction is based on the parameterization of the shower images by extracting geometric and stereoscopic features and machine learning algorithms like random forest or boosted decision trees. In this contribution, we explore deep convolutional neural networks applied directly to the pixelized images of the camera as a promising method for IACT full-event reconstruction and present the performance of the method on observational data using CTLearn, a package for IACT event reconstruction that exploits deep learning.

It has been suggested that the occurrence rate of hot Jupiters (HJs) in open clusters might reach several per cent, significantly higher than that of the field ($\sim$ a per cent). In a stellar cluster, when a planetary system scatters with a stellar binary, it may acquire a companion star which may excite large amplitude von Zeipel-Lidov-Kozai oscillations in the planet's orbital eccentricity, triggering high-eccentricity migration and the formation of an HJ. We quantify the efficiency of this mechanism by modelling the evolution of a gas giant around a solar mass star under the influence of successive scatterings with binary and single stars. We show that the chance that a planet $\in(1,10)$ au becomes an HJ in a Gyr in a cluster of stellar density $n_*=50$ pc$^{-3}$ and binary fraction $f_\mathrm{bin}=0.5$ is about 2\% and an additional 4\% are forced by the companion star into collision with or tidal disruption by the central host. An empirical fit shows that the total percentage of those outcomes asymptotically reaches an upper limit determined solely by $f_\mathrm{bin}$ (e.g., $10\%$ at $f_\mathrm{bin}=0.3$ and 18\% at $f_\mathrm{bin}=1$) on a timescale inversely proportional to $n_*$ ($\sim$ Gyr for $n_*\sim100$ pc$^{-3}$). The ratio of collisions to tidal disruptions is roughly a few, and depends on the tidal model. Therefore, if the giant planet occurrence rate is 10~\%, our mechanism implies an HJ occurrence rate of a few times 0.1~\% in a Gyr and can thus explain a substantial fraction of the observed rate.

A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air-showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 1 January 2004 and 31 December 2020 are given for cosmic rays that have energies in the range 78 EeV to 166 EeV. Details are also given of a further nine very-energetic events that have been used in the calibration procedure adopted to determine the energy of each primary. A sky plot of the arrival directions of the most energetic particles is shown. No interpretations of the data are offered.

This study analyzes 338 new times of visual, CCD, photoelectric and photographic maxima of the classical Cepheid SV Vul. The corresponding observations were made between 1913 and 2022. On this new and large observational basis, the period variations of this star of major astrophysical interest are re-visited. Without contradicting the theory or previous studies, we show that visual observations are important for a long-term monitoring of the period variations of well-selected bright Cepheids. This study establishes the period rate of change of SV Vul at -250 s/yr. Its period is currently shorter than 45 days.

We present PERISTOLE to study the various time delays associated with the pulsar rotation and other general relativistic aspects of binary pulsars. It is made available as an open-source python package which takes some parameters of the double pulsar system as input and outputs the rotational and latitudinal lensing delays along with the geometric and Shapiro delays that arise due to gravitational lensing. This package was intended to provide a way to quickly analyse, evaluate and study the differences between variations of the same systems and also to quantify the consequences that different parameters have over the system. Through this research note, we briefly describe the motivation behind PERISTOLE and showcase its capabilities using the only double pulsar system ever found, J0737-3039.

The Fortran-versions of the CORSIKA air shower simulation code have been at the core of simulations for many astroparticle physics experiments for the last 30 years. Having grown over decades into an ever more complex software, maintainability of CORSIKA has become increasingly difficult, though its performance is still excellent. In 2018, therefore a complete rewrite of CORSIKA has begun in modern modular C++. Today, CORSIKA 8 has reached important milestones with a full-fledged implementation of both the hadronic and electromagnetic cascades, the ability to simulate radio and Cherenkov-light emission from air showers and an unprecedented flexibility to configure simulation media and their geometries. This presentation will discuss the current status of CORSIKA 8, highlight the new possibilities already available, and future prospects of this new air shower simulation framework.

In this study we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian Parameter Estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington-limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124; the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, for its luminosity span several orders of magnitude during its outburst, with a maximum $L_{\rm X}$ that exceeded the Eddington limit by a factor of $\sim 10$. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius and the Alfv\'en radius are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column.

Binary solar system objects are common and range from satellite systems with very large mass ratios $M_1/M_2$ to mass ratios very close to unity. A well-known example of a binary is the Pluto-Charon system. With Charon only eight times less massive than Pluto the question arises as for many other systems, why the mass-ratio is still close to unity. There is much evidence that (binary) planet(esimal) formation happened early, when the protoplanetary gas disk was still around. It is likely that (some of) these binaries grew up together subject to pebble accretion. Here we focus on the question of how the mass arriving in the gravitational influence zone of the binary during pebble accretion, is distributed over the binary components. Does the accretion through time lead to a converging mass ratio, or to a diverging mass ratio? We numerically integrate pebble paths in the same well-known fashion as for a single mass subject to pebble accretion and track what the efficiency of accretion is for the two separate binary components, compared to a single body with the same mass. These numerical simulations are done for a range of binary mass-ratios, mutual separations, Stokes numbers and two orbital distances, 2.5 and 39 au. We find that in the limit where pebbles start to spiral around the primary (this holds for relatively large pebbles), the pebble preferentially collides with the secondary, causing the mass ratio to converge towards unity on Myr timescales. In this regime the total sweep-up efficiency can lower to half that of a pebble-accreting single body because pebbles that are thrown out of the system, after close encounters with the system. The results show that systems such as Pluto-Charon and other larger equal mass binaries could well have co-accreted by means of pebble accretion in the disk phase without producing binaries with highly diverging mass-ratios.

This paper represents an effort to provide robust constraints on cosmology and the galaxy-halo connection using a fully numerical model of small-scale galaxy clustering. We explore two extensions to the standard Halo Occupation Distribution model: assembly bias, whereby halo occupation depends on both halo mass and the larger environment, and velocity bias, whereby galaxy velocities do not perfectly trace the velocity of the dark matter within the halo. Moreover, we incorporate halo mass corrections to account for the impact of baryonic physics on the halo population. We identify an optimal set of clustering measurements to constrain this "decorated" HOD model for both low- and high-luminosity galaxies in SDSS DR7. We find that, for low-luminosity galaxies, a model with both assembly bias and velocity bias provides the best fit to the clustering measurements, with no tension remaining in the fit. In this model we find evidence for both central and satellite galaxy assembly bias at the 99% and 95% confidence levels, respectively. In addition, we find evidence for satellite galaxy velocity bias at the 99.9% confidence level. For high luminosity galaxies, we find no evidence for either assembly bias or velocity bias, but our model exhibits significant tension with SDSS measurements. We find that all of these conclusions still stand when we include the effects of baryonic physics on the halo mass function, suggesting that the tension we find for high luminosity galaxies may be due to a problem with our assumed cosmological model.

Recent observations of the pre-stellar core L1544 and the younger starless core L1498 have revealed that complex organic molecules (COMs) are enhanced in the gas phase toward their outer and intermediate-density shells. Our goal is to determine the level of chemical complexity toward the starless core L1517B, which seems younger than L1498, and compare it with the other two previously studied cores to see if there is a chemical evolution within the cores. We have carried out 3 mm high-sensitivity observations toward two positions in the L1517B starless core: the core's centre and the position where the methanol emission peaks (at a distance of $\sim$5000 au from the core's centre). Our observations reveal that a lower number of COMs and COM precursors are detected in L1517B with respect to L1498 and L1544, and also show lower abundances. Besides methanol, we only detected H$_2$CCO, CH$_3$CHO, CH$_3$CN, CH$_3$NC, HCCCN, and HCCNC. Their measured abundances are $\sim$3 times larger toward the methanol peak than toward the core's centre, mimicking the behaviour found toward the more evolved cores L1544 and L1498. We propose that the differences in the chemical complexity observed between the three studied starless cores are a consequence of their evolution, with L1517B being the less evolved one, followed by L1498 and L1544. Chemical complexity in these cores seems to increase over time, with N-bearing molecules forming first and O-bearing COMs forming at a later stage as a result of the catastrophic depletion of CO.

We present here the first ever mid-infrared spectroscopic time series observation of the transiting exoplanet \object{L 168-9 b} with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope. The data were obtained as part of the MIRI commissioning activities, to characterize the performance of the Low Resolution Spectroscopy (LRS) mode for these challenging observations. To assess the MIRI LRS performance, we performed two independent analyses of the data. We find that with a single transit observation we reached a spectro-photometric precision of $\sim$50 ppm in the 7-8~\micron~range at R=50, consistent with $\sim$25 ppm systematic noise. The derived band averaged transit depth is 524~$\pm$~15~ppm and 547~$\pm$~13~ppm for the two applied analysis methods, respectively, recovering the known transit depth to within 1~$\sigma$. The measured noise in the planet's transmission spectrum is approximately 15-20\% higher than random noise simulations over wavelengths $6.8 \lesssim \lambda \lesssim 11$ $\mu$m. We observed an excess noise at shorter wavelengths, for which possible causes are discussed. This performance was achieved with limited in-flight calibration data, demonstrating the future potential of MIRI for the characterization of exoplanet atmospheres.

The Cosmic Microwave Background (CMB) anisotropies are thought to be statistically isotropic and Gaussian. However, several anomalies are observed, including the CMB Cold Spot, an unexpected cold $\sim 10^{\circ}$ region with $p$-value $\lesssim 0.01$ in standard $\Lambda$CDM. One of the proposed origins of the Cold Spot is an unusually large void on the line of sight, that would generate a cold region through the combination of integrated Sachs-Wolfe and Rees-Sciama effects. In the past decade extensive searches were conducted in large scale structure surveys, both in optical and infrared, in the same area for $z \lesssim 1$ and did find evidence of large voids, but of depth and size able to account for only a fraction of the anomaly. Here we analyze the lensing signal in the Planck CMB data and rule out the hypothesis that the Cold Spot could be due to a large void located anywhere between us and the surface of last scattering. In particular, computing the evidence ratio we find that a model with a large void is disfavored compared to $\Lambda$CDM, with odds 1 : 13 (1 : 20) for SMICA (NILC) maps, compared to the original odds 56 : 1 (21 : 1) using temperature data alone.

Interstellar and cometary ices play an important role in the formation of planetary systems around young stars. Their main constituent is amorphous solid water (ASW). Although ASW is widely studied, vibrational energy dissipation and structural changes due to vibrational excitation are less well understood. The hydrogen-bonding network is likely a crucial component in this. Here we present experimental results on hydrogen-bonding changes in ASW induced by the intense, nearly monochromatic mid-IR free-electron laser (FEL) radiation of the FELIX-2 beamline at the HFML-FELIX facility at the Radboud University in Nijmegen, the Netherlands. Structural changes in ASW are monitored by reflection-absorption infrared spectroscopy and depend on the irradiation history of the ice. The experiments show that FEL irradiation can induce changes in the local neighborhood of the excited molecules due to energy transfer. Molecular Dynamics simulations confirm this picture: vibrationally excited molecules can reorient for a more optimal tetrahedral surrounding without breaking existing hydrogen bonds. The vibrational energy can transfer through the hydrogen-bonding network to water molecules that have the same vibrational frequency. We hence expect a reduced energy dissipation in amorphous material with respect to crystalline material due to the inhomogeneity in vibrational frequencies as well as the presence of specific hydrogen-bonding defect sites which can also hamper the energy transfer.

Several planetary systems are known to host multiple giant planets. However, when two giant planets are accreting from the same disk, it is unclear what effect the presence of the second planet has on the gas accretion process of both planets. In this paper we perform long-term 2D isothermal hydrodynamical simulations (over more than 0.5 Myrs) with the FARGO-2D1D code, considering two non-migrating planets accreting from the same gaseous disk. We find that the evolution of the planets' mass ratio depends on gap formation. However, in all cases, when the planets start accreting at the same time, they end up with very similar masses (0.9 $<m_{p,out}/m_{p,in}<$ 1.1 after 0.5 Myrs). Delaying the onset of accretion of one planet allows the planets' mass ratio to reach larger values initially, but they quickly converge to similar masses afterward (0.8 $<m_{p,out}/m_{p,in}<$ 2 in $10^5$ yrs). In order to reproduce the more diverse observed mass ratios of exoplanets, the planets must start accreting gas at different times, and their accretion must be stopped quickly after the beginning of runaway gas accretion (less than 0.5 Myrs), for example via disk dispersal. The evolution of the planets' mass ratio can have an important impact on the dynamics of the system and may constrain the formation history of Jupiter and Saturn.

We have recently hit the milestone of 5,000 exoplanets discovered. In stark contrast with the Solar System, most of the exoplanets we know to date orbit extremely close to their host stars, causing them to lose copious amounts of gas through atmospheric escape at some stage in their lives. In some planets, this process can be so dramatic that they shrink in timescales of a few million to billions of years, imprinting features in the demographics of transiting exoplanets. Depending on the transit geometry, ionizing conditions, and atmospheric properties, a planetary outflow can be observed using transmission spectroscopy in the ultraviolet, optical or near-infrared. In this review, we will discuss the main techniques to observe evaporating exoplanets and their results. To date, we have evidence that at least 28 exoplanets are currently losing their atmospheres, and the literature has reported at least 42 non-detections.

We present the first rest-frame optical spectrum of a high-redshift quasar observed with JWST/NIRCam in Wide Field Slitless (WFSS) mode. The observed quasar, J0100+2802, is the most luminous quasar known at $z>6$. We measure the mass of the central supermassive black hole (SMBH) by means of the rest-frame optical H$\beta$ emission line, and find consistent mass measurements of the quasar's SMBH of $M_\bullet\approx10^{10}\,M_\odot$ when compared to the estimates based on the properties of rest-frame UV emission lines CIV and MgII, which are accessible from ground-based observatories. To this end, we also present a newly reduced rest-frame UV spectrum of the quasar observed with X-Shooter/VLT and FIRE/Magellan for a total of 16.8 hours. We readdress the question whether this ultra-luminous quasar could be effected by strong gravitational lensing making use of the diffraction limited NIRCam images in three different wide band filters (F115W, F200W, F356W), which improves the achieved spatial resolution compared to previous images taken with the Hubble Space Telescope by a factor of two. We do not find any evidence for a foreground deflecting galaxy, nor for multiple images of the quasar, and determine the probability for magnification due to strong gravitational lensing with image separations below the diffraction limit of $\Delta\theta\lesssim 0.05''$ to be $\lesssim 2.2\times 10^{-3}$. Our observations therefore confirm that this quasar hosts a ten billion solar mass black hole less than $1$ Gyr after the Big Bang, which is challenging to explain with current black hole formation models.

The bulge of M31 is of interest in the context of the nature of galactic bulges and how their structure relates to bulge formation mechanisms and their subsequent evolution. With the UVIT instrument on AstroSat, we have observed the bulge of M31 in five far ultraviolet (FUV) and near ultraviolet (NUV) filters at 1" spatial resolution. Models for the luminosity distribution of the bulge are constructed using the UVIT data and the galaxy image-fitting algorithm GALFIT. We fit the bulge without nuclear region with a Sersic function the five images and find Sersic indices ($\simeq2.1$ to 2.5) similar to previous studies but smaller $R_e$ values ($\simeq0.5$ to 0.6 kpc). When the images include the nuclear region, a multicomponent model is used to find the best-fit. We use an 8-component fit for the FUV 148nm image, which has the highest sensitivity. The other images (169 to 279 nm) are fit with 4-component models. The dust lanes in the bulge region are recovered in the residual image after subtraction of the bright bulge light using the multicomponent model. The dust lanes show that M31's nuclear spiral is visible in absorption at NUV and FUV wavelengths. The bulge images show boxy contours in all five UVIT wavebands, which is confirmed by fitting using GALFIT. The Sersic indices of $\sim$2.1-2.5 are intermediate between expected values for a classical bulge and for a pseudobulge. The boxiness of the bulge provides further evidence that M31's bulge has contributions from a classical bulge and a pseudobulge.

In recent years, a handful of ``dark" binaries have been discovered with a non-luminous compact object. Astrometry and radial velocity measurements of the bright companion allow us to measure the post-supernova orbital elements of such a binary. In this paper, we develop a statistical formalism to use such measurements to infer the pre-supernova orbital elements, and the natal kick imparted by the supernova (SN). We apply this formalism to the recent discovery of an X-ray quiet binary with a black hole, VFTS 243, in the Large Magellanic Cloud. Assuming an isotropic, Maxwellian distribution on natal kicks and using broad agnostic mass priors, we find that kick velocity can be constrained to less than $V_k < 56$ km/s and the dispersion of the kick distribution to $\sigma_k < 66 $ km/s at 90 \% confidence. Assuming that there was a natal kick, we find that at least about $0.6 M_{\odot}$ was lost during the supernova with 90 \% confidence. The pre-SN orbital separation is found to be robustly constrained to be around $0.43$ AU.

We establish constraints on $f(T)$ gravity by considering the possibility of a scenario that supports a phantom crossing of the equation of state parameter $w_{DE}$. After determining the viable parameter space of the model, while checking the impact on the background dynamics, we perform an analysis to obtain constraints on cosmological parameters and determine the viability of this scenario. To this end, we use combined data sets from cosmic chronometers (CC), baryonic acoustic oscillations (BAO), redshift space measurements (RSD), cosmic microwave background (CMB) and Big Bang nucleosynthesis (BBN) priors and Type Ia supernovae (SNe) measurements from the Pantheon set, in which the impact on the absolute magnitude is also considered. It is found that a state where a phantom crossing of $w_{DE}$ happens is favored by data, and while the $f(T)$ model is equivalent to the $\Lambda$CDM one by statistical criteria when CC+SNe+BAO+BBN data is considered, a considerable tension is found when CMB data is added. We also find that the Hubble tension is alleviated in the $f(T)$ model, at the same time that it does not worsen the growth one, indicating a possibility of the scenario as an option to address the current cosmic tensions.

The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $\sigma = 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for.

Context. The present study addresses a key question for our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore how empty of matter local Universe voids are, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained within the CosmicFlows-3 volume. Methods. The cosmic voids matter density profiles are independently obtained using two different methods. They are built on one hand from the galaxy redshift space 2 points-correlation function and, on the other hand, using peculiar velocity gradients from the CosmicFlows-3 dataset. Results. The results are noticeable since when using the redshift survey, all voids show a radial positive gradient of galaxies, while based on the dynamical analysis, only three of these voids display a clear underdensity of matter in their center. Conclusions. It is the first time such a detailed observational analysis of voids is conducted, showing that void emptiness should be derived from dynamical information. Yet, from this limited study, the Hercules void is the best candidate for a local Universe pure "pristine volume" expanding in 3 directions with no dark matter located in that void.

This article publicly releases three-dimensional reconstructions of the local Universe gravitational field below z=0.8 that were computed using the full catalogue CosmicFlows-4 of 56,000 galaxy distances and its sub-sample of 1,008 type Ia supernovae distances. The article also provides some first CF4 measurements of the growth rate of structure using the pairwise correlation of peculiar velocities fsigma8 = 0.44(+/-0.01) and of the bulk flow in the Local Universe of 200+/-88 kms-1 at distance 300 h-1Mpc.

We used 22 $\mu$m (W4) Wide-field Infrared Survey Explorer (WISE) observations of 4420 asteroids to analyze lightcurves and determined spin period estimates for 1929 asteroids. We fit second-order Fourier models at a large number of trial frequencies to the W4 data and analyzed the resulting periodograms. We initially excluded rotational frequencies exceeding 7.57 rotations per day (P < 3.17 hr), which are not sampled adequately by WISE, and periods that exceed twice the WISE observation interval, which is typically 36 hr. Three solutions accurately capture the vast majority of the rotational frequencies in our sample: the best-fit frequency and its mirrors around 3.78 and 7.57 rotations per day. By comparing our solutions to a high-quality control group of 752 asteroid spin periods, we found that one of our solutions is accurate (within 5%) in 88% of the cases. The best-fit, secondary, and tertiary solutions are accurate in 55%, 27%, and 6% of the cases, respectively.

The Earth's turbulent atmosphere results in speckled and blurred images of astronomical objects when observed by ground based visible and near-infrared telescopes. Adaptive optics (AO) systems are employed to reduce these atmospheric effects by using wavefront sensors (WFS) and deformable mirrors. Some AO systems are not fast enough to correct for strong, fast, high turbulence wind layers leading to the wind butterfly effect, or wind-driven halo, reducing contrast capabilities in coronagraphic images. Estimating the effective wind speed of the atmosphere allows us to calculate the atmospheric coherence time. This is not only an important parameter to understand for site characterization but could be used to help remove the wind butterfly in post processing. Here we present a method for estimating the atmospheric effective wind speed from spatio-temporal covariance maps generated from pseudo open-loop (POL) WFS data. POL WFS data is used as it aims to reconstruct the full wavefront information when operating in closed-loop. The covariance maps show how different atmospheric turbulent layers traverse the telescope. Our method successfully recovered the effective wind speed from simulated WFS data generated with the soapy python library. The simulated atmospheric turbulence profiles consist of two turbulent layers of ranging strengths and velocities. The method has also been applied to Gemini Planet Imager (GPI) AO WFS data. This gives insight into how the effective wind speed can affect the wind-driven halo seen in the AO image point spread function. In this paper, we will present results from simulated and GPI WFS data.

The empirical classification of gamma-ray bursts (GRBs) into long and short GRBs based on their durations is already firmly established. This empirical classification is generally linked to the physical classification of GRBs originating from compact binary mergers and GRBs originating from massive star collapses, or Type I and II GRBs, with the majority of short GRBs belonging to Type I and the majority of long GRBs belonging to Type II. However, there is a significant overlap in the duration distributions of long and short GRBs. Furthermore, some intermingled GRBs, i.e., short-duration Type II and long-duration Type I GRBs, have been reported. A multi-wavelength, multi-parameter classification scheme of GRBs is evidently needed. In this paper, we seek to build such a classification scheme with supervised machine learning methods, chiefly XGBoost. We utilize the GRB Big Table and Greiner's GRB catalog and divide the input features into three subgroups: prompt emission, afterglow, and host galaxy. We find that the prompt emission subgroup performs the best in distinguishing between Type I and II GRBs. We also find the most important distinguishing feature in prompt emission to be $T_{90}$, hardness ratio, and fluence. After building the machine learning model, we apply it to the currently unclassified GRBs to predict their probabilities of being either GRB class, and we assign the most probable class of each GRB to be its possible physical class.

Modelling the integrated H I spectra of galaxies has been a difficult task due to their diverse shapes, but more dynamical information is waiting to be explored in Hi line profiles. Based on simple assumptions, we construct a physically motivated model for the integrated Hi spectra: Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization (PANDISC). The model shows great flexibility in reproducing the diverse Hi profiles. We use Monte-Carlo Markov Chain (MCMC) for fitting the model to global H I profiles, producing statistically robust quantitative results. Comparing with several samples of H I data available in the literature , we find the model-fitted results agree with catalogued velocity widths (e.g., W50) down to the lowest S/N. The model is also shown to be useful for applications like the baryonic Tully-Fisher relation (BTFR) and profile-based sample control. By comparing v_r to v_flat , we uncover how the H I width is affected by the structure of the rotation curve, following a trend consistent with the difference in the BTFR slope. We also select a sample of spectra with broad wing-like features suggestive of a population of galaxies with unusual gas dynamics. The PANDISC model bears both promise and limitations for potential use beyond H I lines. Further application on the whole ALFALFA sample will enable us to perform large scale ensemble studies of the H I properties and dynamics in nearby galaxies.

In a novel approach employing implicit likelihood inference (ILI), also known as likelihood-free inference, we calibrate the parameters of cosmological hydrodynamic simulations against observations, which has previously been unfeasible due to the high computational cost of these simulations. For computational efficiency, we train neural networks as emulators on ~1000 cosmological simulations from the CAMELS project to estimate simulated observables, taking as input the cosmological and astrophysical parameters, and use these emulators as surrogates to the cosmological simulations. Using the cosmic star formation rate density (SFRD) and, separately, stellar mass functions (SMFs) at different redshifts, we perform ILI on selected cosmological and astrophysical parameters (Omega_m, sigma_8, stellar wind feedback, and kinetic black hole feedback) and obtain full 6-dimensional posterior distributions. In the performance test, the ILI from the emulated SFRD (SMFs) can recover the target observables with a relative error of 0.17% (0.4%). We find that degeneracies exist between the parameters inferred from the emulated SFRD, confirmed with new full cosmological simulations. We also find that the SMFs can break the degeneracy in the SFRD, which indicates that the SMFs provide complementary constraints for the parameters. Further, we find that the parameter combination inferred from an observationally-inferred SFRD reproduces the target observed SFRD very well, whereas, in the case of the SMFs, the inferred and observed SMFs show significant discrepancies that indicate potential limitations of the current galaxy formation modeling and calibration framework, and/or systematic differences and inconsistencies between observations of the stellar mass function.

We study vertical variations of wave-periods of magnetoacoustic two-fluid waves in the partially ionized lower solar atmosphere, consisting of ion (proton) + electron and neutral (atomic hydrogen) fluids, which are coupled by ion-neutral collisions. The study allows finding the wave period cutoffs and their variations in the solar atmosphere, as well as establishing the role of these cutoffs in determining the wave propagation conditions. The atmosphere is permitted by a uniform vertical magnetic field. We perform numerical simulations in the framework of a one-dimensional (1D), two-fluid model in which plane waves are exited by a harmonic driver in the vertical ion and neutral velocities, operating at the bottom of the solar photosphere. We observe excitation of waves with cutoff wave-periods in addition to waves set directly by the driver. We also see that some waves exited by that driver can reach the solar corona. Despite of its limitations such as the lack of non-adiabatic and non-ideal terms and a simple 1D structure, the developed two-fluid model of the solar atmosphere sheds a new light on the role of cutoffs in setting up the wave propagation conditions in the solar atmosphere and finding periods of waves that may carry their energy from the solar surface to the corona.

Star formation theories have struggled to reproduce binary brown dwarf population demographics (frequency, separation, mass-ratio). Kernel-phase interferometry is sensitive to companions at separations inaccessible to classical imaging, enabling tests of formation at new physical scales below the hydrogen burning limit. We analyze the detections and sensitivity limits from our previous kernel-phase analysis of archival HST/NICMOS surveys of field brown dwarfs. After estimating physical properties of the 105 late M to T dwarfs using Gaia distances and evolutionary models, we use a Bayesian framework to compare these results to a model companion population defined by log-normal separation and power-law mass-ratio distributions. When correcting for Malmquist bias, we find a companion fraction of $F=0.11^{+0.04}_{-0.03}$ and a separation distribution centered at $\rho=2.2^{+1.2}_{-1.0}$ au, smaller and tighter than seen in previous studies. We also find a mass-ratio power-law index which strongly favors equal-mass systems: $\gamma=4.0^{+1.7}_{-1.5}-11^{+4}_{-3}$ depending on the assumed age of the field population ($0.9-3.1$ Gyr). We attribute the change in values to our use of kernel-phase interferometry which enables us to resolve the peak of the semimajor axis distribution with significant sensitivity to low-mass companions. We confirm the previously-seen trends of decreasing binary fraction with decreasing mass and a strong preference for tight and equal-mass systems in the field-age sub-stellar regime; only $0.9^{+1.1}_{-0.6}$ % of systems are wider than 20 au and $<1.0^{+1.4}_{-0.6}$% of systems have a mass-ratio $q<0.6$. We attribute this to turbulent fragmentation setting the initial conditions followed by a brief period of dynamical evolution, removing the widest and lowest-mass companions, before the birth cluster dissolves.

The expansion of the universe on short distance scales is a new frontier to investigate the dark energy. The excess orbital decay in binary pulsars may be related to acceleration by the local cosmic expansion, called the cosmic drag. Modern observations of two independent binaries (PSR B1534+12 and PSR B1913+16) support this interpretation and result in a scale-independent expansion with viscous uniformity, in which binary systems have a smaller expansion rate than the Hubble constant. This paper shows additional evidential binaries (PSR J1012-5307 and PSR J1906+0746), supporting the cosmic drag picture. The total anomaly of the conventional model is about $3.6\,\sigma$ including two evidential binaries reported before. In addition, an observable range of the cosmic drag has been calculated for typical models of both NS-NS binary and NS-WD binary. In this region, six test candidates are listed with predictions of the excess orbital decay.

The hot, neutron-rich, and dense circumstance in core-collapse supernovae provides a source of negatively charged pions that may make up a significant portion of the matter. These abundant thermal pions can play a role to populate light and hidden hypothetical particles. We discuss the dark gauge boson production via reactions involving supernova pions, the rate of which is determined by the isovector nucleon coupling. We take into account the two toy models, the dark photon and the gauged $B-L$ models, that carry the typical distinct isovector nucleon coupling structure in the medium. Pion-induced dark gauge bosons leave an imprint on several observational consequences associated with supernova. Their sizable emissivity and characteristic hard spectral distribution result in the stringent constraints on the dark gauge boson models, in particular at masses above the two electron mass.

Cosmological relaxation of the electroweak scale via Higgs-axion interplay, named as relaxion mechanism, provides a dynamical solution to the Higgs mass hierarchy. In the original proposal by Graham, Kaplan and Rajendran, the relaxion abundance today is too small to explain the dark matter of the universe because of the high suppression of the misalignment angle after inflation. It was then realised by Banerjee, Kim and Perez that reheating effects can displace the relaxion, thus enabling it to account for the dark matter abundance from the misalignment mechanism. However, this scenario is realised in a limited region of parameter space to avoid runaway. We show that in the regime where inflationary fluctuations dominate over the classical slow-roll, the "stochastic misalignment" of the field due to fluctuations can be large. We study the evolution of the relaxion after inflation, including the high-temperature scenario, in which the barriers of the potential shrink and destabilise temporarily the local minimum. We open new regions of parameter space where the relaxion can naturally explain the observed dark matter density in the universe, towards larger coupling, larger mass, larger mixing angle, smaller decay constant, as well as larger scale of inflation.

A GeV scale self-interacting dark matter (SIDM) with a light mediator is a promising scenario to alleviate the small-scale problems of the cold dark matter (CDM) paradigm while being consistent with the latter at large scales, as suggested by astrophysical observations. Conventionally such a scenario leads to an under-abundant SIDM thermal relic due to large annihilation rates of SIDM into mediator particles. In this letter, we propose a minimal realisation of GeV scale SIDM with correct thermal relic, where one of the three singlet fermions, responsible for seesaw origin of light neutrino masses, assists in generating correct SIDM relic density via freeze-out. The singlet fermion, having mass close to but greater than that of SIDM, can lead to efficient annihilation rate into dark matter, thus compensating for the latter's depletion due to annihilation into light mediators. In addition to providing a thermal relic generation mechanism of such light SIDM, there exist exciting discovery prospects at direct search, cosmology and dark photon search experiments. Connection to the origin of light neutrino masses also opens up indirect detection prospects due to long-lived nature of DM.

The Standard Model, extended with three right-handed (RH) neutrinos, is the simplest model that can explain light neutrino masses, the baryon asymmetry of the Universe, and dark matter (DM). Models in which RH neutrinos are light are generally easier to test in experiments. In this work, we show that, even if the RH neutrinos are super-heavy ($M_{i=1,2,3}>10^9$ GeV) -- close to the Grand Unification scale -- the model can be tested thanks to its distinct features on the stochastic Gravitational Wave (GW) background. We consider an early Universe filled with ultralight primordial black holes (PBH) that produce a super-heavy RH neutrino DM via Hawking radiation. The other pair of RH neutrinos generates the baryon asymmetry via thermal leptogenesis, much before the PBHs evaporate. GW interferometers can test this novel spectrum of masses thanks to the GWs induced by the PBH density fluctuations. In a more refined version, wherein a $U(1)$ gauge symmetry breaking dynamically generates the seesaw scale, the PBHs also cause observable spectral distortions on the GWs from the $U(1)$-breaking cosmic strings. Thence, a low-frequency GW feature related to DM genesis and detectable with a pulsar-timing array must correspond to a mid- or high-frequency GW signature related to baryogenesis at interferometer scales.

It is shown that a sub-luminal electromagnetic plasma wave, propagating in phase with a background sub-luminal gravitational wave in a dispersive medium, can undergo parametric amplification. For this phenomena to occur, the dispersive characteristics of the two waves must properly match. The response frequencies of the two waves (medium dependent) must lie within a definite and restrictive range. The combined dynamics is represented by a Whitaker-Hill equation, the quintessential model for parametric instabilities. The exponential growth of the electromagnetic wave is displayed at the resonance; the plasma wave grows at the expense of the background gravitational wave. Different physical scenarios, where the phenomenon can be possible, are discussed.

Variations in the geomagnetic field occur on a vast range of time scales, from milliseconds to millions of years. The advent of satellite measurements has allowed for detailed studies of the short timescale geomagnetic field behaviour, but understanding the long timescale evolution remains challenging due to the sparsity of the paleomagnetic record. This paper introduces a field theory framework for studying magnetic field generation as a result of stochastic fluid motions. By constructing a stochastic kinematic dynamo model, we derive statistical properties of the magnetic field that may be compared to observations from the paleomagnetic record. The fluid velocity is taken to act as a random forcing obeying Gaussian statistics. Using the Martin-Siggia-Rose-Janssen-de Dominicis (MSRJD) formalism, we compute the average magnetic field response function. From this we obtain an estimate for the turbulent contribution to the magnetic diffusivity, and find that it is consistent with results from mean-field dynamo theory. This framework presents much promise for studying the geomagnetic field in a stochastic context.

Bigravity is one of the natural extensions of general relativity and contains an additional massive spin-2 field which can be a good candidate for dark matter. To discuss the production of spin-2 dark matter, we study fixed point solutions of the background equations for axisymmetric Bianchi type-I Universes in two bigravity theories without Boulware-Deser ghost, i.e., Hassan-Rosen bigravity and Minimal Theory of Bigravity. We investigate the local and global stability of the fixed points and classify them. Based on the general analysis, we propose a new scenario where spin-2 dark matter is produced by the transition from an anisotropic fixed point solution to isotropic one. The produced spin-2 dark matter can account for all or a part of dark matter and can be directly detected by laser interferometers in the same way as gravitational waves.

The electroweak axion is identified with the fuzzy dark matter of a mass $m\simeq 10^{-20}$--$10^{-19}\,{\rm eV}$. The model predicts two components of dark matter, one is ultralight and the other is WIMP-like. The Chern-Simons-type interaction between the fuzzy dark matter and photon and the $B+L$ breaking proton decays are predicted.

We investigate production of primordial black holes from first-order electroweak phase transition in the framework of the nearly aligned Higgs effective field theory, in which non-decoupling quantum effects are properly described. Since the mass of such primordial black holes is evaluated to be about $10^{-5}$ of the solar mass, current and future microlensing observations such as Subaru HSC, OGLE, PRIME and Roman Space Telescope may be able to probe the electroweak phase transition. We study parameter regions where primordial black holes can be produced by the first-order electroweak phase transition, and explore their detectability at these observations. Complementarity of primordial black hole observations, gravitational wave observations and collider experiments is also discussed for testing the nature of the electroweak phase transition.

The Hubble tension is analyzed in the framework of quantum cosmological approach. It is found that there arises a new summand in the expression for the total energy density stipulated by the quantum Bohm potential. This additional energy density modifies the expansion history of the early universe and decays faster than radiation in late universe. Similarly to physical models with early dark energy, taking account of this matter source of quantum nature can, in principle, eliminate a discrepancy between the Hubble constant estimates obtained in different approaches. The model been considered allows one to extend the standard cosmology to quantum sector.