Papers for Monday, Mar 21 2022

In antiquity, all of the enduring celestial bodies that were seen to move relative to the background sky of stars were considered planets. During the Copernican revolution, this definition was altered to objects orbiting around the Sun, removing the Sun and Moon but adding the Earth to the list of known planets. The concept of planet is thus not simply a question of nature, origin, composition, mass or size, but historically a concept related to the motion of one body around the other, in a hierarchical configuration. After discussion within the IAU Commission F2 "Exoplanets and the Solar System", the criterion of the star-planet mass ratio has been introduced in the definition of the term "exoplanet", thereby requiring the hierarchical structure seen in our Solar System for an object to be referred to as an exoplanet. Additionally, the planetary mass objects orbiting brown dwarfs, provided they follow the mass ratio criterion, are now considered as exoplanets. Therefore, the current working definition of an exoplanet, as amended in August 2018 by IAU Commission F2 "Exoplanets and the Solar System", reads as follows: - Objects with true masses below the limiting mass for thermonuclear fusion of deuterium (currently calculated to be 13 Jupiter masses for objects of solar metallicity) that orbit stars, brown dwarfs or stellar remnants and that have a mass ratio with the central object below the $L_4$/$L_5$ instability ($M/M_{\rm central}$$<$$2/(25+\sqrt{621}$)$\approx$1/25) are "planets", no matter how they formed. - The minimum mass/size required for an extrasolar object to be considered a planet should be the same as that used in our Solar System, which is a mass sufficient both for self-gravity to overcome rigid body forces and for clearing the neighborhood around the object's orbit. Here we discuss the history and the rationale behind this definition.

We use the SMASH survey to obtain unprecedented deep photometry reaching down to the oldest main sequence turn-offs in the colour-magnitude diagrams (CMDs) of the Small Magellanic Cloud (SMC) and quantitatively derive its star formation history (SFH) using CMD fitting techniques. We identify five distinctive peaks of star formation in the last 3.5 Gyr, at $\sim $3, $\sim$2, $\sim$1.1, $\sim $0.45 Gyr ago, and one presently. We compare these to the SFH of the Large Magellanic Cloud (LMC) finding unequivocal synchronicity, with both galaxies displaying similar periods of enhanced star formation over the past $\sim$3.5 Gyr. The parallelism between their SFHs indicates that tidal interactions between the MCs have recurrently played an important role in their evolution for at least the last $\sim$3.5 Gyr, tidally truncating the SMC and shaping the LMC's spiral arm. We show, for the first time, an SMC-LMC correlated SFH at recent times in which enhancements of star formation are localised in the northern spiral arm of the LMC, and globally across the SMC. These novel findings should be used to constrain not only the orbital history of the MCs but also how star formation should be treated in simulations.

Observations indicate that the convective cores of stars must ingest a substantial amount of material from the overlying radiative zone, but the extent of this mixing and the mechanism which causes it remain uncertain. Recently, Anders et al. (2021) developed a theory of convective penetration and calibrated it with 3D numerical hydrodynamics simulations. Here we employ that theory to predict the extent of convective boundary mixing in early-type main-sequence stars. We find that convective penetration produces enough mixing to explain core masses inferred from asteroseismology and eclipsing binary studies, and matches observed trends in mass and age. While there are remaining uncertainties in the theory, this agreement suggests that most convective boundary mixing in early-type main-sequence stars arises from convective penetration. Finally, we provide a fitting formula for the extent of core convective penetration for main-sequence stars in the mass range from $1.1-60 M_\odot$.

We report on the Swift/XRT observation and classification of eleven blazar candidates at $z>4$. These sources were selected as part of a sample of extremely radio-loud quasars, in order to focus on quasars with jets oriented roughly close to our line-of-sight. Deriving their viewing angles and their jets bulk Lorentz factors was crucial for a strict blazar classification, and it was possible only thanks to X-ray observations. Out of eleven sources, five show strong and hard X-ray fluxes, that allow their blazar classification, two are uncertain, three host relativistic jets that we observe just outside their beaming cone (i.e. are not strictly blazars), while one went undetected by Swift/XRT. Following this approach, we were able to trace the $>10^9M_\odot$ active supermassive black hole population hosted in jetted active galactic nuclei. At $z\geq4$ the massive jetted sources are likely predominant in the overall quasar population: this calls for a deep review of our understanding of the first supermassive black holes formation and evolution. Jets are indeed key actors in fast accretion, and must be searched for across the whole high redshift quasar population. A note of caution must be added: radio-loudness and in general radio features at high redshifts seem not to perfectly reflect high-energy properties. A strong effect due to interaction with CMB radiation is surely in place, that quenches the radio emission with respect to the X-rays, but also more frequent occasions for the jet to be bent seem to play a relevant role in this matter. Classifications and population studies thus must be carefully performed, in order not to be confused by these inconsistencies.

The past 5 years have dramatically changed our view of the disks of gas and dust around young stars. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and extreme adaptive optics systems have revealed that disks are dynamical systems. Most disks contain resolved structures, both in gas and dust, including rings, gaps, spirals, azimuthal dust concentrations, shadows cast by misaligned inner disks, as well as deviations from Keplerian rotation. The origin of these structures and how they relate to the planet formation process remain poorly understood. Spatially resolved kinematic studies offer a new and necessary window to understand and quantify the physical processes (turbulence, winds, radial and meridional flows, stellar multiplicity, instabilities) at play during planet formation and disk evolution. Recent progress, driven mainly by resolved ALMA observations, includes the detection and mass determination of embedded planets, the mapping of the gas flow around the accreting planets, the confirmation of tidal interactions and warped disk geometries, and stringent limits on the turbulent velocities. In this chapter, we will review our current understanding of these dynamical processes and highlight how kinematic mapping provides new ways to observe planet formation in action.

Strong Lensing is a powerful probe of the matter distribution in galaxies and clusters and a relevant tool for cosmography. Analyses of strong gravitational lenses with Deep Learning have become a popular approach due to these astronomical objects' rarity and image complexity. Next-generation surveys will provide more opportunities to derive science from these objects and an increasing data volume to be analyzed. However, finding strong lenses is challenging, as their number densities are orders of magnitude below those of galaxies. Therefore, specific Strong Lensing search algorithms are required to discover the highest number of systems possible with high purity and low false alarm rate. The need for better algorithms has prompted the development of an open community data science competition named Strong Gravitational Lensing Challenge (SGLC). This work presents the Deep Learning strategies and methodology used to design the highest-scoring algorithm in the II SGLC. We discuss the approach used for this dataset, the choice for a suitable architecture, particularly the use of a network with two branches to work with images in different resolutions, and its optimization. We also discuss the detectability limit, the lessons learned, and prospects for defining a tailor-made architecture in a survey in contrast to a general one. Finally, we release the models and discuss the best choice to easily adapt the model to a dataset representing a survey with a different instrument. This work helps to take a step towards efficient, adaptable and accurate analyses of strong lenses with deep learning frameworks.

The chemical composition of the highest end of the ultra-high-energy cosmic ray spectrum is very hard to measure experimentally, and to this day it remains mostly unknown. Since the trajectories of ultra-high-energy cosmic rays are deflected in the magnetic field of the Galaxy by an angle that depends on their atomic number $Z$, it could be possible to indirectly measure $Z$ by quantifying the amount of such magnetic deflections. In this paper we show that, using the angular harmonic cross-correlation between ultra-high-energy cosmic rays and galaxies, we could effectively distinguish different atomic numbers with current data. As an example, we show how, if $Z=1$, the cross-correlation can exclude a $39\%$ fraction of Fe56 nuclei at $2\sigma$ for rays above $100\text{EeV}$.

We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at $z$~2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionisation interstellar (LIS) metal absorption lines and Ly$\alpha$ emission. Outflow velocities are calculated from the centroids of the LIS absorption and/or Ly$\alpha$ emission, as well as the highest-velocity component of the outflow from the blue wings of the LIS absorption lines. Outflow velocities are found to be marginally correlated or independent of galaxy properties, such as star-formation rate (SFR) and star-formation rate surface density ($\Sigma_{\rm SFR}$). Outflow velocity scales with SFR as a power-law with index 0.24, which suggests that the outflows may be primarily driven by mechanical energy generated by supernovae explosions, as opposed to radiation pressure acting on dusty material. On the other hand, outflow velocity and $\Sigma_{\rm SFR}$ are not significantly correlated, which may be due to the limited dynamic range of $\Sigma_{\rm SFR}$ probed by our sample. The relationship between outflow velocity and $\Sigma_{\rm SFR}$ normalised by stellar mass ($\Sigma_{\rm sSFR}$), as a proxy for gravitational potential, suggests that strong outflows (e.g., > 200 km s$^{-1}$) appear ubiquitous above a threshold of log($\Sigma_{\rm sSFR}/\rm{yr}^{-1}\ \rm{kpc}^{-2}$) ~ -11.3, and that above this threshold, outflow velocity uncouples from $\Sigma_{\rm sSFR}$. These results highlight the need for higher resolution spectroscopic data and spatially resolved imaging to test the driving mechanisms of outflows predicted by theory.

Machine learning is a promising tool to reconstruct time-series phenomena, such as variability of active galactic nuclei (AGN), from sparsely-sampled data. Here we use three Continuous Auto-Regressive Moving Average (CARMA) representations of AGN variability -- the Damped Random Walk (DRW) and (over/under-)Damped Harmonic Oscillator (DHO) -- to simulate 10-year AGN light curves as they would appear in the upcoming Vera Rubin Observatory Legacy Survey of Space and Time (LSST), and provide a public tool to generate these for any survey cadence. We investigate the impact on AGN science of five proposed cadence strategies for LSST's primary Wide-Fast-Deep (WFD) survey. We apply for the first time in astronomy a novel Stochastic Recurrent Neural Network (SRNN) algorithm to reconstruct input light curves from the simulated LSST data, and provide a metric to evaluate how well SRNN can help recover the underlying CARMA parameters. We find that the light curve reconstruction is most sensitive to the duration of gaps between observing season, and that of the proposed cadences, those that change the balance between filters, or avoid having long gaps in the {g}-band perform better. Overall, SRNN is a promising means to reconstruct densely sampled AGN light curves and recover the long-term Structure Function of the DRW process (SF$_\infty$) reasonably well. However, we find that for all cadences, CARMA/SRNN models struggle to recover the decorrelation timescale ($\tau$) due to the long gaps in survey observations. This may indicate a major limitation in using LSST WFD data for AGN variability science.

Solar and Heliosphere physics are areas of remarkable data-driven discoveries. Recent advances in high-cadence, high-resolution multiwavelength observations, growing amounts of data from realistic modeling, and operational needs for uninterrupted science-quality data coverage generate the demand for a solar metadata standardization and overall healthy data infrastructure. This white paper is prepared as an effort of the working group "Uniform Semantics and Syntax of Solar Observations and Events" created within the "Towards Integration of Heliophysics Data, Modeling, and Analysis Tools" EarthCube Research Coordination Network (@HDMIEC RCN), with primary objectives to discuss current advances and identify future needs for the solar research cyberinfrastructure. The white paper summarizes presentations and discussions held during the special working group session at the EarthCube Annual Meeting on June 19th, 2020, as well as community contribution gathered during a series of preceding workshops and subsequent RCN working group sessions. The authors provide examples of the current standing of the solar research cyberinfrastructure, and describe the problems related to current data handling approaches. The list of the top-level recommendations agreed by the authors of the current white paper is presented at the beginning of the paper.

The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring formation and evolution scenarios. We include structures from galactic-scale filaments to tenth-of-a-parsec scale filaments, comprising both molecular and atomic structures, from both observational and theoretical perspectives. In addition to the literature overview, we assemble a large amount of catalogue data from different surveys and provide the most comprehensive census of filamentary structures to date. Our census consists of 22 803 filamentary structures, facilitating a holistic perspective and new insights. We use our census to conduct a meta-analysis, leading to a description of filament properties over four orders of magnitudes in length and eight in mass. Our analysis emphasises the hierarchical and dynamical nature of filamentary structures. Filaments do not live in isolation, nor they generally resemble static structures close to equilibrium. We propose that accretion during filament formation and evolution sets some of the key scaling properties of filaments. This highlights the role of accretion during filament formation and evolution and also in setting the initial conditions for star formation. Overall, the study of filamentary structures during the past decade has been observationally driven. While great progress has been made on measuring the basic properties of filaments, our understanding of their formation and evolution is clearly lacking. In this context, we identify a number of directions and questions we consider most pressing for the field.

We present a catalog of spectroscopically measured redshifts over $0 < z < 2$ and emission line fluxes for 1440 galaxies, primarily ($\sim$65\%) from the HALO7D survey, with other spectra from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with $\log(M_{\star}/M_{\odot}) < 9.5$. 810 of these galaxies did not have existing spectroscopic redshifts in our catalogs, including 454 dwarf galaxies. Spectroscopic observations using the DEIMOS spectrograph on the Keck II telescope provide very deep (up to 32 hours exposure, with a median of $\sim$7 hours) optical spectroscopy for this sample of galaxies in the COSMOS, EGS, GOODS-North, and GOODS-South CANDELS fields, and in some areas outside CANDELS. We compare our results to existing spectroscopic and photometric redshifts in these fields, finding only a 1\% rate of discrepancy with other spectroscopic redshifts. We measure a small increase in median photometric redshift error (from 1.0\% to 1.3\%) and catastrophic outlier rate (from 3.5\% to 8\%) with decreasing stellar mass. We obtained successful redshift fits for 75\% of massive galaxies, and demonstrate a similar 70-75\% successful redshift measurement rate in $8.5 < \log(M_{\star}/M_{\odot}) < 9.5$ galaxies, suggesting similar survey sensitivity in this low-mass range. We explore the redshift, mass, and color-magnitude distributions of the catalog galaxies, finding that HALO7D galaxies are representative of CANDELS galaxies in these distributions up to \textit{i}-band magnitudes of 25. The catalogs presented will enable studies of star formation, the mass-metallicity relation, star formation and morphology relations, and other properties of the $z\sim0.7$ dwarf galaxy population.

Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments. The back-reaction of perturbations induced in the disk, both wave-like and non-wavelike, is a change in orbital elements of the planet. The efficiency of orbital migration is a long-standing problem in planet formation theory. We discuss recent progress in planet-disk interactions for different planet masses and disk parameters, in particular the level of turbulence, and progress in modeling observational signatures of embedded planets.

In classical cosmological analysis of large scale structure surveys with 2-pt functions, the parameter measurement precision is limited by several key degeneracies within the cosmology and astrophysics sectors. For cosmic shear, clustering amplitude $\sigma_8$ and matter density $\Omega_m$ roughly follow the $S_8=\sigma_8(\Omega_m/0.3)^{0.5}$ relation. In turn, $S_8$ is highly correlated with the intrinsic galaxy alignment amplitude $A_{\rm{IA}}$. For galaxy clustering, the bias $b_g$ is degenerate with both $\sigma_8$ and $\Omega_m$, as well as the stochasticity $r_g$. Moreover, the redshift evolution of IA and bias can cause further parameter confusion. A tomographic 2-pt probe combination can partially lift these degeneracies. In this work we demonstrate that a deep learning analysis of combined probes of weak gravitational lensing and galaxy clustering, which we call DeepLSS, can effectively break these degeneracies and yield significantly more precise constraints on $\sigma_8$, $\Omega_m$, $A_{\rm{IA}}$, $b_g$, $r_g$, and IA redshift evolution parameter $\eta_{\rm{IA}}$. The most significant gains are in the IA sector: the precision of $A_{\rm{IA}}$ is increased by approximately 8x and is almost perfectly decorrelated from $S_8$. Galaxy bias $b_g$ is improved by 1.5x, stochasticity $r_g$ by 3x, and the redshift evolution $\eta_{\rm{IA}}$ and $\eta_b$ by 1.6x. Breaking these degeneracies leads to a significant gain in constraining power for $\sigma_8$ and $\Omega_m$, with the figure of merit improved by 15x. We give an intuitive explanation for the origin of this information gain using sensitivity maps. These results indicate that the fully numerical, map-based forward modeling approach to cosmological inference with machine learning may play an important role in upcoming LSS surveys. We discuss perspectives and challenges in its practical deployment for a full survey analysis.

The primordial He abundance $Y_P$ is a powerful probe of cosmology. Currently, $Y_P$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($<0.1 Z_\odot$) galaxies (EMPGs) having reliable He/H measurements with HeI$\lambda$10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy and He/H determinations for 10 EMPGs, combining the existing optical data and the Markov chain Monte Carlo algorithm. Adding the existing 3 EMPGs and 51 moderately metal-poor ($0.1-0.4 Z_\odot$) galaxies with reliable He/H estimates, we obtain $Y_P=0.2379^{+0.0031}_{-0.0030}$ by linear regression in the (He/H)-(O/H) plane, where our observations increase the number of EMPGs from 3 to 13 anchoring He/H of the most metal-poor gas in galaxies. Although our $Y_P$ measurement and previous measurements are consistent, our result is slightly (~ 1$\sigma$) smaller due to our EMPGs. Including the existing primordial deuterium $D_P$ constraints, we estimate the effective number of neutrino species to be $N_{eff}=2.41^{+0.19}_{-0.21}$ showing a > 2 $\sigma$ tension with the Standard Model value ($N_{eff}=3.046$), which may be a hint of an asymmetry in electron-neutrino $\nu_e$ and anti-electron neutrino $\bar{\nu}_e$. Allowing the degeneracy parameter of electron-neutrino $\xi_e$ to vary as well as $N_{eff}$ and the baryon-to-photon ratio $\eta$, we obtain $\xi_e$ = $0.05^{+0.03}_{-0.03}$, $N_{eff}=3.22^{+0.33}_{-0.30}$, and $\eta\times10^{10}=6.13^{+0.04}_{-0.04}$ from the $Y_P$ and $D_P$ measurements with a prior of $\eta$ taken from Planck Collaboration et al. (2020). Our constraints suggest a $\nu_e - \bar{\nu}_e$ asymmetry and allow for a high value of $N_{eff}$ within the 1$\sigma$ level, which could mitigate the Hubble tension.

WASP-33 is one of the few $\delta$ Sct stars with a known planetary companion. By analyzing the stellar oscillations, we search for possible star-planet interactions in the pattern of the pulsation. We made use of the Transit and Light Curve Modeller (TLCM) to solve the light curve from the Transiting Exoplanet Survey Satellite (TESS). We include gravity darkening into our analysis. The stellar oscillation pattern of WASP-33 clearly shows signs of several tidally perturbed modes. We find that there are peaks in the frequency spectrum that are at or near the $3$rd, $12$th and $25$th orbital harmonics ($f_{orb} \sim 0.82$ d$^{-1}$). Also, there is a prominent overabundance of pulsational frequencies rightwards of the orbital harmonics, being characteristic of a tidally perturbed stellar pulsation, which is an outcome of star-planet interactions in the misaligned system. There are peaks in both the $\delta$ Sct and $\gamma$ Dor ranges of the Fourier spectrum, implying that WASP-33 is a $\gamma$ Dor -- $\delta$ Sct hybrid pulsator. The transit light curves are best fitted by a gravity darkened stellar model, and the planet parameters are consistent with earlier determinations.

We present a study of 323 photometrically variable young stellar objects that are likely members of the North America and Pelican (NAP) nebulae star forming region. To do so, we utilize over two years of data in the $g$ and $r$ photometric bands from the Zwicky Transient Facility (ZTF). We first investigate periodic variability, finding 46 objects ($\sim$15\% of the sample) with significant periods that phase well, and can be attributed to stellar rotation. We then use the quasi-periodicity (Q) and flux asymmetry (M) variability metrics to assign morphological classifications to the remaining aperiodic light curves. Another $\sim$39\% of the variable star sample beyond the periodic sources are also flux-symmetric, but with a quasi-periodic (moderate $Q$) or stochastic (high $Q$) nature. Concerning flux-asymmetric sources, our analysis reveals $\sim$14\% bursters (high negative $M$) and $\sim$29\% dippers (high positive $M$). We also investigate the relationship between variability slopes in the $g$ vs $g-r$ color-magnitude diagram, and the light curve morphological classes. Burster-type objects have shallow slopes, while dipper-type variables tend to have higher slopes that are consistent with extinction driven variability. Our work is one of the earliest applications of the $Q$ and $M$ metrics to ground-based data. We therefore contrast the $Q$ values of high-cadence and high-precision space-based data, for which these metrics were designed, with $Q$ determinations resulting from degraded space-based lightcurves that have the cadence and photometric precision characteristic of ground-based data.

Close encounters between two bodies in a disk often result in a single orbital deflection. However, within their Jacobi volumes, where the gravitational forces between the two bodies and the central body become competitive, temporary captures with multiple close encounters become possible outcomes: a Jacobi capture. We perform 3-body simulations in order to characterise the dynamics of Jacobi captures in the plane. We find that the phase space structure resembles a Cantor set with a fractal dimension of 0.4. The lifetime distribution decreases exponentially, while the distribution of the closest separation follows a power law with index 0.5. In our first application, we consider the Jacobi capture of the Moon. We demonstrate that both tidal captures and giant impacts are possible outcomes. Their respective 1D cross sections differ within an order of magnitude, evaluated at a heliocentric distance of 1 AU. The impact speed is well approximated by a parabolic encounter, while the impact angles follow that of a uniform beam on a circular target. In our second application, we find that Jacobi captures with gravitational wave dissipation can result in the formation of binary black holes in galactic nuclei. The eccentricity distribution is approximately super-thermal and includes both prograde and retrograde orientations. We estimate a cosmic rate density of 0.083 < R < 14 Gpc^-3 yr^-1. We conclude that dissipative Jacobi captures form an efficient channel for binary formation, which motivates further research into establishing the universality of Jacobi captures across multiple astrophysical scales.

The goal of the search for extraterrestrial intelligence (SETI) is the detection of non-human technosignatures, such as technology-produced emission in radio observations. While many have speculated about the character of such technosignatures, radio SETI fundamentally involves searching for signals that not only have never been detected, but also have a vast range of potential morphologies. Given that we have not yet detected a radio SETI signal, we must make assumptions about their form to develop search algorithms. The lack of positive detections also makes it difficult to test these algorithms' inherent efficacy. To address these challenges, we present Setigen, a Python-based, open-source library for heuristic-based signal synthesis and injection for both spectrograms (dynamic spectra) and raw voltage data. Setigen facilitates the production of synthetic radio observations, interfaces with standard data products used extensively by the Breakthrough Listen project (BL), and focuses on providing a physically-motivated synthesis framework compatible with real observational data and associated search methods. We discuss the core routines of Setigen and present existing and future use cases in the development and evaluation of SETI search algorithms.

Context. There is a population of runaway stars that move at extremely high speeds with respect to their surroundings. The fast motion and the stellar wind of these stars, plus the wind-medium interaction, can lead to particle acceleration and non-thermal radiation. Aims. We characterise the interaction between the winds of fast runaway stars and their environment, in particular to establish their potential as cosmic-ray accelerators and non-thermal emitters. Methods. We model the hydrodynamics of the interaction between the stellar wind and the surrounding material. We self-consistently calculate the injection and transport of relativistic particles in the bow shock using a multi-zone code, and compute their broadband emission from radio to $\gamma$-rays. Results. Both the forward and reverse shocks are favourable sites for particle acceleration, although the radiative efficiency of particles is low and therefore the expected fluxes are in general rather faint. Conclusions. We show that high-sensitivity observations in the radio band can be used to detect the non-thermal radiation associated with bow shocks from hypervelocity and semi-relativistic stars. Hypervelocity stars are expected to be modest sources of sub-TeV cosmic rays, accounting perhaps for a $\sim 0.1$% of that of galactic cosmic rays.

Motivated by the result of Planck+BICEP/Keck recently released, we investigate the consistency of the multi-field inflation models in terms of the spectral index $n_s$ and the tensor-to-scalar ratio $r$. In this study, we focus on double-inflaton models with and without a spectator field. We find that inflaton with a quadratic potential can become viable when three fields with a specific hierarchical mass spectrum are realized such that two fields act as inflatons and the other one is the spectator. We also discuss the conditions to avoid the fine-tuning, by careful study of how the prediction depends on the background trajectory in the inflaton-field space.

The plane-of-the-sky component of interstellar magnetic fields can be traced in two dimensions using polarized dust emission. Its potential to access three-dimensional magnetic fields, including the inclination of the magnetic fields relative to the line-of-sight, is crucial for a variety of astrophysical problems. Based on the statistical features of the polarization fraction and the averaged POS Alfv\'en Mach number $\overline{M_{\rm A}}_{,\bot}$, we present a new method for estimating the inclination angle. The magnetic field fluctuations raised by MHD turbulence are taken into account in our method. By using synthetic dust emission from 3D MHD turbulence simulation, we show that the fluctuations are preferentially perpendicular to the mean magnetic field and magnify the depolarization effect. We analytically and numerically derive that the polarization fraction is characterized by both the inclination angle and magnetic field fluctuations. We propose and demonstrate that the mean inclination angle over a region of interest can be calculated from the polarization fraction in a strongly magnetized reference position, where $\overline{M_{\rm A}}_{,\bot}^2\ll1$. We test the new method in sub-Alfv\'enic, trans-Alfv\'enic, and moderately super-Alfv\'enic situations ($0.4\lesssim M_{\rm A}\lesssim1.2$) and show that it recovers the mean inclination angle well. This method is extended to determine the inclination angle distribution over the sub-regions of a cloud.

Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constraints on the planet formation process. In this chapter, we summarize the new information derived from the exoplanets population and the circumstellar disks observations. We present the new developments in planet formation theory, from dust evolution to the growth of planetary cores by accretion of planetesimals, pebbles, and gas. We review the state-of-the-art models for the formation of diverse planetary systems, including the population synthesis approach which is necessary to compare theoretical model outcomes to the exoplanet population. We emphasize that the planet formation process may not be spatially uniform in the disk and there are preferential locations for the formation of planetesimals and planets. Outside of these locations, a significant fraction of solids is not growing past the pebble-sizes. The reservoir of pebbles plays an important role in the growth of planetary cores in the pebble accretion process. The timescale of the emergence of massive planetary cores is an important aspect of the present models and it is likely that the cores within one disk form at different times. In addition, there is growing evidence that the first planetary cores start forming early, during the circumstellar disk buildup process.

The ultrarelativistic jets triggered by neutrino annihilation processes or Blandford-Znajek (BZ) mechanisms in stellar-mass black hole (BH) hyperaccretion systems are generally considered to power gamma-ray bursts (GRBs). Due to the high accretion rate, the central BHs might grow rapidly on a short timescale, providing a new way to understand "the lower mass gap" problem. In this paper, we use the BH hyperaccretion model to investigate BH mass growth based on observational GRB data. The results show that (i) if the initial BH mass is set as $3~M_\odot$, the neutrino annihilation processes are capable of fueling the BHs to escape the lower mass gap for more than half of long-duration GRBs (LGRBs), while the BZ mechanism is inefficient on triggering BH growths for LGRBs; (ii) the mean BH mass growths in the case of LGRBs without observable supernova (SN) association are much larger than these in the case of LGRBs associated with SNe for both mechanisms, which imply that more massive progenitors or lower SN explosion energies prevail throughout the former cases; (iii) for the short-duration GRBs, the mean BH mass growths are satisfied with the mass supply limitation in the scenario of compact object mergers, but the hyperaccretion processes are unable to rescue BHs from the gap in binary neutron star (NS) mergers or the initial BH mass being $3~M_\odot$ after NS-BH mergers.

We present 21 new long-term variable radio sources found commensally in two years of weekly MeerKAT monitoring of the low-mass X-ray binary GX 339-4. The new sources vary on time scales of weeks to months and have a variety of light curve shapes and spectral index properties. Three of the new variable sources are coincident with multi-wavelength counterparts; and one of these is coincident with an optical source in deep MeerLICHT images. For most sources, we cannot eliminate refractive scintillation of active galactic nuclei as the cause of the variability. These new variable sources represent $2.2\pm0.5$ per cent of the unresolved sources in the field, which is consistent with the 1-2 per cent variability found in past radio variability surveys. However, we expect to find short-term variable sources in the field as well as these 21 new long-term variable sources. We present the radio light curves and spectral index variability of the new variable sources, as well as the absolute astrometry and matches to coincident sources at other wavelengths.

The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks, as well as their spatial distribution, to better understand planet formation. Since PPVI, the field has seen tremendous improvements in observational capabilities, enabling both surveys of large samples of disks and high resolution imaging studies of a few bright disks. Improvements in data quality and sample size have, however, opened up many fundamental questions about properties such as the mass budget of disks, its spatial distribution, and its radial extent. Moreover, the vertical structure of disks has been studied in greater detail with spatially resolved observations, providing new insights on vertical layering and temperature stratification, yet also bringing rise to questions about other properties, such as material transport and viscosity. Each one of these properties - disk mass, surface density distribution, outer radius, vertical extent, temperature structure, and transport - is of fundamental interest as they collectively set the stage for disk evolution and corresponding planet formation theories. In this chapter, we will review our understanding of the fundamental properties of disks including the relevant observational techniques to probe their nature, modelling methods, and the respective caveats. Finally, we discuss the implications for theories of disk evolution and planet formation underlining what new questions have since arisen as our observational facilities have improved.

OB-associations and superbubbles being energetically essential galactic powerhouses are likely to be the important acceleration sites of galactic cosmic rays (CRs). The emission profile of gamma-ray sources related to superbubbles and stellar clusters indicates on continuous particle acceleration by winds of massive stars. One of the most luminous galactic gamma-ray sources is Cygnus Cocoon superbubble, observed by multiple instruments, such as Fermi-LAT, ARGO, and, recently, HAWC. We discuss a model of particle acceleration and transport in a superbubble to explain GeV-TeV gamma-ray spectrum of Cygnus Cocoon, which has a break at the energy of about 1 TeV. It is shown that the gamma rays produced by hadronic interactions of high-energy protons accelerated by an ensemble of shocks from winds of massive stars and supernovae in the Cygnus Cocoon can explain the observations. The proton spectral shape at the highest energies depends on the MHD-fluctuation spectrum in the Cocoon. The viable solutions for Cygnus Cocoon may be applied to some other associations showing similar behaviour. We briefly discuss the similarity and differences of particle acceleration processes in extended superbubbles and compact clusters of young massive stars as represented by Westerlund 1 and 2 gamma-ray sources.

We study the accretion flow from a counter-rotating torus orbiting a central Kerr black hole (BH). We characterize the flow properties at the turning point of the accreting matter flow from the orbiting torus, defined by the condition $u^\phi = 0$ on the flow torodial velocity. The counter-rotating accretion flow and jet-like flow turning point location along BH rotational axis is given. Some properties of the counter-rotating flow thickness and counter-rotating tori energetics are studied.The maximum amount of matter swallowed by the BH from the counter-rotating tori is determined by the background properties. The fast spinning BH energetics depends mostly on BH spin rather than on the properties of the counter-rotating fluids or the tori masses. The turning point is located in a narrow orbital corona (spherical shell), for photons and matter flow constituents, surrounding the BH stationary limit (outer ergosurface), depending on the BH spin-mass ratio and the fluid initial momentum only. The turning corona for jet-like-flow has larger thickness, it is separated from the torus flow turning corona and it is closer to the BH stationary limit. Turning points of matter accreting from torus and from jets are independent explicitly of the details of the accretion and tori model. The turning corona could be observable due to an increase of flow luminosity and temperature. The corona is larger on the BH equatorial plane, where it is the farthest from the central attractor, and narrower on the BH poles.

The atmospheres of brown dwarfs orbiting in close proximity to their parent white dwarf represent some of the most extreme irradiated environments known. Understanding their complex dynamical mechanisms pushes the limits of theoretical and modelling efforts, making them valuable objets to study to test contemporary understanding of irradiated atmospheres. We use the Exo-FMS GCM to simulate the brown dwarfs WD0137-349B, SDSS J141126.20+200911.1B and EPIC212235321B, first coupled to a multi-banded grey radiative-transfer scheme then a spectral correlated-k scheme with high temperature opacity tables. We then post-process the GCM results using gCMCRT to compare to available observational data. Our GCM models predict strongly temperature inverted atmospheres, spanning many decades in pressure due to impact of UV band heating. Post-processing of our models suggest that the day-night contrast is too small in the GCM results. We therefore suggest that the formation of cloud particles as well as atmospheric drag effects such as magnetic drag are important considerations in setting the day-night temperature contrast for these objects.

Using data from the Pantheon SN Ia compilation and the Sloan Digital Sky Survey (SDSS), we propose an estimator for weak lensing convergence incorporating positional and photometric data of foreground galaxies. The correlation between this and the Hubble diagram residuals of the supernovae is $3.7\sigma$, consistent with weak lensing magnification due to dark matter halos centered on galaxies. We additionally constrain the properties of the galactic haloes, such as the mass-to-light ratio $\Gamma$ and radial profile of the halo matter density $\rho(r)$. We derive a new relationship for the additional r.m.s. scatter in magnitudes caused by lensing, finding $\sigma_{\rm lens} = (0.06 \pm 0.017) (d_{\rm M}(z)/ d_{\rm M}(z=1))^{3/2}$ where $d_{\rm M}(z)$ is the comoving distance to redshift $z$. We therefore find that the scatter in flux caused by lensing will be of a similar size as the intrinsic scatter (after standardisation of the magnitude) of modern SN Ia surveys by $z \sim 1.2$, which is closer than generally assumed. We propose a modification of the distance modulus estimator for SN Ia to incorporate lensing, which we anticipate will improve the accuracy of cosmological parameter estimation using high-redshift SN Ia data.

(abridged) Variability is a key property of stars on the asymptotic giant branch (AGB). Their pulsation period is related to the luminosity and mass-loss rate (MLR) of the star. Long-period variables (LPVs) and Mira variables are the most prominent of all types of variability of evolved stars. However, the reddest, most obscured AGB stars are too faint in the optical and have eluded large variability surveys. Our goal is to obtain a sample of LPVs with large MLRs by analysing WISE W1 and W2 light curves (LCs) for about 2000 sources, photometrically selected to include known C-stars with the 11.3 mu silicon carbide dust feature in absorption, and Galactic O-stars with periods longer than 1000 days. Epoch photometry was retrieved from the AllWISE and NEOWISE database and fitted with a sinus curve. Photometry from other variability surveys was also downloaded and fitted. For a subset of 316 of the reddest stars, spectral energy distributions (SEDs) were constructed, and, together with mid-infrared (MIR) spectra when available, fitted with a dust radiative transfer programme in order to derive MLRs. WISE based LCs and fits to the data are presented for all stars. Inspired by a recent paper, a number of non-variable OH/IRs are identified. Based on a selection on amplitude, a sample of about 750 (candidate) LPVs is selected of which 145 have periods over 1000 days, many of them being new. For the subset of the stars with the colours of C-rich extremely red objects (EROs) the fitting of the SEDs (and available MIR spectra) separates them into C- and O-rich objects. The number of Galactic EROs appears to be complete up to about 5~kpc and a total dust return rate in the solar neighbourhood for this class is determined. Based on the EROs in the Magellanic Clouds, a bolometric period luminosity is derived.

In 2010 Jewitt and Li published a paper examining the behavior of comet-asteroid transition object 3200 Phaethon, arguing it was asteroid-like in its behavior throughout most of its orbit, but that near its perihelion, at a distance of only 0.165 AU from the sun, its dayside temperatures would be hot enough to vaporize rock (>1000 K, Hanus et al. 2016). Thus it would act like a "rock comet" as gases produced from evaporating rock were released from the body, in a manner similar to the more familiar sublimation of water ice into vacuum seen for comets coming within ~3 AU of the Sun. In this Note we predict that the same thermal effects that would create "rock comet" behavior with Qgas ~ 10$^{22}$ mol/sec at perihelion would also help greatly bluen Phaethon's surface via preferential thermal alteration and sublimative removal of surface Fe and refractory organics, known reddening and darkening agents. These predictions are testable by searching for signs of spectral bluening of the surfaces of other objects in Phaethon-like small perihelion orbits, and by in situ measurements of Phaethons surface and coma composition near perihelion with the upcoming DESTINY+ mission to Phaethon by JAXA.

The ability to discover new transients via image differencing without direct human intervention is an important task in observational astronomy. For these kind of image classification problems, machine Learning techniques such as Convolutional Neural Networks (CNNs) have shown remarkable success. In this work, we present the results of an automated transient identification on images with CNNs for an extant dataset from the Dark Energy Survey Supernova program (DES-SN), whose main focus was on using Type Ia supernovae for cosmology. By performing an architecture search of CNNs, we identify networks that efficiently select non-artifacts (e.g. supernovae, variable stars, AGN, etc.) from artifacts (image defects, mis-subtractions, etc.), achieving the efficiency of previous work performed with random Forests, without the need to expend any effort in feature identification. The CNNs also help us identify a subset of mislabeled images. Performing a relabeling of the images in this subset, the resulting classification with CNNs is significantly better than previous results.

Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys of multiple star-forming regions cover both spectroscopy of young stars and mm interferometry of their protoplanetary disks. This provides an unprecedented statistical sample of stellar masses and mass accretion rates, as well as disk masses and radii, for almost 1000 young stellar objects within 300 pc from us, while also sampling different evolutionary stages, ages, and environments. At the same time, surveys of debris disks are revealing the bulk properties of this class of more evolved objects. This chapter reviews the statistics of these measured global star and disk properties and discusses their constraints on theoretical models describing global disk evolution. Our comparisons of observations to theoretical model predictions extends beyond the traditional viscous evolution framework to include analytical descriptions of magnetic wind effects. Finally, we discuss how recent observational results can provide a framework for models of planet population synthesis and planet formation.

We present broadband UV/X-ray spectral variability of the changing-look active galactic nucleus NGC 1566 based on simultaneous near-ultraviolet (NUV) and X-ray observations performed by XMM-Newton, Swift, and NuSTAR satellites at five different epochs during the declining phase of the 2018 outburst. We found that the accretion disk, soft X-ray excess, and the X-ray power-law components were extremely variable. Additionally, the X-ray power-law flux was correlated with both the soft excess plus disk and the pure disk fluxes. Our finding shows that at high flux levels the soft X-ray excess and the disk emission both provided the seed photons for thermal Comptonization in the hot corona, whereas at low flux levels where the soft excess was absent, the pure disk emission alone provided the seed photons. The X-ray power-law photon-index was only weakly variable ($\Delta{\Gamma_{hot}}\leq0.06$) and it was not well correlated with the X-ray flux over the declining timescale. On the other hand, we found that the electron temperature of the corona increased from $\sim22$ to $\sim200$ keV with decreasing number of seed photons from June 2018 to August 2019. At the same time, the optical depth of the corona decreased from $\tau_{hot}\sim4$ to $\sim0.7$, and the scattering fraction increased from $\sim1\%$ to $\sim10\%$. These changes suggest structural changes in the hot corona such that it grew in size and became hotter with decreasing accretion rate during the declining phase. The AGN is most likely evolving with decreasing accretion rate towards a state similar to the low/hard state of black hole X-ray binaries.

In this fourth paper of a series on the precision of ages of stellar populations obtained through the full-spectrum fitting technique, we present a first systematic analysis that compare the age, metallicity and reddening of star clusters obtained from resolved and unresolved data (namely colour-magnitude diagrams (CMDs) and integrated-light spectroscopy) using the same sets of isochrones. We investigate the results obtained with both Padova isochrones and MIST isochrones. We find that there generally is a good agreement between the ages derived from CMDs and integrated spectra. However, for metallicity and reddening, the agreement between results from analyses of CMD and integrated spectra is significantly worse. Our results also show that the ages derived with Padova isochrones match those derived using MIST isochrones, both with the full spectrum fitting technique and the CMD fitting method. However, the metallicity derived using Padova isochrones does not match that derived using MIST isochrones using the CMD method. We examine the ability of the full-spectrum fitting technique in detecting age spreads in clusters that feature the extended Main Sequence Turn Off (eMSTO) phenomenon using two-population fits. We find that 3 out of 5 eMSTO clusters in our sample are best fit with one single age, suggesting that eMSTOs do not necessarily translate to detectable age spreads in integrated-light studies.

In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spirals, arcs, and shadows, are seen in all imaged infrared-bright disks to date, and are ubiquitous across spectral types. These substructures are believed to be the result of disk evolution processes, and in particular disk-planet interactions. Since small dust grains that scatter light are tightly bound to the disk's gas, these observations closely trace disk structures predicted by hydrodynamical models and serve as observational tests of the predictions of planet formation theories. We argue that the results of current and next-generation high-contrast imaging surveys will, when combined with complementary data from ALMA, lead to a much deeper understanding of the co-evolution of disks and planets, and the mechanisms by which planets form.

OB associations are low-density groups of young stars that are dispersing from their birth environment into the Galactic field. They are important for understanding the star formation process, early stellar evolution, the properties and distribution of young stars and the processes by which young stellar groups disperse. Recent observations, particularly from Gaia, have shown that associations are highly complex, with a high degree of spatial, kinematic and temporal substructure. The kinematics of associations have shown them to be globally unbound and expanding, with the majority of recent studies revealing evidence for clear expansion patterns in the association subgroups, suggesting the subgroups were more compact in the past. This expansion is often non-isotropic, arguing against a simple explosive expansion, as predicted by some models of residual gas expulsion. The star formation histories of associations are often complex, exhibit moderate age spreads and temporal substructure, but so far have failed to reveal simple patterns of star formation propagation (e.g., triggering). These results have challenged the historical paradigm of the origin of associations as the expanded remnants of dense star clusters and suggests instead that they originate as highly substructured systems without a linear star formation history, but with multiple clumps of stars that have since expanded and begun to overlap, producing the complex systems we observe today. This has wide-ranging consequences for the early formation environments of most stars and planetary systems, including our own Solar System.

Progressive astronomical characterization of planet-forming disks and rocky exoplanets highlight the need for increasing interdisciplinary efforts to understand the birth and life cycle of terrestrial worlds in a unified picture. Here, we review major geophysical and geochemical processes that shape the evolution of rocky planets and their precursor planetesimals during planetary formation and early evolution, and how these map onto the astrophysical timeline and varying accretion environments of planetary growth. The evolution of the coupled core-mantle-atmosphere system of growing protoplanets diverges in thermal, compositional, and structural states to first order, and ultimately shapes key planetary characteristics that can discern planets harboring clement surface conditions from those that do not. Astronomical campaigns seeking to investigate rocky exoplanets will require significant advances in laboratory characterization of planetary materials and time- and spatially-resolved theoretical models of planetary evolution, to extend planetary science beyond the Solar System and constrain the origins and frequency of habitable worlds like our own.

The optical module of the KM3NeT neutrino telescope is an innovative, multi-faceted large area photodetection module. It contains 31 three-inch photomultiplier tubes in a single 0.44 m diameter pressure-resistant glass sphere. The module is a sensory device also comprising calibration instruments and electronics for power, readout and data acquisition. It is capped with a breakout-box with electronics for connection to an electro-optical cable for power and long-distance communication to the onshore control station. The design of the module was qualified for the first time in the deep sea in 2013. Since then, the technology has been further improved to meet requirements of scalability, cost-effectiveness and high reliability. The module features a sub-nanosecond timing accuracy and a dynamic range allowing the measurement of a single photon up to a cascade of thousands of photons, suited for the measurement of the Cherenkov radiation induced in water by secondary particles from interactions of neutrinos with energies in the range of GeV to PeV. A distributed production model has been implemented for the delivery of more than 6000 modules in the coming few years with an average production rate of more than 100 modules per month. In this paper a review is presented of the design of the multi-PMT KM3NeT optical module with a proven effective background suppression and signal recognition and sensitivity to the incoming direction of photons.

In the next decade, two rovers will characterize in situ the mineralogy of rocks on Mars, using for the first time near-infrared reflectance spectrometers: SuperCam onboard the Mars 2020 rover and MicrOmega onboard the ExoMars rover, although this technique is predominantly used in orbit for mineralogical investigations. Until successful completion of sample-return missions from Mars, Martian meteorites are currently the only samples of the red planet available for study in terrestrial laboratories and comparison with in situ data. However, the current spectral database available for these samples does not represent their diversity and consists primarily of spectra acquired on finely crushed samples, albeit grain size is known to greatly affect spectral features. We measured the reflected light of a broad Martian meteorite suite as a means to catalogue and characterize their spectra between 0.4 and 3 microns. These measurements are achieved using a point spectrometer acquiring data comparable to SuperCam, and an imaging spectrometer producing hyperspectral cubes similarly to MicrOmega. Our results indicate that point spectrometry is sufficient to discriminate the different Martian meteorites families, to identify their primary petrology based on band parameters, and to detect their low content in alteration minerals. However, significant spectral mixing occurs in the point measurements, even at spot sizes down to a few millimeters, and imaging spectroscopy is needed to correctly identify the various mineral phases in the meteorites. Bidirectional spectral measurements confirm their non-Lambertian behavior, with backward and suspected forward scattering peaks. With changing observation geometry, the main absorption strengths show variations up to 10-15 percents. All the spectra presented are provided in the supplementary data for further comparison with in situ and orbital measurements.

Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds.

Most stars are born in multiple stellar systems. Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the population statistics spanning system evolution from the protostellar phase through the main-sequence phase and evaluate the influence of the local environment. We describe current models for the origin of stellar multiplicity and review the landscape of numerical simulations and assess their consistency with observations. We review the properties of disks and discuss the impact of multiplicity on planet formation and system architectures. Finally, we summarize open questions and discuss the technical requirements for future observational and theoretical progress.

The assembly and architecture of planetary systems strongly depend on the physical processes governing the evolution and dispersal of protoplanetary disks. Since Protostars and Planets VI, new observations and theoretical insights favor disk winds as being one of those key processes. This chapter provides a comprehensive review of recent observations probing outflowing gas launched over a range of disk radii for a wide range of evolutionary stages, enabling an empirical understanding of how winds evolve. In parallel, we review theoretical advancements in both magnetohydrodynamic and photoevaporative disk wind models and identify predictions that can be confronted with observations. By linking theory and observations we critically assess the role of disk winds in the evolution and dispersal of protoplanetary disks. Finally, we explore the impact of disk winds on planet formation and evolution and highlight theoretical work, observations, and critical tests for future progress.

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

During the third observing run (O3) of the Advanced LIGO and Advanced Virgo detectors, dozens of candidate gravitational-wave (GW) events have been catalogued. A challenge of this observing run has been the rapid identification and public dissemination of compact binary coalescence (CBC) signals, a task carried out by low-latency searches such as PyCBC Live. During the later part of O3, we developed a method of classifying CBC sources, via their probabilities of containing neutron star or black hole components, within PyCBC Live, in order to facilitate immediate follow-up observations by electromagnetic and neutrino observatories. This fast classification uses the chirp mass recovered by the search as input, given the difficulty of measuring the mass ratio with high accuracy for lower-mass binaries. We also use a distance estimate derived from the search output to correct for the bias in chirp mass due to the cosmological redshift. We present results for simulated signals, and for confirmed candidate events identified in low latency over O3.

Hydrocarbon species, and in particular CH$_4$, play a key role in the stratosphere--thermosphere boundary of Jupiter, which occurs around the $\mu$-bar pressure level. Previous analyses of solar occultation, He and Ly-$\alpha$ airglow, and ISO/SWS measurements of the radiance around 3.3 $\mu$m have inferred significantly different methane concentrations. Here we aim to accurately model the CH$_4$ radiance at 3.3 $\mu$m measured by ISO/SWS by using a comprehensive non-local thermodynamic equilibrium model and the most recent collisional rates measured in the laboratory for CH$_4$ to shed new light onto the methane concentration in the upper atmosphere of Jupiter. These emission bands have been shown to present a peak contribution precisely at the $\mu$-bar level, hence directly probing the region of interest. We find that a high CH$_4$ concentration is necessary to explain the data, in contrast with the most recent analyses, and that the observations favour the lower limit of the latest laboratory measurements of the CH$_4$ collisional relaxation rates. Our results provide precise constraints on the composition and dynamics of the lower atmosphere of Jupiter.

We study the production and backreaction of massive vector-like fermions in the background of a classical $SU(2)$ gauge field during axion-driven inflation. We demonstrate all ultraviolet divergences due to the interactions with the fermions can be absorbed by renormalization of the axion wavefunction and the gauge coupling. The effects of the fermion-axion interaction vanish in the massless limit. For very massive fermions, contact interactions are induced between the axion, the gauge field and the gravitational field. In this massive limit, we find the usual axion-gauge field interactions are induced, however, in addition we observe the appearance of axion self-interactions, as well as kinetic braiding of the axion with the Einstein tensor. These new axion derivative interactions present intriguing opportunities for model building and phenomenology.

The presence of a dissipative dark matter (DM) sector may allow for the trapping of a significant DM mass inside stars, either during structure formation or by accretion over their lifetime, influencing stellar behavior well into the Main Sequence stage. Motivated by this scenario, we place an upper bound on the fractional DM mass within current-day Main Sequence stars. Using double-lined spectroscopic binaries (SB2 stars), gravitational masses are extracted and contrasted with luminous masses, derived using a modified mass-luminosity relation which accounts for the effect of DM. High-accuracy mass and luminosity data from a sample of 486 binary stars in the $0.18 < M/M_\odot < 31$ mass range are partitioned into appropriate mass domains and analyzed. A 95% C.L. upper limit of sub-5% is found for the subset of 263 stars in the $1 < M/M_\odot < 2.4$ regime.

For inflation driven by the Higgs field coupled non-minimally to gravity, we study the cutoff energy scale above which perturbation theory breaks down. Employing the metric formulation, we first give an overview of known results and then provide a new way to calculate a lower bound on the cutoff. Our approach neither relies on a gauge choice nor does it require any calculation of amplitudes. Instead, it exploits the fact that the S-matrix is invariant under field redefinitions. In agreement with previous findings, we demonstrate that the cutoff is significantly higher during inflation than in vacuum, which ensures the robustness of semi-classical predictions. Along the way, we generalize our findings to the Palatini formulation and comment on a useful parametrization of the Higgs doublet in both scenarios.

KAGRA uses cryogenics to cool its sapphire test masses down to 20 K to reduce the thermal noise. However, cryocooler vibration and structural resonances of the cryostat couple to test mass and can contaminate the detector sensitivity. We performed vibration analysis of the cooling system at cryogenic temperature to study its impact on detector sensitivity. Our measurement show shield vibration below 1 Hz is not impacted by cryocooler operation or structural resonances and follows ground motion. The noise floor of the shield in 1-100 Hz was observed to be 2-3 order of magnitude larger than seismic motion even without cryocooler operation. The operation of cryocoolers does not change the noise floor, but 2.0 Hz peaks and their harmonics were observed over the entire spectrum (1-100 Hz). These results were used to calculate the coupling of cooling system vibration to the test mass. We conclude that vibration from the cooling system does not limit KAGRA design sensitivity.

The recent observation of $^4$He implies that our universe has a large lepton asymmetry. We consider the Affleck-Dine (AD) mechanism for lepton number generation. In the AD mechanism, non-topological solitons called L-balls are produced, and the generated lepton number is confined in them. The L-balls protect the generated lepton number from being converted to baryon number through the sphaleron processes. We study the formation and evolution of the L-balls and find that the universe with large lepton asymmetry suggested by the recent $^4$He measurement can be realized.

In this work we shall use a bottom-up approach for obtaining viable inflationary Einstein-Gauss-Bonnet models which are also compatible with the GW170817 event. Specifically, we shall use a recently developed theoretical framework in which we shall specify only the tensor-to-scalar ratio, in terms of the $e$-foldings number. Starting from the tensor-to-scalar ratio, we shall reconstruct from it the Einstein-Gauss-Bonnet theory which can yield such a tensor-to-scalar ratio, finding the scalar potential and the Gauss-Bonnet coupling scalar function as functions of the $e$-foldings number. Accordingly, the calculation of the spectral index of the primordial scalar perturbations, and of the tensor spectral index easily is greatly simplified and these observational indices can easily be found. After presenting the general formalism for the bottom-up reconstruction, we exemplify our findings by presenting several Einstein-Gauss-Bonnet models of interest which yield a viable inflationary phenomenology. These models have also an interesting common characteristic, which is a blue tilted tensor spectral index. We also investigate the predicted energy spectrum of the primordial gravitational waves for these Einstein-Gauss-Bonnet models, and as we show, all the models yield a detectable primordial wave energy power spectrum.

This article is an addendum to [1]. We extend our computation of the on-shell scattering amplitudes in an arbitrary inflaton background $\bar{\phi}_1$. Although the effective Einstein frame cutoff for $\bar{\phi}_1>>M_P/\sqrt{\xi}$ turns out to be $\bar{\phi}_1$ or $\xi\bar{\phi}_1$ for the $U(1)$ model, this is not the case for the realistic doublet Higgs model where the effective Einstein frame cutoff turns out to be the standard $M_P/\sqrt{\xi}$ for both the Palatini and metric formulations. Then, as it has been pointed out in [1] the background $\bar{\phi}_1$ is the effective Jordan frame cutoff for both the Palatini and metric formulations.

The weakly interacting massive particles (WIMPs) have been the most popular particle dark matter (DM) candidate for the last several decades, and it is well known that WIMP can be probed via the direct, indirect and collider experiments. However, the direct and indirect signals are highly suppressed in some scalar-mediated DM models, e.g. the lepton portal model with a Majorana DM candidate. As a result, collider searches are considered as the only hope to probe such models. In this white paper, we propose that the gravitational wave (GW) astronomy also serves as a powerful tool to probe such scalar mediated WIMP models via the potential first-order phase transition GW signals. An example for the lepton portal dark matter is provided, showing the complementarity between collider and GW probes.

Interstellar communication transmitters, intended to be discovered and decoded to information bits, are expected to transmit signals that contain message symbols quantized in at least one of the degrees of freedom of the transmitted signal. A hypothesis is proposed that signal quantization, in the form of multiplicative values of one or more signal measurements, may be observable during the reception of hypothetical discoverable interstellar communication signals. In previous work, using single and multiple synchronized radio telescopes, candidate hypothetical interstellar communication signals comprising delta-t delta-f opposite circular polarized pulse pairs have been reported and analyzed (ref. arXiv:2105.03727, arXiv:2106.10168, arXiv:2202.12791). In the latter report, an apparent quantization of delta-f at multiples of 58.575 Hz was observed. In the current work, a machine process has been implemented to further examine anomalous delta-f and delta-t quantization, with results reported in this paper. As in some past work, a 26 foot diameter radio telescope with fixed azimuth and elevation pointing is used to enable a Right Ascension filter to measure signals associated with a celestial direction of interest, relative to other directions, over a 6.3 hour range of Right Ascension. The 5.25 plus or minus 0.15 hour Right Ascension, -7.6 degrees plus or minus 1 degree Declination celestial direction presents repetition and quantization anomalies, during an experiment lasting 157 days, with the first 143 days overlapping the previous experiment.