Papers for Monday, Nov 29 2021

The information content of the minimum spanning tree (MST), used to capture higher-order statistics and other information from the cosmic web, is compared to that of the power spectrum for a $\nu\Lambda$CDM model. The measurements are made in redshift space using haloes from the Quijote simulation of mass $\geq 3.2\times 10^{13}\,h^{-1}{\rm M}_{\odot}$ in a box of length $L_{\rm box}=1\,h^{-1}{\rm Gpc}$. The power spectrum multipoles (monopole and quadrupole) are computed for Fourier modes in the range $0.006 < k < 0.5\, h{\rm Mpc}^{-1}$. For comparison the MST is measured with a minimum length scale of $l_{\min}\simeq13\,h^{-1}{\rm Mpc}$. Combining the MST and power spectrum allows for many of the individual degeneracies to be broken; on its own the MST provides tighter constraints on the sum of neutrino masses $M_{\nu}$, Hubble constant $h$, spectral tilt $n_{\rm s}$, and baryon energy density $\Omega_{\rm b}$ but the power spectrum alone provides tighter constraints on $\Omega_{\rm m}$ and $\sigma_{8}$. The power spectrum on its own gives a standard deviation of $0.25\,{\rm eV}$ on $M_{\nu}$ while the combination of power spectrum and MST gives $0.11\,{\rm eV}$. There is similar improvement of a factor of two for $h$, $n_{\rm s}$, and $\Omega_{\rm b}$. These improvements appear to be driven by the MST's sensitivity to small scale clustering, where the effect of neutrino free-streaming becomes relevant. The MST is shown to be a powerful tool for cosmology and neutrino mass studies, and therefore could play a pivotal role in ongoing and future galaxy redshift surveys (such as DES, DESI, Euclid, and Rubin-LSST).

HR 8799 is a young A5/F0 star hosting four directly imaged giant planets at wide separations ($\sim$16-78 au) which are undergoing orbital motion and have been continuously monitored with adaptive optics imaging since their discovery over a decade ago. We present a dynamical mass of HR 8799 using 130 epochs of relative astrometry of its planets, which include both published measurements and new medium-band 3.1 $\mu$m observations that we acquired with NIRC2 at Keck Observatory. For the purpose of measuring the host star mass, each orbiting planet is treated as a massless particle and is fit with a Keplerian orbit using Markov chain Monte Carlo. We then use a Bayesian framework to combine each independent total mass measurement into a cumulative dynamical mass using all four planets. The dynamical mass of HR 8799 is 1.47$^{+0.12}_{-0.17}$ \Msun assuming a uniform stellar mass prior, or 1.46$^{+0.11}_{-0.15}$ \Msun with a weakly informative prior based on spectroscopy. There is a strong covariance between the planets' eccentricities and the total system mass; when the constraint is limited to low eccentricity solutions of $e<0.1$, which is motivated by dynamical stability, our mass measurement improves to 1.43$^{+0.06}_{-0.07}$ \Msun. Our dynamical mass and other fundamental measured parameters of HR 8799 together with MESA Isochrones & Stellar Tracks grids yields a bulk metallicity most consistent with [Fe/H]$\sim$ -0.25-0.00 dex and an age of 10-23 Myr for the system. This implies hot start masses of 2.7-4.9 \Mjup for HR 8799 b and 4.1-7.0 \Mjup for HR 8799 c, d, and e, assuming they formed at the same time as the host star.

Using the novel TESS 20-sec cadence data, we have discovered an unusual combination of pulsating stars in what we infer to be a binary system. The primary is a standard $\delta$ Scuti star with pulsations over the range 32-41 d$^{-1}$; this is in an inferred wide orbit with a hot subdwarf B (sdB) secondary, which itself has a large-amplitude p-mode pulsation at 524 d$^{-1}$. We establish constraints on the period of the putative binary by using radial velocity measurements of the $\delta$ Scuti star and show that any sdB-mass companion star must have an orbital period greater than $\sim 35$ d. Our identification of this sdB binary serves as an important addition to the relatively small amount of sdB binaries known to have a period longer than a few days. We show that such a binary can be formed through stable, nonconservative mass transfer, without undergoing a common envelope phase.

Accreting planets have been seen at Ha (H alpha), but targeted searches have not been fruitful. For planets, accretion tracers should come from the shock itself, exposing them to extinction by the accreting material. High-resolution (R>5e4) spectrographs at Ha should soon allow studying how the incoming material shapes the line profile. We calculate how much the gas and dust accreting onto a planet reduce the Ha flux from the shock at the planetary surface and how they affect line shapes. We also study the absorption-modified relationship between Ha luminosity and Mdot. We compute the high-resolution radiative transfer of the Ha line using a 1D velocity-density-temperature structure for the inflowing matter in three representative accretion geometries: spherical symmetry, polar inflow, and magnetospheric accretion. For each, we explore wide ranges of Mdot and planet mass M. We use detailed gas opacities and estimate dust opacities. At Mdot<3e-6 MJ/yr, gas extinction is negligible for spherical or polar inflow and at most A_Ha<0.5 mag for magnetospheric accretion. Up to Mdot~3e-4 MJ/yr, the gas has A_Ha<4 mag. This decreases with M. We estimate realistic dust opacities at Ha as ~0.01-10 cm^2/g, i.e., 10-1e4 times lower than in the ISM. Extinction flattens the L_Ha-Mdot relationship, which becomes non-monotonic with a maximum L_Ha~1e-4 LSun near Mdot~1e-4 MJ/yr for M~10 MJ. In magnetospheric accretion, the gas can introduce features in line profiles, but the velocity gradient smears them out in other geometries. For most of parameter space, extinction by the accreting matter should be negligible, simplifying interpretation of observations, especially for planets in gaps. At high Mdot, strong absorption reduces the Ha flux, and some measurements can be interpreted as two Mdot values. Line profiles at R~1e5 can provide complex constraints on the accretion flow's thermal-dynamical structure.

Previous studies have shown that dark matter deficient galaxies (DMDG) such as NGC1052-DF2 (hereafter DF2) can result from tidal stripping. An important question, though, is whether such a stripping scenario can explain DF2's large specific frequency of globular clusters (GCs). After all, tidal stripping and shocking preferentially remove matter from the outskirts. We examine this using idealized, high-resolution simulations of a regular dark matter dominated galaxy that is accreted onto a massive halo. As long as the initial (pre-infall) dark matter halo of the satellite is cored, which is consistent with predictions of cosmological, hydrodynamical simulations, the tidal remnant can be made to resemble DF2 in all its properties, including its GC population. The required orbit has a peri-centre at the 8.3 percentile of the distribution for subhaloes at infall, and thus is not particularly extreme. On this orbit the satellite loses 98.5 (30) percent of its original dark matter (stellar) mass, and thus evolves into a DMDG. The fraction of GCs that is stripped off depends on the initial radial distribution. If, at infall, the median projected radius of the GC population is roughly two times that of the stars, consistent with observations of isolated galaxies, only $\sim 20$ percent of the GCs are stripped off. This is less than for the stars, which is due to dynamical friction counteracting the tidal stirring. We predict that, if indeed DF2 was crafted by strong tides, its stellar outskirts should have a very shallow metallicity gradient.

We present spectral and temporal properties of all the thermonuclear X-ray bursts observed from Aql X-1 by the Neutron Star Interior and Composition Explorer (NICER) between 2017 July and 2021 April. This is the first systematic investigation of a large sample of type I X-ray bursts from Aql X-1 with improved sensitivity at low energies. We detect 22 X-ray bursts including two short recurrence burst events in which the separation was only 451 s and 496 s. We perform time resolved spectroscopy of the bursts using the fixed and scaled background (f_a method) approaches. We show that the use of a scaling factor to the pre-burst emission is the statistically preferred model in about 68% of all the spectra compared to the fixed background approach. Typically the f_a values are clustered around 1-3, but can reach up to 11 in a burst where photospheric radius expansion is observed. Such f_a values indicate a very significant increase in the pre-burst emission especially at around the peak flux moments of the bursts. We show that the use of the f_a factor alters the best fit spectral parameters of the burst emission. Finally, we employed a reflection model instead of scaling the pre-burst emission. We show that reflection models also do fit the spectra and improve the goodness of the fits. In all cases we see that the disc is highly ionized by the burst emission and the fraction of the reprocessed emission to the incident burst flux is typically clustered around 20%.

We present observations of SN 2021csp, the second example of a newly-identified type of supernova (Type Icn) hallmarked by strong, narrow, P Cygni carbon features at early times. The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of -20 within 3 days due to strong interaction between fast SN ejecta (v ~ 30000 km/s) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow line features disappear from the spectrum 10-20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ~60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic supernovae. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss various progenitor scenarios: an ultrastripped progenitor star, a pulsational pair-instability eruption, or a jet-driven fallback supernova from a Wolf-Rayet star. The fallback scenario would naturally explain the similarity between these events and radio-loud fast transients, and suggests a picture in which most stars massive enough to undergo a WR phase collapse directly to black holes at the end of their lives.

Modeling the evolution of progenitors of gravitational-wave merger events in binary stars faces two major uncertainties: the common-envelope phase and supernova kicks. These two processes are critical for the final orbital configuration of double compact-object systems with neutron stars and black holes. Predictive one-dimensional models of common-envelope interaction are lacking and multidimensional simulations are challenged by the vast range of relevant spatial and temporal scales. Here, we present three-dimensional hydrodynamic simulations of the common-envelope interaction of an initially $10\,M_{\odot}$ red supergiant primary star with a black-hole and a neutron-star companion. We show that the high-mass regime is accessible to full ab-initio simulations. Nearly complete envelope ejection is reached assuming that all recombination energy still available at the end of our simulation continues to help unbinding the envelope. In contrast to previous assumptions, we find that the dynamical plunge-in of both companions terminates at orbital separations too wide for gravitational waves to merge the systems in a Hubble time. We discuss the further evolution of the system based on analytical estimates. A subsequent mass-transfer episode from the remaining $3\,M_{\odot}$ core of the supergiant to the compact companion does not shrink the orbit sufficiently either. A neutron-star--neutron-star and neutron-star--black-hole merger is still expected for a fraction of the systems if the supernova kick aligns favorably with the orbital motion. For double neutron star (neutron-star--black-hole) systems we estimate mergers in about $9 \%$ ($1 \%$) of cases while about $77 \%$ ($94 \%$) of binaries are disrupted, i.e., supernova kicks actually enable gravitational-wave mergers in our cases; however, we expect a reduction in predicted gravitational-wave merger events. (abbr.)

Photometric redshift estimation algorithms are often based on representative data from observational campaigns. Data-driven methods of this type are subject to a number of potential deficiencies, such as sample bias and incompleteness. Motivated by these considerations, we propose using physically motivated synthetic spectral energy distributions in redshift estimation. In addition, the synthetic data would have to span a domain in colour-redshift space concordant with that of the targeted observational surveys. With a matched distribution and realistically modelled synthetic data in hand, a suitable regression algorithm can be appropriately trained; we use a mixture density network for this purpose. We also perform a zero-point re-calibration to reduce the systematic differences between noise-free synthetic data and the (unavoidably) noisy observational data sets. This new redshift estimation framework, SYTH-Z, demonstrates superior accuracy over a wide range of redshifts compared to baseline models trained on observational data alone. Approaches using realistic synthetic data sets can therefore greatly mitigate the reliance on expensive spectroscopic follow-up for the next generation of photometric surveys.

Recent imaging of the disc around the V4046$\,$Sgr spectroscopic binary revealed concentric regions of dust rings and gaps. The object's proximity and expected equilibrated state due to its old age (>20$\,$Myr) make it a superb testbed for hydrodynamical studies in direct comparison to observations. We employ two-dimensional hydrodynamical simulations of gas and multiple dust species to test whether the observed structure conforms with the presence of giant planets embedded in the disc. We then perform radiative transfer calculations of sky images, which we filter for the telescope response for comparison with near-infrared and millimetre observations. We find that the existing data are in excellent agreement with a flared disc and the presence of two giant planets, at 9$\,$au and 20$\,$au, respectively. The different ring widths are recovered by diffusion-balanced dust trapping within the gas pressure maxima. In our radiative transfer model, the diffusion in vertical direction is reduced in comparison to the radial value by a factor of five to recover the spectral energy distribution. Further, we report a previously unaddressed, azimuthally-confined intensity decrement on the bright inner ring in the near-infrared scattered light observation. Our model shows that this decrement can be explained by a shadow cast by a circumplanetary disc around the same giant planet that creates the inner cavity in the hydrodynamical simulations. We examine the shape of the intensity indentation and discuss the potential characterisation of a giant planet and its associated disc by its projected shadow in scattered light observations.

The Zwicky Transient Facility (ZTF), a state-of-the-art optical robotic sky survey, registers on the order of a million transient events - such as supernova explosions, changes in brightness of variable sources, or moving object detections - every clear night, and generates associated real-time alerts. We present Alert-Classifying Artificial Intelligence (ACAI), an open-source deep-learning framework for the phenomenological classification of ZTF alerts. ACAI uses a set of five binary classifiers to characterize objects which, in combination with the auxiliary/contextual event information available from alert brokers, provides a powerful tool for alert stream filtering tailored to different science cases, including early identification of supernova-like and anomalous transient events. We report on the performance of ACAI during the first months of deployment in a production setting.

A large fraction of white dwarfs are accreting or have recently accreted rocky material from their planetary systems, thereby polluting their atmospheres with elements heavier than helium. In recent years, the quest for mechanisms that can deliver planetesimals to the immediate vicinity of their central white dwarfs has stimulated a flurry of modelling efforts. The observed time evolution of the accretion rates of white dwarfs through their multi-Gyr lifetime is a crucial test for dynamical models of evolved planetary systems. Recent studies of cool white dwarfs samples have identified a significant decrease of the mass accretion rates of cool, old white dwarfs over Gyr timescales. Here, we revisit those results using updated white dwarf models and larger samples of old polluted H- and He-atmosphere white dwarfs. We find no compelling evidence for a strong decrease of their time-averaged mass accretion rates for cooling times between 1 and 8 Gyrs. Over this period, the mass accretion rates decrease by no more than a factor of the order of 10, which is one order of magnitude smaller than the decay rate found in recent works. Our results require mechanisms that can efficiently and consistently deliver planetesimals inside the Roche radius of white dwarfs over at least 8 Gyrs.

The star NGC2004#115 in the LMC, originally classified as an (SB1) Be spectroscopic binary, bears some morphological resemblance to the Galactic systems LB-1 and HR 6819, both of which are proposed as either Be+black hole (BH) or Be+stripped He-star systems. Two data-sets (ESO/VLT and SALT) of multi-epoch optical spectra of NGC 2004#115, separated by baseline of $\sim$20 years, lead us conclude it is a triple system hosting an inner binary with a period of 2.92 d, eccentricity $\sim$0.0 and mass function $\sim$0.07 $M_\odot$. The inner binary harbours a B-type star (the primary) with projected rotational velocity of 10km/s, and luminosity $\log L/L_\odot$=3.87, contributing $\sim$60% of the V-band light to the system. The secondary is not detected, while the tertiary, which contributes 40% of the light, is tentatively identified as a less luminous B-type star with high projected rotational velocity. No ellipsoidal light variability is detected, with stringent limits being set by MACHO and Gaia data. Assuming the primary to be a main sequence star yields a mass of 8.6$ M_\odot$, while the additional assumption of synchronous rotation constrains the inclination to be almost pole-on with i~9 degrees, implying the secondary is a BH with a mass of $\sim$25 $M_\odot$. A low mass stripped star with similar luminosity is ruled out as a potential solution as its mass implies a Roche radius that is substantially smaller than the stellar radius. The outer period likely exceeds 120 days and, while the disk-like emission is variable (it is almost absent in the SALT dataset), it may be associated with the inner binary rather than the rapidly rotating tertiary. XMM-Newton provides an upper limit of 5x$10^{33}$ ergs/s on the X-ray flux, consistent with, though not constraining of, the system hosting a quiescent B+BH binary. A number of caveats to this scenario are discussed in the paper.

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. As well as creating the $\ell$=1 mode of the CMB sky, this motion affects all astrophysical observations by modulating and aberrating sources across the sky. It can be seen in galaxy clustering, and in principle its time derivative through a dipole-shaped acceleration pattern in quasar positions. Additionally, the dipole modulates the CMB temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect and so these modulated CMB anisotropies can be extracted from the tSZ maps produced by Planck. Unfortunately, this measurement cannot determine if the dipole is due to our motion, but it does provide an independent measure of the dipole and a validation of the y maps. This measurement, and a description of the first-order terms of the CMB dipole, are outlined here.

The development of large constellations of satellites (i.e., so-called megaconstellations or satcons) is poised to increase the number of LEO satellites by more than an order of magnitude in the coming decades. Such a rapid growth of satellite numbers makes the consequences of major fragmentation events ever more problematic. In this study, we investigate the collisional risk to on-orbit infrastructure from kinetic anti-satellite (ASAT) weapon tests, using the 2019 India test as a model. We find that the probability of one or more collisions occurring over the lifetime of ASAT fragments increases significantly in a satcon environment compared with the orbital environment in 2019. For the case of 65,000 satellites in LEO, we find that the chance of one or more satellites being struck by ASAT fragments of size 1 cm or larger is about 30% for a single test. Including sizes down to 3 mm in our models suggests that impacts will occur for any such event. The heavy commercialization of LEO demands a commitment to avoiding debris-generating ASAT tests.

The effect of ionized gas on the propagation of radio signals is known as plasma lensing. Unlike gravitational lensing, plasma lensing causes both magnification and strong de-magnification effects to background sources. We study the cross section of plasma lensing for two density profiles, the Gaussian and power-law models. In general, the cross section increases with the density gradient. Radio sources can be used to measure the free electron density along the line of sight. However, plasma lensing de-magnification causes an underestimate of the electron density. Such a bias increases with the electron density, and can be up to $\sim 15\%$ in the high density region. There is a large probability that high density clumps will be missed due to this bias. The magnification of plasma lensing can also change the luminosity function of the background sources. The number density of sources on both the high and low luminosity ends can be overestimated due to this biasing effect.

We investigate the internal 3D magnetic structure of dense interstellar filaments within NGC 1333 using polarization data at $850 \mu\mathrm{m}$ from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope (JCMT). Theoretical models predict that the magnetic field lines in a filament will tend to be dragged radially inward (i.e., pinched) toward its central axis due to the filament's self-gravity. We study the cross-sectional profiles of the total intensity ($I$) and polarized intensity ($PI$) of dust emission in four segments of filaments unaffected by local star formation, which are expected to retain a pristine magnetic field structure. We find that the filaments' full-width at half-maximum (FWHM) in $PI$ are not the same as those in $I$, with two segments being appreciably narrower in $PI$ (FWHM ratio $\simeq 0.7-0.8$) and one segment being wider (FWHM ratio $\simeq 1.3$). The filament profiles of polarization fraction ($P$) do not show a minimum at the spine of the filament, which is not in line with an anticorrelation between $P$ and $I$ normally seen in molecular clouds and protostellar cores. Dust grain alignment variation with density cannot reproduce the observed $P$ distribution. We demonstrate numerically that the $I$ and $PI$ cross-sectional profiles of filaments in magnetohydrostatic equilibrium will have differing relative widths depending on the viewing angle. The observed variations of FWHM ratios in NGC 1333 are therefore consistent with models of pinched magnetic field structures inside filaments, and especially if they are magnetically near-critical or supercritical.

One widely discussed mechanism to produce highly coherent radio emission of fast radio bursts (FRBs) is coherent emission by bunches, either via curvature radiation or inverse Compton scattering (ICS). It has been suggested that the plasma oscillation effect can significantly suppress coherent emission by bunches, so that the bunching mechanism may not be the dominant mechanism for FRBs. We critically examine two physical conditions for significant plasma suppression and argue that the suppression effect is not important for the coherent ICS mechanism. A suppression factor for curvature radiation $f_{\rm cur}$ is possible if the giant bunch charge is surrounded by a high density plasma so that $\omega_p \gg \omega$ is satisfied, and if the electric field strength along the magnetic field lines ($E_\parallel$) in the emission region is not large enough to separate the surrounding plasma. However, even under these conditions, the derived $f_{\rm cur}$ is about 6 orders of magnitude greater than that derived in previous work, so that the suppression effect is not as significant as predicted before. We conclude that bunched coherent curvature radiation is still a plausible mechanism to power FRB emission, even though a suppression factor of the order $f_{\rm cur} \sim 10^{-3}$ should be considered in future modeling.

In this paper, we present the analysis of the stellar system HIP 101227 to determine the actual number of components in the system, and their properties. We use dynamical modeling and complex spectrophotometric (involving atmospheric modeling) techniques with recent data, to determine the physical properties and orbital solution for the system, respectively, with better accuracy than past studies. Based on our analysis, we found that the system is more consistent with being a quadruple rather than a binary, or a triple system as that suggested by previous studies. The total mass of the system determined from our SED analysis is 3.42 $\pm$ 0.20 $M_{\odot}$, which are distributed almost equally between the four stars. The stars are found to be zero-age main sequence stars; i.e., at the last stage of pre-main sequence,with age less than 200 Myr and spectral types K0. All four stars have very similar physical characteristics, suggesting that fragmentation process is the most likely theory for the formation and evolution of the system.

A key ingredient in any numerical study of supernova explosions is the nuclear network routine that is coupled with the hydrodynamic simulation code. When these studies are performed in more than one dimension, the size of the network is severely limited by computational issues. In this work, we propose a nuclear network, net87, which is close to one hundred nuclei and could be appropriate to simulate supernova explosions in multidimensional studies. One relevant feature is that electron and positron captures on free protons and neutrons have been incorporated to the network. Such addition allows for a better track of both, the neutronized species and the gas pressure. A second important feature is that the reactions are implicitly coupled with the temperature, which enhances the stability in the nuclear statistical equilibrium (NSE) regime. Here we analyze the performance of net87 in light of both: the computational overhead of the algorithm and the outcome in terms of the released nuclear energy and produced yields in typical Type Ia Supernova conditions.

Most massive stars belong to multiple systems, yet the formation process leading to such high multiplicity remains insufficiently understood. To help constrain the different formation scenarios that exist, insights on the low-mass end of the companion mass function of such stars is crucial. However, this is a challenging endeavour as (sub-)solar mass companions at angular separations ($\rho$) below 1" (corresponding to 1000-3000 au in nearby young open clusters and OB associations) are difficult to detect due to the large brightness contrast with the central star. With the Carina High-contrast Imaging Project of massive Stars (CHIPS), we aim to obtain statistically significant constraints on the presence and properties of low-mass companions around massive stars at a previously unreachable observing window ($\Delta \mathrm{mag} \gtrsim 10$ at $\rho \lesssim$ 1"). In this second paper in the series, we focus on the Trumpler 14 cluster, which harbours some of the youngest and most massive O-type stars in the Milky Way. We obtained VLT-SPHERE observations of seven O-type objects in Trumpler 14 using the IRDIFS_EXT mode. These allow us to search for companions at separations larger than 0."15 (approx. 360 au) and down to magnitude contrast $>10 \mathrm{mag}$ in the near-infrared. We used angular and spectral differential imaging along with PSF fitting to detect sources and measure their flux relative to that of the central object. We detected 211 sources with near-infrared magnitude contrast in the range of 2 to 12. The closest companion, at only 0."26, is characterised as a 1.4M$_{\odot}$ stars with an age of 0.6Myr, in excellent agreement with previous age estimates for Tr 14. The mass function peaks at about 0.4M$_{\odot}$ and presents a dearth of stars in the 0.5 to 0.8M$_{\odot}$ mass range compared to previous estimates of the initial mass function in Tr 14.

The [CII] 158$\mu$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$\mu$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

The distribution of ultra-high-energy cosmic-ray arrival directions appears to be nearly isotropic except for a dipole moment of order $6 \times (E/10~\mathrm{EeV})$ per cent. Nonetheless, at the highest energies, as the number of possible candidate sources within the propagation horizon and the magnetic deflections both shrink, smaller-scale anisotropies might be expected to emerge. On the other hand, the flux suppression reduces the statistics available for searching for such anisotropies. In this work, we consider two different lists of candidate sources: a sample of nearby starburst galaxies and the 2MRS catalog tracing stellar mass within $250~\mathrm{Mpc}$. We combine surface-detector data collected at the Pierre Auger Observatory until 2020 and the Telescope Array until 2019, and use them to test models in which UHECRs comprise an isotropic background and a foreground originating from the candidate sources and randomly deflected by magnetic fields. The free parameters of these models are the energy threshold, the signal fraction, and the search angular scale. We find a correlation between the arrival directions of $11.8\%_{-3.1\%}^{+5.0\%}$ of cosmic rays detected with $E \ge 38~\mathrm{EeV}$ by Auger or with $E \gtrsim 49~\mathrm{EeV}$ by TA and the position of nearby starburst galaxies on a ${15.5^\circ}_{-3.2^\circ}^{+5.3^\circ}$ angular scale, with a $4.2\sigma$ post-trial significance, as well as a weaker correlation with the overall galaxy distribution.

We present a machine learning model to classify Active Galactic Nuclei (AGN) and galaxies (AGN-galaxy classifier) and a model to identify type 1 (optically unabsorbed) and type 2 (optically absorbed) AGN (type 1/2 classifier). We test tree-based algorithms, using training samples built from the X-ray Multi-Mirror Mission -Newton (XMM-Newton) catalogue and the Sloan Digital Sky Survey (SDSS), with labels derived from the SDSS survey. The performance was tested making use of simulations and of cross-validation techniques. With a set of features including spectroscopic redshifts and X-ray parameters connected to source properties (e.g. fluxes and extension), as well as features related to X-ray instrumental conditions, the precision and recall for AGN identification are 94 and 93 per cent, while the type 1/2 classifier has a precision of 74 per cent and a recall of 80 per cent for type 2 AGN. The performance obtained with photometric redshifts is very similar to that achieved with spectroscopic redshifts in both test cases, while there is a decrease in performance when excluding redshifts. Our machine learning model trained on X-ray features can accurately identify AGN in extragalactic surveys. The type 1/2 classifier has a valuable performance for type 2 AGN, but its ability to generalise without redshifts is hampered by the limited census of absorbed AGN at high redshift.

NenuFAR, the New Extension in Nancay Upgrading LOFAR, is currently in its early science phase. It is in this context that the Cosmic Filaments and Magnetism Pilot Survey is observing sources with the array as it is still under construction - with 57 (56 core, 1 distant) out of a total planned 102 (96 core, 6 distant) mini-arrays online at the time of observation - to get a first look at the low-frequency sky with NenuFAR. One of its targets is the Coma galaxy cluster: a well-known object, host of the prototype radio halo. It also hosts other features of scientific import, including a radio relic, along with a bridge of emission connecting it with the halo. It is thus a well-studied object. In this paper, we show the first confirmed NenuFAR detection of the radio halo and radio relic of the Coma cluster at 34.4 MHz, with associated intrinsic flux density estimates: we find an integrated flux value of 106.3 +- 3.5 Jy for the radio halo, and 102.0 +- 7.4 Jy for the radio relic. These are upper bound values, as they do not include point-source subtraction. We also give an explanation of the technical difficulties encountered in reducing the data, along with steps taken to resolve them. This will be helpful for other scientific projects which will aim to make use of standalone NenuFAR imaging observations in the future.

The upcoming photometric surveys, such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor unprecedented number of active galactic nuclei (AGN) in a decade long campaign. Motivated by the science goals of LSST, which includes the harnessing of broadband light curves of AGN for photometric reverberation mapping (PhotoRM), we implement the existing formalism to estimate the lagged response of the emission line flux to the continuum variability using only mutli-band photometric light curves. We test the PhotoRM method on a set of 19 artificial light curves simulated using a stochastic model based on the Damped Random Walk process. These light curves are sampled using different observing strategies, including the two proposed by the LSST, in order to compare the accuracy of time-lag retrieval based on different observing cadences. Additionally, we apply the same procedure for time-lag retrieval to the observed photometric light curves of NGC 4395, and compare our results to the existing literature.

We study the vertical perturbations in the galactic disc of the Milky Way-size high-resolution hydrodynamical cosmological simulation named GARROTXA. We detect phase spirals in the vertical projection $Z- V_{Z}$ of disc's stellar particles for the first time in this type of simulations. Qualitatively similar structures were detected in the recent Gaia data, and their origin is still under study. In our model the spiral-like structures in the phase space are present in a wide range of times and locations across the disc. By accounting for an evolving mix of stellar populations, we observe that, as seen in the data, the phase spirals are better observed in the range of younger-intermediate star particles. We measure the intensity of the spiral with a Fourier decomposition and find that these structures appear stronger near satellite pericenters. Current dynamical models of the phase spiral considering a single perturber required a mass at least of the order of 10$^{10}$ M$_\odot$, but all three of our satellites have masses of the order of $\sim$10$^8$ M$_\odot$. We suggest that there are other mechanisms at play which appear naturally in our model such as the physics of gas, collective effect of multiple perturbers, and a dynamically cold population that is continuously renovated by the star formation Complementing collisionless isolated N-body models with the use of fully-cosmological simulations with enough resolution can provide new insights into the nature/origin of the phase spiral.

The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients.

We present the discovery of 32 new double periodic variables (DPVs) located toward the Galactic bulge. We found these objects among the nearly half a million binary stars published by the Optical Gravitational Lensing Experiment project. With this discovery, we increase the number of known DPVs in the Milky Way by a factor of 2. The new set of DPVs contains 31 eclipsing binaries and one ellipsoidal variable star. The orbital periods cover the range from 1.6 to 26 days, while long periods are detected between 47 and 1144 days. Our analysis confirms a known correlation between orbital and long periods that is also observed in similar systems in the Magellanic Clouds.

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim.

The galactic cluster Westerlund 1 contains a rich population of evolved, massive stars, and a high binary fraction has been inferred from previous multiwavelength observations. We use multi-epoch spectroscopy of a large sample of early-type stars to identify new binaries and binary candidates in the cluster. VLT/FLAMES was used to obtain spectra of ~100 OB stars over a 14-month baseline in 2008 and 2009, supplemented with follow-up observations in 2011 and 2013, and we identify 20 new OB I--III binaries, a WN9h binary, and a WC9d binary, greatly increasing the number of directly confirmed binary systems in Westerlund 1, while 12 O9--9.5 Iab--III stars are identified as candidate binaries. The 173.9 day SB1 W1030 represents the first longer-period system identified in the cluster, while the determination of a 53.95 day period for W44/L makes it the first Wolf-Rayet binary in Westerlund 1 with a confirmed orbital period greater than ten days. Our results suggest the binary fraction in the OB population is at least 40%, and may be significantly higher. (Abridged)

Understanding the physics of inflaton condensate fragmentation in the early Universe is crucial as the existence of fragments in the form of non-topological solitons (oscillons or Q-balls) may potentially modify the evolution of the post-inflation Universe. Furthermore, such fragments may evolve into primordial black holes and form dark matter, or emit gravitational waves. Due to the non-perturbative and non-linear nature of the dynamics, most of the studies rely on numerical lattice simulations. Numerical simulations of condensate fragmentation are, however, challenging, and, without knowing where to look in the parameter space, they are likely to be time-consuming as well. In this paper, we provide generic analytical conditions for the perturbations of an inflaton condensate to undergo growth to non-linearity in the cases of both symmetric and asymmetric inflaton potentials. We apply the conditions to various inflation models and demonstrate that our results are in good agreement with explicit numerical simulations. Our analytical conditions are easy to use and may be utilised in order to quickly identify models that may undergo fragmentation and determine conditions under they do so, which can guide subsequent in-depth numerical analyses.

The dynamics in the photosphere is governed by the multi-scale turbulent convection termed as granulation and supergranulation. It is important to derive 3-dimensional velocity vectors to understand the nature of the turbulent convection. However, it is difficult to obtain the velocity component perpendicular to the line-of-sight, which corresponds to the horizontal velocity in disk center observations. We developed a convolutional neural network model with a multi-scale deep learning architecture. The method consists of multiple convolutional kernels with various sizes of the receptive fields, and it performs convolution for spatial and temporal axes. The network is trained with data from three different numerical simulations of turbulent convection, and we introduced a coherence spectrum to assess the horizontal velocity fields that were derived at each spatial scale. The multi-scale deep learning method successfully predicts the horizontal velocities for each convection simulation in terms of the global-correlation-coefficient, which is often used for evaluating the prediction accuracy of the methods. The coherence spectrum reveals the strong dependence of the correlation coefficients on the spatial scales. Although coherence spectra are higher than 0.9 for large-scale structures, they drastically decrease to less than 0.3 for small-scale structures wherein the global-correlation-coefficient indicates a high value of approximately 0.95. We determined that this decrease in the coherence spectrum occurs around the energy injection scales. The accuracy for the small-scale structures is not guaranteed solely by the global-correlation-coefficient. To improve the accuracy in small-scales, it is important to improve the loss function for enhancing the small-scale structures and to utilize other physical quantities related to the non-linear cascade of convective eddies as input data.

We propose a paradigm shift in the data-driven modeling of the instrumental response field of telescopes. By adding a differentiable optical forward model into the modeling framework, we change the data-driven modeling space from the pixels to the wavefront. This allows to transfer a great deal of complexity from the instrumental response into the forward model while being able to adapt to the observations, remaining data-driven. Our framework allows a way forward to building powerful models that are physically motivated, interpretable, and that do not require special calibration data. We show that for a simplified setting of a space telescope, this framework represents a real performance breakthrough compared to existing data-driven approaches with reconstruction errors decreasing 5 fold at observation resolution and more than 10 fold for a 3x super-resolution. We successfully model chromatic variations of the instrument's response only using noisy broad-band in-focus observations.

Software citation has accelerated in astrophysics in the past decade, resulting in the field now having multiple trackable ways to cite computational methods. Yet most software authors do not specify how they would like their code to be cited, while others specify a citation method that is not easily tracked (or tracked at all) by most indexers. Two metadata file formats, codemeta.json and CITATION.cff, developed in 2016 and 2017 respectively, are useful for specifying how software should be cited. In 2020, the Astrophysics Source Code Library (ASCL, ascl.net) undertook a year-long effort to generate and send these software metadata files, specific to each computational method, to code authors for editing and inclusion on their code sites. We wanted to answer the question, "Would sending these files to software authors increase adoption of one, the other, or both of these metadata files?" The answer in this case was no. Furthermore, only 41% of the 135 code sites examined for use of these files had citation information in any form available. The lack of such information creates an obstacle for article authors to provide credit to software creators, thus hindering citation of and recognition for computational contributions to research and the scientists who develop and maintain software.

We find that galaxy angular momenta determined from shapes of spiral galaxies correlate with the initial tidal field in the galaxy's neighborhood. Explicitly, we find that the galaxies tend to be preferentially face on (edge on) oriented when the major (minor) axis of the initial tidal field is aligned with the line of sight. Because tidal-torque theory predicts vanishing of this kind of signal, the observed correlation suggests non-linear origin. This correlation is statistically significant even after considering the number of observables investigated. Our result provides motivation for considering galaxy shape information in the efforts aiming to reconstruct initial conditions in the Universe.

We search for evidence of primordial chirality violation in the galaxy data from the Sloan Digital Sky Survey by comparing how strongly directions of galaxy angular momenta correlate with left and right helical components of a spin vector field constructed from the initial density perturbations. Within uncertainties, galaxy spins correlate with these two helical components identically, which is consistent with Universe without primordial chirality violation. Given current data, it is not yet possible to rule out maximal chiral violation, although the case of vanishing correlation with the right helical component is ruled out at about 3.8$\sigma$.

We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore, by comparing the infrared counterpart's spectral energy distribution to known BeXRBs, we confirm this classification and set a distance of approximately $10-15$ kpc for the source. The source's broadband X-ray spectrum ($1.5-50$ keV) is described by an absorbed power-law with photon index $\Gamma$\,$\sim$\,$0.5$ and cutoff energy at $\sim$\,$13$ keV.

The origin of radio emission in the majority of Active Galactic Nuclei (AGN) is still poorly understood. Various competing mechanisms are likely involved in the production of radio emission and precise diagnostic tools are needed to disentangle them, of which variability is among the most powerful. For the first time, we show evidence for significant radio variability at 5 GHz at milli-arcsecond scales on days to weeks time scales in the highly accreting and extremely radio-quiet (RQ) Narrow Line Seyfert 1 (NLSy1) Mrk110. The simultaneous Swift/XRT light curve indicates stronger soft than hard X-ray variability. The short-term radio variability suggests that the GHz emitting region has a size smaller than ~180 Schwarzschild radii. The high brightness temperature and the radio and X-ray variability rule out a star-formation and a disc wind origin. Synchrotron emission from a low-power jet and/or an outflowing corona is then favoured.

In the era of wide-field surveys like the Zwicky Transient Facility and the Rubin Observatory's Legacy Survey of Space and Time, sparse photometric measurements constitute an increasing percentage of asteroid observations, particularly for asteroids newly discovered in these large surveys. Follow-up observations to supplement these sparse data may be prohibitively expensive in many cases, so to overcome these sampling limitations, we introduce a flexible model based on Gaussian Processes to enable Bayesian parameter inference of asteroid time series data. This model is designed to be flexible and extensible, and can model multiple asteroid properties such as the rotation period, light curve amplitude, changing pulse profile, and magnitude changes due to the phase angle evolution at the same time. Here, we focus on the inference of rotation periods. Based on both simulated light curves and real observations from the Zwicky Transient Facility, we show that the new model reliably infers rotational periods from sparsely sampled light curves, and generally provides well-constrained posterior probability densities for the model parameters. We propose this framework as an intermediate method between fast, but very limited period detection algorithms and much more comprehensive, but computationally expensive shape modeling based on ray-tracing codes.

Investigations of the formation of young stellar objects (YSOs) and planets require the detailed analysis of individual sources as well as statistical analysis of a larger number of objects. The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) project provides such a unique opportunity by establishing a UV spectroscopic library of young high- and low-mass stars in the local universe. Here we analyse optical photometry of the three ULLYSES targets (TXOri, V505Ori, V510Ori) and other YSOs in the $\sigma$Ori cluster taken at the time of the HST observations to provide a reference for those spectra. We identify three populations of YSOs along the line of sight to $\sigma$Ori, separated in parallax and proper motion space. The ULLYSES targets show typical YSO behaviour with pronounced variability and mass accretion rates of the order of 10$^{-8}$M$_\odot$/yr. Optical colours do not agree with standard interstellar reddening and suggest a significant contribution of scattered light. They are also amongst the most variable and strongest accretors in the cluster. V505\,Ori shows variability with a seven day period, indicating an inner disk warp at the co-rotation radius. Uncovering the exact nature of the ULLYSES targets will require improved detailed modelling of the HST spectra in the context of the available photometry, including scattered light contributions as well as non-standard reddening.

Dynamically cold stellar streams are the relics left over from globular cluster dissolution. These relics offer a unique insight into a now fully disrupted population of ancient clusters in our Galaxy. Using a combination of Gaia eDR3 proper motions, optical and near-UV colours we select a sample of likely Red Giant Branch stars from the GD-1 stream for medium-low resolution spectroscopic follow-up. Based on radial velocity and metallicity, we are able to find 14 new members of GD-1, 5 of which are associated with the \emph{spur} and \emph{blob/cocoon} off-stream features. We measured C-abundances to probe for the abundance variations known to exist in globular clusters. These variations are expected to manifest in a subtle way in globular clusters with such low masses ($\sim 10^4 {\rm ~M _{\odot}}$) and metallicities (${\rm [Fe/H]}\sim-2.1 {\rm ~dex}$). We find that the C-abundances of the stars in our sample display a small but significant ($ 3\sigma$ level) spread. Furthermore, we find $\ge 4\sigma$ variations in Al- and Mg-abundances as well as a $\sim2\sigma$ O-abundance variation among the stars in our sample that have been observed by APOGEE. These abundance patterns match the ones found in Galactic globular clusters of similar metallicity. Our results suggest that GD-1 represents the first fully disrupted globular cluster where light-element abundance spreads have been found.

With a 100mx110m off-axis paraboloid dish, the Green Bank Telescope (GBT) is the largest fully steerable radio telescope on Earth. A major challenge facing large ground-based radio telescopes is achieving sufficient pointing accuracy for observing at high frequencies, up to 116 GHz in the case of the GBT. Accurate pointing requires the ability to blindly acquire source locations and perform ad hoc corrections determined by observing nearby calibrator sources in order to obtain a starting position accurate to within a small margin of error of the target's location. The required pointing accuracy is dependent upon the half-power beamwidth, and for the higher-frequency end of GBT observing, this means that pointing must be accurate to within a few arcseconds RMS. The GBT's off-axis design is advantageous in that it eliminates blockage of the dish and reduces sidelobe interference, and there is no evidence that the resulting asymmetric structure adversely affects pointing accuracy. However, factors such as gravitational flexure, thermal deformation, azimuth track tilt and irregularity, and small misalignments and offset errors within the telescope's structure cause pointing inaccuracies. A pointing model was developed for the GBT to correct for these effects. The model utilizes standard geometrical corrections along with metrology data from the GBT's structural temperature sensors and data from measurements of the track levels. In this paper we provide a summary of the GBT's pointing model and associated corrections, as well as a discussion of relevant metrology systems and an analysis of its current nighttime pointing accuracy.

We present the first statistical analysis of complexity changes affecting the magnetic structure of interplanetary coronal mass ejections (ICMEs), with the aim of answering the questions: How frequently do ICMEs undergo magnetic complexity changes during propagation? What are the causes of such changes? Do the in situ properties of ICMEs differ depending on whether they exhibit complexity changes? We consider multi-spacecraft observations of 31 ICMEs by MESSENGER, Venus Express, ACE, and STEREO between 2008 and 2014 while radially aligned. By analyzing their magnetic properties at the inner and outer spacecraft, we identify complexity changes which manifest as fundamental alterations or significant re-orientations of the ICME. Plasma and suprathermal electron data at 1 au, and simulations of the solar wind enable us to reconstruct the propagation scenario for each event, and to identify critical factors controlling their evolution. Results show that ~65% of ICMEs change their complexity between Mercury and 1 au and that interaction with multiple large-scale solar wind structures is the driver of these changes. Furthermore, 71% of ICMEs observed at large radial (>0.4 au) but small longitudinal (<15 degrees) separations exhibit complexity changes, indicating that propagation over large distances strongly affects ICMEs. Results also suggest ICMEs may be magnetically coherent over angular scales of at least 15 degrees, supporting earlier theoretical and observational estimates. This work presents statistical evidence that magnetic complexity changes are consequences of ICME interactions with large-scale solar wind structures, rather than intrinsic to ICME evolution, and that such changes are only partly identifiable from in situ measurements at 1 au.

The nearby V4046 Sgr spectroscopic binary hosts a gas-rich disc known for its wide cavity and dusty ring. We present high resolution ($\sim$20 mas or 1.4 au) ALMA observations of the 1.3mm continuum of V4046 Sgr which, combined with SPHERE--IRDIS polarised images and a well-sampled spectral energy distribution (SED), allow us to propose a physical model using radiative transfer (RT) predictions. The ALMA data reveal a thin ring at a radius of 13.15$\pm$0.42 au (Ring13), with a radial width of 2.46$\pm$0.56 au. Ring13 is surrounded by a $\sim$10 au-wide gap, and it is flanked by a mm-bright outer ring (Ring24) with a sharp inner edge at 24 au. Between 25 and $\sim$35 au the brightness of Ring24 is relatively flat and then breaks into a steep tail that reaches out to $\sim$60 au. In addition, central emission is detected close to the star which we interpret as a tight circumbinary ring made of dust grains with a lower size limit of 0.8 mm at 1.1 au. In order to reproduce the SED, the model also requires an inner ring at $\sim$5 au (Ring5) composed mainly of small dust grains, hiding under the IRDIS coronagraph, and surrounding the inner circumbinary disc. The surprisingly thin Ring13 is nonetheless roughly 10 times wider than its expected vertical extent. The strong near-far disc asymmetry at 1.65 $\mu$m points at a very forward-scattering phase function and requires grain radii of no less than 0.4 $\mu$m.

We present a tool for the comparison and validation of the integration packages suitable for Solar System dynamics. iCompare, written in Python, compares the ephemeris prediction accuracy of a suite of commonly-used integration packages (JPL/HORIZONS, OpenOrb, OrbFit at present). It integrates a set of test particles with orbits picked to explore both usual and unusual regions in Solar System phase space and compares the computed to reference ephemerides. The results are visualized in an intuitive dashboard. This allows for the assessment of integrator suitability as a function of population, as well as monitoring their performance from version to version (a capability needed for the Rubin Observatory's software pipeline construction efforts). We provide the code on GitHub with a readily runnable version in Binder (https://github.com/dirac-institute/iCompare).

We study two intriguing disintegrating exoplanets, Kepler-1520b and K2-22b, and attempt to constrain the size of the underlying objects. These two planets are being disintegrated by their host stars, spewing dust and debris pulled from their surface into tails that trail and precede the exoplanet in its orbit, making it difficult to discern the true nature of the object. We attempted to peer through the dust cloud to put a constraint on the maximum radii of these exoplanets. While previous studies have done this in the past by selecting shallow transit events, we attempt a new statistical approach to model the intrinsic astrophysical and photon noise distributions simultaneously. We assume that the lightcurve flux distribution is distributed as a convolution of a Gaussian photon noise component and a Raleigh astrophysical component. The Raleigh curve has a finite flux maximum, which we fit with a Hamiltonian Markov Chain. With these methods, a more accurate flux maximum may be estimated, producing a more accurate and better final value for the size of these exoplanets. To determine statistical significance, we used the python package PyMC3 to find the posterior distribution for our data with Gaussian, Rayleigh, and joint function curves and plotting it against our collected flux. After completing this analysis, we were unable to constrain the radii of the exoplanets, as the forward scattering by dust dominates over dust extinction. However, this does mean that we were able better able to constrain the astrophysical variability and its maximum with our analysis.

Formed in protoplanetary disks around young stars, giant planets can leave observational features such as spirals and gaps in their natal disks through planet-disk interactions. Although such features can indicate the existence of giant planets, protoplanetary disk signals can overwhelm the innate luminosity of planets. Therefore, in order to image planets that are embedded in disks, it is necessary to remove the contamination from the disks to reveal the planets possibly hiding within their natal environments. We observe and directly model the detected disk in the Keck/NIRC2 vortex coronagraph $L'$-band observations of the single-armed protoplanetary disk around HD 34282. Despite a non-detection of companions for HD 34282, this direct disk modeling improves planet detection sensitivity by up to a factor of 2 in flux ratio and ${\sim}10 M_{\rm Jupiter}$ in mass. This suggests that performing disk modeling can improve directly imaged planet detection limits in systems with visible scattered light disks, and can help to better constrain the occurrence rates of self-luminous planets in these systems.

Magnetic monopoles have a long history of theoretical predictions and experimental searches, carrying direct implications for fundamental concepts such as electric charge quantization. We analyze in detail for the first time magnetic monopole production from collisions of cosmic rays bombarding the atmosphere. This source of monopoles is independent of cosmology, has been active throughout Earth's history, and supplies an irreducible monopole flux for all terrestrial experiments. Using results for robust atmospheric collider flux of monopoles, we systematically establish direct comparisons of previous ambient monopole searches with monopole searches at particle colliders and set leading limits on magnetic monopole production in the $\sim 5-100$ TeV mass-range.

The first Global e-Competition on Astronomy and Astrophysics was held on-line in September-October 2020 as a replacement for the International Olympiad on Astronomy and Astrophysics, which was postponed due to the ongoing pandemic caused by COVID-19. Despite the short time available for organisation, the competition was run successfully with 325 students from over 42 countries participating with no major issues. The feedback from the participants was positive and reflects the ways in which such events can boost interest in astronomy and astronomy education. With on-line activities set to be more prevalent in the future, we present an overview of the competition process and some of the lessons learned in hindsight as a guide for other event organisers.

We analyze families of hybrid equations of state of cold QCD matter, which combine input from gauge/gravity duality and from various ab initio methods for nuclear matter at low density, and predict that all neutron stars are fully hadronic without quark matter cores. We focus on constraints from recent measurements by the NICER telescope on the radius and mass of the millisecond pulsar PSR J0740+6620. These results are found to be consistent with our approach: they set only mild constraints on the hybrid equations of state, and favor the most natural models which are relatively stiff at low density. Adding an upper bound on the maximal mass of neutron stars, as suggested by the analysis of the GW170817 neutron star merger event, tightens the constraints considerably. We discuss updated predictions on observables such as the transition density and latent heat of the nuclear to quark matter transition as well as the masses, radii, and tidal deformabilities of neutron stars.

We introduce the Broadband Reflector Experiment for Axion Detection (BREAD) conceptual design and science program. This haloscope plans to search for bosonic dark matter across the [10$^{-3}$, 1] eV ([0.24, 240] THz) mass range. BREAD proposes a cylindrical metal barrel to convert dark matter into photons, which a novel parabolic reflector design focuses onto a photosensor. This unique geometry enables enclosure in standard cryostats and high-field solenoids, overcoming limitations of current dish antennas. A pilot 0.7 m$^{2}$ barrel experiment planned at Fermilab is projected to surpass existing dark photon coupling constraints by over a decade with one-day runtime. Axion sensitivity requires $<10^{-20}$ W/$\sqrt{\textrm{Hz}}$ sensor noise equivalent power with a 10 T solenoid and 10 m$^{2}$ barrel. We project BREAD sensitivity for various sensor technologies and discuss future prospects.

Flavor violating axion couplings can be in action before recombination, and they can fill the early universe with an additional radiation component. Working within a model-independent framework, we consider an effective field theory for the axion field and quantify axion production. Current cosmological data exclude already a fraction of the available parameter space, and the bounds will improve significantly with future CMB-S4 surveys. Remarkably, we find that future cosmological bounds will be comparable or even stronger than the ones obtained in our terrestrial laboratories.

The generalized hybrid metric-Palatini gravity is a theory of gravitation that has an action composed of a Lagrangian $f(R,\cal R)$, where $f$ is a function of the metric Ricci scalar $R$ and a new Ricci scalar $\cal R$ formed from a Palatini connection, plus a matter Lagrangian. This theory can be rewritten by trading the new geometric degrees of freedom of $f(R,\cal R)$ into two scalar fields, $\varphi$ and $\psi$, yielding an equivalent scalar-tensor theory. Given a spacetime theory, the next step is to find solutions. To construct solutions it is often necessary to know the junction conditions between two regions at a separation hypersurface $\Sigma$, with each region being an independent solution. The junction conditions for the generalized hybrid metric-Palatini gravity are found here, in the geometric and in the scalar-tensor representations, and in addition, for each representation, the junction conditions for a matching with a thin-shell and for a smooth matching at $\Sigma$ are worked out. These junction conditions are applied to three configurations, a star, a quasistar with a black hole, and a wormhole. The star has a Minkowski interior, a thin shell at the interface with all the energy conditions being satisfied, and a Schwarzschild exterior with mass $M$, and for this theory the matching can only be performed at the shell radius given by $r_\Sigma=\frac{9M}4$, the Buchdahl radius in general relativity. The quasistar with a black hole has an interior Schwarzschild black hole surrounded by a thick shell that matches smoothly to a mass $M$ Schwarzschild exterior at the light ring, and with the energy conditions being satisfied everywhere. The wormhole has an interior that contains the throat, a thin shell at the interface, and a Schwarzschild-AdS exterior with mass $M$ and negative cosmological constant $\Lambda$, with the null energy condition being obeyed.

We find an exact solution of the equations of motion of a two-field cosmological model, which realizes multi-field dark energy. The latter is characterized by field-space trajectories with turning rates that are always large. We study a class of two-field models and show that it is possible to have such trajectories, giving accelerated space-time expansion, even when the scalar potential preserves the rotational invariance of the field-space metric. For the case of Poincar\'e-disk field space, we derive the form of the scalar potential compatible with such background solutions and, furthermore, we find the exact solutions analytically. Their field-space trajectories are spirals inward, toward the center of the Poincar\'e disk. Interestingly, the functional form of the relevant scalar potential is compatible with a certain hidden symmetry, although the latter is broken by the presence of a constant term.

In both particle physics and cosmic ray muon applications, a high-resolution muon momentum measurement capability plays a significant role not only in providing valuable information on the properties of subatomic particles but also in improving the utilizability of muons. Currently, muon momentum is estimated by reconstructing the muon path using a strong magnetic field and muon trackers. Alternatively, time-of-flight or multiple Coulomb scattering techniques are less frequently applied, especially when there is a need to avoid using a magnetic field. However, the measurement resolution is much lower than that of magnetic spectrometers, approximately 20% in the muon momentum range of 0.5 to 4.5 GeV/c whereas it is nearly 10% or less when using magnets and trackers. Here, we propose a different paradigm to estimate muon momentum that utilizes multi-layer pressurized gas Cherenkov radiators. Using the fact that the gas refractive index varies with pressure and temperature, we can optimize the muon Cherenkov threshold momentum for which a muon signal will be detected. By analyzing the optical signals from Cherenkov radiation, we show that the actual muon momentum can be estimated with a minimum resolution of +-0.05 GeV/c for a large number of radiators over the range of 0.1 to 10.0 GeV/c. The results also show that our spectrometer correctly classifies the muon momentum (~87% classification rate) in the momentum range of 0.1 to 10.0 GeV/c. We anticipate our new spectrometer will to provide an alternative substitute for the bulky magnets without degrading measurement resolution. Furthermore, we expect it will significantly improve the quality of imaging or reduce the scanning time in cosmic muon applications by being incorporated with existing instruments.

We review the origins, motivations, and implications for cosmology and black holes, of our proposal that "dark energy" is not a quantum vacuum energy, but rather arises from a Weyl scaling invariant nonderivative component of the gravitational action.

Primordial black holes, which could have been formed in the very early Universe due to the collapse of large curvature fluctuations, are nowadays one of the most attractive and fascinating research areas in cosmology for their possible theoretical and observational implications. This review article presents the current results and developments on the conditions for primordial black hole formation from the collapse of curvature fluctuations in spherical symmetry on a Friedman-Lemaitre-Robertson-Walker background and its numerical simulation. We review the appropriate formalism for the conditions of primordial black hole formation, and we detail a numerical implementation. We then focus on different results regarding the threshold and the black hole mass using different sets of curvature fluctuations. Finally, we present the current state of analytical estimations for the primordial black hole formation threshold, contrasted with numerical simulations.