Papers for Monday, Oct 24 2022

Upcoming cosmological surveys will provide unprecedented amount of data, which will require innovative statistical methods to maximize the scientific exploitation. Standard cosmological analyses based on abundances, two-point and higher-order statistics of cosmic tracers have been widely used to investigate the properties of the cosmic web and Large Scale Structure. However, these statistics can only exploit a subset of the entire information content available. Our goal is thus to implement new data analysis techniques based on machine learning to extract cosmological information through forward modelling, by directly exploiting the spatial coordinates and other observed properties of galaxies and galaxy clusters. Specifically, we investigated a new representation of large-scale structure data in the form of graphs. This data format can be directly fed to Graph Neural Networks, a recently proposed class of supervised Deep Learning algorithms. We tested the method on dark matter halo catalogues in different cosmologies, finding promising results. In particular, the method can discriminate among different dark energy models with high accuracy, through both binary classification ($99\%-$accuracy) and multi-class classification ($97\%-$accuracy). Moreover, it provides constraints on the value of $w_0$, through regression, with high precision.

As one of the most commonly observed disk substructures, dust rings from high-resolution disk surveys appear to have different radial widths. Recent observations on PDS 70 and AB Aur reveal not only planets in the disk, but also the accompanying wide dust rings. We use three-dimensional dust-and-gas disk simulations to study whether gap-opening planets are responsible for the large ring width in disk observations. We find that gap-opening planets can widen rings of dust trapped at the pressure bump via planetary perturbations, even with the mid-plane dust-to-gas ratio approaching order unity and with the dust back-reaction accounted for. We show that the planet-related widening effect of dust rings can be quantified using diffusion-advection theory, and provide a generalized criterion for an equilibrated dust ring width in three-dimensional disk models. We also suggest that the ring width can be estimated using the gas turbulent viscosity $\alpha_{\rm turb}$, but with cautions about the Schmidt number greater than order unity.

We present a systematic numerical relativity study of the impact of different treatment of microphysics and grid resolution in binary neutron star mergers. We consider series of simulations at multiple resolutions comparing hydrodynamics, neutrino leakage scheme, leakage augmented with the M0 scheme and the more consistent M1 transport scheme. Additionally, we consider the impact of a sub-grid scheme for turbulent viscosity. We find that viscosity helps to stabilise the remnant against gravitational collapse but grid resolution has a larger impact than microphysics on the remnant's stability. The gravitational wave (GW) energy correlates with the maximum remnant density, that can be thus inferred from GW observations. M1 simulations shows the emergence of a neutrino trapped gas that locally decreases the temperature a few percent when compared to the other simulation series. This out-of-thermodynamics equilibrium effect does not alter the GW emission at the typical resolutions considered for mergers. Different microphysics treatments impact significantly mass, geometry and composition of the remnant's disc and ejecta. M1 simulations show systematically larger proton fractions. The different ejecta compositions reflect into the nucleosynthesis yields, that are robust only if both neutrino emission and absorption are simulated. Synthetic kilonova light curves calculated by means of spherically-symmetric radiation-hydrodynamics evolutions up to 15 days post-merger are mostly sensitive to ejecta's mass and composition; they can be reliably predicted only including the various ejecta components. We conclude that advanced microphysics in combination with resolutions higher than current standards appear essential for robust long-term evolutions and astrophysical predictions.

We develop a new method to infer the temporal, geometric, and energetic properties of galaxy outflows, by combining stellar spectral modeling to infer starburst ages, and absorption lines to measure velocities. If winds are accelerated with time during a starburst event, then these two measurements enable us to solve for the wind radius, similarly to length scales and the Hubble parameter in Big Bang cosmology. This wind radius is the vital, but hard-to-constrain parameter in wind physics. We demonstrate the method using spectra of 87 starburst galaxies at z=0.05-0.44, finding that winds accelerate throughout the starburst phase and grow to typical radii of ~1 kpc in ~10 Myr. Mass flow rates increase rapidly with time, and the mass-loading factor exceeds unity at about 10 Myr - while still being accelerated, the gas will likely unbind from the local potential and enrich the circumgalactic medium. We model the mechanical energy available from stellar winds and supernovae, and estimate that a negligible amount is accounted for in the cool outflow at early times. However, the energy deposition increases rapidly and ~10% of the budget is accounted for in the cool flow at 10 Myr, similar to some recent hydrodynamical simulations. We discuss how this model can be developed, especially for high-redshift galaxies.

Halo Models of large scale structure provide powerful and indispensable tools for phenomenological understanding of the clustering of matter in the Universe. While the halo model builds structures out of the superposition of haloes, defining halo profiles in their outskirts - beyond their virial radii - becomes increasingly ambiguous, as one cannot assign matter to individual haloes in a clear way. In this paper, we tackle this question by using numerical N-body simulations to find a systematic definition of mean halo profile that can be extended to large radii. These profiles must be compensated and are the key ingredients for the computation of cosmological correlation functions in the Amended Halo Model. The latter, introduced in our earlier work (arXiv:1912.04872), provides a more physically accurate phenomenological description of nonlinear structure formation, which respects conservation laws on large scales. Here, we extend this model from the matter auto-power spectrum to the halo-matter cross-power spectra, using N-body simulations. We find that the (dimensionless) compensated halo profile, $r^3 \times \rho(r)/M_{200c}$ has a near-universal maximum in the small range $0.03-0.04$ around the virial radius, $r \simeq r_{\rm 200c}$, independent of the halo mass. The profiles cross zero into negative values in the halo outskirts, beyond 2-3$\times r_{\rm 200c}$, consistent with our previous results. We provide preliminary fitting functions for compensated Navarro-Frenk-White (NFW) profiles, and this can be used to compute more physical observables in the Amended Halo Model of large scale structure.

We present the results of the broadband X-ray spectral analysis of simultaneous NuSTAR and XMM-Newton observations of four nearby Compton-thick active galactic nuclei (AGN) candidates selected from the Swift-Burst Alert Telescope (BAT) 150-month catalog. This work is part of a larger effort to identify and characterize all Compton-thick (NH >= 10^24 cm^-2) AGN in the local Universe (z < 0.05). We used three physically motivated models -- MYTorus, borus02, and UXClumpy -- to fit and characterize these sources. Of the four candidates analyzed, 2MASX J02051994-0233055 was found to be an unobscured (NH < 10^22 cm^-2) AGN, 2MASX J04075215-6116126 and IC 2227 to be Compton-thin (10^22 cm^-2 < NH < 10^24 cm^-2) AGN, and one, ESO 362-8, was confirmed to be a Compton-thick AGN. Additionally, every source was found to have a statistically significant difference between their line-of-sight and average torus hydrogen column density, further supporting the idea that the obscuring material in AGN is inhomogeneous. Furthermore, half of the sources in our sample (2MASX J02051994-0233055 and 2MASX J04075215-6116126) exhibited significant luminosity variation in the last decade, suggesting that this might be a common feature of AGN.

Fuzzy dark matter (FDM) is a proposed modification for the standard cold dark matter (CDM) model motivated by small-scale discrepancies in low-mass galaxies. Composed of ultra-light (mass $\sim 10^{-22}$ eV) axions with kpc-scale de Broglie wavelengths, this is one of a class of candidates that predicts that the first collapsed objects form in relatively massive dark matter halos. This implies that the formation history of the first stars and galaxies would be very different, potentially placing strong constraints on such models. Here we numerically simulate the formation of the first stars in an FDM cosmology, following the collapse in a representative volume all the way down to primordial protostar formation including a primordial non-equilibrium chemical network and cooling for the first time. We find two novel results: first, the large-scale collapse results in a very thin and flat gas ``pancake"; second, despite the very different cosmology, this pancake fragments until it forms protostellar objects indistinguishable from those in CDM. Combined, these results indicate that the first generation of stars in this model are also likely to be massive and, because of the sheet morphology, do not self-regulate, resulting in a massive Pop III starburst. We estimate the total number of first stars forming in this extended structure to be $10^4$ over 20 Myr using a simple model to account for the ionizing feedback from the stars, and should be observable with JWST. These predictions provide a potential smoking gun signature of FDM and similar dark matter candidates.

Magnetars are a special subset of the isolated neutron star family, with X-ray and radio emission mainly powered by the decay of their immense magnetic fields. Many attributes of magnetars remain poorly understood: spin-down glitches or the sudden reductions in the star's angular momentum, radio bursts reminiscent of extra-galactic Fast Radio Bursts (FRBs), and transient pulsed radio emission lasting months to years. Here we unveil the detection of a large spin-down glitch event ($|\Delta\nu/\nu| = 5.8_{-1.6}^{+2.6}\times10^{-6}$) from the magnetar SGR~1935+2154 on 2020 October 5 (+/- 1 day). We find no change to the source persistent surface thermal or magnetospheric X-ray behavior, nor is there evidence of strong X-ray bursting activity. Yet, in the subsequent days, the magnetar emitted three FRB-like radio bursts followed by a month long episode of pulsed radio emission. Given the rarity of spin-down glitches and radio signals from magnetars, their approximate synchronicity suggests an association, providing pivotal clues to their origin and triggering mechanisms, with ramifications to the broader magnetar and FRB populations. We postulate that impulsive crustal plasma shedding close to the magnetic pole generates a wind that combs out magnetic field lines, rapidly reducing the star's angular momentum, while temporarily altering the magnetospheric field geometry to permit the pair creation needed to precipitate radio emission.

Context. With the detection of thousands of exoplanets, characterising their dynamical evolution in detail represents a key step in the understanding of their formation. Studying the dissipation of tides occurring both in the host star and in the planets is of great relevance in order to investigate the distribution of the angular momentum occurring among the objects populating the system and to studying the evolution of the orbital parameters. From a theoretical point of view, the dissipation of tides throughout a body may be studied by relying on the so-called phase or time-lag equilibrium tides model in which the reduced tidal quality factor Q'p, or equivalently the product between the love number and the time lag (k2DeltaT), describe how efficiently tides are dissipated within the perturbed body. Constraining these factors by looking at the current configuration of the exoplanetary system is extremely challenging, and simulations accounting for the evolution of the system as a whole might help to shed some light on the mechanisms governing this process. Aims. We aim to constrain the tidal dissipation factors of hot-Jupiter-like planets by studying the orbital evolution of Kepler-91b. Methods. We firstly carried out a detailed asteroseismc characterisation of Kepler-91 and computed a dedicated stellar model using both classical and astereoseismic constraints. We then coupled the evolution of the star to the one of the planets by means of our orbital evolution code and studied the evolution of the system by accounting for tides dissipated both in the planet and in the host star. Results. We found that the maximum value for k2DeltaT (or equivalently the minimum value for Q'p) determining the efficiency of equilibrium tides dissipation occurring within Kepler-91b is 0.4 pm 0.25 s (4.5+5.8 * 10^5).

While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $\delta\rho/\rho$, pressure, $\delta P/P$, X-ray surface brightness, Sunyaev-Zeldovich (SZ) y-parameter, and the characteristic Mach number of gas motions, $M_{\rm 1d}$. In relaxed clusters, accounting for halo ellipticities reduces $\delta\rho/\rho$ or $\delta P/P$ by a factor of up to 2 within $r_{500c}$. We confirm a strong linear correlation between $\delta\rho/\rho$ (or $\delta P/P$) and $M_{\rm 1d}$ in relaxed clusters, with the proportionality coefficient $\eta \approx 1$. For unrelaxed clusters, the correlation is less strong and has a larger $\eta\approx 1.3\pm 0.5$ ($1.5\pm0.5$) for $\delta\rho/\rho$ ($\delta P/P$). Examination of the power-law scaling of $M_{\rm 1d}$ with $\delta\rho/\rho$ shows that it is almost linear for relaxed clusters, while for the unrelaxed ones, it is closer to $\delta\rho/\rho\propto M_{\rm 1d}^2$, supporting an increasing role of non-linear terms and compressive modes. In agreement with previous studies, we observe a strong correlation of $M_{\rm 1d}$ with radius. Correcting for these correlations leaves a residual scatter in $M_{\rm 1d}$ of $\sim 4 (7)$ per cent for relaxed (perturbed) clusters. Hydrostatic mass bias correlates with $M_{\rm 1d}$ as strongly as with $\delta\rho/\rho$ in relaxed clusters. The residual scatters after correcting for derived trends is $\sim 6-7$ per cent. These predictions can be verified with existing X-ray and SZ observations of galaxy clusters combined with forthcoming velocity measurements with X-ray microcalorimeters.

The neutral hydrogen 21 cm line is an excellent tracer of the atomic interstellar medium in the cold and the warm phases. Combined 21 cm emission and absorption observations are very useful to study the properties of the gas over a wide range of density and temperature. In this work, we have used 21 cm absorption spectra from recent interferometric surveys, along with the corresponding emission spectra from earlier single dish surveys to study the properties of the atomic gas in the Milky Way. In particular, we focus on a comparison of properties between lines of sight through the gas disk in the Galactic plane and high Galactic latitude lines of sight through more diffuse gas. As expected, the analysis shows a lower average temperature for the gas in the Galactic plane compared to that along the high latitude lines of sight. The gas in the plane also has a higher molecular fraction, showing a sharp transition and flattening in the dust - gas correlation. On the other hand, the observed correlation between 21 cm brightness temperature and optical depth indicates some intrinsic difference in spin temperature distribution and a fraction of gas in the Galactic plane having intermediate optical depth (for 0.02 < $\tau$ < 0.2) but higher spin temperature, compared to that of the diffuse gas at high latitude with the same optical depth. This may be due to a small fraction of cold gas with slightly higher temperature and lower density present on the Galactic plane.

We analyze spectra of a gravitationally lensed galaxy, known as the Sunburst Arc, that is leaking ionizing photons, also known as the Lyman continuum (LyC). Magnification from gravitational lensing permits the galaxy to be spatially resolved into one region that leaks ionizing photons, and several that do not. Rest-frame ultraviolet and optical spectra from Magellan target ten different regions along the lensed Arc, including six multiple images of the LyC leaking region, as well as four regions that do not show LyC emission. The rest-frame optical spectra of the ionizing photon emitting regions reveal a blue-shifted ($\Delta V$=27 km s$^{-1}$) broad emission component (FWHM=327 km s$^{-1}$) comprising 55% of the total [OIII] line flux, in addition to a narrow component (FWHM = 112 km s$^{-1}$), suggesting the presence of strong highly ionized gas outflows. This is consistent with the high-velocity ionized outflow inferred from the rest-frame UV spectra. In contrast, the broad emission component is less prominent in the non-leaking regions, comprising $\sim$26% of total [OIII] line flux. The high ionization absorption lines are prominent in both leaker and non-leaker but low ionization absorption lines are very weak in the leaker, suggesting that the line of sight gas is highly ionized in the leaker. Analyses of stellar wind features reveal that the stellar population of the LyC leaking regions is considerably younger ($\sim$3 Myr) than the non-leaking regions ($\sim$12 Myr), highlighting that stellar feedback from young stars may play an important role in ionizing photon escape.

Using ion measurements from Ultra-Low-Energy Isotope Spectrometer (ULEIS) observations onboard Advanced Composition Explorer (ACE) and Solar Isotope Spectrometer (SIS) observations onboard the Solar Terrestrial Observatory (STEREO)-A and STEREO-B spacecraft, we have identified 854 3He-rich time periods between 1997 September and 2021 March. We include all event types with observed 3He enhancements such as corotating interaction regions (CIRs), gradual solar energetic particle (SEP) events, interplanetary shocks, and impulsive SEP events. We employ two different mass separation techniques to obtain 3He, 4He, Fe, and O fluences for each event, and we determine the 3He/4He and Fe/O abundance ratios between 0.32 to 0.45 MeV/nucleon and 0.64 to 1.28 MeV/nucleon. We find a clear correlation in the 3He/4He and Fe/O abundance ratios between both energy ranges. We find two distinct trends in the 3He/4He vs. Fe/O relation. For low 3He/4He values, there is a positive linear correlation between 3He/4He and Fe/O. However, at 3He/4He ~ 0.3, Fe/O appears to reach a limit and the correlation weakens significantly. We provide a live catalogue of 3He rich time periods that includes the robust determination of the onset and end times of the 3He enhancements in SEP-associated periods for different types of events observed my multiple spacecraft. This catalogue is available for public use. New releases will follow after major additions such as adding new periods from new missions (e.g., Parker Solar Probe and Solar Orbiter), identifying event types (impulsive SEP events, etc.), or adding new parameters such as remote observations detailing characteristics of the active regions.

The spectacular outbursts of energy associated with supernovae (SNe) have long motivated research into their potentially hazardous effects on Earth and analogous environments. Much of this research has focused primarily on the atmospheric damage associated with the prompt arrival of ionizing photons within days or months of the initial outburst, and the high-energy cosmic rays that arrive thousands of years after the explosion. In this study, we turn the focus to persistent X-ray emission, arising in certain SNe that have interactions with a dense circumstellar medium, and observed months and/or years after the initial outburst. The sustained high X-ray luminosity leads to large doses of ionizing radiation out to formidable distances. We provide an assessment of the threat posed by these X-ray luminous SNe by analyzing the collective X-ray observations from Chandra, Swift-XRT, XMM-Newton, NuSTAR, and others. We find that this threat is particularly acute for SNe showing evidence of strong circumstellar interaction, such as Type IIn explosions, which have significantly larger ranges of influence than previously expected, and lethal consequences up to $\sim$ 50 pc away. Furthermore, X-ray bright SNe could pose a substantial and distinct threat to terrestrial biospheres, and tighten the Galactic habitable zone. We urge follow-up X-ray observations of interacting SNe for months and years after the explosion to shed light on the physical nature of the emission and its full time evolution, and to clarify the danger that these events pose for life in our Galaxy and other star-forming regions.

We simulate the formation and collapse of prestellar cores at few-AU resolution in a set of radiation-magneto-hydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming-in onto individual massive prestellar cores within the GMC. We identify two distinct modes of fragmentation: "quasi-spherical" and "filamentary". In both modes the fragments eventually become embedded in a quasi-steady accretion disk or toroid with radii ~ 500-5000 AU and opening angles $H/R \sim 0.5-1$. The disks/toroids are Toomre stable but the accreted pre-existing fragments are found orbiting the outer disk, appearing as disk fragmentation. Each core converts nearly 100 percent of the gas mass into a few massive stars forming near the disk center. Large and massive disks around high-mass stars are supported by magnetic pressure in the outer disk, at radii >200-1000 AU, and turbulent pressure in the inner disk. The most massive core accretes several times more mass than its initial mass, forming a (proto)star cluster of 8 massive stars enshrouded by a toroid, suggesting a competitive accretion scenario for ultra-high-mass star formation. We also find that the HII regions produced by a single massive star remain trapped in the dense circumstellar disks for a few hundred kiloyears, while the dynamic motions of massive stars in wide binaries or multiple systems displace the stars from the densest parts of the disk, allowing UV radiation to escape producing steady or pulsating bipolar HII regions.

Hierarchical Bayesian inference is an essential tool for studying the population properties of compact binaries with gravitational waves. The basic premise is to infer the unknown prior distribution of binary black hole and/or neutron star parameters such component masses, spin vectors, and redshift. These distributions shed light on the fate of massive stars, how and where binaries are assembled, and the evolution of the Universe over cosmic time. Hierarchical analyses model the binary black hole population using a prior distribution conditioned on hyper-parameters, which are inferred from the data. However, a misspecified model can lead to faulty astrophysical inferences. In this paper we answer the question: given some data, which prior distribution -- from the set of all possible prior distributions -- produces the largest possible population likelihood? This distribution (which is not a true prior) is $\pistroke$ (pronounced ``pi stroke''), and the associated \textit{maximum population likelihood} is $\Lstroke$ (pronounced ``L stroke''). We postulate that the structure of $\pistroke$ is a linear superposition of delta functions. We show how $\pistroke$ and $\Lstroke$ can be used for model exploration/criticism. We apply this $\Lstroke$ formalism to study the population of binary black hole mergers observed in LIGO--Virgo--KAGRA's third Gravitational-Wave Transient Catalog. Based on our results, we discuss possible improvements for gravitational-wave population models.

We present analyses of improved photometric and spectroscopic observations for two detached eclipsing binaries at the turnoff of the open cluster NGC 752: the 1.01 day binary DS And and the 15.53 d BD $+$37 410. For DS And, we find $M_1 = 1.692\pm0.004\pm0.010 M_\odot$, $R_1 = 2.185\pm0.004\pm0.008 R_\odot$, $M_2 = 1.184\pm0.001\pm0.003 M_\odot$, and $R_2 = 1.200\pm0.003\pm0.005 R_\odot$. We either confirm or newly identify unusual characteristics of both stars in the binary: the primary star is found to be slightly hotter than the main sequence turn off and there is a more substantial discrepancy in its luminosity compared to models (model luminosities are too large by about 40%), while the secondary star is oversized and cooler compared to other main sequence stars in the same cluster. The evidence points to non-standard evolution for both stars, but most plausible paths cannot explain the low luminosity of the primary star. BD $+$37 410 only has one eclipse per cycle, but extensive spectroscopic observations and the TESS light curve constrain the stellar masses well: $M_1 = 1.717\pm0.011 M_\odot$ and $M_2 = 1.175\pm0.005 M_\odot$. The radius of the main sequence primary star near $2.9R_\odot$ definitively requires large convective core overshooting ($> 0.2$ pressure scale heights) in models for its mass, and multiple lines of evidence point toward an age of $1.61\pm0.03\pm0.05$ Gyr (statistical and systematic uncertainties). Because NGC 752 is currently undergoing the transition from non-degenerate to degenerate He ignition of its red clump stars, BD $+$37 410 A directly constrains the star mass where this transition occurs.

A recent study shows that gravitational scattering of dark matter, in the form of massive objects with mass $m \sim 10^3-10^4M_{\odot}$, could provide a possible solution to alleviate the small-scale structure problems of cold dark matter. The scattering cross section is velocity-dependent so that this scenario can explain why self-interaction of dark matter is significant in dwarf galaxies, but not in massive galaxies and galaxy clusters. In this Letter, we show that this kind of dark massive objects could be made of sterile neutrinos with a possible rest mass range $m_{\nu} \sim 7.6$ keV $-$ 71 MeV. This mass range generally satisfies most of the current observational constraints. The entire structure of the sterile neutrino halos can be simply predicted from standard physics.

We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope SZ survey (SPT-SZ) and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array (ALMA). This work builds on the previous average rest-frame spectrum, given in Spilker et al. 2014 for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective dust masses. The spectrum spans $1.9$$<$z$<$$6.9$ and covers rest-frame frequencies of 240$-$800 GHz. We detect multiple bright line features from $^{12}$CO, $[$CI$]$, and H$_2$O, as well as fainter molecular transitions from $^{13}$CO, HCN, HCO$^+$, HNC, CN, and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright $^{12}$CO emission lines, and contain either warmer and more excited dense gas tracers, or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L$_{\mathrm{IR}}$$>$$10^{12}$L$_\odot$), and will inform studies of chemical evolution before the peak epoch of star formation at $z=2-3$.

Planet engulfment can be inferred from enhancement of refractory elements in the photosphere of the engulfing star following accretion of rocky planetary material. Such refractory enrichments are subject to stellar interior mixing processes, namely thermohaline mixing induced by an inverse mean-molecular-weight gradient between the convective envelope and radiative core. Using MESA stellar models, we quantified the strength and duration of engulfment signatures following planet engulfment. We found that thermohaline mixing dominates during the first $\sim$5$-$45 Myr post-engulfment, weakening signatures by a factor of $\sim$2 before giving way to depletion via gravitational settling on longer timescales. Solar metallicity stars in the 0.5-1.2 $M_{\odot}$ mass range have observable signature timescales of $\sim$1 Myr$-$8 Gyr, depending on the engulfing star mass and amount of material engulfed. Early type stars exhibit larger initial refractory enhancements but more rapid depletion. Solar-like stars ($M$ = 0.9$-$1.1 $M_{\odot}$) maintain observable signatures ($>$0.05 dex) over timescales of $\sim$20 Myr$-$1.7 Gyr for nominal 10 $M_{\oplus}$ engulfment events, with longer-lived signatures occurring for low-metallicity and/or hotter stars (1 $M_{\odot}$, $\sim$2$-$3 Gyr). Engulfment events occurring well after the zero-age main sequence produce larger signals due to suppression of thermohaline mixing by gravitational settling of helium (1 $M_{\odot}$, $\sim$1.5 Gyr). These results indicate that it may be difficult to observe engulfment signatures in solar-like stars that are several Gyr old.

Strong gravitational lenses allow us to peer into the farthest reaches of space by bending the light from a background object around a massive object in the foreground. Unfortunately, these lenses are extremely rare, and manually finding them in astronomy surveys is difficult and time-consuming. We are thus tasked with finding them in an automated fashion with few if any, known lenses to form positive samples. To assist us with training, we can simulate realistic lenses within our survey images to form positive samples. Naively training a ResNet model with these simulated lenses results in a poor precision for the desired high recall, because the simulations contain artifacts that are learned by the model. In this work, we develop a lens detection method that combines simulation, data augmentation, semi-supervised learning, and GANs to improve this performance by an order of magnitude. We perform ablation studies and examine how performance scales with the number of non-lenses and simulated lenses. These findings allow researchers to go into a survey mostly ``blind" and still classify strong gravitational lenses with high precision and recall.

We propose a novel subhalo abundance matching (SHAM) model that uses the virial mass of the main progenitor of each (sub)halo $M_{\rm prog}$ as a proxy of the galaxy stellar mass $M_*$ at the time of observation. This $M_{\rm prog}$ model predicts the two-point correlation functions depending on the choice of the epoch $z_{\rm prog}$ at which $M_\mathrm{prog}$ is measured. With $z_{\rm prog}$ as a fitting parameter, we apply the $M_{\rm prog}$ model to the latest observed angular correlation functions of galaxy samples at $z\simeq0.4$ with varying stellar mass thresholds from $M_{*,~{\rm lim}}/(h^{-2}M_\odot)=10^{11}$ to $10^{8.6}$. The $M_{\rm prog}$ model can reproduce the observations very well over the entire scale range of $0.1\textrm{--}10~h^{-1}{\rm Mpc}$. We find that, for the samples of $10^9\leq M_{*,~{\rm lim}}/(h^{-2}M_\odot)\leq10^{10.2}$, the correlation functions predicted by the widely-used $V_{\rm peak}$ model lacks amplitudes at the scale of $\lesssim1~h^{-1}{\rm Mpc}$, demonstrating the high capability of our $M_\mathrm{prog}$ model to explain observed clustering measurements. The best-fit $z_{\rm prog}$ parameter is highest ($z_{\rm prog}\simeq5$) for intermediate mass galaxies at $M_*\simeq10^{9.7}~h^{-2}M_\odot$, and becomes smaller towards $z_\mathrm{prog}\simeq1$ for both lower- and higher-mass galaxies. We interpret these trends as reflecting the downsizing in the in-situ star formation in lower-mass galaxies and the larger contribution of ex-situ stellar mass growth in higher-mass galaxies.

It is commonly accepted that radio-loud active galactic nuclei are hosted exclusively by giant elliptical galaxies. We analyze high-resolution optical Hubble Space Telescope images of a sample of radio galaxies with extended double-lobed structures associated with disk-like optical counterparts. After systematically evaluating the probability of chance alignment between the radio lobes and the optical counterparts, we obtain a sample of 18 objects likely to have genuine associations. The host galaxies have unambiguous late-type morphologies, including spiral arms, large-scale dust lanes among the edge-on systems, and exceptionally weak bulges, as judged by the low global concentrations, small global S\'{e}rsic indices, and low bulge-to-total light ratios (median $B/T = 0.13$). With a median S\'{e}rsic index of 1.4 and low effective surface brightnesses, the bulges are consistent with being pseudo bulges. The majority of the hosts have unusually large stellar masses (median $M_* = 1.3\times 10^{11}\, M_\odot$) and red optical colors (median $g-r = 0.69\,$mag), consistent with massive, quiescent galaxies on the red sequence. We suggest that black hole mass (stellar mass) plays a fundamental role in launching large-scale radio jets, and that the rarity of extended radio lobes in late-type galaxies is the consequence of the steep stellar mass function at the high-mass end. The disk radio galaxies have mostly Fanaroff-Riley type II morphologies yet lower radio power than sources of a similar type traditionally hosted by ellipticals. The radio jets show no preferential alignment with the minor axis of the galactic bulge or disk, apart from a possible mild tendency for alignment among the most disk-dominated systems.

Contemporary radiation-magnetohydrodynamic simulation of the AGNs predicts presence of the hot and strong accretion disk wind, which gets unstable far from the central region and turns into gas clumps. These inner-wind and outer clumps may be actually observed as the ultrafast outflows (UFOs) and the clumpy absorbers, respectively. We may call this picture as the "hot inner and clumpy outer wind model". Observationally, it is challenging to place constraints on the origin of the UFOs and clumpy absorbers due to complicated spectral variations. We developed a novel method, "spectral-ratio model fitting", to resolve parameter degeneracy of the clumpy absorbers and other spectral components. In this method, the parameters of the absorber in the line of sight are estimated from the ratio of the partially absorbed spectrum to the non-absorbed one. We applied this method to the narrow-line Seyfert 1 galaxy IRAS 13224-3809 observed by XMM-Newton in 2016 for 1.5 Ms, where the source showed extreme spectral variability and complex absorption features. As a result, we found that the soft spectral variation is mostly explained by a change of the covering fraction of the mildly-ionized clumpy absorbers, and that these absorbers are outflowing with such a high velocity that is comparable to that of the UFO (~ 0.2-0.3 c). This result implies that the formation of the clumpy absorbers and the UFO shares the same origin, supporting the "hot inner and clumpy outer wind model".

We present optical photometric and spectroscopic analysis of a Type Iax supernova (SN) 2020rea situated at the brighter luminosity end of Type Iax supernovae (SNe). The light curve decline rate of SN~2020rea is $\Delta$m$_{15}$(g) = 1.31$\pm$0.08 mag which is similar to SNe 2012Z and 2005hk. Modelling the pseudo bolometric light curve with a radiation diffusion model yields a mass of $^{56}$Ni of 0.13$\pm$0.01 M$_{\odot}$ and an ejecta mass of 0.77$^{+0.11}_{-0.21}$ M$_{\odot}$. Spectral features of SN~2020rea during the photospheric phase show good resemblance with SN 2012Z. TARDIS modelling of the early spectra of SN~2020rea reveals a dominance of Iron Group Elements (IGEs). The photospheric velocity of the Si {\sc II} line around maximum for SN~2020rea is $\sim$ 6500 km s$^{-1}$ which is less than the measured velocity of the Fe {\sc II} line and indicates significant mixing. The observed physical properties of SN~2020rea match with the predictions of pure deflagration model of a Chandrasekhar mass C-O white dwarf. The metallicity of the host galaxy around the SN region is 12+log(O/H) = 8.56$\pm$0.18 dex which is similar to that of SN 2012Z.

Hot Jupiters (HJs) are giant planets with orbital periods of the order of a few days with semimajor axis within $\sim$0.1 au. Several theories have been invoked in order to explain the origin of this type of planets, one of them being the high-eccentricity migration. This migration can occur through different high-eccentricity mechanisms. Our investigation focused on six different kinds of high-eccentricity mechanisms, namely, direct dispersion, coplanar, Kozai-Lidov, secular chaos, E1 and E2 mechanisms. We investigated the efficiency of these mechanisms for the production of HJ candidates in multi-planet systems initially tightly-packed in the semimajor axis, considering a large set of numerical simulations of the exact equations of motion in the context of the N-body problem. In particular, we analyzed the sensitivity of our results to the initial number of planets, the initial semimajor axis of the innermost planetary orbit, the initial configuration of planetary masses, and to the inclusion of general relativity effects. We found that the E1 mechanism is the most efficient in producing HJ candidates both in simulations with and without the contribution of general relativity, followed by the Kozai-Lidov and E2 mechanisms. Our results also revealed that, except for the initial equal planetary mass configuration, the E1 mechanism was notably efficient in the other initial planetary mass configurations considered in this work. Finally, we investigated the production of HJ candidates with prograde, retrograde, and alternating orbits. According to our statistical analysis, the Kozai-Lidov mechanism has the highest probability of significantly exciting the orbital inclinations of the HJ candidates.

Observations show that individual protoplanetary disk lifetimes vary from \mbox{$<$ 1 Myr} to $\gg$ 20 Myr. The disk lifetime distribution is currently unknown. For the example of a Gaussian distribution of the disk lifetime, I suggest a simple method for deducing such a disk lifetimes distribution. The median disk lifetimes inferred with this method is also shown.

A linear Ogorodnikov-Milne model is applied to study the three-dimensional kinematics of classical Cepheids in the Milky Way. A sample of 832 classical Cepheids from Mr'oz et al. (2019) with distances, line-of-sight velocities, and proper motions from the Gaia DR2 catalogue is used. The Cepheid space velocities have been freed from the differential Galactic rotation found by us previously based on a nonlinear rotation model. Based on a complete Ogorodnikov-Milne model, involving the line-of-sight velocities and proper motions of stars, we have estimated the angular velocity of rotation around the Galactic $y$ axis, $\Omega_y=0.64\pm0.17$~km s$^{-1}$ kpc$^{-1}$. We think that this rotation is associated with the warp of the Galactic thin disk. Our calculations using only the proper motions of Cepheids under the assumption of no deformations due to the disk warp have shown the presence of a residual rotation around the $y$ axis with an angular velocity $\Omega_y=0.54\pm0.15$~km s$^{-1}$ kpc$^{-1}$ and the presence of a positive rotation around the $x$ axis with an angular velocity $\Omega_x=0.33\pm0.10$~km s$^{-1}$ kpc$^{-1}$.

Nearby M dwarf systems currently offer the most favorable opportunities for spectroscopic investigations of terrestrial exoplanet atmospheres. The LTT~1445 system is a hierarchical triple of M dwarfs with two known planets orbiting the primary star, LTT~1445A. We observe four transits of the terrestrial world LTT~1445Ab ($R=1.3$ R$_\oplus$, $M=2.9$ M$_\oplus$) at low resolution with Magellan II/LDSS3C. We use the combined flux of the LTT~1445BC pair as a comparison star, marking the first time that an M dwarf is used to remove telluric variability from time-series observations of another M dwarf. We find H$\alpha$ in emission from both LTT~1445B and C, as well as a flare in one of the data sets from LTT~1445C. These contaminated data are removed from the analysis. We construct a broadband transit light curve of LTT~1445Ab from 620--1020 nm. Binned to 3-minute time bins we achieve an rms of 43 ppm for the combined broadband light curve. We construct a transmission spectrum with 20 spectrophotometric bins each spanning 20 nm and compare it to models of clear, 1$\times$ solar composition atmospheres. We rule out this atmospheric case with a surface pressure of 10 bars to $3.4 \sigma$ confidence, and with a surface pressure of 1 bar to $2.9 \sigma$ confidence. Upcoming secondary eclipse observations of LTT~1445Ab with JWST will further probe the cases of a high mean molecular weight atmosphere, a hazy or cloudy atmosphere, or no atmosphere at all on this terrestrial world.

We present core-collapse supernovae simulations including nuclear reaction networks which impact explosion dynamics and nucleosynthesis. The different composition treatment can lead to changes in the neutrino heating in the vicinity of the shock, by modifying the amount of nucleons and thus the $\mathrm{\nu}$-opacity of the region. This reduces the ram pressure outside the shock and allows an easier expansion. The energy released by the nuclear reactions during collapse also slows down the accretion, and aids the shock expansion. In addition, nuclear energy generation in the post-shocked matter produces more energetic explosions, up to $20\,\%$. Nucleosynthesis is affected due to the different dynamic evolution of the explosion. Our results indicate that the energy generation from nuclear reactions helps to sustain late outflows from the vicinity of the proto-neutron star (PNS), synthesizing more neutron-rich species. Furthermore, we show that there are systematic discrepancies between the ejecta calculated with in-situ and ex-situ reaction networks. The mass fractions of some Ca, Ti, Cr, and Fe isotopes are consistently under-produced in post-processing calculations, leading to different nucleosynthesis paths. Therefore, large in-situ nuclear reaction networks are needed for a more accurate nucleosynthesis.

We study the formation and evolution of elliptical galaxies and how they suppress star formation and maintain it quenched. A one-zone chemical model which follows in detail the time evolution of gas mass and its chemical abundances during the active and passive evolution, is adopted. The model includes both gas infall and outflow as well as detailed stellar nucleosynthesis. Elliptical galaxies with different infall masses, following a down-sizing in star formation scenario, are considered. In the chemical evolution simulation we include a novel calculation of the feedback processes. We include heating by stellar wind, core-collapse SNe, Type Ia SNe (usually not highlighted in galaxy formation simulations) and AGN feedback. The AGN feedback is a novelty in this kind of models and is computed by considering a Bondi-Eddington limited accretion onto the central supermassive black hole. We successfully reproduce several observational features, such as the [$\alpha$/Fe] ratios increasing with galaxy mass, mass-metallicity, $\rm M_{BH}-\sigma$ and $\rm M_{BH}-M_{*}$ relations. Moreover, we show that stellar feedback and in particular Type Ia SNe, has a main role in maintaining quenched the star formation after the occurrence of the main galactic wind, especially in low-mass ellipticals. For larger systems, the contribution from AGN to thermal energy of gas appears to be necessary. However, the effect of the AGN on the development of the main galactic wind is negligible, unless an unreasonable high AGN efficiency or an extremely low stellar feedback are assumed. We emphasize the important role played by Type Ia SNe in the energy budget of early-type galaxies.

Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescope of the recently discovered extreme debris disc around TYC 4209-1322-1. By combining these observations with other time-domain optical and mid-infrared data, we explore the disc variability of the last four decades with particular emphasis on the last 12 yr. During the latter interval the disc showed substantial changes, the most significant was the brightening and subsequent fading between 2014 and 2018 as outlined in WISE data. The Spitzer light curves outline the fading phase and a subsequent new brightening of the disc after 2018, revealing an additional flux modulation with a period of ~39 days on top of the long-term trend. We found that all these variations can be interpreted as the outcome of a giant collision that happened at an orbital radius of ~0.3 au sometime in 2014. Our analysis implies that a collision on a similar scale could have taken place around 2010, too. The fact that the disc was already peculiarly dust rich 40 yr ago, as implied by IRAS data, suggests that these dust production events belong to a chain of large impacts triggered by an earlier even more catastrophic collision.

NGC~6946, also known as the `Fireworks' galaxy, is an unusual galaxy that hosts a total of 225 supernova remnant (SNR) candidates, including 147 optically identified with high [SII]/Ha line ratios. In addition, this galaxy shows prominent HI holes, which were analyzed in previous studies. Indeed, the connection between SNRs and HI holes together with their physical implications in the surrounding gas is worth of attention. This paper explores the connection between the SNRs and the HI holes, including an analysis of their physical link to observational optical properties inside and around the rims of the holes, using new integral field unit (IFU) data from the Metal-THINGS survey. We present an analysis combining previously identified HI holes, SNRs candidates, and new integral field unit (IFU) data from Metal-THINGS of the spiral galaxy NGC 6946. We analyze the distributions of the oxygen abundance, star formation rate surface density, extinction, ionization, diffuse ionized gas, and the Baldwin-Phillips-Terlevich classification throughout the galaxy. By analyzing in detail the optical properties of the 121 previously identify HI holes in NGC 6946, we find that the SNRs are concentrated at the rims of the HI holes. Furthermore, our IFU data shows that the star formation rate and extinction are enhanced at the rims of the holes. To a lesser degree, the oxygen abundance and ionization parameter show hints of enhancement on the rims of the holes. Altogether, this provides evidence of induced star formation taking place at the rims of the holes, whose origin can be explained by the expansion of superbubbles created by multiple supernova explosions in large stellar clusters dozens of Myr ago.

Impact-induced erosion of the Earth's early crust during accretion of terrestrial bodies can significantly modify the primordial chemical composition of the Bulk Silicate Earth (BSE, that is, the composition of the crust added to the present-day mantle). In particular, it can be particularly efficient in altering the abundances of elements having a strong affinity for silicate melts (i.e. incompatible elements) as the early differentiated crust was preferentially enriched in those. Here, we further develop an erosion model (EROD) to quantify the effects of collisional erosion on the final composition of the BSE. Results are compared to the present-day BSE composition models and constraints on Earth's accretion processes are provided. The evolution of the BSE chemical composition resulting from crustal stripping is computed for entire accretion histories of about 50 Earth analogs in the context of the Grand Tack model. The chosen chemical elements span a wide range of incompatibility degrees. We find that a maximum loss of 40wt% can be expected for the most incompatible lithophile elements such as Rb, Th or U in the BSE when the crust is formed from low partial melting rates. Accordingly, depending on both the exact nature of the crust-forming processes during accretion and the accretion history itself, Refractory Lithophile Elements (RLE) may not be in chondritic relative proportions in the BSE. In that case, current BSE estimates may need to be corrected as a function of the geochemical incompatibility of these elements. Alternatively, if RLE are indeed in chondritic relative proportions in the BSE, accretion scenarios that are efficient in affecting the BSE chemical composition should be questioned.

We investigate the star formation history of the IC1396 region by studying its kinematics and completing the population census. We use multiwavelength data, combining optical spectroscopy (to identify and classify new members), near-infrared photometry (to trace shocks, jets, and outflows and the interactions between the cluster members and the cloud), along with Gaia EDR3 to identify new potential members in the multidimensional proper motion/parallax space. The revised Gaia EDR3 distance is 925$\pm$73 pc, slightly closer than previously obtained with DR2. The Gaia data reveal four distinct subclusters in the region. These subclusters are consistent in distance but display differences in proper motion. This, with their age differences, hints towards a complex and varied star formation history. Gaia data also unveil the intermediate-mass objects that tend to evade spectroscopic and disk surveys. Our analysis allows us to identify 334 new members. We estimate an average age of $\sim$4 Myr, confirming previous age estimates. With the new members added to our study, we estimate a disk fraction of 28\%, lower than previous values, due to our method detecting mainly new, diskless intermediate-mass stars. We find age differences between the subclusters, which evidences a complex star formation history with different episodes of star formation.

The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-modality which quickly dissolves and only shows in young stars. Unlike other explanations, such as a 1:1 orbital resonance with the Galactic bar, this also accounts for the observed steepening of the stellar warp at $L_\mathrm{gap}$, since the adjacent spiral arms represent different parts of the warped gas disc. In this scenario the gap is clearly present only in young stars, as observed, while most purely stellar dynamical origins would affect all disc populations, including older disc stars.

We present the first set of rate coefficients for the rotational excitation of the 7 lowest levels of hydrogen fluoride (HF) induced by collision with water molecules, the dominant collider in cometary comas, in the 5-150 K temperature range. The calculations are performed with a quantum statistical approach from an accurate rigid rotor ab initio interaction potential. Rate coefficients for excitation of HF by electron-impact are also computed, within the Born approximation, in the 10-10,000 K temperature range. These rate coefficients are then used in a simplified non-local thermodynamic equilibrium (non-LTE) model of a cometary coma that also includes solar radiative pumping and radiative decay. We investigate the range of H2O densities that lead to non-LTE populations of the rotational levels of HF. We show that to describe the excitation of HF in comets, considering collisions with both water molecules and electrons is needed as a result of the large dipole of HF.

We use the SOUX sample of $\sim$700 AGN to form average optical-UV-X-rays SEDs on a 2D grid of $M_{\mathrm{BH}}$ and $L_{2500}$. We compare these with the predictions of a new AGN SED model, QSOSED, which includes prescriptions for both hot and warm Comptonisation regions as well as an outer standard disc. This predicts the overall SED fairly well for 7.5<log($M_{\mathrm{BH}}/M_{\mathrm{\odot}}$)<9.0 over a wide range in $L/L_{\mathrm{Edd}}$, but at higher masses the outer disc spectra in the model are far too cool to match the data. We create optical-UV composites from the entire SDSS sample and use these to show that the mismatch is due to there being no significant change in spectral shape of the optical-UV continuum across several decades of $M_{\mathrm{BH}}$ at constant luminosity. We show for the first time that this cannot be matched by standard disc models with high black hole spin. These apparently fit, but are not self-consistent as they do not include the General Relativistic effects for the emission to reach the observer. At high spin, increased gravitational redshift compensates for almost all of the higher temperature emission from the smaller inner disc radii. The data do not match the predictions made by any current accretion flow model. Either the disc is completely covered by a warm Comptonisation layer whose properties change systematically with $L/L_{\mathrm{Edd}}$, or the accretion flow structure is fundamentally different to that of the standard disc models.

The atmospheres of ultra hot Jupiters (UHJs) are prime targets for the detection of molecules and atoms at both low and high spectral resolution. We study the atmospheres of the UHJs WASP-121b and WASP-189b by performing 3D general circulation models (GCMs) of these planets using high temperature correlated-k opacity schemes with ultra-violet (UV) absorbing species included. The GCM results are then post-processed at low and high spectral resolutions and compared to available data. The high resolution results are cross-correlated with molecular and atomic templates to produce mock molecular detections. Our GCM models produce similar temperature-pressure (T-p) structure trends to previous 1D radiative-convective equilibrium models of UHJs. Furthermore, the inclusion of UV opacities greatly shapes the thermal and dynamical properties of the high-altitude, low-pressure regions of the UHJ atmospheres, with sharp T-p inversions due to the absorption of UV light. This suggests that optical wavelength, high-resolution observations probe a dynamically distinct upper atmospheric region, rather than the deeper jet forming layers.

Some icy moons and small bodies in the solar system are believed to host subsurface liquid water oceans. The interaction of these saline, electrically conductive oceans with time-varying external magnetic fields generates induced magnetic fields. Magnetometry observations of these induced fields in turn enable the detection and characterization of these oceans. We present a framework for characterizing the interiors of icy moons using multi-frequency induction and Bayesian inference applied to magnetometry measurements anticipated from the upcoming Europa Clipper mission. Using simulated data from the Europa Clipper Magnetometer (ECM), our approach can accurately retrieve a wide range of plausible internal structures for Europa. In particular, the ocean conductivity is recovered to within ${\pm}50\%$ for all internal structure scenarios considered and the ocean thickness can be retrieved to within ${\pm}25~\mathrm{km}$ for five out of seven scenarios. Characterization of the ice shell thickness to ${\pm}50\%$ is possible for six of seven scenarios. Our recovery of the ice shell thickness is highly contingent on accurate modeling of magnetic fields arising from the interaction of Europa with the ambient magnetospheric plasma, while the ocean thickness is more modestly affected and the ocean conductivity retrieval is largely unchanged. Furthermore, we find that the addition of a priori constraints (e.g., static gravity measurements) can yield improved ocean characterization compared to magnetometry alone, suggesting that multi-instrument techniques can play a key role in revealing the interiors of Europa and other ocean worlds.

The SVOM satellite, to be launched at the end of 2023, is primarily devoted to the multi-wavelength observation of gamma-ray bursts and other higher-energy transients. Thanks to its onboard Microchannel X-ray Telescope and Visible-band Telescope, it is also very well adapted to the electromagnetic follow-up of gravitational wave events. We discuss the SVOM follow-up strategy for gravitational wave trigger candidates provided by LIGO-Virgo-KAGRA. In particular, we make use of recent developments of galaxy catalogs adapted to the horizon of gravitational wave detectors to optimise the chance of counterpart discovery. We also take into account constraints specific to the SVOM platform. Finally, we implement the production of the SVOM observation plan following a gravitational wave alert and quantify the efficiency of several optimisations introduced in this work.

Dynamical evolution within planetary systems can cause planets to be engulfed by their host stars. Following engulfment, the stellar photosphere abundance pattern will reflect accretion of rocky material from planets. Multi-star systems are excellent environments to search for such abundance trends because stellar companions form from the same natal gas cloud and are thus expected to share primordial chemical compositions to within 0.03$-$0.05 dex. Abundance measurements have occasionally yielded rocky enhancements, but few observations targeted known planetary systems. To address this gap, we carried out a Keck-HIRES survey of 36 multi-star systems where at least one star is a known planet host. We found that only HAT-P-4 exhibits an abundance pattern suggestive of engulfment, but is more likely primordial based on its large projected separation (30,000 $\pm$ 140 AU) that exceeds typical turbulence scales in molecular clouds. To understand the lack of engulfment detections among our systems, we quantified the strength and duration of refractory enrichments in stellar photospheres using MESA stellar models. We found that observable signatures from 10 $M_{\oplus}$ engulfment events last for $\sim$90 Myr in 1 $M_{\odot}$ stars. Signatures are largest and longest lived for 1.1$-$1.2 $M_{\odot}$ stars, but are no longer observable $\sim$2 Gyr post-engulfment. This indicates that engulfment will rarely be detected in systems that are several Gyr old.

Low-mass neutron stars are directly associated with the nuclear saturation parameters because their central density is definitely low. We have already found a suitable combination of nuclear saturation parameters for expressing the neutron star mass and gravitational redshift, i.e., $\eta\equiv (K_0L^2)^{1/3}$ with the incompressibility for symmetric nuclear matter, $K_0$, and the density-dependent nuclear symmetry energy, $L$. In this study, we newly find another suitable combination given by $\eta_\tau\equiv (-K_\tau L^5)^{1/6}$ with the isospin dependence of incompressibility for asymmetric nuclear matter, $K_\tau$, and derive the empirical relations for the neutron star mass and gravitational redshift as a function of $\eta_\tau$ and the normalized central number density. With these empirical relations, one can evaluate the mass and gravitational redshift of the neutron star, whose central number density is less than threefold the saturation density, within $\sim 10\%$ accuracy, and the radius within a few \% accuracies. In addition, we discuss the neutron star mass and radius constraints from the terrestrial experiments, using the empirical relations, together with those from the astronomical observations. Furthermore, we find a tight correlation between $\eta_\tau$ and $\eta$. With this correlation, we derive the constraint on $K_\tau$ as $-348\le K_\tau\le -237$ MeV, assuming that $L=60\pm 20$ and $K_0=240\pm 20$ MeV.

A major conundrum of particle physics is what mass ordering (MO) follow neutrinos. The construction of next-generation neutrino detectors of unprecedented size, sensitivity and budget is underway and an answer is expected in the next decade through the combined study of reactor, atmospheric and accelerator neutrinos. In this Letter, the potential of an additional pathway is pointed out. Due to the MSW effect the flavor content of the neutrino flux from a Core-Collapse Supernovae (CCSNe) is highly dependent on the true neutrino MO. To exploit this feature, an analysis strategy robust to systematic uncertainties is identified for the first time and, by means of it, is shown that for a paradigmatic galactic CCSN a MO separation similar to 5~$\sigma$ can be achieved.

We revisit the framework of axion-like inflation in view of the possibility that the coupling of the inflaton to a non-Abelian topological charge density could lead to the generation of a rapidly thermalizing heat bath. Both dispersive (mass) and absorptive (friction) effects are included. For phenomenologically viable parameters, the system remains in a weak regime of warm inflation (thermal friction $\ll$ Hubble rate). For tensor perturbations we derive an interpolating formula that incorporates both vacuum and thermal production. The latter yields a model-independent frequency shape $\sim f_0^3$ in the LISA window, whose coefficient allows to measure the maximal shear viscosity of the thermal epoch. It is a challenge, however, to find models where the coefficient is large enough to be observable.

The spatial structure of the lowest 0$_1^+$, 0$_2^+$, 2$_1^+$ and 2$_2^+$ states of the $^{12}$C nucleus is studied within the 3$\alpha$ model with the Buck, Friedrich, and Wheatley $\alpha \alpha$ potential with Pauli forbidden states in the $S$ and $D$ waves. The Pauli forbidden states in the three-body system are treated by the exact orthogonalization method. The largest contributions to the ground and excited 2$_1^+$ bound states energies come from the partial waves $(\lambda, \ell)=(2,2)$ and $(\lambda, \ell)=(4,4)$. In contrast to the bound states, for the Hoyle resonance 0$_2^+$ and its analog state 2$_2^+$, dominant contributions come from the $(\lambda, \ell)=(0,0)$ and $(\lambda, \ell)=(2,2)$ configurations, respectively. The estimated probability density functions for the $^{12}$C(0$_1^+$) ground and 2$_1^+$ excited bound states show mostly a triangular structure, where the $\alpha$ particles move at a distance of about 2.5 fm from each other. However, the spatial structure of the Hoyle resonance and its analog state have a strongly different structure, like $^8$Be + $\alpha$. In the Hoyle state, the last $\alpha$ particle moves far from the doublet at the distance between $R=3.0$ fm and $R=5.0$ fm. In the Hoyle analog 2$_2^+$ state the two alpha particles move at a distance of about 15 fm, but the last $\alpha$ particle can move far from the doublet at the distance up to $R=30.0$ fm.

We analyze the phase transition in improved holographic QCD to obtain an estimate of the gravitational wave signal emitted in the confinement transition of a pure SU(3) Yang-Mills dark sector. We derive the effective action from holography and show that the energy budget and duration of the phase transition can be calculated with minor errors. These are used as input to obtain a prediction of the gravitational wave signal. To our knowledge, this is the first computation of the gravitational wave signal in a holographic model designated to match lattice data on the thermal properties of pure Yang-Mills.

In this work, we explore how modified gravity theories based on the non-metricity scalar, known as $f(Q)$ gravity, affects the propagation of gravitational waves from inspiraling of binary systems. We discuss forecast constraints on $f(Q)$ gravity by considering standard siren events in two contexts: i) simulated sources of gravitational waves as black hole - neutron star binary systems, emitting in the frequency band of the third-generation detector represented by the Einstein Telescope (ET); ii) three standard siren mock catalogs based on the merger of massive black hole binaries that are expected to be observed in the operating frequency band of the Laser Interferometer Space Antenna (LISA). We find that, within the ET sensitivity, it will be possible to test deviations from general relativity at $<3\%$ accuracy in the redshift range $0<z<4$, while the main free parameter of the theory is globally constrained at 1.6\% accuracy within the same range. In light of LISA's forecasts, in the best scenario, we find that the main free parameter of the theory will be constrained at 1.6\% accuracy up to high redshifts. Therefore, we conclude that future gravitational wave observations by ET and LISA will provide a unique way to test, with good accuracy, the nature of gravity up to very large cosmic distances.

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning-Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece-Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece-Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

We calculate the Big Bang Nucleosynthesis abundances for helium-4 and deuterium for a range of neutron lifetimes, $\tau_n = 840 - 1050$ s, using the state-of-the-art Python package \textsc{PRyMordial}. We show the results for two different nuclear reaction rates, calculated by NACRE II [1] and the PRIMAT [2] collaborations.

The DAMIC-M (DArk Matter In CCDs at Modane) experiment employs thick, fully depleted silicon charged-coupled devices (CCDs) to search for dark matter particles with a target exposure of 1 kg-year. A novel skipper readout implemented in the CCDs provides single electron resolution through multiple non-destructive measurements of the individual pixel charge, pushing the detection threshold to the eV-scale. DAMIC-M will advance by several orders of magnitude the exploration of the dark matter particle hypothesis, in particular of candidates pertaining to the so-called "hidden sector." A prototype, the Low Background Chamber (LBC), with 20g of low background Skipper CCDs, has been recently installed at Laboratoire Souterrain de Modane and is currently taking data. We will report the status of the DAMIC-M experiment and first results obtained with LBC commissioning data.

Nuclear data is critical for many modern applications from stockpile stewardship to cutting edge scientific research. Central to these pursuits is a robust pipeline for nuclear modeling as well as data assimilation and dissemination. We summarize a small portion of the ongoing nuclear data efforts at Los Alamos for medium mass to heavy nuclei. We begin with an overview of the NEXUS framework and show how one of its modules can be used for model parameter optimization using Bayesian techniques. The mathematical framework affords the combination of different measured data in determining model parameters and their associated correlations. It also has the advantage of being able to quantify outliers in data. We exemplify the power of this procedure by highlighting the recently evaluated 239-Pu cross section. We further showcase the success of our tools and pipeline by covering the insight gained from incorporating the latest nuclear modeling and data in astrophysical simulations as part of the Fission In R-process Elements (FIRE) collaboration.

We study numerically the scalar wave emission by a non-spherical oscillation of neutron stars in a scalar-tensor theory of gravity with kinetic screening, considering both the monopole and quadrupole mode emission. In agreement with previous results in the literature, we find that the monopole is always suppressed by the screening effect, regardless of the size of the screening radius, $r_{\rm sc}$. For the quadrupole mode, however, our analysis shows that the suppression only occurs for screening radius larger than the wavelength of scalar waves, $\lambda_{\rm wave}$, but not for $r_{\rm sc} < \lambda_{\rm wave}$. This demonstrates that to fully understand the nature of this theory, it is necessary to study other more complex systems, such as neutron star binaries, considering a wide range of $r_{\rm sc}$ values.