Papers for Monday, Apr 03 2023

`ronswanson` provides a simple-to-use framework for building so-called table or template models for `astromodels`the modeling package for multi-messenger astrophysical data-analysis framework, `3ML`. With `astromodels` and `3ML` one can build the interpolation table of a physical model result of an expensive computer simulation. This then enables efficient reevaluation of the model while, for example, fitting it to a dataset. While `3ML` and `astromodels` provide factories for building table models, the construction of pipelines for models that must be run on high-performance computing (HPC) systems can be cumbersome. `ronswanson` removes this complexity with a simple, reproducible templating system. Users can easily prototype their pipeline on multi-core workstations and then switch to a multi-node HPC system. `ronswanson` automatically generates the required `Python` and `SLURM` scripts to scale the execution of `3ML` with `astromodel`'s table models on an HPC system.

There should be about 10,000 stellar hierarchical systems within 100 pc with primary stars more massive than 0.5 Msun, and a similar amount of less massive hierarchies. A list of 8000 candidate multiples is derived from wide binaries found in the Gaia Catalog of Nearby Stars where one or both components have excessive astrometric noise or other indicators of inner subsystems. A subset of 1243 southern candidates were observed with high angular resolution at the 4.1 m telescope, and 503 new pairs with separations from 0.03" to 1" were resolved. These data allow estimation of the inner mass ratios and periods and help to quantify the ability of Gaia to detect close pairs. Another 621 hierarchies with known inner periods come from the Gaia catalog of astrometric and spectroscopic orbits. These two non-overlapping groups, combined with existing ground-based data, bring the total number of known nearby hierarchies to 2754, reaching a completeness of ~22% for stars above 0.5 Msun. Distributions of their periods and mass ratios are briefly discussed, and the prospects of further observations are outlined.

The newfound ability to detect SO2 in exoplanet atmospheres presents an opportunity to measure sulfur abundances and so directly test between competing modes of planet formation. In contrast to carbon and oxygen, whose dominant molecules are frequently observed, sulfur is much less volatile and resides almost exclusively in solid form in protoplanetary disks. This dichotomy leads different models of planet formation to predict different compositions of gas giant planets. Whereas planetesimal-based models predict roughly stellar C/S and O/S ratios, pebble accretion models more often predict superstellar ratios. To explore the detectability of SO2 in transmission spectra and its ability to diagnose planet formation, we present a grid of atmospheric photochemical models and corresponding synthetic spectra for WASP-39b (where SO2 has been detected). Our 3D grid contains 11^3 models (spanning 1--100x the solar abundance ratio of C, O, and S) for thermal profiles corresponding to the morning and evening terminators, as well as mean terminator transmission spectra. Our models show that for a WASP-39b-like O/H and C/H enhancement of ~10x Solar, SO2 can only be seen for C/S and O/S <~1.5, and that WASP-39b's reported SO2 abundance of 1--10 ppm may be more consistent with planetesimal accretion than with pebble accretion models (although some pebble models also manage to predict similarly low ratios). More extreme C/S and O/S ratios may be detectable in higher-metallicity atmospheres, suggesting that smaller and more metal-rich gas and ice giants may be particularly interesting targets for testing planet formation models. Future studies should explore the dependence of SO2 on a wider array of planetary and stellar parameters, both for the prototypical SO2 planet WASP-39b, as well as for other hot Jupiters and smaller gas giants.

We present the first results on the spatial distribution of dust attenuation at $1.0<z<2.4$ traced by the Balmer Decrement, H$\alpha$/H$\beta$, in emission-line galaxies using deep JWST NIRISS slitless spectroscopy from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). H$\alpha$ and H$\beta$ emission line maps of emission-line galaxies are extracted and stacked in bins of stellar mass for two grism redshift bins, $1.0<z_{grism}<1.7$ and $1.7<z_{grism}<2.4$. Surface brightness profiles for the Balmer Decrement are measured and radial profiles of the dust attenuation towards H$\alpha$, $A_{\mathrm{H}\alpha}$, are derived. In both redshift bins, the integrated Balmer Decrement increases with stellar mass. Lower mass ($7.6\leqslant$Log($M_{*}$/M$_{\odot}$)$<10.0$) galaxies have centrally concentrated, negative dust attenuation profiles whereas higher mass galaxies ($10.0\leqslant$Log($M_{*}$/M$_{\odot}$)$<11.1$) have flat dust attenuation profiles. The total dust obscuration is mild, with on average $0.07\pm0.07$ and $0.14\pm0.07$ mag in the low and high redshift bins respectively. We model the typical light profiles of star-forming galaxies at these redshifts and stellar masses with GALFIT and apply both uniform and radially varying dust attenuation corrections based on our integrated Balmer Decrements and radial dust attenuation profiles. If these galaxies were observed with typical JWST NIRSpec slit spectroscopy ($0.2\times0.5^{\prime\prime}$ shutters), on average, H$\alpha$ star formation rates (SFRs) measured after slit-loss corrections assuming uniform dust attenuation will overestimate the total SFR by $6\pm21 \%$ and $26\pm9 \%$ at $1.0\leqslant z < 1.7$ and $1.7\leqslant z < 2.4$ respectively.

Sunskirting asteroid (3200) Phaethon has been repeatedly observed in STEREO HI1 imagery to anomalously brighten and produce an antisunward tail for a few days near each perihelion passage, phenomena previously attributed to the ejection of micron-sized dust grains. Color imaging by the SOHO LASCO coronagraphs during the 2022 May apparition indicate that the observed brightening and tail development instead capture the release of sodium atoms, which resonantly fluoresce at the 589.0/589.6 nm D lines. While HI1's design bandpass nominally excludes the D lines, filter degradation has substantially increased its D line sensitivity, as quantified by the brightness of Mercury's sodium tail in HI1 imagery. Furthermore, the expected fluorescence efficiency and acceleration of sodium atoms under solar radiation readily reproduce both the photometric and morphological behaviors observed by LASCO and HI1 during the 2022 apparition and the 17 earlier apparitions since 1997. This finding connects Phaethon to the broader population of sunskirting and sungrazing comets observed by SOHO, which often also exhibit bright sodium emission with minimal visible dust, but distinguishes it from other sunskirting asteroids without detectable sodium production under comparable solar heating. These differences may reflect variations in the degree of sodium depletion of near-surface material, and thus the extent and/or timing of any past or present resurfacing activity.

The lower limit on the chicken density function (CDF) of the observable Universe was recently determined to be approximately 10$^{-21}$ chickens pc$^{-3}$. For over a year, however, the scientific community has struggled to determine the upper limit to the CDF. Here we aim to determine a reasonable upper limit to the CDF using multiple observational constraints. We take a holistic approach to considering the effects of a high CDF in various domains, including the Solar System, interstellar medium, and effects on the cosmic microwave background. We find the most restrictive upper limit from the domains considered to be 10$^{23}$ pc$^{-3}$, which ruffles the feathers of long-standing astrophysics theory.

How many gravitational-wave observations from compact object mergers have we seen to date? This seemingly simple question has a surprisingly complex answer that even ChatGPT struggles to answer. To shed light on this, we present a database with the literature's answers to this question. We find values spanning 67-100 for the number of detections from double compact object mergers to date, emphasizing that the exact number of detections is uncertain and depends on the chosen data analysis pipeline and underlying assumptions. We also review the number of gravitational-wave detections expected in the coming decades with future observing runs, finding values up to millions of detections per year in the era of Cosmic Explorer and Einstein Telescope. We present a publicly available code to visualize the detection numbers, highlighting the exponential growth in gravitational-wave observations in the coming decades and the exciting prospects of gravitational-wave astrophysics. See this http URL We plan to keep this database up-to-date and welcome comments and suggestions for additional references.

We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$\alpha$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$\alpha$ forest absorption signal for various surveys, including eBOSS, DESI, and PFS. We show that CO$\times$Ly-$\alpha$ forest can significantly enhance the detection signal-to-noise ratio of CO, with a $200$ to $300 \%$ improvement when cross-correlated with the forest observed in the Prime Focus Spectrograph (PFS) survey and a $50$ to $75\%$ enhancement for the currently available eBOSS or the upcoming DESI observations. We compare to the signal-to-noise improvements expected for a galaxy survey and show that CO$\times$Ly-$\alpha$ is competitive with even a spectroscopic galaxy survey in raw signal-to-noise. Furthermore, our study suggests that the clustering of CO emission is tightly constrained by CO$\times$Ly-$\alpha$ forest, due to the increased signal-to-noise ratio and the simplicity of Ly-$\alpha$ absorption power spectrum modeling. Any foreground contamination or systematics are expected not to be shared between LIM surveys and Ly-$\alpha$ forest observations; this provides an unbiased inference. Our findings highlight the potential benefits of utilizing the Ly-$\alpha$ forest to aid in the initial detection of signals in line intensity experiments. For example, we also estimate that [CII]$\times$Ly-$\alpha$ forest measurements from EXCLAIM and DESI/eBOSS, respectively, should have a larger signal-to-noise ratio than planned [CII]$\times$quasar observations by about an order of magnitude. Our results can be readily applied to actual data thanks to the observed quasar spectra in eBOSS Stripe 82, which overlaps with several LIM surveys.

The immense popularity of spheroid-based scored events (colloquially ``football games'') motivates the desire to better understand the underlying mechanisms affecting their outcomes. By construction of these events, participants must distinguish the spheroidal ball from not only the background, but also their team and enemy players, which are marked by self-assigned linear combinations of specific frequencies of electromagnetic magnetic radiation, known as uniform color. We investigate chromatic effects on the outcome of such events. We do this by finding the correlation between the color contrast and the success of several key spheroidal ball match tactics. We perform this analysis for the 2020 NFL regular season, focusing on moves in which uniform colors may be a factor in performance. We conduct a primary analysis using each team's cumulative results over the season, but in doing so neglect non-uniformity in the chosen uniform color per individual. We then conduct a secondary analysis of the performance per game of a single team, the Seattle Seahawks, which exhibited large uniform color variability for the 2020 NFL regular season. In this work, tackles and completions are considered. The Pearson correlation coefficient is then calculated for both tactics. We find little evidence of chromaticity effects, with correlation values of $r_t=-0.0885\pm 0.1819$ and $r_c-0.0292\pm0.1825 $, respectively, for the primary analysis.

We present ALMA observations of a diffuse gas tracer, CO(J = 1-0), and a warmer and denser gas tracer, CO(J = 3-2), in the overlapping region of interacting galaxies NGC 4567/4568, which are in the early stage of interaction. To comprehend the impact of galaxy interactions on molecular gas properties, we focus on interacting galaxies during the early stage and study their molecular gas properties. In this study, we investigate the physical states of a filamentary molecular structure at the overlapping region, which was previously reported. Utilising new higher-resolution CO(J = 1-0) data, we identify molecular clouds within overlapping and disc regions. Although the molecular clouds in the filament have a factor of two higher an average virial parameter (0.56+-0.14) than that in the overlapping region (0.28+-0.12) and in the disc region (0.26+-0.16), all identified molecular clouds are gravitationally bound. These clouds in the filament also have a larger velocity dispersion than that in the overlapping region, suggesting that molecular gas and/or atomic gas with different velocities collide there. We calculate the ratio of the integrated intensity of CO(J = 3-2) and CO(J = 1-0) (= R3-2/1-0) on the molecular cloud scale. The maximum R3-2/1-0 is 0.17+-0.04 for all identified clouds. The R3-2/1-0 of the molecular clouds in the filament is lower than that of the surrounding area. This result contradicts the predictions of previous numerical simulations, which suggested that the molecular gas on the collision front of galaxies is compressed and becomes denser. Our results imply that NGC 4567/4568 is in a very early stage of interaction; otherwise, the molecular clouds in the filament would not yet fulfil the conditions necessary to trigger star formation.

Missions like NASA's Imaging X-ray Polarimetry Explorer (IXPE) are poised to provide an unprecedented view of the Universe in polarized X-rays. Polarization probes physical anisotropies, a fact exploited by particle physicists to look for the anisotropic $a\boldsymbol{E}\cdot\boldsymbol{B}$ operator in the axion-like-particle (ALP) Lagrangian. Such studies have typically focused on polarization in the radio and microwaves, through local or cosmic birefringence effects. To such polarization studies we add X-rays emanating from magnetars -- a class of neutron stars with near-critical strength magnetic fields -- that are important targets for IXPE. ALPs produced in the neutron star core convert to X-rays in the magnetosphere; such X-rays are polarized along the direction parallel to the dipolar magnetic field at the point of conversion. We develop the full theoretical formalism for ALP-induced polarization in the presence of dipolar magnetic fields. For uncorrelated photon and ALP production mechanisms, we completely disentangle the ALP contributions to the Stokes parameters in terms of the ALP intensity, the ALP-to-photon conversion probability, and the ALP-induced birefringence. In the proper limit, our results demonstrate that the inclusion of ALPs suppresses the observed degree of circular polarization compared to its pure astrophysical value. Our results can also be used to impose limits on ALP couplings with IXPE polarization data from magnetars 4U 0142+61 and 1RXS J170849.0-400910, the subject of upcoming work.

Sun-like dwarf stars in the solar neighborhood reflect ages, an ``average'' chemical evolution, and departures from that average. We show the chemical, and kinematic properties of four groups of Sunlike dwarfs form a continuum related to age. We plot [Fe/H] vs. age, as well as kinematical values for the four groups. The vertical (negative) scatter in [Fe/H] increases with age in a systematic way: as the age increases, [Fe/H] decreases. The sets of Solar and metal-poor stars in the solar neighborhood are related by distributions in [Fe/H] vs. age, as well as in Galactic position (XYZ) and velocity space (UVW). Among the samples there are no clusters of points that set one sample apart from the others. The distributions vary slowly from one set to the next, suggesting a mixture of stellar populations. A plot in Energy vs angular momentum phase space, with coordinate origin moved to the Galactic center, highlights different aspects of the kinematics of the four groups of stars. We finally compare the kinematic properties of these four groups with those of two sets of ultra metal-poor stars.

Very metal-poor stars ([Fe/H] < -2) in the Milky Way are fossil records of early chemical evolution and the assembly and structure of the Galaxy. However, they are rare and hard to find. Gaia DR3 has provided over 200 million low-resolution (R = 50) XP spectra, which provides an opportunity to greatly increase the number of candidate metal-poor stars. In this work, we utilise the XGBoost classification algorithm to identify about 188,000 very metal-poor star candidates. Compared to past work, we increase the candidate metal-poor sample by about an order of magnitude, with comparable or better purity than past studies. Firstly, we develop three classifiers for bright stars (BP < 16). They are classifier-T (for Turn-off stars), classifier-GC (for Giant stars with high completeness), and classifier-GP (for Giant stars with high purity) with expected purity of 47%/47%/74% and completeness of 40%/94%/65% respectively. These three classifiers obtained a total of 11,000/116,000/45,000 bright metal-poor candidates. We apply model-T and model-GP on faint stars (BP > 16) and obtain 13,000/48,500 additional metal-poor candidates with purity 40%/50%, respectively. We make our metal-poor star catalogs publicly available, for further exploration of the metal-poor Milky Way.

Humans like to party, and New Year celebrations are a great way to do that. However New Years celebrations that rely on an orbital year don't line up with those that use a Lunar Calendar, as there are currently 12.368 synodic months (moonths) in a year. There is cyclostratigraphic, paleontological, and tidal rhythmite data that reveal that over billions of years the interplay of angular momentum between the Sun, Earth and Moon has changed the rate of rotation of Earth, and at the same time evolved the orbit of the Moon, and therefore the length of a Lunar month. Using a subset of this data and referencing literature models of the Moon's orbital evolution, we create our own simple model to determine "True Happy New Years", time periods when there were an integer number of lunar synodic months in an Earth orbital year. This would allow modern calendars to pick a shared New Year's Day, and party accordingly. We then predict the next True Happy New Year to be in 252 million years, and offer suggestions to begin the party planning process early, so that we as a planet may be ready.

The lunar crater record features $\sim 50$ basins. The radiometric dating of Apollo samples indicates that the Imbrium basin formed relatively late -- from the planet formation perspective -- some $\simeq 3.9$ Ga. Here we develop a dynamical model for impactors in the inner solar system to provide context for the interpretation of the lunar crater record. The contribution of cometary impactors is found to be insignificant. Asteroids produced most large impacts on the terrestrial worlds in the last $\simeq 3$ Gyr. The great majority of early impactors were rocky planetesimals left behind at $\sim 0.5$--1.5 au after the terrestrial planet accretion. The population of terrestrial planetesimals was reduced by disruptive collisions in the first $t \sim 20$ Myr after the gas disk dispersal. We estimate that there were $\sim 4 \times 10^5$ diameter $d>10$ km bodies when the Moon formed (total planetesimal mass $\sim 0.015$ $M_{\rm Earth}$ at $t \sim 50$ Myr). The early bombardment of the Moon was intense. To accommodate $\sim 50$ known basins, the lunar basins that formed before $\simeq 4.35$--4.41 Ga must have been erased. The late formation of Imbrium occurs with a $\sim 15$--35\% probability in our model. About 20 $d>10$-km bodies were expected to hit the Earth between 2.5 and 3.5 Ga, which is comparable to the number of known spherule beds in the late Archean. We discuss implications of our model for the lunar/Martian crater chronologies, Late Veneer, and noble gases in the Earth atmosphere.

Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to estimate the selection function for different sub-samples of stars in the Gaia catalogue. Comparing the number of stars in a given sub-sample to those in the overall Gaia catalogue, provides an estimate of the sub-sample membership probability, as a function of sky position, magnitude and colour. This estimate must differentiate the stochastic absence of sub-sample stars from selection effects. When multiplied with the overall Gaia catalogue selection function this provides the total selection function of the sub-sample. We present the method by estimating the selection function of the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. The developed method estimates the selection function of the stars present in a sub-sample of Gaia data, given that the sub-sample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package.

We present a systematic search for radio counterparts of novae using the Australian Square Kilometer Array Pathfinder (ASKAP). Our search used the Rapid ASKAP Continuum Survey, which covered the entire sky south of declination $+41^{\circ}$ ($\sim34,000$ square degrees) at a central frequency of 887.5 MHz, the Variables and Slow Transients Pilot Survey, which covered $\sim5,000$ square degrees per epoch (887.5 MHz), and other ASKAP pilot surveys, which covered $\sim200-2000$ square degrees with 2-12 hour integration times. We crossmatched radio sources found in these surveys over a two-year period, from April 2019 to August 2021, with 440 previously identified optical novae, and found radio counterparts for four novae: V5668 Sgr, V1369 Cen, YZ Ret, and RR Tel. Follow-up observations with the Australian Telescope Compact Array confirm the ejecta thinning across all observed bands with spectral analysis indicative of synchrotron emission in V1369 Cen and YZ Ret. Our light-curve fit with the Hubble Flow model yields a value of $1.65\pm 0.17 \times 10^{-4} \rm \:M_\odot$ for the mass ejected in V1369 Cen. We also derive a peak surface brightness temperature of $250\pm80$ K for YZ Ret. Using Hubble Flow model simulated radio lightcurves for novae, we demonstrate that with a 5$\sigma$ sensitivity limit of 1.5 mJy in 15-min survey observations, we can detect radio emission up to a distance of 4 kpc if ejecta mass is in the range $10^{-3}\rm \:M_\odot$, and upto 1 kpc if ejecta mass is in the range $10^{-5}-10^{-3}\rm \:M_\odot$. Our study highlights ASKAP's ability to contribute to future radio observations for novae within a distance of 1 kpc hosted on white dwarfs with masses $0.4-1.25\:\rm M_\odot$ , and within a distance of 4 kpc hosted on white dwarfs with masses $0.4-1.0\:\rm M_\odot$.

We carry out hydrodynamical simulations to study the eccentricity growth of a 1-30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dominate from 1:3 up to 3:5. The maximum values of eccentricity reached in our simulations are 0.65-0.75. We calculate the eccentricity growth rate for different planet masses and disc parameters and derive analytical dependencies on these parameters. We observe that the growth rate is proportional to both the planet's mass and the characteristic disc mass for a wide range of parameters. In a separate set of simulations, we derived the width of the 1:3 Lindblad resonance.

The inference of cosmological quantities requires accurate and large hydrodynamical cosmological simulations. Unfortunately, their computational time can take millions of CPU hours for a modest coverage in cosmological scales ($\approx (100 {h^{-1}}\,\text{Mpc})^3)$). The possibility to generate large quantities of mock Lyman-$\alpha$ observations opens up the possibility of much better control on covariance matrices estimate for cosmological parameters inference, and on the impact of systematics due to baryonic effects. We present a machine learning approach to emulate the hydrodynamical simulation of intergalactic medium physics for the Lyman-$\alpha$ forest called LyAl-Net. The main goal of this work is to provide highly efficient and cheap simulations retaining interpretation abilities about the gas field level, and as a tool for other cosmological exploration. We use a neural network based on the U-net architecture, a variant of convolutional neural networks, to predict the neutral hydrogen physical properties, density, and temperature. We train the LyAl-Net model with the Horizon-noAGN simulation, though using only 9% of the volume. We also explore the resilience of the model through tests of a transfer learning framework using cosmological simulations containing different baryonic feedback. We test our results by analysing one and two-point statistics of emulated fields in different scenarios, as well as their stochastic properties. The ensemble average of the emulated Lyman-$\alpha$ forest absorption as a function of redshift lies within 2.5% of one derived from the full hydrodynamical simulation. The computation of individual fields from the dark matter density agrees well with regular physical regimes of cosmological fields. The results tested on IllustrisTNG100 showed a drastic improvement in the Lyman-$\alpha$ forest flux without arbitrary rescaling.

Nontrivial dark sector physics continues to be an interesting avenue in our quest to the nature of dark matter. In this paper, we study the cosmological signatures of mass-varying dark matter where its mass changes from zero to a nonzero value in the early Universe. We compute the changes in various observables, such as, the matter and the cosmic microwave background anisotropy power spectrum. We explain the origin of the effects and point out a qualitative similarity between this model and a warm dark matter cosmology with no sudden mass transition. We also do a simple frequentist analysis of the linear matter power spectrum to estimate the constraint on the parameters of this model from latest cosmological observation data.

Superoutbursts in WZ Sge-type dwarf novae (DNe) are characterized by both early superhumps and ordinary superhumps originating from the 2:1 and 3:1 resonances, respectively. However, some WZ Sge-type DNe show a superoutburst lacking early superhumps; it is not well established how these differ from superoutbursts with an early superhump phase. We report time-resolved photometric observations of the WZ Sge-type DN V627 Peg during its 2021 superoutburst. The detection of ordinary superhumps before the superoutburst peak highlights that this 2021 superoutburst of V627 Peg, like that {in} 2014, did not feature an early superhump phase. The duration of stage B superhumps was slightly longer in the 2010 superoutburst accompanying early superhumps than that in the 2014 and 2021 superoutbursts which lacked early superhumps. This result suggests that an accretion disk experiencing the 2:1 resonance may have a larger mass at the inner part of the disk and hence take more time for the inner disk to become eccentric. The presence of a precursor outburst in the 2021 superoutburst suggests that the maximum disk radius should be smaller than that of the 2014 superoutburst, even though the duration of quiescence was longer than that before the 2021 superoutburst. This could be accomplished if the 2021 superoutburst was triggered as an inside-out outburst or if the mass transfer rate in quiescence changes by a factor of two, suggesting that the outburst mechanism and quiescence state of WZ Sge-type DNe may have more variety than ever thought.

In this chapter we review our knowledge of our galaxy's cometary population outside our Oort Cloud - exocomets and Interstellar Objects (ISOs). We start with a brief overview of planetary system formation, viewed as a general process around stars. We then take a more detailed look at the creation and structure of exocometary belts, as revealed by the unprecedented combination of theoretical and observational advances in recent years. The existence and characteristics of individual exocomets orbiting other stars is summarized, before looking at the mechanisms by which they may be ejected into interstellar space. The discovery of the first two ISOs is then described, along with the surprising differences in their observed characteristics. We end by looking ahead to what advances may take place in the next decade.

The dark matter content of globular clusters, highly compact gravity-bound stellar systems, is unknown. It is also generally unknow*able*, due to their mass-to-light ratios typically ranging between 1$-$3 in solar units, accommodating a dynamical mass of dark matter at best comparable to the stellar mass. That said, recent claims in the literature assume densities of dark matter around 1000 GeV/cm$^3$ to set constraints on its capture and annihilation in white dwarfs residing in the globular cluster M4, and to study a number of other effects of dark matter on compact stars. Motivated by these studies, we use measurements of stellar kinematics and luminosities in M4 to look for a dark matter component via a spherical Jeans analysis; we find no evidence for it, and set the first empirical limits on M4's dark matter distribution. Our density upper limits, a few $\times \ 10^4$ GeV/cm$^3$ at 1 parsec from the center of M4, do not negate the claims (nor confirm them), but do preclude the use of M4 for setting limits on non-annihilating dark matter kinetically heating white dwarfs, which require at least $10^5$ GeV/cm$^3$ densities. The non-robust nature of globular clusters as dynamical systems, combined with evidence showing that they may originate from molecular gas clouds in the absence of dark matter, make them unsuitable as laboratories to unveil dark matter's microscopic nature in current or planned observations.

To simulate the heated exterior of a meteorite, we impact a granular bed with melted tin. The morphology of tin remnant and crater is found to be sensitive to the temperature and solidification of tin. By employing deep learning and convolutional neural network, we can quantify and map the complex impact patterns onto network systems based on feature maps and Grad-CAM results. This gives us unprecedented details on how the projectile deforms and interacts with the granules, which information can be used to trace the development of different remnant shapes. Furthermore, full dynamics of granular system is revealed by the use of Particle Image Velocimetry. Kinetic energy, temperature and diameter of the projectile are used to build phase diagrams for the morphology of both crater and tin remnant. In addition to successfully reproducing key features of simple and complex craters, we are able to detect a possible artifact when compiling crater data from field studies. The depth of craters from high-energy impacts in our work is found to be independent of their width. However, when mixing data from different energy, temperature and diameter of projectile, a bogus power-law relationship appears between them. Like other controlled laboratory researches, our conclusions have the potential to benefit the study of paint in industry and asteroid sampling missions on the surface of celestial bodies.

Context: With the data releases from the astrometric space mission Gaia, the exploration of the structure of the Milky Way has developed in unprecedented detail and unveiled many previously unknown structures in the Galactic disc and halo. One such feature is the phase spiral where the stars in the Galactic disc form a spiral density pattern in the $Z-V_Z$ plane. Aims: We aim to characterize the shape, rotation, amplitude, and metallicity of the phase spiral in the outer disc of the Milky Way. This will allow us to better understand which physical processes caused the phase spiral and can give further clues to the Milky Way's past and the events that contributed to its current state. Methods: We use Gaia data release 3 (DR3) to get full position and velocity data on approximately 31.5 million stars, and metallicity for a subset of them. We then compute the angular momenta of the stars and develop a model to characterise the phase spiral in terms of amplitude and rotation at different locations in the disc. Results: We find that the rotation angle of the phase spiral changes with Galactic azimuth and Galactocentric radius, making the phase spiral appear to rotate about $3^\circ$ per degree in Galactic azimuth. Furthermore, we find that the phase spiral in the $2200 - 2400$ kpc km s$^{-1}$ range of angular momentum is particularly strong compared to the phase spiral that can be observed in the solar neighbourhood. The metallicity of the phase spiral appears to match that of the Milky Way disc field stars. Conclusions: We created a new model capable of fitting several key parameters of the phase spiral. We have been able to determine the rotation rate of the phase spiral and found a peak in the phase spiral amplitude which manifests as a very clear phase spiral when using only stars with similar angular momentum.

Over the past few years, wide-field time-domain surveys like ZTF and OGLE have led to discoveries of various types of interesting short-period stellar variables, such as ultracompact eclipsing binary white dwarfs, rapidly rotating magnetised white dwarfs (WDs), transitional cataclysmic variables between hydrogen-rich and helium accretion, and blue large-amplitude pulsators (BLAPs), which greatly enrich our understandings of stellar physics under some extreme conditions. In this paper, we report the first-two-year discoveries of short-period variables (i.e., P<2 hr) by the Tsinghua University-Ma Huateng Telescopes for Survey (TMTS). TMTS is a multi-tube telescope system with a field of view up to 18 deg^2, which started to monitor the LAMOST sky areas since 2020 and generated uninterrupted minute-cadence light curves for about ten million sources within 2 years. Adopting the Lomb-Scargle periodogram with period-dependent thresholds for the maximum powers, we identify over 1 100 sources that exhibit a variation period shorter than 2 hr. Compiling the light curves with the Gaia magnitudes and colours, LAMOST spectral parameters, VSX classifications, and archived observations from other prevailing time-domain survey missions, we identified 1 076 as delta Scuti stars, which allows us study their populations and physical properties in the short-period regime. The other 31 sources include BLAPs, subdwarf B variables (sdBVs), pulsating WDs, ultracompact/short-period eclipsing/ellipsoidal binaries, cataclysmic variables below the period gap, etc., which are highly interesting and worthy of follow-up investigations.

We investigate the role of dense environments in suppressing star formation by studying $\rm \log_{10}(M_\star/M_\odot) > 9.7$ star-forming galaxies in nine clusters from the Local Cluster Survey ($0.0137 < z < 0.0433$) and a large comparison field sample drawn from the Sloan Digital Sky Survey. We compare the star-formation rate (SFR) versus stellar mass relation as a function of environment and morphology. After carefully controlling for mass, we find that in all environments, the degree of SFR suppression increases with increasing bulge-to-total (B/T) ratio. In addition, the SFRs of cluster and infall galaxies at a fixed mass are more suppressed than their field counterparts at all values of B/T. These results suggest a quenching mechanism that is linked to bulge growth that operates in all environments and an additional mechanism that further reduces the SFRs of galaxies in dense environments. We limit the sample to $B/T < 0.3$ galaxies to control for the trends with morphology and find that the excess population of cluster galaxies with suppressed SFRs persists. We model the timescale associated with the decline of SFRs in dense environments and find that the observed SFRs of the cluster core galaxies are consistent with a range of models including: a mechanism that acts slowly and continuously over a long (2-5 Gyr) timescale, and a more rapid ($<1$ Gyr) quenching event that occurs after a delay period of 1-6 Gyr. Quenching may therefore start immediately after galaxies enter clusters.

We investigate the physical nature of Active Galactic Nuclei using machine learning tools. We show that the redshift $z$, the bolometric luminosity $L_{\rm Bol}$, the central mass of the supermassive black hole $M_{\rm BH}$, the Eddington ratio $\lambda_{\rm Edd}$ as well as the AGN class (obscured or unobscured) can be reconstructed through multi-wavelength photometric observations only. A Support Vector Regression (SVR) ML-model is trained on 7616 of spectroscopically observed AGN from the SPIDERS-AGN survey, previously cross-matched with soft X-ray observations (from ROSAT or XMM), WISE mid-infrared photometry, and optical photometry from SDSS $ugriz$ filters. We build a catalogue of 21364 AGN to be reconstructed with the trained SVR: for 9944 sources, we found archival redshift measurements. All AGN are classified as either Type 1/2 using a Random Forest (RF) algorithm on a subset of known sources. All known photometric measurement uncertainties are incorporated using a simulation-based approach. We present the reconstructed catalogue of 21364 AGN with redshifts ranging from $ 0 < z < 2.5$. $z$ estimations are made for 11420 new sources, with an outlier rate within 10%. Type 1/2 AGN can be identified with respective efficiencies of 88% and 93%: the estimated classification of all sources is given in the dataset. $L_{\rm Bol}$, $M_{\rm BH}$, and $\lambda_{\rm Edd}$ values are given for 16907 new sources with their estimated error. These results have been made publicly available. The release of this catalogue will advance AGN studies by presenting key parameters of the accretion history of 6 dex in luminosity over a wide range of $z$. Similar applications of ML techniques using photometric data only will be essential in the future, with large datasets from eROSITA, JSWT and the VRO poised to be released in the next decade.

Coronal jets are eruptions identified by a collimated, sometimes twisted spire. They are small-scale energetic events compared with flares. Using multi-wavelength observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) and a magnetogram from Hinode/Spectro-Polarimeter (Hinode/SP), we study the formation and evolution of a jet occurring on 2019 March 22 in the active region NOAA 12736. A zero-$\beta$ magnetohydrodynamic (MHD) simulation is conducted to probe the initiation mechanisms and appearance of helical motion during this jet event. As the simulation reveals, there are two pairs of field lines at the jet base, indicating two distinct magnetic structures. One structure outlines a flux rope lying low above the photosphere in the north of a bald patch region and the other structure shows a null point high in the corona in the south. The untwisting motions of the observed flux rope was recovered by adding an anomalous (artificial) resistivity in the simulation. A reconnection occurs at the bald patch in the flux rope structure, which is moving upwards and simultaneously encounters the field lines of the null point structure. The interaction of the two structures results in the jet while the twist of the flux rope is transferred to the jet by the reconnected field lines. The rotational motion of the flux rope is proposed to be an underlying trigger of this process and responsible for helical motions in the jet spire.

We have conducted a revised analysis of the first-order phase transition that is associated with symmetry breaking in a classically scale-invariant model that has been extended with a new $SU(2)$ gauge group. By incorporating recent developments in the understanding of supercooled phase transitions, we were able to calculate all of its features and significantly limit the parameter space. We were also able to predict the gravitational wave spectra generated during this phase transition and found that this model is well-testable with LISA. Additionally, we have made predictions regarding the relic dark matter abundance. Our predictions are consistent with observations but only within a narrow part of the parameter space. We have placed significant constraints on the supercool dark matter scenario by improving the description of percolation and reheating after the phase transition, as well as including the running of couplings. Finally, we have also analyzed the renormalization-scale dependence of our results.

The infall and merger scenario of massive clusters in the Milky Way's potential well, as one of the Milky Way formation mechanisms, is reexamined to understand how the stars of the merging clusters are redistributed during and after the merger process using, for the first time, simulations with a high resolution concentrated in the 300 pc around the Galactic center. We adopted simulations developed in the framework of the "Modelling the Evolution of Galactic Nuclei" (MEGaN) project. We compared the evolution of representative clusters in the mass and concentration basis in the vicinity of a supermassive black hole. We used the spatial distribution, density profile, and the $50\%$ Lagrange radius (half mass radius) as indicators along the complete simulation to study the evolutionary shape in physical and velocity space and the final fate of these representative clusters. We detect that the least massive clusters are quickly (<10 Myr) destroyed. Instead, the most massive clusters have a long evolution, showing variations in the morphology, especially after each passage close to the supermassive black hole. The deformation of the clusters depends on the concentration, with general deformations for the least concentrated clusters and outer strains for the more concentrated ones. At the end of the simulation, a dense concentration of stars belonging to the clusters is formed. The particles that belong to the most massive and most concentrated clusters are concentrated in the innermost regions, meaning that the most massive and concentrated clusters contribute with a more significant fraction of particles to the final concentration, which suggests that the population of stars of the nuclear star cluster formed through this mechanism comes from massive clusters rather than low-mass globular clusters.

The code HARM\_COOL, a conservative scheme for relativistic magnetohydrodynamics, is being developed in our group and works with a tabulated equation of state of dense matter. This EOS can be chosen and used during dynamical simulation, instead of the simple ideal gas one. In this case, the inversion scheme between the conserved and primitive variables is not a trivial task. In principle, the code needs to solve numerically five coupled non-linear equations at every time-step. The 5-D recovery schemes were originally implemented in HARM and worked accurately for a simple polytropic EOS which has an analytic form. Our current simulations support the composition-dependent EOS, formulated in terms of rest-mass density, temperature and electron fraction. In this proceeding, I discuss and compare several recovery schemes that have been included in our code. I also present and discuss their convergence tests. Finally, I show set of preliminary results of a numerical simulation, addressed to the post-merger system formed after the binary neutron stars (BNS) coalescence.

We investigate the question of the nature of compact stars, considering they may be neutron stars or hybrid stars containing a quark core, within the present constraints given by gravitational waves, radio-astronomy, X-ray emissions from millisecond pulsars and nuclear physics. A Bayesian framework is used to combine together all these constraints and to predict tidal deformabilities and radii for a 1.4~M$_\odot$ compact star. We find that present gravitation wave and radio-astronomy data favors asy-stiff EoS compatible with nuclear physics and that GW170817 waveform is best described for binary hybrid stars. In addition, this data favors stiff quark matter, independently of the nuclear EoS. Combining this result with constraints from X-ray observation supports the existence of canonical $1.4$~M$_\odot$ mass hybrid star, with a radius predicted to be $R_{1.4}=12.02(8)$~km.

The TNG300-1 run of the IllustrisTNG simulations includes 1697 clusters of galaxies with $M_{200c}>10^{14}$M$_\odot$ covering the redshift range $0.01-1.04$. We build mock spectroscopic redshift catalogues of simulated galaxies within these clusters and apply the caustic technique to estimate the cumulative cluster mass profiles. We compute the total true cumulative mass profile from the 3D simulation data and calculate the ratio of caustic mass to total 3D mass, $\mathcal{F}_\beta$, as a function of cluster-centric distance and identify the radial range where $\mathcal{F}_\beta$ is roughly constant. The filling factor, $\mathcal{F}_\beta=0.41\pm 0.08$, is constant on a plateau that covers a wide cluster-centric distance range, $(0.6-4.2)R_{200c}$. This calibration is insensitive to redshift. The calibrated caustic mass profiles are unbiased, with an average uncertainty of $23\%$. At $R_{200c}$, the average $M^C/M^{3D}=1.03\pm 0.22$; at $2R_{200c}$, the average $M^C/M^{3D}=1.02\pm 0.23$. Simulated galaxies are unbiased tracers of the mass distribution. IllustrisTNG is a broad statistical platform for application of the caustic technique to large samples of clusters with spectroscopic redshifts for $\gtrsim 200$ members in each system. These observations will allow extensive comparisons with weak lensing masses and will complement other techniques for measuring the growth rate of structure in the universe.

The nonlinear outcome of plasma instabilities range from a gentle reconfiguration of the initial state to an explosive one, and a non-disruptive outcome in between can nevertheless still be a route to turbulence. The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, all these possibilities can occur, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. Here we derive several general identities from the evolution equation for entropy that reveals the mechanism by which TI saturates: whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change sign. This temporal event implies that the gas pressure force reverses direction, slowing the continued growth of the condensation. For isobaric evolution, this `pressure reversal' occurs nearly simultaneously for every fluid element in the condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave-mediated as a result of the pressure reversal occurring at different times for different locations in the cloud core. We show that both the contraction rebound event and the subsequent pulsation behavior follow analytically from an analysis of the new identities.

Terminal velocity reached by bubble walls in first order phase transitions is an important parameter determining both primordial gravitational-wave spectrum and production of baryon asymmetry in models of electroweak baryogenesis. We developed a numerical code to study the real-time evolution of expanding bubbles and investigate how their walls reach stationary states. Our results agree with profiles obtained within the so-called bag model with very good accuracy, however, not all such solutions are stable and realised in dynamical systems. Depending on the exact shape of the potential there is always a range of wall velocities where no steady state solutions exist. This behaviour in deflagrations was explained by hydrodynamical obstruction where solutions that would heat the plasma outside the wall above the critical temperature and cause local symmetry restoration are forbidden. For even more affected hybrid solutions causes are less straight forward, however, we provide a simple numerical fit allowing one to verify if a solution with a given velocity is allowed simply by computing the ratio of the nucleation temperature to the critical one for the potential in question.

I solve the problem of nomenclature of planets, stars, and moons, and in doing so repair two of the IAU's blunders. Drawing and improving upon foundational work by Chen & Kipping, I describe a single, physics-based taxonomy that christens all objects in hydrostatic equilibrium as "stars," a category that contains several subcategories based on the relevant pressure terms in the equation of state. I also acknowledge dynamical considerations, which allow me to describe a single designation scheme for all "stars" following the Washington Multiplicity Catalog convention. Under this unified scheme, what we used to call "Planet Earth" is now the moon rock "star" Sun Da.

We observe neither life beyond Earth, nor moons around exoplanets, despite the prevalence of Earth-like planets across the galaxy. We suggest Moonfall as a possible mechanism to explain both simultaneously.

The massive elliptical galaxy M87 has been the subject of several supermassive black hole mass measurements from stellar dynamics, gas dynamics, and recently the black hole shadow by the Event Horizon Telescope (EHT). This uniquely positions M87 as a benchmark for alternative black hole mass determination methods. Here we use stellar kinematics extracted from integral-field spectroscopy observations with Adaptive Optics (AO) using MUSE and OASIS. We exploit our high-resolution integral field spectroscopy to spectrally decompose the central AGN from the stars. We derive an accurate inner stellar-density profile and find it is flatter than previously assumed. We also use the spectrally-extracted AGN as a reference to accurately determine the observed MUSE and OASIS AO PSF. We then perform Jeans Anisotropic Modelling (JAM), with a new flexible spatially-variable anisotropy, and measure the anisotropy profile, stellar mass-to-light variations, inner dark matter fraction, and black hole mass. Our preferred black hole mass is $M_{\rm BH}=(8.7\pm1.2) \times 10^9 \ M_\odot $. However, using the inner stellar density from previous studies, we find a preferred black hole mass of $M_{\rm BH} = (5.5^{+0.5}_{-0.3}) \times 10^9 \ M_\odot $, consistent with previous work. We conduct numerous systematic tests of the kinematics and model assumptions and conclude that uncertainties in the black hole mass of M87 from previous determinations may have been underestimated and further analyses are needed.

We revive a cross-platform focal-plane wavefront sensing and control algorithm originally released in 1980 and show that it can provide significant contrast improvements over conventional control methods on coronagraphic instruments. Its simplicity makes it applicable to various coronagraph models and we demonstrate it on a classical Lyot coronagraph and a phase-apodized pupil Lyot coronagraph, both in simulation and in laboratory experiments. Surprisingly, it had been forgotten for decades, but we present its unbeatable advantages considering the increase in computational power in the last 40 years. We consider it a major game changer in the planning for future, space-based high-contrast imaging missions and recommend it be intensively revisited by all readers.

A complete demographic of active galactic nuclei (AGN) is essential to understand the evolution of the Universe. Optical surveys estimate the population of AGN in the local Universe to be of $\sim$ 4%. However, these results could be biased towards bright sources, not affected by the host galaxy attenuation. An alternative method for detecting these objects is through the X-ray emission. In this work, we aim to complement the AGN population of the optical CALIFA survey (941 sources), by using X-ray data from Chandra, which provides the best spatial resolution to date, essential to isolate the nuclear emission from the host galaxy. We study a total of 138 sources with available data. We find 34 new bonafide AGN and 23 AGN candidates, which could increase the AGN population to 7-10\% among the CALIFA survey. X-rays are particularly useful for low-luminosity AGN since they are excluded by the criterion of large equivalent width of the $H\alpha$ emission line when applied to optical selections. Indeed, placing such a restrictive criteria might cause a loss of up to 70% of AGN sources. X-ray detected sources are preferentially located in the right side of the [$OIII$]/$H\beta$ versus [$NII$]/$H\alpha$ diagram, suggesting that this diagram might be the most reliable at classifying AGN sources. Our results support the idea that multi-wavelength studies are the best way to obtain a complete AGN population.

Far away observers can in principle bound from below the dimensionless maximum-density parameter $\Lambda\equiv4\pi R^2\rho_{\text{max}}$ of a compact star by measuring the gravitational redshift factor $z\equiv\nu_{\text{e}}/\nu_{\infty}-1$ of photons that were emitted from the {\it surface} of the star: $\Lambda\geq{3\over2}[1-(1+z)^{-2}]$ [here $R$ is the radius of the star and $\{\nu_{\text{e}},\nu_{\infty}\}$ are respectively the frequency of the emitted light as measured at the location of the emission and by asymptotic observers]. However, if photons that were created somewhere {\it inside} the star can make their way out and reach the asymptotic observers, then the measured redshift parameter $z$ may not determine uniquely the surface properties of the star, thus making the above bound unreliable. In the present compact paper we prove that in these cases, in which the creation depth of a detected photon is not known to the far away observers, the empirically measured redshift parameter can still be used to set a (weaker) lower bound on the dimensionless density parameter of the observed star: $\Lambda\geq{3\over2}[1-(1+z)^{-2/3}]$.

Astronomical images can be fascinating to the general public, but the interaction is typically limited to contemplation. Numerical simulations of astronomical systems do permit a closer interaction, but are generally unknown outside the research community. We are developing "Protoplanet Express", a video game based on hydrodynamical simulations of protoplanetary discs. In the game, the player visits several discs, finds its relevant features and learns about them. Here we present the current version of the game, discuss its reception, and consider its further development.

Accurately predicting the z-component of the interplanetary magnetic field, particularly during the passage of an interplanetary coronal mass ejection (ICME), is a crucial objective for space weather predictions. Currently, only a handful of techniques have been proposed and they remain limited in scope and accuracy. Recently, a robust machine learning (ML) technique was developed for predicting the minimum value of Bz within ICMEs based on a set of 42 'features', that is, variables calculated from measured quantities upstream of the ICME and within its sheath region. In this study, we investigate these so-called explanatory variables in more detail, focusing on those that were (1) statistically significant; and (2) most important. We find that number density and magnetic field strength accounted for a large proportion of the variability. These features capture the degree to which the ICME compresses the ambient solar wind ahead. Intuitively, this makes sense: Energy made available to CMEs as they erupt is partitioned into magnetic and kinetic energy. Thus, more powerful CMEs are launched with larger flux-rope fields (larger Bz), at greater speeds, resulting in more sheath compression (increased number density and total field strength).

We consider a wide class of four-dimensional effective field theories in which gravity is coupled to multiple four-forms and their dual scalar fields, with membrane sources charged under the corresponding three-form potentials. Four-form flux, quantised in units of the membrane charges, generically generates a landscape of vacua with a range of values for the cosmological constant that is scanned through membrane nucleation. We list various ways in which the landscape can be made sufficiently dense to be compatible with observations of the current vacuum without running into the empty universe problem. Further, we establish the general criteria required to ensure the absolute stability of the Minkowski vacuum under membrane nucleation and the longevity of those vacua that are parametrically close by. This selects the current vacuum on probabilistic grounds and can even be applied in the classic model of Bousso and Polchinski, albeit with some mild violation of the membrane weak gravity conjecture. We present other models where the membrane weak gravity conjecture is not violated but where the same probabilistic methods can be used to tackle the cosmological constant problem.

In what will likely be a litany of generative-model-themed arXiv submissions celebrating April the 1st, we evaluate the capacity of state-of-the-art transformer models to create a paper detailing the detection of a Pulsar Wind Nebula with a non-existent Imaging Atmospheric Cherenkov Telescope (IACT) Array. We do this to evaluate the ability of such models to interpret astronomical observations and sources based on language information alone, and to assess potential means by which fraudulently generated scientific papers could be identified during peer review (given that reliable generative model watermarking has yet to be deployed for these tools). We conclude that our jobs as astronomers are safe for the time being. From this point on, prompts given to ChatGPT and Stable Diffusion are shown in orange, text generated by ChatGPT is shown in black, whereas analysis by the (human) authors is in blue.

Taking axion inflation as an example, we estimate the maximal temperature ($T_{\rm max}^{ }$) that can be reached in the post-inflationary universe, as a function of the confinement scale of a non-Abelian dark sector ($\Lambda_{\rm IR}^{ }$). Below a certain threshold $\Lambda_{\rm IR}^{ } < \Lambda_{\rm 0}^{ } \sim 2\times 10^{-8}_{ } m_{\rm pl}^{ }$, the system heats up to $T_{\rm max}^{ } \sim \Lambda_{\rm 0}^{ } > T_{\rm c}^{ }$, and a first-order thermal phase transition takes place. On the other hand, if $\Lambda_{\rm IR}^{ } > \Lambda_{\rm 0}^{ }$, then $T_{\rm max}^{ } \sim \Lambda_{\rm IR}^{ } < T_{\rm c}^{ }$: very high temperatures can be reached, but there is no phase transition. If the inflaton thermalizes during heating-up (which we find to be unlikely), or if the plasma includes light degrees of freedom, then heat capacity and entropy density are larger, and $T_{\rm max}^{ }$ is lowered towards $\Lambda_{\rm 0}^{ }$. The heating-up dynamics generates a gravitational wave background. Its contribution to $N^{ }_{\rm eff}$ at GHz frequencies, the presence of a monotonic $\sim f_{\rm 0}^3$ shape at $(10^{-4}_{ } - 10^2_{ })\,$Hz frequencies, and the frequency domain of peaked features that may originate via first-order phase transitions, are discussed.

The free propagator of a massless mode in an expanding universe can be written as a sum of two terms, a lightcone and a tail part. The latter describes a subluminal (time-like) signal. We show that the inflationary gravitational wave background, influencing cosmic microwave background polarization, and routinely used for constraining inflationary models through the so-called $r$ ratio, originates exclusively from the tail part.

Spin precession is one of the key physical effects that could unveil the origin of the compact binaries detected by ground- and space-based gravitational-wave (GW) detectors, and shed light on their possible formation channels. Efficiently and accurately modeling the GW signals emitted by these systems is crucial to extract their properties. Here, we present SEOBNRv5PHM, a multipolar precessing-spin waveform model within the effective-one-body (EOB) formalism for the full signal (i.e. inspiral, merger and ringdown) of binary black holes (BBHs). In the non-precessing limit, the model reduces to SEOBNRv5HM, which is calibrated to $442$ numerical-relativity (NR) simulations, 13 waveforms from BH perturbation theory, and non-spinning energy flux from second-order gravitational self-force theory. We remark that SEOBNRv5PHM is not calibrated to precessing-spin NR waveforms from the Simulating eXtreme Spacetimes Collaboration. We validate SEOBNRv5PHM by computing the unfaithfulness against 1543 precessing-spin NR waveforms, and find that for 99.8% (84.4%) of the cases, the maximum value, in the total mass range 20-300 $M_\odot$, is below 3% (1%). These numbers reduce to 95.3% (60.8%) when using the previous version of the SEOBNR family, SEOBNRv4PHM, and to 78.2% (38.3%) when using the state-of-the-art frequency-domain multipolar precessing-spin phenomenological IMRPhenomXPHM model. Due to much better computational efficiency of SEOBNRv5PHM compared to SEOBNRv4PHM, we are also able to perform extensive Bayesian parameter estimation on synthetic signals and GW events observed by LIGO-Virgo detectors. We show that SEOBNRv5PHM can be used as a standard tool for inference analyses to extract astrophysical and cosmological information of large catalogues of BBHs.

We propose a new method to detect sub-GeV dark matter, through their scatterings from free leptons and the resulting kinematic shifts. Specially, such an experiment can detect dark matter interacting solely with muons. The experiment proposed here is to directly probe muon-philic dark matter, in a model-independent way. Its complementarity with the muon on target proposal, is similar to, e.g. XENON/PandaX and ATLAS/CMS on dark matter searches. Moreover, our proposal can work better for relatively heavy dark matter such as in the sub-GeV region. We start with a small device of a size around 0.1 to 1 meter, using atmospheric muons to set up a prototype. Within only one year of operation, the sensitivity on cross section of dark matter scattering with muons can already reach $\sigma_D\sim 10^{-19 (-20,\,-18)}\rm{cm}^{2}$ for a dark mater $\rm{M_D}=100\, (10,\,1000)$ MeV. We can then interface the device with a high intensity muon beam of $10^{12}$/bunch. Within one year, the sensitivity can reach $\sigma_D\sim 10^{-27 (-28,\,-26)}\rm{cm}^{2}$ for $\rm{M_D}=100\, (10,\,1000)$ MeV.

The structure of a vortex in the inner crust of a pulsar is calculated microscopically in the Wigner-Seitz cell approximation, simulating the conditions of the inner crust of a cold, non-accreting neutron star, in which a lattice of nuclei coexists with a sea of superfluid neutrons. The calculation is based on the axially deformed Hartree-Fock-Bogolyubov framework, using effective interactions. The present work extends and improves previous studies in four ways: i) it allows for the axial deformation of protons induced by the large deformation of neutrons due to the appearance of vortices; ii) it includes the effect of Coulomb exchange; iii) considers the possible effects of the screening of the pairing interaction; and iv) it improves the numerical treatment. We also demonstrate that the binding energy of the nucleus-vortex system can be used as a proxy to the pinning energy of a vortex and discuss in which conditions this applies. From our results, we can estimate the mesoscopic pinning forces per unit length acting on vortices. We obtain values ranging between $10^{14}$ to $10^{16}$ dyn/cm, consistent with previous findings.

Spacecraft observations showed that electron heat conduction in the solar wind is probably regulated by whistler waves, whose origin and efficiency in electron heat flux suppression is actively investigated. In this paper, we present Particle-In-Cell simulations of a combined whistler heat flux and temperature anisotropy instability that can operate in the solar wind. The simulations are performed in a uniform plasma and initialized with core and halo electron populations typical of the solar wind. We demonstrate that the instability produces whistler waves propagating both along (anti-sunward) and opposite (sunward) to the electron heat flux. The saturated amplitudes of both sunward and anti-sunward whistler waves are strongly correlated with their {\it initial} linear growth rates, $B_{w}/B_0\sim (\gamma/\omega_{ce})^{\nu}$, where for typical electron betas we have $0.6\lesssim \nu\lesssim 0.9$. The correlations of whistler wave amplitudes and spectral widths with plasma parameters (electron beta and temperature anisotropy) revealed in the simulations are consistent with those observed in the solar wind. The efficiency of electron heat flux suppression is positively correlated with the saturated amplitude of sunward whistler waves. The electron heat flux can be suppressed by 10--60% provided that the saturated amplitude of sunward whistler waves exceeds about 1% of background magnetic field. Other experimental applications of the presented results are discussed.