Papers for Wednesday, Mar 16 2022

We measure the colour evolution and quenching time-scales of $z=1.0-1.8$ galaxies across the green valley. We derive rest-frame $NUVrK$ colours and select blue-cloud, green-valley and red-sequence galaxies from the spectral energy distribution modelling of CANDELS GOODS-South and UDS multi-band photometry. Separately, we constrain the star-formation history (SFH) parameters (ages, $\tau$) of these galaxies by fitting their deep archival HST grism spectroscopy. We derive the galaxy colour-age relation and show that only rapidly evolving galaxies with characteristic delayed-$\tau$ SFH time-scales of $<0.5$ Gyr reach the red sequence at these redshifts, after a period of accelerated colour evolution across the green valley. These results indicate that the stellar mass build-up of these galaxies stays minimal after leaving the blue cloud and entering the green valley (i.e., it may represent $\lesssim 5\%$ of the galaxies' final, quiescent masses). Visual inspection of age-sensitive features in the stacked spectra also supports the view that these galaxies follow a quenching sequence along the blue-cloud $\rightarrow$ green-valley $\rightarrow$ red-sequence track. For this rapidly evolving population, we measure a green-valley crossing time-scale of $0.99^{+0.42}_{-0.25}$ Gyr and a crossing rate at the bottom of the green valley of $0.82^{+0.27}_{-0.25}$ mag/Gyr. Based on these time-scales, we estimate that the number density of massive ($M_\star>10^{10} M_\odot$) red-sequence galaxies doubles every Gyr at these redshifts, in remarkable agreement with the evolution of the quiescent galaxy stellar mass function. These results offer a new approach to measuring galaxy quenching over time and represent a pathfinder study for future JWST, Euclid, and Roman Space Telescope programs.

Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild-unstable but Ledoux-stable due to a composition gradient. Over many convective overturn timescales the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.

Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial AGN disk. We find numbers of order $10^{2-4}$, confined to the inner $r_{\rm cap} \sim 3000 r_s \sim 0.03\rm pc$. These stars can profoundly alter the chemistry of AGN disks, enriching them in helium and depleting them in hydrogen, both by order-unity amounts. We further consider mergers between these stars and other disk objects, suggesting that star-star mergers result in rapid mass loss from the remnant to restore an equilibrium mass, while star-compact object mergers may result in exotic outcomes and even host binary black hole mergers within themselves. Finally, we examine how these stars react as the disk dissipates towards the end of its life, and find that they may return mass to the disk fast enough to extend its lifetime by a factor of several and/or may drive powerful outflows from the disk. Post-AGN, these stars rapidly lose mass and form a population of stellar mass black holes around $10M_{\odot}$. Due to the complex and uncertain interactions between embedded stars and the disk, their plausible ubiquity, and their order unity impact on disk structure and evolution, they must be included in realistic disk models.

We investigate the impact of higher-order gravitational lens properties and properties of the background source on our approach to directly infer local lens properties from observables in multiple images of strong gravitationally lensed extended, static background sources developed in papers I to VI. As the degeneracy between local lens and source properties only allows to determine relative local lens properties between the multiple image positions, we cannot distinguish common scalings and distortions caused by lensing from intrinsic source characteristics. The consequences of this degeneracy for lens modelling and our approach and ways to break it are detailed here. We also set up quantitative measures around the critical curve to find clear limits on the validity of the approximation that source properties are negligible to infer local lens properties at critical points. The impact of the source on the local lens properties depends on the reduced shear at the image position and the amplitude and orientation of the source ellipticity, as we derive in this paper. Similarly, we investigate the role of third-order lens properties (flexion), in two galaxy-cluster simulations and in the Lenstool-reconstruction of the galaxy-cluster lens CL0024. In all three cases, we find that flexion is negligible in over 90% of all pixels of the lensing region for our current imprecision of local lens properties of about 10%. Decreasing the imprecision to 2%, higher-order terms start to play a role, especially in regions with shear components close to zero.

As long as astronomers have searched for exoplanets, the intrinsic variability of host stars has interfered with the ability to reliably detect and confirm exoplanets. One particular source of false positives is the presence of stellar magnetic or chromospheric activity that can mimic the radial-velocity reflex motion of a planet. Here we present the results of a photometric data analysis for the known planet hosting star, BD-06 1339, observed by the Transiting Exoplanet Survey Satellite (TESS) during Sector 6 at 2 minute cadence. We discuss evidence that suggests the observed 3.9 day periodic radial velocity signature may be caused by stellar activity rather than a planetary companion, since variability detected in the photometric data are consistent with the periodic signal. We conclude that the previously reported planetary signature is likely the result of a false positive signal resulting from stellar activity, and discuss the need for more data to confirm this conclusion.

We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics.

[abridged] Probing the faint end of the number counts at mm wavelengths is important to identify the origin of the extragalactic background light in this regime. Aided by strong gravitational lensing, ALMA observations towards massive galaxy clusters have opened a window to disentangle this origin, allowing to resolve sub-mJy dusty star-forming galaxies. We aim to derive number counts at 1.1 mm down to flux densities fainter than 0.1 mJy, based on ALMA observations towards five Hubble Frontier Fields (FF) galaxy clusters, following a statistical approach to correct for lensing effects. We created a source catalog that includes 29 ALMA 1.1 mm continuum detections down to a 4.5sigma significance. We derived source intrinsic flux densities using public lensing models. We folded the uncertainties in both magnifications and source redshifts into the number counts through Monte Carlo simulations. We derive cumulative number counts over two orders of magnitude down to 0.01 mJy after correction for lensing effects. Cosmic variance estimates are all exceeded by uncertainties in our median combined cumulative counts that come from both our Monte Carlo simulations and Poisson statistics. Our number counts are consistent to 1sigma with most of recent ALMA estimates and galaxy evolution models. However, below 0.1 mJy, they are lower by 0.4 dex compared to all deep ALMA studies but ASPECS-LP. Importantly, the flattening found for our cumulative counts extends further to 0.01 mJy. Our results bring further support in line of the flattening of the number counts reported previously by us and ASPECS-LP, which has been interpreted by a recent galaxy evolution model as a measurement of the "knee" of the infrared luminosity function at high redshift. Our estimates of the contribution to the EBL in the FFs suggest that we may be resolving most of the EBL at 1.1mm down to 0.01 mJy.

The diffuse light that spreads through groups and clusters of galaxies is made of free-floating stars not bound to any galaxy. This is known as the intracluster light (ICL) and holds important clues for understanding the evolution of these large structures. The study of this light has gained traction in the past 20 years thanks to technological and data processing advances that have permitted us to reach unprecedented observational depths. This progress has led to ground-breaking results in the field, such as pinpointing the origin of the ICL and its potential to map dark matter in clusters of galaxies. We now enter an era of deep and wide surveys that promise to uncover the faint Universe as never seen before, adding to our growing understanding of the properties of the ICL and, consequently, of the formation of the largest gravitationally bound structures in the Universe. The goal of this Review is to summarize the most recent results on ICL, synthesizing the current knowledge in the field and providing a global perspective that may benefit future ICL studies.

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but the best way to implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modeling choices for SMBH accretion and feedback in a suite of $\sim500$ cosmological zoom-in simulations across a wide range of halo mass (10^10-10^13 Msun). Within the suite, we vary the numerical schemes for BH accretion and feedback, the accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models that satisfy the observational constraints, and for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many successful models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, over-massive BHs tend to over-quench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass halos, but by feedback efficiency in low-mass halos. The amount of star formation suppression by SMBH feedback in low-mass halos is determined primarily by the time-integrated feedback energy. For massive galaxies, the "responsiveness" of a model (i.e. how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching.

Dusty quasars might be in a young stage of galaxy evolution with prominent quasar feedback. We study a population of luminous, extremely red quasars at redshifts $z\sim$~2--4 that has a suite of extreme spectral properties that may all be related to exceptionally powerful quasar-driven outflows. We present Keck/KCWI observations of the reddest known ERQ, at $z=$\,2.3184, with an extremely fast [\ion{O}{III}]~$\lambda$5007 outflow at $\sim$6000~km~s$^{-1}$. The Ly$\alpha$ halo spans $\sim$100~kpc. The halo is kinematically quiet, with velocity dispersion $\sim$300~km~s$^{-1}$ and no broadening above the dark matter circular velocity down to the spatial resolution $\sim$6~kpc from the quasar. We detect spatially-resolved \ion{He}{II}~$\lambda$1640 and \ion{C}{IV}~$\lambda$1549 emissions with kinematics similar to the Ly$\alpha$ halo and a narrow component in the [\ion{O}{III}]~$\lambda$5007. Quasar reddening acts as a coronagraph allowing views of the innermost halo. A narrow Ly$\alpha$ spike in the quasar spectrum is inner halo emission, confirming the broad \ion{C}{IV}~$\lambda$1549 in the unresolved quasar is blueshifted by $2242$~km~s$^{-1}$ relative to the halo frame. We propose the inner halo is dominated by past/moderate-speed outflow and the outer halo dominated by inflow. The high central concentration of the surface brightness and the circularly symmetric morphology of the inner halo are consistent with the ERQ being young. The \ion{He}{II}~$\lambda$1640/Ly$\alpha$ ratio of the inner halo and the asymmetry level of the overall halo are dissimilar to Type II quasars, implying central ionizing photons can escape along the line of sight. While quasar dominates the ionization, we find no evidence of mechanical quasar feedback on circumgalactic scales.

The number of potentially habitable planets continues to increase, but we lack the time and resources to characterize all of them. With $\sim$30 known potentially habitable planets and an ever-growing number of candidate and confirmed planets, a robust statistical framework for prioritizing characterization of these planets is desirable. Using the $\sim$2 Gy it took life on Earth to make a detectable impact on the atmosphere as a benchmark, we use a Bayesian statistical method to determine the probability that a given radius around a star has been continuously habitable for 2 Gy. We perform this analysis on 9 potentially habitable exoplanets with planetary radii $<$1.8 R$_\oplus$ and/or planetary masses $<$10 M$_\oplus$ around 9 low-mass host stars ($\sim$0.5-1.1 M$_\odot$) with measured stellar mass and metallicity, as well as Venus, Earth, and Mars. Ages for the host stars are generated by the analysis. The technique is also used to provide age estimates for 2768 low-mass stars (0.5-1.3 M$_\odot$) in the TESS Continuous Viewing Zones.

Ambipolar diffusion is a process occurring in partially ionised astrophysical systems that imparts a complicated mathematical and physical nature to Ohm's law. The numerical codes that solve the magnetohydrodynamic (MHD) equations have to be able to deal with the singularities that are naturally created in the system by the ambipolar diffusion term. The global aim is to calculate a set of theoretical self-similar solutions to the nonlinear diffusion equation with cylindrical symmetry that can be used as tests for MHD codes which include the ambipolar diffusion term. First, following the general methods developed in the applied mathematics literature, we obtained the theoretical solutions as eigenfunctions of a nonlinear ordinary differential equation. Phase-plane techniques were used to integrate through the singularities at the locations of the nulls, which correspond to infinitely sharp current sheets. In the second half of the paper, we consider the use of these solutions as tests for MHD codes. To that end, we used the Bifrost code, thereby testing the capabilities of these solutions as tests as well as (inversely) the accuracy of Bifrost's recently developed ambipolar diffusion module. The obtained solutions are shown to constitute a demanding, but nonetheless viable, test for MHD codes that incorporate ambipolar diffusion. The Bifrost code is able to reproduce the theoretical solutions with sufficient accuracy up to very advanced diffusive times. Using the code, we also explored the asymptotic properties of our theoretical solutions in time when initially perturbed with either small or finite perturbations. The functions obtained in this paper are relevant as physical solutions and also as tests for general MHD codes. They provide a more stringent and general test than the simple Zeldovich-Kompaneets-Barenblatt-Pattle solution.

The small-scale moving intensity enhancements remotely observed in the extreme ultraviolet images of the solar active regions, which we refer to as active region moving campfires (ARMCs), are related to local plasma temperature and/or density enhancements. Their dynamics is driven by the physical processes in the entire coronal plasma. Our previous study of ARMCs indicates that they have characteristic velocities at around the background sound speed. The main goal of our work is to carry out a simultaneous analysis of EUV images from two observational missions, SDO/AIA and Hi-C 2.1. The aims of the performed cross-validating analysis of both SDO/AIA and Hi-C 2.1 data were to reveal how the observed moving features are distributed over the studied active region, AR12712, test the existence of different groups of ARMCs with distinct physical characteristics. We use the statistical model of intensity centroid convergence and tracking that was developed in our previous paper. Furthermore, a Gaussian mixture model fit of the observed complex of moving ARMCs is elaborated to reveal the existence of distinct ARMC groups and to study the physical characteristics of these different groups. In data from the 171\AA, 193\AA\ and 211\AA\ channels of SDO/AIA, we identified several groups of ARMCs with respect to both blob intensity and velocity profiles. The existence of such groups is confirmed by the cross-validation of the 172\AA\ data sets from Hi-C 2.1. The ARMCs studied in this paper have characteristic velocities in the range of the typical sound speeds in coronal loops. Hence, these moving objects differ from the well-known rapid Alfv\'enic velocity jets from magnetic reconnection sites. This is also proven by the fact that ARMCs propagate along the active region magnetic structure (strands). The nature of the discovered statistical grouping of the ARMC events is not known.

The total masses of galaxy clusters characterize many aspects of astrophysics and the underlying cosmology. It is crucial to obtain reliable and accurate mass estimates for numerous galaxy clusters over a wide range of redshifts and mass scales. We present a transfer-learning approach to estimate cluster masses using the ugriz-band images in the SDSS Data Release 12. The target masses are derived from X-ray or SZ measurements that are only available for a small subset of the clusters. We designed a semi-supervised deep learning model consisting of two convolutional neural networks. In the first network, a feature extractor is trained to classify the SDSS photometric bands. The second network takes the previously trained features as inputs to estimate their total masses. The training and testing processes in this work depend purely on real observational data. Our algorithm reaches a mean absolute error (MAE) of 0.232 dex on average and 0.214 dex for the best fold. The performance is comparable to that given by redMaPPer, 0.192 dex. We have further applied a joint integrated gradient and class activation mapping method to interpret such a two-step neural network. The performance of our algorithm is likely to improve as the size of training dataset increases. This proof-of-concept experiment demonstrates the potential of deep learning in maximizing the scientific return of the current and future large cluster surveys.

We develop a new method for analytical inversion of binned exoplanet transit spectra and for retrieval of planet parameters. The method has a geometrical interpretation and treats each observed spectrum as a single vector $\vec r$ in the multidimensional spectral space of observed bin values. We decompose the observed $\vec{r}$ into a wavelength-independent component $\vec{r}_\parallel$ corresponding to the spectral mean across all observed bins, and a transverse component $\vec{r}_\perp$ which is wavelength-dependent and contains the relevant information about the atmospheric chemistry. The method allows us to extract, without any prior assumptions or additional information, the relative mass (or volume) mixing ratios of the absorbers in the atmosphere, the scale height to stellar radius ratio, $H/R_S$, and the atmospheric temperature. The method is illustrated and validated with several examples of increasing complexity.

We explore atmospheric escape from close-in exoplanets with the highest mass loss rates. First, we locate the transition from stellar X-ray and UV-driven escape to rapid Roche lobe overflow, which occurs once the 10-100 nbar pressure level in the atmosphere reaches the Roche lobe. Planets enter this regime when the ratio of the substellar radius to the polar radius along the visible surface pressure level, that aligns with a surface of constant Roche potential, is X/Z~$\gtrsim$~1.2 for Jovian planets (Mp~$\gtrsim$~100 M$_{\Earth}$) and X/Z~$\gtrsim$~1.02 for sub-Jovian planets ($M_p \approx$~10--100 M$_{\Earth}$). Around a sun-like star, this regime applies to orbital periods of less than two days for planets with radii of about 3--14 R$_{\Earth}$. Our results agree with the properties of known transiting planets and can explain parts of the sub-Jovian desert in the population of known exoplanets. Second, we present detailed numerical simulations of atmospheric escape from a planet like Uranus or Neptune orbiting close to a sun-like star that support the results above and point to interesting qualitative differences between hot Jupiters and sub-Jovian planets. We find that hot Neptunes with solar metallicity hydrogen and helium envelopes have relatively more extended upper atmospheres than typical hot Jupiters, with a lower ionization fraction and higher abundances of escaping molecules. This is consistent with existing ultraviolet transit observations of warm Neptunes and it might provide a way to use future observations and models to distinguish solar metallicity atmospheres from higher metallicity atmospheres.

We revisit the eccentric neutron star (NS)-white dwarf (WD) binary model for the periodic activity of fast radio burst (FRB) sources, by including the effects of gravitational wave (GW) radiation. In this model, the WD fills its Roche lobe at the periastron and mass transfer occurs from the WD to the NS. The accreted materials can be fragmented and arrive at the NS episodically, resulting in multiple bursts through curvature radiation. Consequently, the WD may be kicked away owing to the conservation of angular momentum. To initiate the next mass transfer, the WD has to refill its Roche lobe through GW radiation. In this scenario, whether the periodic activity can show up relies on three timescales, i.e., the orbital period $P_{\rm orb}$, the timescale $T_{\rm GW}$ for the Roche lobe to be refilled, and the time span $T_{\rm frag}$ for all the episodic events corresponding to each mass transfer process. Only when the two conditions $T_{\rm GW} \gtrsim P_{\rm orb}$ and $T_{\rm frag} < P_{\rm orb}$ are both satisfied, the periodic activity will manifest itself and the period should be equal to $P_{\rm orb}$. In this spirit, the periodic activity is more likely to show up for relatively long periods ($P_{\rm orb} \gtrsim$ several days). Thus, it is reasonable that FRBs 180916 and 121102, the only two sources having been claimed to manifest periodic activity, both correspond to relatively long periods.

V476 Cyg (Nova Cyg 1920) is a bright, fast nova reaching a photographic magnitude of 2.0. Using the Zwicky Transient Facility (ZTF) public database, I found that this nova is currently a dwarf nova with a cycle length of ~24 d. Compared to other classical novae currently in dwarf nova-type states, outbursts of V476 Cyg are rapidly rising and short with durations of a few days. Based on the AAVSO observations, this nova was probably already in the dwarf nova-type phase in 2016, 96 years after the nova eruption. I found a possible orbital period of 0.1018002(6) d using the ZTF data, which would place the object in the period gap. This supposed short orbital period appears to explain the features and faint absolute magnitudes of the observed dwarf nova outbursts. If this period is confirmed, V476 Cyg is a classical nova with the shortest orbital period with distinct dwarf nova outbursts and in which a nova eruption was recorded in the modern era. I also compared with the outburst properties with V446 Her (Nova Her 1960), which currently shows SS Cyg-type outbursts. The transition to the dwarf nova-phase in V476 Cyg occurred much earlier (~100 yr) than what has been supposed (~1000 yr) for classical novae below the period gap. V476 Cyg would not only provide an ideal laboratory of the behavior of an irradiated accretion disk in which tidal instability is expected to work, but also an ideal laboratory of the effect of a massive white dwarf on dwarf nova outbursts.

The nearby, massive, runaway star Zeta Ophiuchi has a large bow shock detected in optical and infrared, and, uniquely among runaway O stars, diffuse X-ray emission is detected from the shocked stellar wind. Previous investigation of NGC 7635 found that numerical simulations over-predict the X-ray emission by ~10 times. Here we make the first detailed computational investigation of the bow shock of Zeta Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare the detected X-ray emission with simulated emission maps. We re-analysed archival {\it Chandra} observations of the thermal diffuse X-ray emission from the shocked wind region of the bow shock, finding total unabsorbed X-ray flux (0.3-2 keV band) corresponding to a diffuse luminosity of $L_\mathrm{X}=2.33~(0.79-3.45)\times10^{29}$ergs$^{-1}$. 3D MHD simulations were used to model the interaction of the star's wind with a uniform ISM using a range of stellar and ISM parameters motivated by observational constraints. Synthetic infrared, Ha, soft X-ray, emission measure, and radio 6\,GHz emission maps were generated from three simulations, for comparison with relevant observations. Simulations where the space velocity of Zeta Ophiuchi has a significant radial velocity produce infrared emission maps with opening angle of the bow shock in better agreement with observations than for the case where motion is fully in the plane of the sky. The simulation with the highest pressure has the closest match, with flux level within a factor of 2 of the observational lower limit, and emission weighted temperature of $\log_{10}(T_\mathrm{A}/\mathrm{K})=6.4$, although the morphology of the diffuse emission appears somewhat different. Observed X-ray emission is a filled bubble brightest near the star whereas simulations predict brightening towards the contact discontinuity as density increases.

We analyze two binary systems containing giant stars, V723 Mon ("the Unicorn") and 2M04123153+6738486 ("the Giraffe"). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modeling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\,giant} = 3,800\,\rm K$ and $4,000\,\rm K$, $R_{\rm giant}= 22.5\,R_{\odot}$ and $25\,R_{\odot}$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\,M_{\odot}$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\,2} = 5,800\,\rm K$ and $5,150\,\rm K$, $R_2= 8.3\,R_{\odot}$ and $9\,R_{\odot}$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8\,M_{\odot}$ and $\approx 1.8\,M_{\odot}$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1-2\,M_{\odot}$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely-observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs.

We applied the image-based approach with a convolutional neural network model to the sample of low-redshifts galaxies with $-24^{m}<M_{r}<-19.4^{m}$ from the SDSS DR9. We divided it into two subsamples, SDSS DR9 galaxy dataset and Galaxy Zoo 2 (GZ2) dataset, considering them as the inference and training datasets, respectively. As a result, we created the morphological catalog of 315782 galaxies at 0.02<z<0.1, where morphological five classes and 34 detailed features (bar, rings, number of spiral arms, mergers, etc.) were first defined for 216148 galaxies (inference dataset) by the image-based CNN classifier. For the rest of galaxies the initial morphological classification was re-assigned as in the GZ2 project. Our method shows the promising performance of morphological classification attaining more 93 % of accuracy for five classes morphology prediction except the cigar-shaped (75 %) and completely rounded (83 %) galaxies. Main results are presented in the catalog of 19468 completely rounded, 27321 rounded in-between, 3235 cigar-shaped, 4099 edge-on, 18615 spiral, and 72738 general low-redshift galaxies of the studied SDSS sample. As for the classification of galaxies by their detailed structural morphological features, our CNN model gives the accuracy in the range of 92-99 % depending on features, a number of galaxies with the given feature in the inference dataset, and the galaxy image quality. We demonstrate that implication of the CNN model with adversarial validation and adversarial image data augmentation improves classification of smaller and fainter SDSS galaxies with $m_{r}$ <17.7.

ixpeobssim is a simulation and analysis framework, based on the Python programming language and the associated scientific ecosystem, specifically developed for the Imaging X-ray Polarimetry Explorer (IXPE). Given a source model and the response functions of the telescopes, it is designed to produce realistic simulated observations, in the form of event lists in FITS format, containing a strict super-set of the information provided by standard IXPE level-2 files. The core ixpeobssim simulation capabilities are complemented by a full suite of post-processing applications, allowing for the implementation of complex, polarization-aware analysis pipelines, and facilitating the inter-operation with the standard visualization and analysis tools traditionally in use by the X-ray community. We emphasize that, although a significant part of the framework is specific to IXPE, the modular nature of the underlying implementation makes it potentially straightforward to adapt it to different missions with polarization capabilities.

Gamma Cygni is a young supernova remnant located in the Cygnus region. MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes detected TeV emission (MAGIC J2019+408) to the north-west of this remnant, about 5 arcmin from its border. We want to identify the radio sources within the region encompassing gamma Cygni and MAGIC J2019+408 to shed light on their nature and investigate if these radio sources could be potential contributors to gamma-ray emission. We carried out a detailed study of the data we obtained with a survey of the Cygnus region at 325 and 610 MHz with the Giant Metrewave Radio Telescope (GMRT). We detected several radio sources in the region where the radio and the TeV emission overlap, as well as several areas of enhanced radio emission. In particular, two of these areas of diffuse enhanced emission may correspond to the supernova remnant interacting with a high density region, which seems to be the best candidate for the MAGIC source. Another two radio sources, which may or may not contribute to the gamma rays, are also spatially coincident with the emission peak of the MAGIC TeV source. One of them displays a rather peculiar extended morphology whose nature is completely unknown. We have identified the radio sources overlapping gamma Cygni and MAGIC J2019+408 and have shown that their potential gamma-ray contribution is likely not dominant. In addition, some of the studied sources show peculiar physical characteristics that deserve deeper multi-wavelength observations.

Andromeda (M31) is the nearest giant spiral galaxy to our Milky Way (MW) and the most massive member of our Local Group. We obtain a magnitude-limited sample of M31 disc PNe with chemical abundance measurements through the direct detection of the [OIII] 4363 \r{A} line. This leads to $205$ and $200$ PNe with oxygen and argon abundance measurements respectively. We find that high- and low-extinction M31 disc PNe have statistically distinct argon and oxygen abundance distributions. In the radial range $2-30$ kpc, the older low-extinction disc PNe are metal-poorer on average with a slightly positive radial oxygen abundance gradient ($0.006 \pm 0.003$ dex/kpc) and slightly negative radial argon abundance gradient ($-0.005 \pm 0.003$ dex/kpc), while the younger high-extinction disc PNe are metal-richer on average with steeper radial abundance gradients for both oxygen ($-0.013 \pm 0.006$ dex/kpc) and argon ($-0.018 \pm 0.006$ dex/kpc), similar to the gradients measured for M31 HII regions. These abundance gradients are consistent with a major merger in the M31 disc, with the majority of the low-extinction PNe being the older pre-merger disc stars in the thicker disc, and the majority of the high-extinction PNe being younger stars in the thin disc, formed during and after the merger event. The chemical abundance of the M31 thicker disc has been radially homogenized because of the major merger. Accounting for differences in disc-scale lengths, we find the radial oxygen abundance gradient of the M31 thicker disc is amongst the most positive values observed till date in spiral galaxies, flatter than that of the MW thick disc which has a negative radial gradient. On the other hand, the thin discs of the MW and M31 have remarkably similar negative oxygen abundance gradient values.

The Doppler effect is commonly used to infer the velocity difference between stars based on the relative shifts in the rest-frame wavelengths of their spectral features. In wide binaries, the difference in gravitational redshift from the surfaces of the constituent stars with distinct compactness dominates at separations above 0.01 pc. I suggest that this effect became apparent for wide pairs in the Gaia eDR3 catalogue but incorrectly interpreted as a possible modification of Newtonian gravity in the internal kinematics of very wide binaries.

We present the first neutron star merger simulations performed with the newly developed Numerical Relativity code SPHINCS_BSSN. This code evolves the spacetime on a mesh using the BSSN formulation, but matter is evolved via Lagrangian particles according to a high-accuracy version of general-relativistic Smooth Particle Hydrodynamics (SPH). Our code contains a number of new methodological elements compared to other Numerical Relativity codes. The main focus here is on the new elements that were introduced to model neutron star mergers. These include a) a refinement (fixed in time) of the spacetime-mesh, b) corresponding changes in the particle--mesh mapping algorithm and c) a novel way to construct SPH initial data for binary systems via the recently developed "Artificial Pressure Method." This latter method makes use of the spectral initial data produced by the library LORENE, and is implemented in a new code called SPHINCS_ID. While our main focus is on introducing these new methodological elements and documenting the current status of SPHINCS_BSSN, we also show as a first application a set of neutron star merger simulations employing "soft" ($\Gamma=2.00$) and "stiff" ($\Gamma=2.75$) polytropic equations of state.

Thanks to tremendous experimental efforts, galactic cosmic-ray fluxes are being measured up to the unprecedented per cent precision level. The logarithmic slope of these fluxes is a crucial quantity that promises us information on the diffusion properties and the primary or secondary nature of the different species. However, these measured slopes are sometimes interpreted in the pure diffusive regime, guiding to misleading conclusions. In this paper, we have studied the propagation of galactic cosmic rays by computing the fluxes of species between H and Fe using the USINE code and considering all the relevant physical processes and an updated set of cross-section data. We show that the slope of the well-studied secondary-to-primary B/C ratio is distinctly different from the diffusion coefficient slope, by an offset of about 0.2 in the rigidity range in which the AMS-02 data reach their best precision (several tens of GV). Furthermore, we have demonstrated that none of the species from H to Fe follows the expectations of the pure-diffusive regime. We argue that these differences arise from propagation processes such as fragmentation, convection, and reacceleration, which cannot be neglected. On this basis, we also provide predictions for the spectral slope of elemental fluxes not yet analysed by the AMS collaboration.

Construction of the precise shape of an asteroid is critical for spacecraft operations as the gravitational potential is determined by spatial mass distribution. The typical approach to shape determination requires a prolonged mapping phase of the mission over which extensive measurements are collected and transmitted for Earth-based processing. This paper presents a set of approaches to explore an unknown asteroid with onboard calculations, and to land on its surface area selected in an optimal fashion. The main motivation is to avoid the extended period of mapping or preliminary ground observations that are commonly required in spacecraft missions around asteroids. First, range measurements from the spacecraft to the surface are used to incrementally correct an initial shape estimate according to the Bayesian framework. Then, an optimal guidance scheme is proposed to control the vantage point of the range sensor to construct a complete 3D model of the asteroid shape. This shape model is then used in a nonlinear controller to track a desired trajectory about the asteroid. Finally, a multi resolution approach is presented to construct a higher fidelity shape representation in a specified location while avoiding the inherent burdens of a uniformly high resolution mesh. This approach enables for an accurate shape determination around a potential landing site. We demonstrate this approach using several radar shape models of asteroids and provide a full dynamical simulation about asteroid 4769 Castalia.

The noise of gravitational-wave (GW) interferometers limits their sensitivity and impacts the data quality, hindering the detection of GW signals from astrophysical sources. For transient searches, the most problematic are transient noise artifacts, known as glitches, that happen at a rate around $ 1 \text{ min}^{-1}$, and can mimic GW signals. Because of this, there is a need for better modeling and inclusion of glitches in large-scale studies, such as stress testing the pipelines. In this proof-of concept work we employ Generative Adversarial Networks (GAN), a state-of-the-art Deep Learning algorithm inspired by Game Theory, to learn the underlying distribution of blip glitches and to generate artificial populations. We reconstruct the glitch in the time-domain, providing a smooth input that the GAN can learn. With this methodology, we can create distributions of $\sim 10^{3}$ glitches from Hanford and Livingston detectors in less than one second. Furthermore, we employ several metrics to measure the performance of our methodology and the quality of its generations. This investigation will be extended in the future to different glitch classes with the final goal of creating an open-source interface for mock data generation.

The recent detection of a Fast Radio Burst (FRB) from a Galactic magnetar secured the fact that neutron stars (NSs) with super-strong magnetic fields are capable of producing these extremely bright coherent radio bursts. One of the leading mechanisms to explain the origin of such coherent radio emission is the curvature radiation process within the dipolar magnetic field structure. It has, however, already been demonstrated that magnetars likely have a more complex magnetic field topology. Here we critically investigate curvature radio emission in the presence of inclined dipolar and quadrupolar ("quadrudipolar") magnetic fields and show that such field structures differ in their angular characteristics from a purely dipolar case. We analytically show that the shape of open field lines can be modified significantly depending on both the ratio of quadrupole to dipole field strength and their inclination angle at the NS surface. This creates multiple points along each magnetic field line that coincides with the observer's line of sight and may explain the complex spectral and temporal structure of the observed FRBs. We also find that in quadrudipole, the radio beam can take a wider angular range and the beam width can be wider than in pure dipole. This may explain why the pulse width of the transient radio pulsation from magnetars is as large as that of ordinary radio pulsars.

Evidence for the presence of short-lived radioactive isotopes when the Solar System formed is preserved in meteorites, providing insights into the conditions at the birth of our Sun. A low-mass core-collapse supernova had been postulated as a candidate for the origin of $^{10}$Be, reinforcing the idea that a supernova triggered the formation of the Solar System. We present a detailed study of the production of $^{10}$Be by the $\nu$ process in supernovae, which is very sensitive to the reaction rate of the major destruction channel, $^{10}$Be(p,$\alpha$)$^{7}$Li. With data from recent nuclear experiments that show the presence of a resonant state in $^{11}$B at $\approx$~193 keV, we derive new values for the $^{10}$Be(p,$\alpha$)$^{7}$Li reaction rate which are significantly higher than previous estimates. We show that, with the new $^{10}$Be(p,$\alpha$)$^{7}$Li reaction rate, a low mass CCSN is unlikely to produce enough $^{10}$Be to explain the observed $^{10}$Be/ $^{9}$Be ratio in meteorites, even for a wide range of neutrino spectra considered in our models. These findings point towards spallation reactions induced by solar energetic particles in the early Solar System as the origin of $^{10}$Be.

LIGO has detected several binary black hole (BBH) merger events that may have originated in the accretion disks of Active Galactic Nuclei (AGN). These events require individual black hole masses that fall within the pair instability supernova mass gap, and therefore these black holes may have been grown from hierarchical mergers. AGN disks are a prime environment for hierarchical mergers, and thus a potential location for the progenitors of BBH gravitational wave events. Understanding how a BBH embedded in an AGN disk interacts with the surrounding environment is thus crucial for determining if this interaction can lead to its merger. However, there are few high fidelity simulations of this process, and almost all are two-dimensional. We present the results from 3D, high-resolution, local shearing-box simulations of an embedded BBH interacting with an AGN disk. In these first simulations of their kind, we focus on determining the mass accretion rate and the orbital evolution rate at different BBH separations. We find that circular, equal-mass BBHs with separations greater than 10% of their Hill radius contract while accreting at a super-Eddington rate. At smaller separations, however, our 3D simulations find that BBHs expand their orbits. This result suggests that it may be difficult for an AGN disk to push a BBH to merger, but we discuss several mechanisms, including MHD turbulence and radiative and mechanical feedback, that could alleviate this difficulty.

We aim to analyze the chemical evolution of the Small Magellanic Cloud adding 12 additional clusters to our existing sample having accurate and homogeneously derived metallicities. We are particularly interested in seeing if there is any correlation between age and metallicity for the different structural components to which the clusters belong. Spectroscopic metallicities of red giant stars are derived from the measurement of the equivalent width of the near-IR calcium triplet lines. Cluster membership analysis was carried out using criteria that include radial velocities, metallicities, proper motions and distance from the cluster center. The mean cluster radial velocity and metallicity were determined with a typical error of 2.1 km/s and 0.03 dex, respectively. We added this information to that available in the literature for other clusters studied with the same method, compiling a final sample of 48 clusters with metallicities homogeneously determined. Clusters of the final sample are distributed in an area of ~ 70 deg^2 and cover an age range from 0.4 Gyr to 10.5 Gyr. The metallicity distribution of our new cluster sample shows a lower probability of being bimodal than suggested in previous studies. The separate chemical analysis of clusters in the six components (Main Body, Counter-Bridge, West Halo, Wing/Bridge, Northern Bridge and Southern Bridge) shows that only clusters belonging to the Northern Bridge appear to trace a V-Shape, showing a clear inversion of the metallicity gradient in the outer regions. There is a suggestion of a metallicity gradient in the West Halo, similar to that previously found for field stars. It presents, however, a very large uncertainty. Also, clusters belonging to the West Halo, Wing/Bridge and Southern Bridge exhibit a well-defined age-metallicity relation with relatively little scatter in abundance at fixed age compared to other regions.

The origin of Uranus and Neptune is still unknown. In particular, it has been challenging for planet formation models to form the planets in their current radial distances within the expected lifetime of the solar nebula. In this paper, we simulate the in-situ formation of Uranus and Neptune via pebble accretion and show that both planets can form within ~ 3 Myr at their current locations, and have final compositions that are consistent with the heavy-element to H-He ratios predicted by structure models. We find that Uranus and Neptune could have been formed at their current locations. In several cases a few earth masses of heavy elements are missing, suggesting that Uranus and/or Neptune may have accreted $\sim$1 -- 3 M Earth's mass of heavy elements after their formation via planetesimal accretion and/or giant impacts.

The quasi-molecular mechanism of recombination, recently suggested by Kereselidze et al., is a non-standard process where an electron and two neighboring protons in the early universe directly form an ionized hydrogen molecule in a highly excited state, which then descends to lower levels or dissociates. It has been suggested that the increased binding energy due to the participation of a second proton may lead to an earlier cosmic recombination that alleviates the Hubble tension. Revisiting the quasi-molecular channel of recombination in more details, we find that the original work significantly overestimated the probability of finding a pair of adjacent protons in the relevant epoch ($z\sim $ a few thousand). Our new estimation suggests that the quasi-molecular mechanism of recombination cannot be the primary cause of the Hubble tension.

Interacting supernovae (SNe) IIn and Ibn show narrow emission lines and have always been a mysterious and unsolved genre in SNe physics. We present a comprehensive analysis of the temporal and spectroscopic behaviour of a group of interacting SNe~IIn and Ibn. We choose SNe~2012ab, 2020cui, 2020rc and 2019uo as representative members of these SN sub-types to probe the nature of explosion. Our study reveals that SNe~IIn are heterogeneous, bright depicting multi-staged temporal evolution while SNe~Ibn are moreover homogeneous, comparatively fainter than SNe~IIn and short lived, but limited in sample to firmly constrain the homogeneity. The spectroscopic features display a great diversity in H$\alpha$ and He profiles for both SNe~IIn and Ibn. The representative SN~Ibn also show flash ionisation signatures of CIII and NIII. Modelling of H$\alpha$ reveals that SNe~IIn have in general an asymmetric CSM which interacts with SN ejecta resulting in diversity in H$\alpha$ profiles.

The Cherenkov Telescope Array (CTA) is the next-generation ground-based very-high-energy gamma-ray observatory. The Large-Sized Telescope (LST) of CTA is designed to detect gamma rays between 20 GeV and a few TeV with a 23-meter diameter mirror. We have developed the focal plane camera of the first LST, which has 1855 photomultiplier tubes (PMTs) and the readout system which samples a PMT waveform at GHz with switched capacitor arrays, Domino Ring Sampler ver4 (DRS4). To measure the precise pulse charge and arrival time of Cherenkov signals, we developed a method to calibrate the output voltage of DRS4 and the sampling time interval, as well as an analysis method to correct the spike noise of DRS4. Since the first LST was inaugurated in 2018, we have performed the commissioning tests and calibrated the camera. We characterised the camera in terms of the charge pedestal under various conditions of the night sky background, the charge resolution of each pixel, the charge uniformity of the whole camera, and the time resolutions with a test pulse and calibration laser.

Heating by the central star is one of the key factors determining the physical structure of protoplanetary disks. Due to the large optical thickness in the radial direction, disk midplane regions are heated by the infrared radiation from the disk surface (atmosphere), which in turn is directly heated by the star. It was previously shown that interception of the stellar radiation by inhomogeneities on the disk surface can cause perturbations that propagate towards the star. In this work, we investigate the occurrence of such waves within a detailed 1+1D numerical model of the protoplanetary disk. We confirm the previous findings that in the disk, that is optically thick to its own radiation, the surface perturbations indeed occur and propagate towards the star. However, contrary to some analytical predictions, the thermal waves in sufficiently massive disks affect only the upper layers without significant fluctuations of temperature in the midplane. Our results indicate the need to study this instability within more consistent hydrodynamic models.

The extended source effect on microlensing magnification is non-negligible and must be taken into account for in an analysis of microlensing. However, the evaluation of the extended source magnification is numerically expensive because it includes the two-dimensional integral over source profile. Various studies have developed methods to reduce this integral down to the one-dimensional-integral or integral-free form, which adopt some approximations or depend on the exact form of the source profile, e.g. disk, linear/quadratic limb-darkening profile. In this paper, we develop a new method to evaluate the extended source magnification based on fast Fourier transformation (FFT), which does not adopt any approximations and is applicable to any source profiles. Our implementation of the FFT based method enables the fast evaluation of the extended source magnification as fast as $\sim1$ msec (CPU time on a laptop) and guarantees an accuracy better than 0.3%. The FFT based method can be used for the template fitting to a huge data set of light curves from the existing and upcoming surveys.

Context. While synthesis of organic molecules in molecular clouds or protoplanetary disks is complex, observations of interstellar grains, analyses of carbonaceous chondrites, and UV photochemistry experiments are rapidly developing and providing constraints on and clues to the complex organic molecule synthesis in space. It motivates us to construct a theoretical synthesis model. Aims. We develop a new code to simulate global reaction sequences of organic molecules to apply it for sugar synthesis by intermittent UV irradiation on the surface of icy particles in a protoplanetary disk. Here we show the first results of our new simulation. Methods. We apply a Monte Carlo method to select reaction sequences from all possible reactions, using the graph-theoretic matrix model for chemical reactions and modeling reactions on the icy particles during UV irradiation. Results. We here obtain the results consistent with the organic molecules in carbonaceous chondrites and obtained by the experiments, however, through a different pathway from the conventional formose reactions previously suggested. During UV irradiation, loosely-bonded O-rich large molecules are continuously created and destroyed. After UV irradiation is turned off, the ribose abundance rapidly increases, through the decomposition of the large molecules with break-ups of O-O bonds and replacements of C-OH by C-H to reach O/C = 1 for sugars. The sugar abundance is regulated mostly by the total atomic ratio H/O of starting materials, but not by their specific molecule forms. Deoxyribose is simultaneously synthesized, and most of the molecules end up with complex C-rich molecules.

NGC 4258 is one of the most important anchors for calibrating the Cepheid period--luminosity relations (PLRs) owing to its accurate distance measured from water maser motions. We expand on previous efforts and carry out a new Cepheid search in this system using the Hubble Space Telescope (HST). We discover and measure a sample of 669 Cepheids in four new and two archival NGC 4258 fields, doubling the number of known Cepheids in this galaxy and obtaining an absolute calibration of their optical PLRs. We determine a Wesenheit PLR of $-2.574(\pm0.034) -3.294(\pm0.042) \log P$, consistent with an independent Large Magellanic Cloud (LMC) calibration at the level of $0.032\pm0.044$~mag in its zeropoint, after accounting for a metallicity dependence of $-0.20\pm0.05$~mag\,dex$^{-1}$ (Riess et al. 2006). Our determination of the PLR slope also agrees with the LMC-based value within their uncertainties. We attempt to characterize the metallicity effect of Cepheid PLRs using only the NGC 4258 sample, but a relatively narrow span of abundances limits our sensitivity and yields a Wesenheit zero-point dependence of $-0.07 \pm 0.21$ mag\,dex$^{-1}$. The Cepheid measurements presented in this study have been used as part of the data to derive the Hubble constant in a companion paper by the SH0ES team.

We propose a generalization of our previous KDE (kernel density estimation) method for estimating luminosity functions (LFs). This new upgrade further extend the application scope of our KDE method, making it a very flexible approach which is suitable to deal with most of bivariate LF calculation problems. From the mathematical point of view, usually the LF calculation can be abstracted as a density estimation problem in the bounded domain of $\{Z_1<z<Z_2,~ L>f_{\mathrm{lim}}(z) \}$. We use the transformation-reflection KDE method ($\hat{\phi}$) to solve the problem, and introduce an approximate method ($\hat{\phi}_{\mathrm{1}}$) based on one-dimensional KDE to deal with the small sample size case. In practical applications, the different versions of LF estimators can be flexibly chosen according to the Kolmogorov-Smirnov test criterion. Based on 200 simulated samples, we find that for both cases of dividing or not dividing redshift bins, especially for the latter, our method performs significantly better than the traditional binning method $\hat{\phi}_{\mathrm{bin}}$. Moreover, with the increase of sample size $n$, our LF estimator converges to the true LF remarkably faster than $\hat{\phi}_{\mathrm{bin}}$. To implement our method, we have developed a public, open-source Python Toolkit, called \texttt{kdeLF}. With the support of \texttt{kdeLF}, our KDE method is expected to be a competitive alternative to existing nonparametric estimators, due to its high accuracy and excellent stability. \texttt{kdeLF} is available at \url{this http URL} with extensive documentation available at \url{this http URL}.

The similarity of Venus and Earth in bulk properties make Venus an appealing target for future colonization. Several proposals have been put forward for colonizing and even terraforming Venus despite the extreme conditions on the planet's surface. Such a terraforming project would face large challenges centered around removing Venus's massive carbon dioxide atmosphere and replacing it with a habitable environment. I review past proposals and propose a new method for terraforming Venus by building an artificial surface in the much more hospitable upper atmosphere where the temperature and pressure are both Earth-like. Such a surface could be built with locally produced materials and would float above the bulk of the atmosphere using nitrogen as a lifting gas. This would allow the engineering of a breathable atmosphere above the surface and would remove the need to import or export extreme amounts of mass, except for comparatively modest quantities of water. The engineering, logistical, and energy requirements of this method are surveyed. I find that such a terraforming project could be completed in a minimum of 200 years in a best-case scenario, comparable to other proposals, with significantly lower resource costs.

The purpose of this work is to characterize the diffuse Galactic polarized synchrotron. We present EE, BB, and EB power spectra estimated cross-correlating Planck and WMAP polarization frequency maps at 23 and 30 GHz, for a set of six sky regions covering from 30% to 94% of the sky. The EE and BB angular power spectra show a steep decay of the spectral amplitude as a function of multipole, approximated by a power law with power indices around -2.9 for both components. The B/E ratio is about 0.22. The EB cross-component is compatible with zero at 1$\sigma$, with upper constraint on the EB/EE ratio of 1.2% at the 2$\sigma$ level. The recovered SED, in the frequency range 23-30 GHz, shows E and B power-law spectral indices compatible between themselves with a value of about -3.

A supernova remnant (SNR) candidate SRGe~J0023+3625 = G116.6-26.1 was recently discovered in the \textit{SRG}/eROSITA all-sky X-ray survey. This large ($\sim 4$ deg in diameter) SNR candidate lacks prominent counterparts in other bands. Here we report detection of radio emission from G116.6-26.1 in the LOFAR Two-metre Sky Survey (LoTTS-DR2). Radio images show a shell-like structure coincident with the X-ray boundary of the SNR. The measured surface brightness of radio emission from this SNR is very low. Extrapolation of the observed surface brightness to 1~GHz places G116.6-26.1 well below other objects in the $\Sigma-D$ diagram. We argue that the detected radio flux might be consistent with the minimal level expected in the van der Laan adiabatic compression model, provided that the volume emissivity of the halo gas in the LOFAR band is $\sim 10^{-42}\,{\rm Wm^{-3}Hz^{-1} sr^{-1}}$. If true, this SNR can be considered as a prototypical example of an evolved SNR in the Milky Way halo. In the X-ray and radio bands, such SNRs can be used as probes of thermal and non-thermal components constituting the Milky Way halo.

We investigate the formation of the large scale structures in the present accelerated era in $f(R)$ gravity background. This is done by considering the linear growth of matter perturbations at low redshift $z<1$. The effect of $f(R)$ alters the behaviour of the matter density perturbations from the matter dominated universe to the late-time accelerated universe which is encoded in the Newtonian gravitational constant as $G\rightarrow G_{eff}$. The modified gravitational constant ($G_{eff}$) depends on the form of $f(R)$. The late-time accelerated expansion affects the formation of large scale structures by slowing down the growth of matter density. On the other hand, $f(R)$ increases the growth rate of the matter density perturbations. We have found that the source term in $f(R)$ background, $G_{eff}\Omega_m$ overcomes the accelerated expansion and the effect of accelerated expansion suppresses the formation of the large scale structures in the asymptotic future.

At fixed stellar mass $M_*$, the effective radius $R_{\rm e}$ of massive satellite early-type galaxies (ETGs) in galaxy clusters is, on average, larger at lower redshift. We study theoretically this size evolution using the state-of-the-art cosmological simulation IllustrisTNG100: we sampled $75$ simulated satellite ETGs at redshift $z=0$ with $M_* \ge 10^{10.4} M_{\odot}$ belonging to the two most massive ($\approx 10^{14.6} M_{\odot} $) haloes of the simulation. We traced back in time the two clusters' main progenitors and we selected their satellite ETGs at $z>0$ with the same criterion adopted at $z=0$. The $R_{\rm e}-M_*$ relation of the simulated cluster satellite ETGs, which is robustly measured out to $z=0.85$, evolves similarly to the observed relation over the redshift range $0\lesssim z \lesssim 0.85$. In the simulation the main drivers of this evolution are the acquisition of new galaxies ("newcomers") by the clusters and the transformation of member galaxies located at large clustercentric distance ("suburbanites") at $z=0.85$, which end up being massive satellite ETGs at $z=0$. Though several physical processes contribute to change the population of satellite ETGs in the considered redshift interval, the shape of the stellar mass function of the simulated cluster ETGs is not significantly different at $z=0.85$ and at $z=0$, consistent with observations.

This talk has the goal of introducing two upcoming exciting avenues of discovery in precision Cosmology: spectral distortions (SDs) and gravitational waves (GWs). The former signals offer a clear window into the state of the primordial plasma at times prior to recombination, and thus sheds light on small-scale primordial perturbations, dark matter decay and black hole formation to just mention a few scenarios. The latter signals, which are already offering new insights into the landscape of black hole and neutron star mergers, will reveal intricate dynamics of the early universe including primordial black holes, inflationary potentials, and even reheating dynamics. An elegant link is drawn between these two future observations since primordial GW backgrounds will source SDs, a coupling which offers unique insight to over six decades of GW frequencies. More importantly the SD visibility window bridges the gap between astrophysical high- and cosmological low-frequency measurements. This means SDs will not only complement other GW observations, but will be the sole probe of physical processes at certain scales.

WFI J161953.3+031909 was considered to be an eclipsing novalike object in the period gap. Using the Zwicky Transient Facility (ZTF) public database, I found that this object is the first eclipsing Z Cam star in the period gap, and is also most likely an ER UMa star with supercycles of 60-80 d. The longest outburst (most likely a superoutburst) comprised 35-46% of the supercycle in the extreme case. If superhumps are confirmed, this becomes the second object showing both ER UMa and Z Cam states after NY Ser. These objects have anomalously high mass-transfer rates despite that they are in the period gap. I refined the orbital period to be 0.099419808(8) d. We can expect to learn from WFI J161953.3+031909 using eclipses what is actually happening in the ER UMa-type and Z Cam-type disk. I also provide a revised classification of an eclipsing IW And star for BMAM-V383 = IPHAS J200822.55+300341.6 by detecting an IW And-type standstill.

The Stokes profiles of Fe I lines in the photosphere of the Sun are calculated within the Unno-Beckers-Landi-Dagl`Innocenti theory. Estimates of the magnetic strengthening of the lines were obtained. The changes in the Stokes profiles depending on the excitation potential, wavelength, equivalent width, Lande factor, micro-macroturbulent velocities, radial velocity, damping constant, atmospheric model, magnetic field strength and direction are considered. The graphically presented variations of the Stokes profiles make it possible to determine the initial values of the input parameters for solving the problems of magnetic field vector reconstruction by the inversion method. The presented dependencies of the magnetic strengthening on the line parameters will help to correctly select magnetically sensitive lines for the investigation of sunspots, flux tubes, plages, and other magnetic features.

The relation between properties of galaxies and dark matter halos they reside in can be valuable to understand structure formation and evolution, particularly the baryonic-to-halo mass ratio (BHMR) and its dynamic evolution. We first review some unique properties of self-gravitating collisionless dark matter flow (SG-CFD), followed by their implications in derivation of BHMR. To maximize system entropy, the long-range interaction requires a broad size of halos to be formed. These halos facilitate an inverse mass and energy cascade from small to large mass scales that involves a constant rate of energy transfer $\epsilon_u$=$-4.6\times10^{-7}m^2/s^3$. Considering galaxies and halos with a total baryonic mass $m_b$, halo mass $m_h$, halo virial size $r_h$, and flat rotation velocity $v_f$, the baryonic-to-halo mass relation can be analytically derived by combining the baryonic Tully-Fisher relation and constant rate of energy cascade $\epsilon_u$. We found a maximum BHMR ratio ~0.076 for halos with a critical mass $m_{hc}$~$10^{12}M_{sun}$ at z=0. That ratio is much lower for both smaller and larger halos such that two regimes can be identified: i) for incompressible small halos with mass $m_h$<$m_{hc}$, we have $\epsilon_u$$\propto$$v_f/r_h$, $v_f$$\propto$$r_h$, and $m_b$$\propto$$m_h^{4/3}$; ii) for large halos with mass $m_h$>$m_{hc}$, we have $\epsilon_u$$\propto$$v_f^3/r_h$, $v_f$$\propto$$r_h^{1/3}$, and $m_b$$\propto$$m_h^{4/9}$. Combined with proposed double-$\lambda$ halo mass function, the average BHMR ratio for all halos (~0.024 at z=0) can be analytically derived, along with its redshift evolution. The fraction of baryons in galaxy is ~7.6% at z=0 and increases with time $\propto$$t^{1/3}$. Most (~92.4%) of the baryons are not in galaxies. The SPARC (Spitzer Photometry & Accurate Rotation Curves) data with ~175 late-type galaxies were used for the derivation and comparison.

We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune the parameters of a popular star formation and feedback model within Enzo to match simulations. The parameters that are adjusted include the star formation efficiency, coupling of thermal energy from stellar feedback, and volume into which the energy is deposited. This number translates to a factor of more than three improvements over manual calibration. Despite using fewer simulations, we obtain a better agreement to the observed baryon makeup of a Milky-Way (MW) sized halo. Switching to a different strategy, we improve the consistency of the recommended parameters from the machine. Given the success of the calibration, we then apply the technique to reconcile metal transport between grid-based and particle-based simulation codes using an isolated galaxy. It is an improvement over manual exploration while hinting at a less known relation between the diffusion coefficient and the metal mass in the halo region. The exploration and calibration of the parameters of the sub-grid models with a machine learning approach is concluded to be versatile and directly applicable to different problems.

Numerical simulations are at the center of predicting the space debris environment of the upcoming decades. In light of debris generating events, such as continued anti-satellite weapon tests and planned mega-constellations, accurate predictions of the space debris environment are critical to ensure the long-term sustainability of critical satellite orbits. Given the computational complexity of accurate long-term trajectory propagation for a large number of particles, numerical models usually rely on Monte-Carlo approaches for stochastic conjunction assessment. On the other hand, deterministic methods bear the promise of higher accuracy and can serve to validate stochastic approaches. However, they pose a substantial challenge to computational feasibility. In this work, we present the architecture and proof of concept results for a numerical simulation capable of modeling the long term debris evolution over decades with a deterministic conjunction tracking model. For the simulation, we developed an efficient propagator in modern C++ accounting for Earth's gravitational anomalies, solar radiation pressure, and atmospheric drag. We utilized AutoPAS, a sophisticated particle container, which automatically selects the most efficient data structures and algorithms. We present results from a simulation of 16 024 particles in low-Earth orbit over 20 years. Overall, conjunctions are tracked for predicted collisions and close encounters to allow a detailed study of both. We analyze the runtime and computational cost of the simulation in detail. In summary, the obtained results show that modern computational tools finally enable deterministic conjunction tracking and can serve to validate prior results and build higher-fidelity numerical simulations of the long-term debris environment.

The physical origin of fast radio bursts (FRBs) remains unclear. Finding multiwavelength counterparts of FRBs can be a breakthrough in understanding their nature. In this work, we perform a systematic search for astronomical transients (ATs) whose positions are consistent with FRBs. We find an unclassified optical transient AT2020hur ($\alpha=01^{\mathrm{h}} 58^{\mathrm{m}} 00.750^{\mathrm{s}} \pm 1$ arcsec, $\delta=65^{\circ} 43^{\prime} 00.30^{\prime \prime} \pm 1$ arcsec) that is spatially coincident with the repeating FRB 180916B ($\alpha=01^{\mathrm{h}} 58^{\mathrm{m}} 00.7502^{\mathrm{s}} \pm 2.3$ mas, $\delta=65^{\circ} 43^{\prime} 00.3152^{\prime \prime} \pm 2.3$ mas). The chance possibility for the AT2020hur-FRB 180916B association is about 0.04%, which corresponds to a significance of $3.5\sigma$. We develop a giant flare afterglow model to fit AT2020hur. Although the giant flare afterglow model can interpret the observations of AT2020hur, the derived kinetic energy of such a GF is at least three orders of magnitude larger than that of the typical GF, and there is a lot of fine tuning and coincidences required for this model. Another possible explanation is that AT2020hur might consist of two or more optical flares originating from the FRB source, e.g. fast optical bursts produced by the inverse Compton scattering of FRB emission. Besides, AT2020hur is located in one of the activity windows of FRB 180916B, which is an independent support for the association. This coincidence may be due to the reason that the optical counterpart is subject to the same periodic modulation as FRB 180916B, as implied by the prompt FRB counterparts. Future simultaneous observations of FRBs and their optical counterparts may help reveal their physical origin.

Solar coronal mass ejections are well known to expand as they propagate through the heliosphere. Despite this, their cross-sections are usually modeled as static plasma columns within the magnetohydrodynamics (MHD) framework. We test the validity of this approach using in-situ plasma data from 151 magnetic clouds (MCs) observed by the WIND spacecraft and 45 observed by the Helios spacecrafts. We find that the most probable cross-section expansion speeds for the WIND events are only $\approx 0.06$ times the Alfv\'en speed inside the MCs while the most probable cross-section expansion speeds for the Helios events is $\approx 0.03$. MC cross-sections can thus be considered to be nearly static over an Alfv\'en crossing timescale. Using estimates of electrical conductivity arising from Coulomb collisions, we find that the Lundquist number inside MCs is high ($\approx 10^{13}$), suggesting that the MHD description is well justified. The Joule heating rates using our conductivity estimates are several orders of magnitude lower than the requirement for plasma heating inside MCs near the Earth. While the (low) heating rates we compute are consistent with the MHD description, the discrepancy with the heating requirement points to possible departures from MHD and the need for a better understanding of plasma heating in MCs.

In recent years we have seen an important increase in the number of protonated molecules detected in cold dense clouds. Here we report the detection in TMC-1 of HCCNCH+, the protonated form of HCCNC, which is a metastable isomer of HC3N. This is the first protonated form of a metastable isomer detected in a cold dense cloud. The detection was based on observations carried out with the Yebes 40m and IRAM 30m telescopes, which revealed four harmonically related lines. We derive a rotational constant B = 4664.431891 +/- 0.000692 MHz and a centrifugal distortion constant D = 519.14 +/- 4.14 Hz. From a high-level ab initio screening of potential carriers we confidently assign the series of lines to the ion HCCNCH+. We derive a column density of (3.0 +/- 0.5)e10 cm-2 for HCCNCH+, which results in a HCCNCH+/HCCNC abundance ratio of 0.010 +/- 0.002. This value is well reproduced by a state-of-the-art chemical model, which however is subject to important uncertainties regarding the chemistry of HCCNCH+. The observational and theoretical status of protonated molecules in cold dense clouds indicate that there exists a global trend in which protonated-to-neutral abundance ratios MH+/M increase with increasing proton affinity of the neutral M, although if one restricts to species M with high proton affinities (>700 kJ/mol), MH+/M ratios fall in the range 0.001-0.1, with no apparent correlation with proton affinity. We suggest various protonated molecules that are good candidates for detection in cold dense clouds in the near future.

Using numerical simulations, we investigate the gravitational evolution of filamentary molecular cloud structures and their condensation into dense protostellar cores. One possible process is the so called 'edge effect', the pile-up of matter at the end of the filament due to self-gravity. This effect is predicted by theory but only rarely observed. To get a better understanding of the underlying processes we used a simple analytic approach to describe the collapse and the corresponding collapse time. We identify a model of two distinct phases: The first phase is free fall dominated, due to the self-gravity of the filament. In the second phase, after the turning point, the collapse is balanced by the ram pressure, produced by the inside material of the filament, which leads to a constant collapse velocity. This approach reproduces the established collapse time of uniform density filaments and agrees well with our hydrodynamic simulations. In addition, we investigate the influence of different radial density profiles on the collapse. We find that the deviations compared to the uniform filament are less than 10%. Therefore, the analytic collapse model of the uniform density filament is an excellent general approach.

KAGRA, the kilometer-scale underground gravitational-wave detector, is located at Kamioka, Japan. In April 2020, an astrophysics observation was performed at the KAGRA detector in combination with the GEO 600 detector; this observation operation is called O3GK. The optical configuration in O3GK is based on a power recycled Fabry-P\'{e}rot Michelson interferometer; all the mirrors were set at room temperature. The duty factor of the operation was approximately 53%, and the strain sensitivity was $3\times10^{-22}~/\sqrt{\rm{Hz}}$ at 250 Hz. In addition, the binary-neutron-star (BNS) inspiral range was approximately 0.6 Mpc. The contributions of various noise sources to the sensitivity of O3GK were investigated to understand how the observation range could be improved; this study is called a "noise budget". According to our noise budget, the measured sensitivity could be approximated by adding up the effect of each noise. The sensitivity was dominated by noise from the sensors used for local controls of the vibration isolation systems, acoustic noise, shot noise, and laser frequency noise. Further, other noise sources that did not limit the sensitivity were investigated. This paper provides a detailed account of the KAGRA detector in O3GK including interferometer configuration, status, and noise budget. In addition, strategies for future sensitivity improvements such as hardware upgrades, are discussed.

A currently unsolved question in supernova research is the origin of stripped-envelope supernovae (SESNe). Such SNe lack spectral signatures of hydrogen (Type Ib), or hydrogen and helium (Type Ic), indicating that the outer stellar layers have been stripped during their evolution. The mechanism for this is not well understood, and to disentangle the different scenarios determination of nucleosynthesis yields from observed spectra can be attempted. However, the interpretation of observations depends on the adopted spectral models. A previously missing ingredient in these is the inclusion of molecular effects, which can be significant. We aim to investigate how the molecular chemistry in stripped-envelope supernovae affects physical conditions and optical spectra and produces ro-vibrational emission in the mid-infrared (MIR). We also aim to assess the diagnostic potential of observations of such MIR emission with JWST. We couple a chemical kinetic network including carbon, oxygen, silicon, and sulfur-bearing molecules into the NLTE spectral synthesis code SUMO. We let four species - CO, SiO, SiS and SO - participate in the NLTE cooling of the gas to achieve self-consistency between the molecule formation and the temperature. We apply the new framework to model the spectrum of a Type Ic supernova in the 100-600d time range. Molecules are predicted to form in SESN ejecta in significant quantities (typical mass $10^{-3}$ $M_\odot$) throughout the 100-600d interval. The impact on the temperature and optical emission depends on the density of the oxygen zones and varies with epoch. For example, the [O I] 6300, 6364 feature can be quenched by molecules from 200 to 450d depending on density. The MIR predictions show strong emission in the fundamental bands of CO, SiO, and SiS, and in the CO and SiO overtones.

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

We present systematic and uniform analysis of NuSTAR data with 10-78 keV S/N > 50, of a sample of 60 SWIFT BAT selected AGNs, 10 of which are radio-loud. We measure their high energy cutoff Ecut or coronal temperature Te using three different spectral models to fit their NuSTAR spectra, and show a threshold in NuSTAR spectral S/N is essential for such measurements. High energy spectral breaks are detected in the majority of the sample, and for the rest strong constraints to Ecut or Te are obtained. Strikingly, we find extraordinarily large Ecut lower limits ( > 400 keV, up to > 800 keV) in 10 radio-quiet sources, whereas none in the radio-loud sample. Consequently and surprisingly, we find significantly larger mean Ecut/Te of radio-quiet sources compared with radio-loud ones. The reliability of these measurements are carefully inspected and verified with simulations. We find a strong positive correlation between Ecut and photon index {\Gamma}, which can not be attributed to the parameter degeneracy. The strong dependence of Ecut on {\Gamma}, which could fully account for the discrepancy of Ecut distribution between radio-loud and radio-quiet sources, indicates the X-ray coronae in AGNs with steeper hard X-ray spectra have on average higher temperature and thus smaller opacity. However, no prominent correlation is found between Ecut and {\lambda}edd. In the l-{\Theta} diagram, we find a considerable fraction of sources lie beyond the boundaries of forbidden regions due to runaway pair production, posing (stronger) challenges to various (flat) coronal geometries.

The relation between X-ray and UV/optical variability in AGNs has been explored in many individual sources, however a large sample study is yet absent. Through matching the XMM-Newton serendipitous X-ray and UV source catalogs with SDSS quasars, we build a sample of 802 epoch-pairs of 525 quasars showing clear variability in logFx-logFuv space. After correcting for the effect of photometric noise, we find 35.6\pm2.1% of the epoch-pairs show asynchronous variability between X-ray and UV (brightening in one band but dimming in the other). This indicates only in 28.8\pm4.2% of the epoch-pairs the X-ray and UV variability are intrinsically coordinated. The variability synchronicity exhibits no dependence on physical parameters of quasars or the time lag of the epoch-pairs, except for stronger variability tends to have stronger synchronicity. Switches between synchronous and asynchronous variability are also seen in individual sources. The poor coordination clearly contradicts both the X-ray reprocessing model and the accretion rate variation model for AGN variability. The ratios of the observed X-ray variability amplitude to that in UV span a broad range and peak at ~ 2. The dominant fraction of the ratios appear too small to be attributed to X-ray reprocessing, and too large for accretion rate variation. The inhomogeneous disk model which incorporates both X-ray and UV/optical variability in AGNs is favored by the observed stochastic relation between X-ray and UV variations.

In the disks of spiral galaxies, diffuse soft X-ray emission is known to be strongly correlated with star-forming regions. However, this emission is not simply from a thermal-equilibrium plasma and its origin remains greatly unclear. In this work, we present an X-ray spectroscopic analysis of the emission from a star-forming enhanced region off the nucleus of the spiral galaxy M51, named as the northern hot spot (NHS) hereafter. We analyze the high spectral resolution data from XMM-Newton/RGS observations and unambiguously detect an abnormally high G-ratio ($3.2^{+6.9}_{-1.5}$) of the O VII K$\alpha$ triplet emission from the NHS, which is spatially confirmed by oxygen emission line maps from the RGS data. A physical model consisting of a thermal plasma and its charge exchange (CX) with cold gas gives an excellent description of the high G-ratio and the entire RGS spectra. Complementary to the RGS data, the Chandra/ACIS-S data are also analyzed because its high spatial resolution is beneficial for decomposing the diffuse and discrete source contributions. Its imaging spectra of the diffuse emission either from the NHS or other parts of the M51 disk can be well interpreted by the same physical model. Though the surface brightness of the diffuse emission from the NHS is about twice that from other parts of the disk, the hot plasma has the similar temperature of ~0.34 keV and approximately solar metallicity. For both parts of the disk, the CX contributes ~50% to the diffuse soft X-ray emission in the 7-30 angstrom band, and has an interface area of about five times the geometric area of the M51 disk. Our results suggest that the hot plasma can interact with the cold interstellar medium efficiently, and have strong implications for our understanding of soft X-ray emission as a tracer of the stellar feedback in active star-forming galaxies.

We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-$C_\ell$ method and offer a view complementary to that of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, such as baryons and intrinsic alignments (IA), making this analysis an important cross-check. In the context of $\Lambda$CDM, and using the same fiducial model as in the DES Y3 real space analysis, we find ${S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to ${S_8 = 0.784\pm 0.026 }$ when including shear ratios. This constraint is within expected statistical fluctuations from the real space analysis, and in agreement with DES~Y3 analyses of non-Gaussian statistics, but favors a slightly higher value of $S_8$, which reduces the tension with the Planck cosmic microwave background 2018 results from $2.3\sigma$ in the real space analysis to $1.5\sigma$ in this work. We explore less conservative IA models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}{h{\rm Mpc}^{-1}}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, which is found to be about 20\% lower than predicted by Planck 2018, as reflected by the $1.5\sigma$ lower $S_8$ value.

Supernova (SN) feedback plays a crucial role in simulations of galaxy formation. Because blastwaves from individual SNe occur on scales that remain unresolved in modern cosmological simulations, SN feedback must be implemented as a subgrid model. Differences in the manner in which SN energy is coupled to the local interstellar medium and in which excessive radiative losses are prevented have resulted in a zoo of models used by different groups. However, the importance of the selection of resolution elements around young stellar particles for SN feedback has largely been overlooked. In this work, we examine various selection methods using the smoothed particle hydrodynamics code SWIFT. We run a suite of isolated disk galaxy simulations of a Milky Way-mass galaxy and small cosmological volumes, all with the thermal stochastic SN feedback model used in the EAGLE simulations. We complement the original mass-weighted neighbour selection with a novel algorithm guaranteeing that the SN energy distribution is as close to isotropic as possible. Additionally, we consider algorithms where the energy is injected into the closest, least dense, or most dense neighbour. We show that different neighbour-selection strategies cause significant variations in star formation rates, gas densities, wind mass loading factors, and morphology. The isotropic method results in more efficient feedback and better numerical convergence than the conventional mass-weighted selection. We conclude that the manner in which the feedback energy is distributed among the resolution elements surrounding a feedback event is as important as changing the amount of energy by factors of a few.

The goal of this paper is to investigate the optimal parameters space for using the reloaded Rossiter-McLaughlin technique to detect differential rotation (DR) and centre-to-limb convective variations. We simulated a star-planet system with and without convective effects to map the optimal regions of the parameter space for retrieving the injected differential rotation. Our simulations explored all possible ranges of projected obliquity (spin-orbit angle), stellar inclination, and impact parameter, as well as differences in instrumental configuration, stellar magnitude, and exposure time. We find that DR is more easily retrieved at low-impact parameters, corresponding to system configurations in which the transiting planet crosses the largest number of stellar latitudes. The main hot-spots for detection (i.e. areas in which DR detectability is high) are $120^{\rm{o}}<|\lambda|<180^{\rm{o}}$ for $i_*<90^{\rm{o}}$ and $|\lambda|<60^{\rm{o}}$ for $i_*>90^{\rm{o}}$ on average, and they tend to shrink as the impact parameter increases. Additionally, in contrast to the crucial impact of brightness, we identify that exposure time has a negligible impact on the difficulty of detecting DR as the increase in signal-to-noise ratio (S/N) at longer exposure times is counteracted by the degraded sampling rate. We determine that an ESPRESSO-like setup of instrumental configuration and sensitivity might retrieve DR up to $V = 12$, compared to $V = 10$ for HARPS. We reach no clear conclusion about limb-dependent convective effects and the possible confusion with DR; preliminary results suggest, however, that under certain circumstances, while it seems that one effect could be mistaken for the other, the accuracy of the fit (in particular of $\alpha$) does not hold up under additional scrutiny.

We study the bispectrum of large-scale structure in the EFTofLSS including corrections up to two-loop. We derive an analytic result for the double-hard limit of the two-loop correction, and show that the UV-sensitivity can be absorbed by the same four EFT operators that renormalize the one-loop bispectrum. For the single-hard region, we employ a simplified treatment, introducing one extra EFT parameter. We compare our results to N-body simulations, and show that going from one- to two-loop extends the wavenumber range with percent-level agreement from $k \simeq 0.08$ to $0.15~h/\mathrm{Mpc}$.

We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest density planet orbiting an M dwarf (M0V). It orbits a solar-metallicity M dwarf discovered using TESS photometry and confirmed with precise radial velocities (RV) from HPF and NEID. With a planetary radius of $12.0^{+0.4}_{-0.5}$ $R_{\oplus}$ and mass of $85.3^{+8.8}_{-8.7}$ $M_{\oplus}$, not only does this object add to the small sample of gas giants ($\sim 10$) around M dwarfs, but also, its low density ($\rho =$ $0.27^{+0.05}_{-0.04}$ $\textrm{g~cm}^{-3}$) provides an opportunity to test theories of planet formation. We present two hypotheses to explain its low density; first, we posit that the low metallicity of its stellar host ($\sim$ 0.3 dex lower than the median metallicity of M dwarfs hosting gas giants) could have played a role in the delayed formation of a solid core massive enough to initiate runaway accretion. Second, we use our eccentricity estimate of $0.14 \pm 0.06$ to estimate the significance of tidal heating for inflating the radius of TOI-3757 b. The low density and large scale height of TOI-3757 b makes it an excellent target for transmission spectroscopy studies of atmospheric escape and composition (TSM $\sim$ 190). We use HPF to perform transmission spectroscopy of TOI-3757 b using the helium line. Doing this, we place an upper limit of 6.9 \% (with 90\% confidence) on the maximum depth of the absorption from the metastable transition of He at $\sim$ 10830 \AA.

Observations of the redshifted 21-cm signal emitted by neutral hydrogen represent a promising probe of large-scale structure in the universe. However, cosmological 21-cm signal is challenging to observe due to astrophysical foregrounds which are several orders of magnitude brighter. Traditional linear foreground removal methods can optimally remove foregrounds for a known telescope response but are sensitive to telescope systematic errors such as antenna gain and delay errors, leaving foreground contamination in the recovered signal. Non-linear methods such as principal component analysis, on the other hand, have been used successfully for foreground removal, but they lead to signal loss that is difficult to characterize and requires careful analysis. In this paper, we present a systematics-robust foreground removal technique which combines both linear and non-linear methods. We first obtain signal and foreground estimates using a linear filter. Under the assumption that the signal estimate is contaminated by foreground residuals induced by parameterizable systematic effects, we infer the systematics-induced contamination by cross-correlating the initial signal and foreground estimates. Correcting for the inferred error, we are able to subtract foreground contamination from the linearly filtered signal up to the first order in the amplitude of the telescope systematics. In simulations of an interferometric 21-cm survey, our algorithm removes foreground leakage induced by complex gain errors by one to two orders of magnitude in the power spectrum. Our technique thus eases the requirements on telescope characterization for modern and next-generation 21-cm cosmology experiments.

The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold clues about the cause for the accelerated expansion of the universe and possibly the nature of dark matter. During the next decade, LSST will be able to confirm or dispute if tensions seen today in cosmological data are due to new physics. New and unexpected phenomena could confirm or disrupt our current understanding of the universe. Findings from LSST will guide the path forward post-LSST. The Rubin Observatory will still be a uniquely powerful facility even then, capable of revealing further insights into the physics of the dark universe. These could be obtained via innovative observing strategies, e.g., targeting new probes at shorter timescales than with LSST, or via modest instrumental changes, e.g., new filters, or through an entirely new instrument for the focal plane. This White Paper highlights some of the opportunities in each scenario from Rubin observations after LSST.

In this paper, we propose a self-consistent test for a Hubble constant estimate using galaxy cluster and type Ia supernovae (SNe Ia) observations. The approach consists, in a first step, of obtaining the observational value of the galaxy cluster scaling-relation $Y_{SZE}D_{A}^{2}/C_{XZS}Y_X = C $ by combining the X-Ray and SZ observations of galaxy clusters at low redshifts ($z < 0.1$) from the first {\it Planck mission} all-sky data set ($0.044 \leq z \leq 0.444$), along with SNe Ia observations and making use of the cosmic distance duality relation validity. Then, by considering a flat $\Lambda$CDM model for $D_A$, the constant $C$ from the first step and the Planck prior on $\Omega_M$ parameter, we obtain $H_0$ by using the galaxy cluster data with $z>0.1$. As a result, we obtain a low $H_0$ estimate, $H_0=63.583^{+6.496}_{-5.685}$ km/s/Mpc, in full agreement with the latest results from {\it Planck mission}.

There are indications that the magnetic field evolution in galaxies might be massively shaped by tidal interactions and mergers between galaxies. The details of the connection between the evolution of magnetic fields and that of their host galaxies is still a field of research. We use a combined approach of magnetohydrodynamics for the baryons and an N-body scheme for the dark matter to investigate magnetic field amplification and evolution in interacting galaxies. We find that, for two colliding equal-mass galaxies and for varying initial relative spatial orientations, magnetic fields are amplified during interactions, yet cannot be sustained. Furthermore, we find clues for an active mean-field dynamo.

Characterizing accurately the polarized dust emission from our Galaxy will be decisive for the quest for the Cosmic Microwave Background (CMB) primordial $B$-modes. The incomplete modeling of its potentially complex spectral properties could lead to biases in the CMB polarization analyses and to a spurious detection of the tensor-to-scalar ratio $r$. Variations of the dust properties along and between lines of sight lead to unavoidable distortions of the spectral energy distribution (SED) that can not be easily anticipated by standard component separation methods. This issue can be tackled using a moment expansion of the dust SED, an innovative parametrization method imposing minimal assumptions on the sky complexity. In the recent work [Vacher \emph{et al.} (2022)]\cite{Vacher_2022}, we apply this formalism to the $B$-mode cross-angular power spectra computed from simulated \lb{} polarization data at frequencies between 100 and 402\,GHz, containing CMB, dust and instrumental noise. Thanks to the moment expansion, we can measure an unbiased value of the tensor-to-scalar ratio with a dispersion compatible with the target values aimed by the instrument.

Optical navigation on a CubeSat must rely on the best extraction of the directions of some beacons from on-board images. We present an experiment on OPS-SAT, a CubeSat of the European Space Agency (ESA), that will characterize an onboard algorithm to this aim, named Angle-based Correlation (AbC). OPS-SAT is a 3-unit CubeSat with an Attitude Determination and Control System (ADCS) and an imager that have proved their reliability with typical performance at CubeSat scale. We selected a few star-fields that all present enough visible stars within a 10 ? field of view. When our experiment is run, OPS-SAT is pointed to the most convenient star-field at that time. There, the star-field is imaged and subwindows are extracted from the image around the expected location of each star, based on the attitude-quaternion reported by the ADCS. The AbC reconstructs the absolute direction of the central body, in principle unknown, which is the pointed known star in the experiment. The method intensively uses the quaternion algebra. The beacon location is first consolidated in the field of view with the AbC. Then, the field of view is finely positioned against the sky, again with the AbC. A covariance is associated with the found beacon direction. Our experiment with OPS-SAT manages the pointing and the imager, and processes the taken images. Then, it downloads the on-board computed absolute directions and their covariances, to be compared with the actual directions. After a campaign of intensive use of the experiment, the statistical performance of the algorithm will be established and compared to the on-board computed covariances. As a bonus, an assessment of OPS-SAT's inertial pointing stability will be available. The AbC can theoretically get rid of the Attitude Control Error (ACE) of the platform and of the Attitude Knowledge Error (AKE) estimated by the ADCS, and potentially converge to (...)

Next-generation tests of fundamental physics and cosmology using large scale structure require measurements over large volumes of the Universe, including high redshifts inaccessible to present-day surveys. Line intensity mapping, an emerging technique that detects the integrated emission of atomic and molecular lines without resolving sources, can efficiently map cosmic structure over a wide range of redshifts. Observations at millimeter wavelengths detect far-IR emission lines such as CO/[CII], and take advantage of observational and analysis techniques developed by CMB experiments. These measurements can provide constraints with unprecedented precision on the physics of inflation, neutrino masses, light relativistic species, dark energy and modified gravity, and dark matter, among many other science goals. In this white paper we forecast the sensitivity requirements for future ground-based mm-wave intensity mapping experiments to enable transformational cosmological constraints. We outline a staged experimental program to steadily improve sensitivity, and describe the necessary investments in developing detector technology and analysis techniques.

The AMS-02 and HELIX experiments should provide soon $^{10}$Be/$^9$Be cosmic-ray data of unprecedented precision. We propose an analytical formula to quickly and accurately determine $L$ from these data. Our formula is validated against the full calculation performed with the propagation code USINE. We compare the constraints on $L$ set by Be/B and $^{10}$Be/$^9$Be, relying on updated sets of production cross sections. The best-fit $L$ from AMS-02 Be/B data is shifted from 5 kpc to 3.8 kpc when using the updated cross sections. We obtain consistent results from the Be/B analysis with USINE, $L=3.8^{+2.8}_{-1.6}$ kpc (data and cross-section uncertainties), and from the analysis of $^{10}$Be/$^9$Be data with the simplified formula, $L=4.7\pm0.6$ (data uncertainties) $\pm2$ (cross-section uncertainties) kpc. The analytical formula indicates that improvements on $L$ brought by future data will be limited by production cross-section uncertainties, unless either $^{10}$Be/$^9$Be measurements are extended up to several tens of GeV/n or nuclear data for the production of $^{10}$Be and $^{9}$Be are improved; new data for the production cross section of $^{16}$O into Be isotopes above a few GeV/n are especially desired.

Theoretical works have shown that off-plane motions of bars can heat stars in the vertical direction during buckling but is not clear how do they affect the rest of components of the Stellar Velocity Ellipsoid (SVE). We study the 2D spatial distribution of the vertical, $\sigma_{z}$, azimuthal, $\sigma_{\phi}$ and radial, $\sigma_{r}$ velocity dispersions in the inner regions of Auriga galaxies, a set of high-resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the influence of the bar on the stellar kinematics. $\sigma_{z}$ and $\sigma_{\phi}$ maps exhibit non-axisymmetric features that closely match the bar light distribution with low $\sigma$ regions along the bar major axis and high values in the perpendicular direction. On the other hand, $\sigma_{r}$ velocity dispersion maps present more axisymmetric distributions. We show that isophotal profile differences best capture the impact of the bar on the three SVE components providing strong correlations with bar morphology proxies although there is no relation with individual $\sigma$. Time evolution analysis shows that these differences are a consequence of the bar formation and that they tightly coevolve with the strength of the bar. We discuss the presence of different behaviours of $\sigma_{z}$ and its connection with observations. This work helps us understand the intrinsic $\sigma$ distribution and motivates the use of isophotal profiles as a mean to quantify the effect of bars.

Joint studies of imaging and spectroscopic samples, informed by theory and simulations, offer the potential for comprehensive tests of the cosmological model over redshifts z<1.5. Spectroscopic galaxy samples at these redshifts can be increased beyond the planned Dark Energy Spectroscopic Instrument (DESI) program by at least an order of magnitude, thus offering significantly more constraining power for these joint studies. Spectroscopic observations of these galaxies in the latter half of the 2020's and beyond would leverage the theory and simulation effort in this regime. In turn, these high density observations will allow enhanced tests of dark energy, physics beyond the standard model, and neutrino masses that will greatly exceed what is currently possible. Here, we present a coordinated program of simulations, theoretical modeling, and future spectroscopy that would enable precise cosmological studies in the accelerating epoch where the effects of dark energy are most apparent.

We present an investigation of the internal kinematic properties of M79 (NGC 1904). Our study is based on radial velocity measurements obtained from the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters for more than 1700 individual stars distributed between $\sim 0.3^{\prime\prime}$ and $770^{\prime\prime}$ ($\sim14$ three-dimensional half-mass radii), from the center. Our analysis reveals the presence of ordered line-of-sight rotation with a rotation axis almost aligned along the East-West direction and a velocity peak of $1.5$ km s$^{-1}$ at $\sim 70^{\prime\prime}$ from the rotation axis. The velocity dispersion profile is well described by the same King model that best fits the projected density distribution, with a constant central plateau at $\sigma_0\sim 6$ km s$^{-1}$. To investigate the cluster rotation in the plane of the sky, we have analyzed the proper motions provided by the Gaia EDR3, finding a signature of rotation with a maximum amplitude of $\sim 2.0$ km s$^{-1}$ at $\sim 80^{\prime\prime}$ from the cluster center. Analyzing the three-dimensional velocity distribution, for a sub-sample of 130 stars, we confirm the presence of systemic rotation and find a rotation axis inclination angle of $37${\deg} with respect to the line-of-sight. As a final result, the comparison of the observed rotation curves with the results of a representative N-body simulation of a rotating star cluster shows that the present-day kinematic properties of NGC 1904 are consistent with those of a dynamically old system that has lost a significant fraction of its initial angular momentum.

Space-based photometric missions widely use statistical validation tools for vetting transiting planetary candidates, particularly when other traditional methods of planet confirmation are unviable. In this paper, we refute the planetary nature of three previously validated planets -- Kepler-854~b, Kepler-840~b, and Kepler-699~b -- and possibly a fourth, Kepler-747~b, using updated stellar parameters from Gaia and phase-curve analysis. In all four cases, the inferred physical radii rule out their planetary nature given the stellar radiation the companions receive. For Kepler-854~b, the mass derived from the host star's ellipsoidal variation, which had not been part of the original vetting procedure, similarly points to a non-planetary value. To contextualize our understanding of the phase curve for stellar mass companions in particular and extend our understanding of high-order harmonics, we examine Kepler eclipsing binaries with periods between 1.5 and 10 days. Using a sample of 20 systems, we report a strong power-law relation between the second cosine harmonic of the phase-curve signal and the higher cosine harmonics, which supports the hypothesis that those signals arise from the tidal interaction between the binary components. We find that the ratio between the second and third-harmonic amplitudes is $2.24 \pm 0.48$, in good agreement with the expected value of 2.4 from the classical formalism for the ellipsoidal distortion.

In recent times, large organic molecules of exceptional complexity have been found in diverse regions of the interstellar medium. In this context, we aim to provide accurate frequencies of the ground vibrational state of two key aliphatic aldehydes, n-butanal and its branched-chain isomer, i-butanal, to enable their eventual detection in the interstellar medium. We employ a frequency modulation millimeter-wave absorption spectrometer to measure the rotational features of n- and i-butanal. We use the spectral line survey ReMoCA performed with the Atacama Large Millimeter/submillimeter Array to search for n- and i-butanal toward the star-forming region Sgr B2(N). We also search for both aldehydes toward the molecular cloud G+0.693-0.027 with IRAM 30 m and Yebes 40 m observations. Several thousand rotational transitions belonging to the lowest-energy conformers have been assigned in the laboratory spectra up to 325 GHz. A precise set of the relevant rotational spectroscopic constants has been determined for each structure. We report non-detections of n- and i-butanal toward both sources, Sgr B2(N1S) and G+0.693-0.027. We find that n- and i-butanal are at least 2-6 and 6-18 times less abundant than acetaldehyde toward Sgr B2(N1S), respectively, and that n-butanal is at least 63 times less abundant than acetaldehyde toward G+0.693-0.027. Comparison with astrochemical models indicates good agreement between observed and simulated abundances (where available). Grain-surface chemistry appears sufficient to reproduce aldehyde ratios in G+0.693-0.027; gas-phase production may play a more active role in Sgr B2(N1S). Our astronomical results indicate that the family of interstellar aldehydes in the Galactic center region is characterized by a drop of one order of magnitude in abundance at each incrementation in the level of molecular complexity.

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time.

Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass.

We study the theory and phenomenology of massive spin-2 fields during the inflation with nonzero background chemical potential, and extend the cosmological collider physics to tensor modes. We identify a unique dimension-5 and parity-violating chemical potential operator for massive spin-2 fields, which leads to a ghost-free linear theory propagating one scalar mode and two tensor modes. The chemical potential greatly boosts the production of one tensor mode even for very heavy spin-2 particles, and thereby leads to large and distinct cosmological collider signals for massive spin-2 particles. The large signals show up at the tree-level in both the curvature trispectrum and the tensor-curvature mixed bispectrum.

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

We outline the unique opportunities and challenges in the search for "ultraheavy" dark matter candidates with masses between roughly $10~{\rm TeV}$ and the Planck scale $m_{\rm pl} \approx 10^{16}~{\rm TeV}$. This mass range presents a wide and relatively unexplored dark matter parameter space, with a rich space of possible models and cosmic histories. We emphasize that both current detectors and new, targeted search techniques, via both direct and indirect detection, are poised to contribute to searches for ultraheavy particle dark matter in the coming decade. We highlight the need for new developments in this space, including new analyses of current and imminent direct and indirect experiments targeting ultraheavy dark matter and development of new, ultra-sensitive detector technologies like next-generation liquid noble detectors, neutrino experiments, and specialized quantum sensing techniques.

Magnetic reconnection plays an important role in converting energy while modifying field topology. This process takes place in varied plasma environments in which the transport of magnetic flux is intrinsic. Identifying active magnetic reconnection sites in in-situ observations is challenging. A new technique, Magnetic Flux Transport (MFT) analysis, has been developed recently and proven in numerical simulation for identifying active reconnection efficiently and accurately. In this study, we examine the MFT process in 37 previously reported electron diffusion region (EDR)/reconnection-line crossing events at the dayside magnetopause and in the magnetotail and turbulent magnetosheath using Magnetospheric Multiscale measurements. The coexisting inward and outward MFT flows at an X-point provides a signature that magnetic field lines become disconnected and reconnected. The application of MFT analysis to in-situ observations demonstrates that MFT can successfully identify active reconnection sites under complex varied conditions, including asymmetric and turbulent upstream conditions. It also provides a higher rate of identification than plasma outflow jets alone. MFT can be applied to in situ measurements from both single- and multi-spacecraft missions and laboratory experiments.

Theoretical investigations into the evolution of the early universe are an essential part of particle physics that allow us to identify viable extensions to the Standard Model as well as motivated parameter space that can be probed by various experiments and observations. In this white paper, we review particle physics models of the early universe. First, we outline various models that explain two essential ingredients of the early universe (dark matter and baryon asymmetry) and those that seek to address current observational anomalies. We then discuss dynamics of the early universe in models of neutrino masses, axions, and several solutions to the electroweak hierarchy problem. Finally, we review solutions to naturalness problems of the Standard Model that employ cosmological dynamics.

This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dark matter searches using photons -- radio and gamma rays -- to cosmic rays, neutrinos and gravitational waves. Examples include astrophysical sources of cosmic messengers and their interactions which can mimic dark matter signatures. In turn, these depend on commensurate studies in understanding the cosmic environment -- gas distributions, magnetic field configurations -- as well as relevant nuclear astrophysics. We also cover the astrophysics governing celestial bodies and galaxies used to probe dark matter, from black holes to dwarf galaxies. Finally, we cover astrophysical backgrounds related to probing the dark matter distribution and kinematics, which impact a wide range of dark matter studies. In the future, the rise of multi-messenger astronomy, and novel analysis methods to exploit it for dark matter, will offer various strategic ways to continue to enhance our understanding of astrophysical backgrounds to deliver improved sensitivity to dark matter.

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online.

Intriguing signals with excesses over expected backgrounds have been observed in many astrophysical and terrestrial settings, which could potentially have a dark matter origin. Astrophysical excesses include the Galactic Center GeV gamma-ray excess detected by the Fermi Gamma-Ray Space Telescope, the AMS antiproton and positron excesses, and the 511 and 3.5 keV X-ray lines. Direct detection excesses include the DAMA/LIBRA annual modulation signal, the XENON1T excess, and low-threshold excesses in solid state detectors. We discuss avenues to resolve these excesses, with actions the field can take over the next several years.

This white paper discusses the current landscape and prospects for experiments sensitive to particle dark matter processes producing photons and cosmic rays. Much of the gamma-ray sky remains unexplored on a level of sensitivity that would enable the discovery of a dark matter signal. Currently operating GeV-TeV observatories, such as Fermi-LAT, atmospheric Cherenkov telescopes, and water Cherenkov detector arrays continue to target several promising dark matter-rich environments within and beyond the Galaxy. Soon, several new experiments will continue to explore, with increased sensitivity, especially extended targets in the sky. This paper reviews the several near-term and longer-term plans for gamma-ray observatories, from MeV energies up to hundreds of TeV. Similarly, the X-ray sky has been and continues to be monitored by decade-old observatories. Upcoming telescopes will further bolster searches and allow new discovery space for lines from, e.g., sterile neutrinos and axion-photon conversion. Furthermore, this overview discusses currently operating cosmic-ray probes and the landscape of future experiments that will clarify existing persistent anomalies in cosmic radiation and spearhead possible new discoveries. Finally, the article closes with a discussion of necessary cross section measurements that need to be conducted at colliders to reduce substantial uncertainties in interpreting photon and cosmic-ray measurements in space.

Many parity violating gravity models suffer from the ghost instability problem. In this paper we consider a symmetric teleparallel gravity model which extends the general relativity equivalent model by several parity violating interactions between the gravitational field and a scalar field. These interactions exclude higher derivatives and are quadratic in the non-metricity tensor. Through investigations on the linear cosmological perturbations, our results show that in general this model suffers from the difficulty caused by the existence of ghost mode in the vector perturbations. However, in cases where a special condition is imposed on the coefficients of these parity violating interactions, this model can be ghost free.

As the number of gravitational wave observations has increased in recent years, the variety of sources has broadened. Here we investigate whether it is possible for the current generation of detectors to distinguish between very short-lived gravitational wave signals from mergers between high-mass black holes, and the signal produced by a close encounter between two black holes which results in gravitational capture, and ultimately a merger. We compare the posterior probability distributions produced by analysing simulated signals from both types of progenitor events, both under ideal and realistic scenarios. We show that while, under ideal conditions it is possible to distinguish both progenitors, under more realistic conditions they are indistinguishable. This has important implications for the interpretation of such short signals, and we therefore advocate that these signals be the focus of additional investigation even when satisfactory results have been achieved from standard analyses.

The state of the space environment plays a significant role for forecasting of geomagnetic storms produced by disturbances of the solar wind (SW). Coronal mass ejections (CMEs) passing through the heliosphere often have a prolonged (up to several days) trail with declining speed, which affects propagation of the subsequent SW streams. We studied the CME and the post-eruption plasma flows behind the CME rear in the event on 2010 August 18 observed in quadrature by several space-based instruments. Observations of the eruption in the corona with EUV telescopes and coronagraphs revealed several discrete outflows followed by a continuous structureless post-eruption stream. The interplanetary coronal mass ejection (ICME), associated with this CME, was registered by PLAsma and SupraThermal Ion Composition (PLASTIC) instrument aboard the Solar TErrestrial RElations Observatory (STEREO-A) between August 20, 16:14 UT and August 21, 13:14 UT, after which the SW disturbance was present over 3 days. Kinematic consideration with the use of the gravitational and Drag-Based models has shown that the discrete plasma flows can be associated with the ICME, whereas the post-eruption outflow was arrived in the declining part of the SW transient. We simulated the Fe-ion charge distributions of the ICME and post-CME parts of SW using the plasma temperature and density in the ejection region derived from the Differential Emission Measure analysis. The results demonstrate that in the studied event the post-ICME trailing region was associated with the post-eruption flow from the corona, rather then with the ambient SW entrained by the CME.

We describe the mapping at high density of topological structure of baryonic matter to a nuclear effective field theory that implements hidden symmetries emergent from strong nuclear correlations. The theory so constructed is found to be consistent with no conflicts with the presently available observations in both normal nuclear matter and compact-star matter. The hidden symmetries involved are "local flavor symmetry" of the vector mesons identified to be (Seiberg-)dual to the gluons of QCD and hidden "quantum scale symmetry" with an IR fixed point with a "genuine dilaton (GD)" characterized by non-vanishing pion and dilaton decay constants. Both the skyrmion topology for $N_f \geq 2$ baryons and the fractional quantum Hall (FQH) droplet topology for $N_f=1$ baryons are unified in the "homogeneous/hidden" Wess-Zumino term in the hidden local symmetry (HLS) Lagrangian. The possible indispensable role of the FQH droplets in going beyond the density regime of compact stars approaching scale-chiral restoration is explored by moving toward the limit where both the dilaton and the pion go massless.

Experimental searches for new, "fifth" forces are attracting a lot of attention because they allow to test theoretical extensions to the standard model. Here, we report a new experimental search for possible fifth forces, specifically spin-and-velocity dependent forces, by using a K-Rb-$^{21}$Ne co-magnetometer and a tungsten ring featuring a high nucleon density. Taking advantage of the high sensitivity of the co-magnetometer, the pseudomagnetic field from the fifth force is measured to be $<7$\,aT. This sets new limits on coupling constants for the neutron-nucleon and proton-nucleon interactions in the range of $\ge 0.1$ m. The coupling constant limits are established to be $|g_V^n|<6.6\times 10^{-11}$ and $|g_V^p|<3.0\times 10^{-10}$, which are more than one order of magnitude tighter than astronomical and cosmological limits on the coupling between the new gauge boson such as Z$'$ and standard model particles.

The nature of dark matter is one of the major puzzles of fundamental physics, integral to the understanding of our universe across almost every epoch. The search for dark matter takes place at different energy scales, and use data ranging from particle colliders to astrophysical surveys. We focus here on CMB-S4, a future ground-based Cosmic Microwave Background (CMB) experiment, which is expected to provide exquisite measurements of the CMB temperature and polarization anisotropies. These measurements (on their own and in combination with other surveys) will allow for new means to shed light on the nature of dark matter.

There seems to exist agreement about the fact that inflation squeezes the quantum state of cosmological perturbations and entangles modes with wavenumbers $\vec k$ and $-\vec k$. Paradoxically, this result has been used to justify both the classicality as well as the quantumness of the primordial perturbations at the end of inflation. We reexamine this question and point out that the definition of two-mode squeezing of the modes $\vec k$ and $-\vec k$ used in previous work rests on choices that are only justified for systems with time-independent Hamiltonians and finitely many degrees of freedom. We argue that for quantum fields propagating on generic time-dependent Friedmann-Lema\^itre-Robertson-Walker backgrounds, the notion of squeezed states is subject to ambiguities, which go hand in hand with the ambiguity in the definition of particles. In other words, we argue that the question "does the cosmic expansion squeeze and entangle modes with wavenumbers $\vec k$ and $-\vec k$?" contains the same ambiguity as the question "does the cosmic expansion create particles?". When additional symmetries are present, like in the (quasi) de Sitter-like spacetimes used in inflationary models, one can resolve the ambiguities, and we find that the answer to the question in the title turns out to be in the negative. We further argue that this fact does not make the state of cosmological perturbations any less quantum, at least when deviations from Gaussianity can be neglected.

We consider conformal vector models which could play the role of a cosmological dark radiation component. We analyse the propagation of gravitational waves in the presence of this vector background and find a suppression in the tensor transfer function at large scales. We also find that although the cosmological background metric is isotropic, anisotropies are imprinted in the tensor power spectrum. In addition, the presence of the background vector fields induces a net polarization of the gravitational wave background and, for certain configurations of the vector field, a linear to circular polarization conversion. We also show that this kind of effects are also present for vector models with more general potential terms.

The absence of clear signals from particle dark matter in direct detection experiments motivates new approaches in disparate regions of viable parameter space. In this Snowmass white paper, we outline the Windchime project, a program to build a large array of quantum-enhanced mechanical sensors. The ultimate aim is to build a detector capable of searching for Planck mass-scale dark matter purely through its gravitational coupling to ordinary matter. In the shorter term, we aim to search for a number of other physics targets, especially some ultralight dark matter candidates. Here, we discuss the basic design, open R&D challenges and opportunities, current experimental efforts, and both short- and long-term physics targets of the Windchime project.

Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities.

This white paper provides a comprehensive review of our present understanding of experimental neutrino anomalies that remain unresolved, charting the progress achieved over the last decade at the experimental and phenomenological level, and sets the stage for future programmatic prospects in addressing those anomalies. It is purposed to serve as a guiding and motivational "encyclopedic" reference, with emphasis on needs and options for future exploration that may lead to the ultimate resolution of the anomalies. We see the main experimental, analysis, and theory-driven thrusts that will be essential to achieving this goal being: 1) Cover all anomaly sectors -- given the unresolved nature of all four canonical anomalies, it is imperative to support all pillars of a diverse experimental portfolio, source, reactor, decay-at-rest, decay-in-flight, and other methods/sources, to provide complementary probes of and increased precision for new physics explanations; 2) Pursue diverse signatures -- it is imperative that experiments make design and analysis choices that maximize sensitivity to as broad an array of these potential new physics signatures as possible; 3) Deepen theoretical engagement -- priority in the theory community should be placed on development of standard and beyond standard models relevant to all four short-baseline anomalies and the development of tools for efficient tests of these models with existing and future experimental datasets; 4) Openly share data -- Fluid communication between the experimental and theory communities will be required, which implies that both experimental data releases and theoretical calculations should be publicly available; and 5) Apply robust analysis techniques -- Appropriate statistical treatment is crucial to assess the compatibility of data sets within the context of any given model.

The Snowmass Community Survey was designed by the Snowmass Early Career (SEC) Survey Core Initiative team between April 2020 and June 2021, and released to the community on June 28, 2021. It aims to be a comprehensive assessment of the state of the high-energy particle and astrophysics (HEPA) community, if not the field, though the Snowmass process is largely based within the United States. Among other topics, some of the central foci of the Survey were to gather demographic, career, physics outlook, and workplace culture data on a large segment of the Snowmass community. With nearly $1500$ total interactions with the Survey, the SEC Survey team hopes the findings and discussions within this report will be of service to the community over the next decade. Some conclusions should reinforce the aspects of HEPA which are already functional and productive, while others should strengthen arguments for cultural and policy changes within the field.

The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics.

Coherent elastic neutrino-nucleus scattering (CE$\nu$NS) is a process in which neutrinos scatter on a nucleus which acts as a single particle. Though the total cross section is large by neutrino standards, CE$\nu$NS has long proven difficult to detect, since the deposited energy into the nucleus is $\sim$ keV. In 2017, the COHERENT collaboration announced the detection of CE$\nu$NS using a stopped-pion source with CsI detectors, followed up the detection of CE$\nu$NS using an Ar target. The detection of CE$\nu$NS has spawned a flurry of activities in high-energy physics, inspiring new constraints on beyond the Standard Model (BSM) physics, and new experimental methods. The CE$\nu$NS process has important implications for not only high-energy physics, but also astrophysics, nuclear physics, and beyond. This whitepaper discusses the scientific importance of CE$\nu$NS, highlighting how present experiments such as COHERENT are informing theory, and also how future experiments will provide a wealth of information across the aforementioned fields of physics.