Papers for Wednesday, Jul 13 2022

We use very high resolution cosmological zoom simulations to follow the early evolution of twelve first-generation haloes formed from gaussian initial conditions with scale-free power spectra truncated on small scales by a gaussian. Initial collapse occurs with a diverse range of sheet- or filament-like caustic morphologies, but in almost all cases it gives rise to a numerically converged density cusp with $\rho = Ar^{-3/2}$ and total mass comparable to that of the corresponding peak in the initial linear density field. The constant $A$ can be estimated to within about 10 per cent from the properties of this peak. This outcome agrees with earlier work on the first haloes in cold and warm dark matter universes. Within central cusps, the velocity dispersion is close to isotropic, and equidensity surfaces tend to align with those of the main body of the halo at larger radii. As haloes grow, their cusps are often (but not always) overlaid with additional material at intermediate radii to produce profiles more similar to the Einasto or NFW forms typical of more massive haloes. Nevertheless, to the extent that we can resolve them, cusps survive at the smallest radii. Major mergers can disturb them, but the effect is quite weak in the cases that we study. The cusps extend down to the resolution limits of our simulations, which are typically a factor of several larger than the cores that would be produced by phase-space conservation if the initial power spectrum cutoff arises from free streaming.

Fuzzy dark matter (FDM) is a well motivated candidate for dark matter (DM) as its tiny mass and large de-Broglie wavelength suppress small-scale matter fluctuations, thereby solving some of the small-scale discrepancies in $\Lambda$CDM. Although it has been ruled out as the single component of DM by several observables, there is still a region in the FDM parameter space (the "FDM window", $10^{-25}\,\mathrm{eV}\lesssim m_\mathrm{FDM}\lesssim10^{-23}\,\mathrm{eV}$) where FDM is allowed to comprise a large portion of the total DM. In this work, for the first time, we study the signature of FDM (comprised of ultra-light axions) in fractions less than unity on the 21cm signal and its detectability by 21cm interferometers such as HERA, taking into account the degeneracy with both astrophysical and cosmological parameters, using a new pipeline that combines modified versions of the CAMB and 21cmFAST codes. Our forecasts imply that HERA in its design performance will be sensitive to FDM fractions as small as 1% in the FDM window, and improve over existing bounds for other masses by up to an order of magnitude.

The disks of Active Galactic Nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long Gamma-Ray Bursts (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN disks is key to their proper identification and to better constrain the star and compact object population in AGN disks. Some of these transients would appear unusual because they take place in much higher densities than the interstellar medium. Others, which are the subject of this paper, would additionally be modified by radiation diffusion since they are generated within optically thick regions of the accretion disks. Here we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black-hole masses and disk locations. We find that the radiation from radio to UV and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disk. In addition, photon diffusion can significantly delay the emergence of the emission peak, turning a beamed, fast transient into a slow, isotropic, and dimmer one. These would appear as broadband-correlated AGN variability with dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN disks as well as the disk structure.

Gaia DR3 provides >181, 000 radial velocity solutions for single-lined spectroscopic binaries (SB1s) which can be used to search for non-interacting compact object+star binary candidates by selecting systems with large mass functions. We selected 234 such systems and identified 114 systems with good RV solutions in DR3. We used light curves from ASAS-SN and TESS to identify and remove 31 eclipsing binaries to produce a catalog of 80 compact object+star candidates, including 38 ellipsoidal variables. The positions of these candidates on Gaia and 2MASS CMDs suggest that many of these systems are binaries with luminous companions. We compared the periods and eccentricities of detached eclipsing binaries in Gaia DR3 and ASAS-SN, and found that ~11% and ~60% of the binaries had different periods and eccentricities. We also compared RV solutions for 311 binaries in both Gaia DR3 and the Ninth Catalog of Spectroscopic Binary Orbits (SB9), and found similar results. We do not identify any strong candidates for non-interacting compact object+star binaries.

Non-thermal emission from relativistic Cosmic Ray (CR) electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure formation shocks. Emission caused by CR protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We introduce CRESCENDO: Cosmic Ray Evolution with SpeCtral Electrons aND prOtons, an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR (re-)acceleration at intra-cluster shocks, based on results of recent PIC simulations, adiabatic changes and radiative losses of electrons. We show its performance in test cases as well as idealized galaxy cluster (GC) simulations. We apply the model to an idealized GC merger following best-fit parameters for CIZA J2242.4+5301-1 and study CR injection, radio relic morphology, spectral steepening and synchrotron emission.

Measurements of the quantity of radio-mode feedback injected by an active galactic nucleus into the cluster environment have mostly relied on X-ray observations, which reveal cavities in the intracluster medium excavated by the radio lobes. However, the sensitivity required to accurately constrain the dimensions of these cavities has proven to be a major limiting factor and is the main bottleneck on high-redshift measurements. We describe a hybrid method based on a combination of X-ray and radio observations, which aims to enhance our ability to study radio-mode feedback. In this paper, we present one of the first samples of galaxy clusters observed with the International LOFAR Telescope (ILT) at 144 MHz and use this sample to test the hybrid method at lower frequencies than before. By comparing our measurements with results found in literature based on the traditional method using only X-ray observations, we find that the hybrid method provides consistent results to the traditional method. In addition, we find that the correlation between the traditional method and the hybrid method improves as the X-ray cavities are more clearly defined. This suggests that using radio lobes as proxies for cavities may help to circumvent systematic uncertainties in the cavity volume measurements. Encouraged by the high volume of unique ILT observations successfully processed, this hybrid method enables radio-mode feedback to be studied at high redshifts for the first time even for large samples of clusters.

The Local Volume Mapper (LVM) project is one of three surveys that form the Sloan Digital Sky Survey V. It will map the interstellar gas emission in a large fraction of the southern sky using wide-field integral field spectroscopy. Four 16-cm telescopes in siderostat configuration feed the integral field units (IFUs). A reliable acquisition and guiding (A&G) strategy will help ensure that we meet our science goals. Each of the telescopes hosts commercial CMOS cameras used for A&G. In this work, we present our validation of the camera performance. Our tests show that the cameras have a readout noise of around 5.6e- and a dark current of 21e-/s, when operated at the ideal gain setting and at an ambient temperature of 20{\deg}C. To ensure their performance at a high-altitude observing site, such as the Las Campanas Observatory, we studied the thermal behaviour of the cameras at different ambient pressures and with different passive cooling solutions. Using the measured properties, we calculated the brightness limit for guiding exposures. With a 5 s exposure time, we reach a depth of around 16.5 Gaia gmag with a signal-to-noise ratio (SNR)>5. Using Gaia Early Data Release 3, we verified that there are sufficient guide stars for each of the around 25000 survey pointings. For accurate acquisition, we also need to know the focal plane geometry. We present an approach that combines on-chip astrometry and using a point source microscope to measure the relative positions of the IFU lenslets and the individual CMOS pixels to around 2 $\mu$m accuracy.

How quasars powered by supermassive black holes (SMBHs) formed less than a billion years after the Big Bang is still one of the outstanding problems in astrophysics 20 years after their discovery$^{1-4}$. Cosmological simulations suggest that rare cold flows converging on primordial haloes in low-shear environments could have created these quasars if they were 10$^4$ - 10$^5$ M$_{\odot}$ at birth but could not resolve their formation$^{5-8}$. Semianalytical studies of the progenitor halo of a primordial quasar found that it favours the formation of such seeds but could not verify if one actually appeared$^9$. Here we show that a halo at the rare convergence of strong, cold accretion flows creates massive BH seeds without the need for UV backgrounds, supersonic streaming motions, or even atomic cooling. Cold flows drive violent, supersonic turbulence in the halo that prevents star formation until it reaches a mass that triggers sudden, catastrophic baryon collapse that forms 31,000 and 40,000 M$_{\odot}$ stars. This simple, robust process ensures that haloes capable of forming quasars by z $>$ 6 produce massive seeds. The first quasars were thus a natural consequence of structure formation in cold dark matter cosmologies, not exotic, finely-tuned environments as previously thought$^{10-14}$.

The large opacity fluctuations in the $z > 5.5$ Ly$\alpha$ forest may indicate inhomogeneous progress of reionization. To explain the observed large scatter of the effective Ly$\alpha$ optical depth ($\tau_{\rm eff}$) of the intergalactic medium (IGM), fluctuation of UV background ($\Gamma$ model) or the IGM gas temperature ($T$ model) have been proposed, which predict opposite correlations between $\tau_{\rm eff}$ and galaxy density. In order to address which model can explain the large scatter of $\tau_{\rm eff}$, we search for Ly$\alpha$ emitters (LAEs) around two (J1137+3549 and J1602+4228) quasar sightlines with $\tau_{\rm eff}\sim3$ and J1630+4012 sightline with $\tau_{\rm eff}\sim5.5$. Using a narrowband imaging with Subaru/Hyper Suprime-Cam, we draw LAE density maps to explore their spatial distributions. Overdensities are found within 20 $h^{-1}$Mpc of the quasar sightlines in the low $\tau_{\rm eff}$ regions, while a deficit of LAEs is found in the high $\tau_{\rm eff}$ region. Although the $\tau_{\rm eff}$ of the three quasar sightlines are neither high nor low enough to clearly distinguish the two models, these observed $\tau_{\rm eff}$-galaxy density relations all consistently support the $\Gamma$ model rather than the $T$ model in the three fields, along with the previous studies. The observed overdensities near the low $\tau_{\rm eff}$ sightlines may suggest that the relic temperature fluctuation does not affect reionization that much. Otherwise, these overdensities could be attributed to other factors besides the reionization process, such as the nature of LAEs as poor tracers of underlying large-scale structures.

Using dark matter only N-body cosmological simulations, we measure the pericentre distance of dark matter halos on their first infall into group and cluster halos. We find that the pericentre distance (R$_{\rm{peri}}$) is an important parameter as it significantly affects the strength of tidal mass loss in dense environments, and likely other environmental mechanisms as well. We examine what determines the R$_{\rm{peri}}$ value and find that, for most infallers, the dominant parameter is V$_{\rm{\perp}}$, the tangential component of the orbital velocity as the halo enters the group/cluster halo for the first time. This means that the strength of tidal stripping acting near the cores of groups/clusters are strongly influenced by the external peculiar velocity field of the large scale structure surrounding them, which differs between clusters, and is sensitive to the mass ratio of infaller to host. We find that filament feeding also partially contributes to feeding in low V$_{\rm{\perp}}$ halos. Dynamical friction can also play a role in reducing R$_{\rm{peri}}$ but this is only significant for those few relatively massive infallers ($>$10\% of the mass of their host). These results highlight how the response of galaxies to dense environments will sensitively depend on dynamics inherited from far outside those dense environments.

Theoretical models of protoplanetary disks including stellar irradiation often show a spontaneous amplification of scale height perturbations, produced by the enhanced absorption of starlight in enlarged regions. In turn, such regions cast shadows on adjacent zones that consequently cool down and shrink, eventually leading to an alternating pattern of overheated and shadowed regions. Previous investigations have proposed this to be a real self-sustained process, the so-called self-shadowing or thermal wave instability, which could naturally form frequently observed disk structures such as rings and gaps, and even potentially enhance the formation of planetesimals. All of these, however, have assumed in one way or another vertical hydrostatic equilibrium and instantaneous radiative diffusion throughout the disk. In this work we present the first study of the stability of accretion disks to self-shadowing that relaxes these assumptions, relying instead on radiation-hydrodynamical simulations. We first construct hydrostatic disk configurations by means of an iterative procedure and show that the formation of a pattern of enlarged and shadowed regions is a direct consequence of assuming instantaneous radiative diffusion. We then let these solutions evolve in time, which leads to a fast damping of the initial shadowing features in layers close to the disk surface. These thermally relaxed layers grow towards the midplane until all temperature extrema in the radial direction are erased in the entire disk. Our results suggest that radiative cooling and gas advection at the disk surface prevent a self-shadowing instability from forming, by damping temperature perturbations before these reach lower, optically thick regions.

We introduce a novel machine learning framework for estimating the Bayesian posteriors of morphological parameters for arbitrarily large numbers of galaxies. The Galaxy Morphology Posterior Estimation Network (GaMPEN) estimates values and uncertainties for a galaxy's bulge-to-total light ratio ($L_B/L_T$), effective radius ($R_e$), and flux ($F$). To estimate posteriors, GaMPEN uses the Monte Carlo Dropout technique and incorporates the full covariance matrix between the output parameters in its loss function. GaMPEN also uses a Spatial Transformer Network (STN) to automatically crop input galaxy frames to an optimal size before determining their morphology. This will allow it to be applied to new data without prior knowledge of galaxy size. Training and testing GaMPEN on galaxies simulated to match $z < 0.25$ galaxies in Hyper Suprime-Cam Wide $g$-band images, we demonstrate that GaMPEN achieves typical errors of $0.1$ in $L_B/L_T$, $0.17$ arcsec ($\sim 7\%$) in $R_e$, and $6.3\times10^4$ nJy ($\sim 1\%$) in $F$. GaMPEN's predicted uncertainties are well-calibrated and accurate ($<5\%$ deviation) -- for regions of the parameter space with high residuals, GaMPEN correctly predicts correspondingly large uncertainties. We also demonstrate that we can apply categorical labels (i.e., classifications such as "highly bulge-dominated") to predictions in regions with high residuals and verify that those labels are $\gtrsim 97\%$ accurate. To the best of our knowledge, GaMPEN is the first machine learning framework for determining joint posterior distributions of multiple morphological parameters and is also the first application of an STN to optical imaging in astronomy.

We present measurements of black hole masses and Eddington ratios for a sample of 38 bright (M$_{1450}$ < -24.4 mag) quasars at 5.8 < z < 7.5, derived from VLT/X-shooter near-IR spectroscopy of their broad CIV and MgII emission lines. The black hole masses (on average M$_{BH}$ ~ 4.6 x 10$^9$ M$_{\odot}$) and accretion rates (with Eddington ratios ranging between 0.1 and 1.0) are broadly consistent with that of similarly luminous 0.3 < z < 2.3 quasars, but there is evidence for a mild increase in the median Eddington ratio going towards z > 6. Combined with deep ALMA observations of the [CII] 158 $\mu$m line from the quasar host galaxies and VLT/MUSE investigations of the extended Ly$\alpha$ halos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations, z > 5.7 black holes appear to be over-massive with respect to their host galaxies, and their accretion properties do not change with host galaxy morphology. Under the assumption that the kinematics of the T ~ 10$^4$ K gas, traced by the extended Ly$\alpha$ halos, are dominated by the gravitational potential of the dark matter halo, we report a similar relation between the black hole mass and circular velocity to the one reported for z ~ 0 galaxies. These results paint a picture where the first supermassive black holes reside in massive halos at z > 6 and lead the first stages of galaxy formation by rapidly growing in mass with a duty cycle of order unity. However, this duty cycle needs to drastically drop towards lower redshifts, while the host galaxies continue forming stars at a rate of hundreds of solar masses per year, sustained by the large reservoirs of cool gas surrounding them.

Coronal mass ejections (CMEs) are a prominent contributor to solar system space weather and might have impacted the Sun's early angular momentum evolution. A signal diagnostic of CMEs on the Sun is coronal dimming: a drop in coronal emission, tied to the mass of the CME, that is the direct result of removing emitting plasma from the corona. We present the results of a coronal dimming analysis of Fe XII 1349 A and Fe XXI 1354 A emission from $\epsilon$ Eridani ($\epsilon$ Eri), a young K2 dwarf, with archival far-ultraviolet observations by the Hubble Space Telescope's Cosmic Origins Spectrograph. Following a flare in February 2015, $\epsilon$ Eri's Fe XXI emission declined by $81\pm5$%. Although enticing, a scant 3.8 min of preflare observations allows for the possibility that the Fe XXI decline was the decay of an earlier, unseen flare. Dimming nondetections following each of three prominent flares constrain the possible mass of ejected Fe XII-emitting (1 MK) plasma to less than a few $\times10^{15}$ g. This implies that CMEs ejecting this much or more 1 MK plasma occur less than a few times per day on $\epsilon$ Eri. On the Sun, $10^{15}$ g CMEs occur once every few days. For $\epsilon$ Eri, the mass loss rate due to CME-ejected 1 MK plasma could be $<0.6$ $\dot{M}_\odot$, well below the star's estimated 30 $\dot{M}_\odot$ mass loss rate (wind + CMEs). The order-of-magnitude formalism we developed for these mass estimates can be broadly applied to coronal dimming observations of any star.

We describe Keck-telescope spectrophotometry and imaging of the companion of the ``black widow" pulsar PSR~J0952$-$0607, the fastest known spinning neutron star (NS) in the disk of the Milky Way. The companion is very faint at minimum brightness, presenting observational challenges, but we have measured multicolor light curves and obtained radial velocities over the illuminated ``day" half of the orbit. The model fits indicate system inclination $i=59.8\pm 1.9^\circ$ and a pulsar mass $M_{NS} = 2.35\pm 0.17 M_\odot$, the largest well-measured mass found to date. Modeling uncertainties are small, since the heating is not extreme; the companion lies well within its Roche lobe and a simple direct-heating model provides the best fit. If the NS started at a typical pulsar birth mass, nearly $1 M_\odot$ has been accreted; this may be connected with the especially low intrinsic dipole surface field, estimated at $6\times 10^7$G. Joined with reanalysis of other black widow and redback pulsars, we find that the minimum value for the maximum NS mass is $M_{\rm max} > 2.19 M_\odot$$(2.09 M_\odot)$ at $1\sigma$$(3\sigma)$ confidence. This is $\sim 0.15 M_\odot$ heavier than the lower limit on $M_{\rm max}$ implied by the white-dwarf--pulsar binaries measured via radio Shapiro-delay techniques.

The SOSS mode of the NIRISS instrument is poised to be one of the workhorse modes for exoplanet atmosphere observations with the newly launched James Webb Space Telescope. One of the challenges of the SOSS mode, however, is the physical overlap of the first two diffraction orders of the G700XD grism on the detector. Recently, the ATOCA algorithm was developed and implemented as an option in the official JWST pipeline, as a method to extract SOSS spectra by decontaminating the detector -- that is, separating the first and second orders. Here, we present APPleSOSS (A Producer of ProfiLEs for SOSS), which generates the spatial profiles for each diffraction order upon which ATOCA relies. We validate APPleSOSS using simulated SOSS time series observations of WASP-52b, and compare it to ATOCA extractions using two other spatial profiles (a best and worst case scenario on-sky), as well as a simple box extraction performed without taking into account the order contamination. We demonstrate that APPleSOSS traces retain a high degree of fidelity to the true underlying spatial profiles, and therefore yield accurate extracted spectra. We further confirm that the effects of the order contamination for relative measurements (e.g., exoplanet transmission or emission observations) is small -- the transmission spectrum obtained from each of our four tests, including the contaminated box extraction, deviates by $\lesssim$0.1$\sigma$ from the atmosphere model input into our noiseless simulations. We further confirm via a retrieval analysis that the atmosphere parameters (metallicity and C/O) obtained from each transmission spectrum are consistent at the 1$\sigma$ level with the true underlying values.

The surfaces of neutron stars are likely sources of strongly polarized soft X rays due to the presence of strong magnetic fields. Scattering transport in the surface layers is critical to the determination of the emergent anisotropy of light intensity, and is strongly influenced by the complicated interplay between linear and circular polarization information. We have developed a magnetic Thomson scattering simulation to model the outer layers of fully-ionized atmospheres in such compact objects. Here we summarize emergent intensities and polarizations from extended atmospheric simulations, spanning considerable ranges of magnetic colatitudes. General relativistic propagation of light from the surface to infinity is fully included. The net polarization degrees are moderate and not very small when summing over a variety of field directions. These results provide an important foundation for observations of magnetars to be acquired by NASA's new IXPE X-ray polarimeter and future X-ray polarimetry missions.

It is argued that the new morphological and spectral information gleaned from the recently published LoFAR Two meter Sky Survey data release 2 (LoTSS-2 at 144 MHz) observations of the cluster Abell 980 (A980), in combination with its existing GMRT and VLA observations at higher frequencies, provide the much-needed evidence to strengthen the proposal that the cluster's radio emission comes mainly from two double radio sources, both produced by the brightest cluster galaxy (BCG) in two major episodes of jet activity. The two radio lobes left from the previous activity have become diffuse and developed an ultra-steep radio spectrum while rising buoyantly through the confining hot intra-cluster medium (ICM) and, concomitantly, the host galaxy has drifted to the cluster centre and entered a new active phase manifested by a coinciding younger double radio source. The new observational results and arguments presented here bolster the case that the old and young double radio sources in A980 conjointly represent a `double-double' radio galaxy whose two lobe-pairs have lost colinearity due to the (lateral) drift of their parent galaxy, making this system by far the most plausible case of a `Detached-Double-Double Radio Galaxy' (dDDRG).

We present robust, model-marginalized limits on both the total neutrino mass ($\sum m_\nu$) and abundance ($N_{\rm eff}$) to minimize the role of parameterizations, priors and models when extracting neutrino properties from cosmology. The cosmological observations we consider are CMB temperature fluctuation and polarization measurements, Supernovae Ia luminosity distances, BAO observations and determinations of the growth rate parameter from the Data Release 16 of the Sloan Digital Sky Survey IV. The degenerate neutrino mass spectrum (which implies $\sum m_\nu>0$) is weakly (moderately) preferred over the normal and inverted hierarchy possibilities, which imply the priors $\sum m_\nu>0.06$ and $\sum m_\nu>0.1$ eV respectively. Concerning the underlying cosmological model, the $\Lambda$CDM minimal scenario is almost always strongly preferred over the possible extensions explored here. The most constraining $95\%$ CL bound on the total neutrino mass in the $\Lambda$CDM+$\sum m_\nu$ picture is $\sum m_\nu< 0.087$ eV. The parameter $N_{\rm eff}$ is restricted to $3.08\pm 0.17$ ($68\%$ CL) in the $\Lambda$CDM+$N_{\rm eff}$ model. These limits barely change when considering the $\Lambda$CDM+$\sum m_\nu$+$N_{\rm eff}$ scenario. Given the robustness and the strong constraining power of the cosmological measurements employed here, the model-marginalized posteriors obtained considering a large spectra of non-minimal cosmologies are very close to the previous bounds, obtained within the $\Lambda$CDM framework in the degenerate neutrino mass spectrum. Future cosmological measurements may improve the current Bayesian evidence favouring the degenerate neutrino mass spectra, challenging therefore the consistency between cosmological neutrino mass bounds and oscillation neutrino measurements, and potentially suggesting a more complicated cosmological model and/or neutrino sector.

(Abridged) We present a cloud decomposition of $^{12}$CO (2--1) observations of the merger and nearest luminous infrared galaxy, NGC 3256. 185 spatially and spectrally resolved clouds are identified across the central $\approx$ 130 kpc$^{2}$ at 90 pc resolution and completeness is estimated. We compare our cloud catalogue from NGC 3256 to ten galaxies observed in the PHANGS-ALMA survey. Distributions in NGC 3256 of cloud velocity dispersions, luminosities, CO-estimated masses, mass surface densities, virial masses, virial parameters, size-linewidth coefficients, and internal turbulent pressures are significantly higher than in the PHANGS-ALMA galaxies. Cloud radii are slightly larger in NGC 3256 and free-fall times are shorter. The distribution of cloud eccentricities in NGC 3256 is indistinguishable from many PHANGS-ALMA galaxies, possibly because the dynamical state of clouds in NGC 3256 is similar to that of nearby spiral galaxies. However, the narrower distribution of virial parameters in NGC 3256 may reflect a narrower range of dynamical states than in PHANGS-ALMA galaxies. No clear picture of cloud alignment is detected, despite the large eccentricities. Correlations between cloud properties point to high external pressures in NGC 3256 keeping clouds bound and collapsing given such high velocity dispersions and star-formation rates. A fit to the cloud mass function gives a high-mass power-law slope of $-2.75^{+0.07}_{-0.01}$, near the average from PHANGS-ALMA galaxies. We also compare our results to a pixel-based analysis of these observations and find molecular-gas properties agree qualitatively, though peak brightness temperatures are somewhat higher and virial parameters and free-fall times are somewhat lower in this cloud-based analysis.

We present near-infrared spectra of three low-luminosity protostars and one background star in the Perseus molecular cloud, acquired using the Infrared Camera (IRC) onboard the \textit{AKARI} space telescope. For the comparison with different star-forming environments, we also present spectra of the massive protostar AFGL 7009S, where the protostellar envelope is heated significantly, and the low-mass protostar RNO 91, which is suspected to be undergoing an episodic burst. We detected ice absorption features of \ch{H2O}, \ch{CO2}, and \ch{CO} at all spectra around the wavelengths of 3.05, 4.27, and 4.67 $\mu$m, respectively. At least two low-luminosity protostars, we also detected the \ch{XCN} ice feature at 4.62 $\mu$m. The presence of the crystalline \ch{H2O} ice and \ch{XCN} ice components indicates that the low-luminosity protostars experienced a hot phase via accretion bursts during the past mass accretion process. We compared the ice abundances of low-luminosity protostars with those of the embedded low-mass protostars and the dense molecular clouds and cores, suggesting that their ice abundances reflect the strength of prior bursts and the timescale after the last burst.

We present results from three-dimensional hybrid-kinetic simulations of Alfv\'enic turbulence in a high-beta, collisionless plasma. The key feature of such turbulence is the interplay between local wave-wave interactions between the fluctuations in the cascade and the non-local wave-particle interactions associated with kinetic micro-instabilities driven by anisotropy in the thermal pressure (namely, firehose, mirror, and ion-cyclotron). We present theoretical estimates for, and calculate directly from the simulations, the effective collisionality and plasma viscosity in pressure-anisotropic high-beta turbulence, demonstrating that, for strong Alfv\'enic turbulence, the effective parallel viscous scale is comparable to the driving scale of the cascade. Most of the cascade energy (80-90%) is dissipated as ion heating through a combination of Landau damping and anisotropic viscous heating. The kinetic-energy spectrum of the turbulence has a slope steeper than $-5/3$ due to the anisotropic viscous stress. The magnetic-energy spectrum is shallower than $-5/3$ near the ion-Larmor scale due to fluctuations produced by the firehose instability. Our results have implications for models of particle heating in low-luminosity accretion onto supermassive black holes, the effective viscosity of the intracluster medium, and the interpretation of near-Earth solar-wind observations.

After a successful launch, the James Webb Space Telescope is preparing to undertake one of its principal missions, the characterization of the atmospheres of exoplanets. The Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph (NIRISS) is the only observing mode that has been specifically designed for this objective. It features a wide simultaneous spectral range (0.6--2.8\,\micron) through two spectral diffraction orders. However, due to mechanical constraints, these two orders overlap slightly over a short range, potentially introducing a ``contamination'' signal in the extracted spectrum. We show that for a typical box extraction, this contaminating signal amounts to 1\% or less over the 1.6--2.8\,\micron\ range (order 1), and up to 1\% over the 0.85--0.95\,\micron\ range (order 2). For observations of exoplanet atmospheres (transits, eclipses or phase curves) where only temporal variations in flux matter, the contamination signal typically biases the results by order of 1\% of the planetary atmosphere spectral features strength. To address this problem, we developed the Algorithm to Treat Order ContAmination (ATOCA). By constructing a linear model of each pixel on the detector, treating the underlying incident spectrum as a free variable, ATOCA is able to perform a simultaneous extraction of both orders. We show that, given appropriate estimates of the spatial trace profiles, the throughputs, the wavelength solutions, as well as the spectral resolution kernels for each order, it is possible to obtain an extracted spectrum accurate to within 10\,ppm over the full spectral range.

We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological inference using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensitivity of modular graph structure to the underlying cosmology in the noise-free limit, b) show that networks automatically combine mass and clustering information through comparisons to traditional statistics, c) demonstrate that graph neural networks can still extract information when catalogues are subject to noisy survey cuts, and d) illustrate how nonlinear IMNN summaries can be used as asymptotically optimal compressed statistics for Bayesian implicit likelihood inference. We reduce the area of joint $\Omega_m, \sigma_8$ parameter constraints with small ($\sim$100 object) halo catalogues by a factor of 42 over the two-point correlation function, and demonstrate that the networks automatically combine mass and clustering information. This work utilises a new IMNN implementation over graph data in Jax, which can take advantage of either numerical or auto-differentiability. We also show that graph IMNNs successfully compress simulations far from the fiducial model at which the network is fitted, indicating a promising alternative to $n$-point statistics in catalogue-based analyses.

The surfaces of airless bodies such as asteroids are exposed to many phenomena that can alter their physical properties. Bennu, the target of the OSIRIS-REx mission, has demonstrated how complex the surface of a small body can be. In 2019 November, the potentially hazardous asteroid 2015 JD1 experienced a close approach of 0.0331 au from the Earth. We present results of the physical characterization of 2015 JD1 based on ground-based radar, spectroscopy, and photometric observations acquired during 2019 November. Radar polarimetry measurements from the Arecibo Observatory indicate a morphologically complex surface. The delay-Doppler images reveal a contact-binary asteroid with an estimated visible extent of ~150 m. Our observations suggest that 2015 JD1 is an E-type asteroid with a surface composition similar to aubrites, a class of differentiated enstatite meteorites. The dynamical properties of 2015 JD1 suggest it came from the $\nu_6$ resonance with Jupiter, and spectral comparison with major E-type bodies suggest that it may have been derived from a parental body similar to the progenitor of the E-type (64) Angelina. Significantly, we find rotational spectral variation across the surface of 2015 JD1 from red to blue spectral slope. Our compositional analysis suggests that the spectral slope variation could be due to the lack of iron and sulfides in one area of the 2015 JD1 and/or differences in grain sizes.

Where the cosmic baryons lie in and around galactic dark matter halos is only weakly constrained. We develop a method to quickly paint on models for their distribution. Our approach uses the statistical advantages of $N$-body simulations, while painting on the profile of gas around individual halos in ways that can be motivated by semi-analytic models or zoom-in hydrodynamic simulations of galaxies. Possible applications of the algorithm include extragalactic dispersion measures to fast radio bursts (FRBs), the Sunyaev-Zeldovich effect, baryonic effects on weak lensing, and cosmic metal enrichment. As an initial application, we use this tool to investigate how the baryonic profile of foreground galactic-mass halos affects the statistics of the dispersion measure (DM) towards cosmological FRBs. We show that the distribution of DM is sensitive to the distribution of baryons in galactic halos, with viable gas profile models having significantly different probability distributions for DM to a given redshift. We also investigate the requirements to statistically measure the circumgalactic electron profile for FRB analyses that stack DM with impact parameter to foreground galaxies, quantifying the size of the contaminating "two-halo" term from correlated systems and the number of FRBs for a high significance detection. Publicly available Python modules implement our CGMBrush algorithm.

We present photoionization modeling of galaxy populations at $z\sim0$, $2$, and $> 6$ to bridge optical and far-infrared (FIR) emission-line diagrams. We collect galaxies with measurements of optical and/or FIR ([O III] 88 $\mathrm{\mu m}$ and [C II] 158 $\mathrm{\mu m}$) emission line fluxes and plot them on the [O III]$\lambda5007/\mathrm{H\beta}$--[N II]$\lambda6585/\mathrm{H\alpha}$ (BPT) and L([O III]88)/SFR--L([C II]158)/SFR diagrams, where SFR is the star-formation rate and L([O III]88) and L([C II]158) are the FIR line luminosities. We aim to explain the galaxy distributions on the two diagrams with photoionization models that employ three nebular parameters: the ionization parameter $U$, hydrogen density $n_\text{H}$, and gaseous metallicity $Z_\text{gas}$. Our models successfully reproduce the nebular parameters of local galaxies, and then predict the distributions of the $z\sim0$, $2$, and $> 6$ galaxies on the diagrams. The predicted distributions illustrate the redshift evolution on all the diagrams; e.g., [O III]$/\mathrm{H\beta}$ and [O III]88/[C II]158 ratios continuously decrease from $z > 6$ to $0$. Specifically, the $z > 6$ galaxies exhibit $\sim\!0.5$ dex higher $U$ than low-redshift galaxies at a given $Z_\text{gas}$ and show predicted flat distributions on the BPT diagram at $\log{\mathrm{[O III]/H\beta}} = 0.5$-$0.8$. We find that some of the $z > 6$ galaxies exhibit high L([O III]88)/SFR ratios. To explain these high ratios, our photoionization models require a low stellar-to-gaseous metallicity ratio or bursty/increasing star-formation history at $z > 6$. The James Webb Space Telescope will test the predictions and scenarios for the $z > 6$ galaxies proposed by our photoionization modeling.

Ram pressure stripping (RPS) is a process that removes the interstellar medium (ISM) quickly, playing a vital role in galaxy evolution. Previous RPS studies have treated the ISM as single-phase or lack the resolution and physical processes to properly capture the full multiphase ISM. To improve this simplification, we introduce an inflowing, hot intracluster medium (ICM) into a self-consistently modeled ISM in a local patch of star-forming galactic disks using the TIGRESS framework. Our simulations reveal that the workings of RPS are not only direct acceleration of the ISM by ICM ram pressure but also mixing-driven momentum transfer involving significant phase transition and radiative cooling. The hot ICM passes through the low-density channels of the porous, multiphase ISM, shreds the cool ISM, and creates mixing layers. The ICM momentum is transferred through the mixing layers while populating the intermediate temperature gas and radiating thermal energy away. The mixed gas extends beyond galactic disks and forms stripped tails that cool back unless the ICM fluxes are large enough to prevent cooling until they escape the simulation domain. The mixing-driven momentum transfer predicts that the more ICM mixes in, the faster the ISM moves, resulting in the anti-correlation of outflow velocity and gas metallicity of the stripped ISM. The compression of the ISM disks due to the ICM ram pressure enhances star formation rates up to 50% compared to the model without ICM. With the ICM ram pressure higher than the disk anchoring pressure, star formation is quenched within ~100 Myr.

Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $\lambda$ = 450 $\mu$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $\lambda = 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $\mu$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed.

As a typical data-driven method, deep learning becomes a natural choice for analysing astronomical data. In this study, we built a deep convolutional neural network to estimate basic stellar parameters $T\rm{_{eff}}$, log g, metallicity ([M/H] and [Fe/H]) and [$\alpha$/M] along with 10 individual elemental abundances ([C/Fe], [N/Fe], [O/Fe], [Mg/Fe], [Al/Fe], [Si/Fe], [Cl/Fe], [Ca/Fe], [Mn/Fe], [Ni/Fe]). The neural network is trained using common stars between the APOGEE survey and the LAMOST survey. We used low-resolution spectra from LAMOST survey as input, and measurements from APOGEE as labels. For stellar spectra with the signal-to-noise ratio in g band larger than 10 in the test set, the mean absolute error (MAE) is 29 K for $T\rm{_{eff}}$, 0.07 dex for log g, 0.03 dex for both [Fe/H] and [M/H], and 0.02 dex for [$\alpha$/M]. The MAE of most elements is between 0.02 dex and 0.04 dex. The trained neural network was applied to 1,210,145 giants, including sub-giants, from LAMOST DR8 within the range of stellar parameters 3500 K < $T\rm{_{eff}}$ < 5500 K, 0.0 dex < log g < 4.0 dex, -2.5 dex < [Fe/H] < 0.5 dex. The distribution of our results in the chemical spaces is highly consistent with APOGEE labels and stellar parameters show consistency with external high-resolution measurements from GALAH. The results in this study allows us to further studies based on LAMOST data and deepen our understanding of the accretion and evolution history of the Milky Way. The electronic version of the value added catalog is available at this http URL .

The ultracool M-dwarf star TRAPPIST-1 is surrounded by seven planets configured in a resonant chain. Transit-timing variations have shown that the planets are caught in multiple three-body resonances and that their orbits are slightly eccentric, probably caused by resonant forcing. The current values of the eccentricities could be a remnant from their formation. Here we run numerical simulations using fictitious forces of trapping the fully-grown planets in resonances as they migrated in the gas disc, followed by numerical simulations detailing their tidal evolution. For a reduced disc scale height $h\sim 0.03$--0.05, the eccentricities of the planets upon capture in resonance are higher than their current values by factors of a few. We show that the current eccentricities and spacing of planets d to h are natural outcomes of coupled tidal evolution wherein the planets simultaneously damp their eccentricities and separate due to their resonant interaction. We further show that the planets evolve along a set of equilibrium curves in semimajor axis--eccentricity phase space that are defined by the resonances, and that conserve angular momentum. As such, the current 8:5--5:3--(3:2)$^2$--4:3--3:2 resonant configuration cannot be reproduced from a primordial (3:2)$^4$--4:3--3:2 resonant configuration from tidal dissipation in the planets alone. We use our simulations to constrain the long-term tidal parameters $k_2/Q$ for planets b to e, which are in the range $10^{-3}$ to $10^{-2}$, and show that these are mostly consistent with those obtained from interior modelling following reasonable assumptions.

Nonrestricted hierarchical three-body configurations are common in various scales of astrophysical systems. Dynamical structures of the quadrupole-order resonance (the von Zeipel-Lidov-Kozai resonance) and the octupole-order resonance (the apsidal resonance) under the nonrestricted hierarchical planetary systems are investigated in this work by taking advantage of perturbative treatments. Under the quadrupole-order Hamiltonian model, the distribution of libration and circulation regions as well as the distribution of flipping region are analytically explored in the parameter space spanned by the conserved quantities. The fundamental frequencies of system are produced and then the nominal location of octupole-order resonance is identified. From the viewpoint of perturbative theory, the quadrupole-order Hamiltonian determines the unperturbed dynamical model and the octupole-order Hamiltonian plays an role of perturbation to the quadrupole-order dynamics. The resonant Hamiltonian for octupole-order resonances is formulated by means of averaging theory, giving rise to a new constant of motion. Phase portraits are produced to analyse dynamical structures of octupole-order resonance, including resonant centres, saddle points, dynamical separatrices and islands of libration. By analysing phase portraits, it is found that there are four branches of libration centre and eight libration zones in the considered space. Applications to orbit flips show that there are five flipping regions.

The recent Gaia Third Data Release contains a homogeneous analysis of millions of high-quality Radial Velocity Spectrometer (RVS) stellar spectra by the GSP-Spec module. This lead to the estimation of millions of individual chemical abundances and allows us to chemically map the Milky Way. Among the published GSP-Spec abundances, three heavy-elements produced by neutron-captures in stellar interiors can be found: Ce, Zr and Nd. We use a sample of about 30,000 LTE Ce abundances, selected after applying different combinations of GSP-Spec flags. Thanks to the Gaia DR3 astrometric data and radial velocities, we explore the cerium content in the Milky Way and, in particular, in its halo and disc components. The high quality of the Ce GSP-Spec abundances is quantified thanks to literature comparisons. We found a rather flat [Ce/Fe] versus [M/H] trend. We also found a flat radial gradient in the disc derived from field stars and, independently, from about 50 open clusters, in agreement with previous studies. The [Ce/Fe] vertical gradient has also been estimated. We also report an increasing [Ce/Ca] vs [Ca/H] in the disc, illustrating the late contribution of AGB with respect to SN II. Our cerium abundances in the disc, including the young massive population, are well reproduced by a new three-infall chemical evolution model. Among the halo population, the M 4 globular cluster is found to be enriched in cerium. Moreover, eleven stars with cerium abundances belonging to the Thamnos, Helmi Stream and Gaia-Sausage-Enceladus accreted systems were identified from chemo-dynamical diagnostics. We found that the Helmi Stream could be slightly underabundant in cerium, compared to the two other systems. This work illustrates the high quality of the GSP-Spec chemical abundances, that significantly contributes to unveil the heavy elements evolution history of the Milky Way.

Quantitative and analytical analysis of modulation process of the collimator is a great challenge, and is also of great value to the design and development of Fourier transform imaging telescopes. The Hard X-ray Imager (HXI), as one of the three payloads onboard the Advanced Space-based Solar Observatory(ASO-S) mission, adopts modulating Fourier-Transformation imaging technique and will be used to explore mechanism of energy release and transmission in solar flare activities. In this paper, a mathematical model is developed to analyze the modulation function under a simplified condition first. Then its behavior under six degrees of freedom is calculated after adding the rotation matrix and translation change to the model. In addition, unparalleled light and extended sources also are considered so that our model can be used to analyze the X-ray beam experiment. Next, applied to the practical HXI conditions, the model has been confirmed not only by Geant4 simulations but also by some verification experiments. Furthermore, how this model will help to improve the image reconstruction process after the launch of ASO-S is also presented.

We search for compact objects in binaries based on Gaia DR3. A sample of ten targets is derived under the conditions: low temperature ($T_{\rm eff} < 6000$ K), high radial velocity variation ($\Delta V_{\rm r} > 200$ km s$^{-1}$), high mass function ($f(M_2) > 1 M_\odot$), and ellipsoidal-like light curves. Two targets have LAMOST spectroscopic observations, one of which is a double-lined spectroscopic binary. The observational data of seven targets are not self-consistent, since their photometric periods are even shorter than the theoretical minimum orbital periods calculated by the stellar parameters from Gaia DR3. According to the gathered data, two targets may contain compact objects and are worth follow-up observations. This work may serve as an example to demonstrate the feasibility of searching for compact objects in the massive Gaia data.

The angular momentum of the Earth-Moon system was initially dominated by Earth's rotation with a short solar day of around 4 hours duration. Since then, Earth gradually transferred angular momentum through tidal friction to the orbit of the Moon, resulting in an increasing orbital radius and a deceleration of Earth's rotation. Geologic observations of tidal deposits can be used to verify and constrain models of lunar orbital evolution. In this work we reexamine the oldest tidal record suitable for analysis from the Moodies Group, South Africa with an age of 3.22 billion years. Time frequency analysis of the series of thicknesses of the sandstone layers yields a periodicity at 15.0 layers, taking into account the possibility of missing laminae. Assuming a mixed tidal system, the duration of two neap-spring-neap cycles was 30.0 lunar days for dominant semidiurnal or 30.0 sidereal days for dominant diurnal tides. We derive the relationship between this observation and the past Earth-Moon distance and re-visit related published work. We find that the Earth-Moon distance 3.22 billion years ago was around 70% of today's value. The Archean solar day was around 13 hours long with around 700 solar days per year. The ratio of solar to lunar tide-raising torque controls the leakage of angular momentum from the Earth-Moon system, but deviation from the assumed ratio of 0.211 results in only moderate changes. A duration of a postulated 21-hour atmospheric resonance shorter than 200 million years would be consistent with our observation; it would significantly alter Earth-Moon distance.

We present a near-UV space telescope on a ~70kg micro-satellite with a moderately fast repointing capability and a near real-time alert communication system that has been proposed in response to a call for an ambitious Czech national mission. The mission, which has recently been approved for Phase 0, A, and B1 study shall measure the brightness evolution of kilonovae, resulting from mergers of neutron stars in the near-UV band and thus it shall distinguish between different explosion scenarios. Between the observations of transient sources, the satellite shall perform observations of other targets of interest, a large part of which will be chosen in open competition.

The Bullet cluster (1E 0657-55.8) is a massive merging system at redshift $z$=0.296, known to host a powerful radio halo and a relic. Here we present high fidelity MeerKAT L-band (0.9-1.7 GHz) observations of the Bullet cluster, in which we trace a larger extent of both the radio halo and relic in comparison to previous studies. The size of the recovered radio halo is 1.6 Mpc $\times$ 1.3 Mpc and the largest linear size of the relic is ~988 kpc. We detect a new decrement feature on the southern outskirts of the halo emission, where a region appears to have a lower surface brightness in comparison to its surroundings. The fainter extension on the outskirts of the radio halo suggests a new wave of merger activity. An in-band spectral index map of the radio halo reveals a hint of radial steepening towards the edges of the diffuse source, likely due to synchrotron electron ageing. The integrated spectral index of the radio halo is 1.1$\pm$0.2. We perform a radio-X-ray surface brightness point-to-point analysis, which reveals a linear correlation for the radio halo. This indicates that the radio halo emission is produced by primary re-acceleration mechanisms. Finally, we derive a radio Mach number of M_R= 5.4$\pm$0.9 for the relic shock region, which is higher than the Mach number inferred by earlier analyses based on X-ray data. Discrepancies between radio and X-ray Mach numbers have been observed for multiple systems, with studies suggesting that this is due to various factors, including relic orientation.

The environment of S187, a nearby H II region (1.4$\pm$0.3 kpc), is analyzed. A surrounding shell has been studied in the H I line, molecular lines, and also in infrared and radio continua. We report the first evidence of a clumpy HI environment in its photodissociation region. A background radio galaxy enables the estimation of the properties of cold atomic gas. The estimated atomic mass fraction of the shell is $\sim$260~M$_{\odot}$, the median spin temperature is $\sim$50~K, the shell size is $\sim$4 pc with typical wall width around 0.2 pc. The atomic shell consists of $\sim$100 fragments. The fragment sizes correlate with mass with a power-law index of 2.39-2.50. The S187 shell has a complex kinematical structure, including the expanding quasi spherical layer, molecular envelope, an atomic sub-bubble inside the shell and two dense cores (S187~SE and S187~NE) at different stages of evolution. The atomic sub-bubble inside the shell is young, contains a Class II young stellar object and OH maser in the centre and the associated YSOs in the walls of the bubble. S187~SE and S187~NE have similar masses ($\sim$1200~M$_\odot$ and $\sim$900~M$_\odot$, respectively). S187~SE is embedded into the atomic shell and has a number of associated objects including high mass protostars, outflows, maser sources and other indicators of ongoing star formation. No YSOs inside S187~NE were detected, but indications of compression and heating by the H II region exist.

We present a new scheme to compensate for the small-scales approximations resulting from Particle-Mesh (PM) schemes for cosmological N-body simulations. This kind of simulations are fast and low computational cost realizations of the large scale structures, but lack resolution on small scales. To improve their accuracy, we introduce an additional effective force within the differential equations of the simulation, parameterized by a Fourier-space Neural Network acting on the PM-estimated gravitational potential. We compare the results for the matter power spectrum obtained to the ones obtained by the PGD scheme (Potential gradient descent scheme). We notice a similar improvement in term of power spectrum, but we find that our approach outperforms PGD for the cross-correlation coefficients, and is more robust to changes in simulation settings (different resolutions, different cosmologies).

We report the results of the photodynamical analyses of four compact, tight triple stellar systems, KICs 6964043, 5653126, 5731312, 8023317, based largely on $Kepler$ and $TESS$ data. All systems display remarkable eclipse timing and eclipse depth variations, the latter implying a non-aligned outer orbit. Moreover, KIC 6964043 is also a triply eclipsing system. We combined photometry, ETV curves, and archival spectral energy distribution data to obtain the astrophysical parameters of the constituent stars and the orbital elements with substantial precision. KICs 6964043 and 5653126 were found to be nearly flat with mutual inclinations $i_{mut}=4.1$ deg and $12.3$ deg, respectively, while KICs 5731312, 8023317 ($i_{mut}=39.4$ deg and $55.7$ deg, respectively) are found to lie in the high $i_{mut}$ regime of the von Zeipel-Kozai-Lidov (ZKL) theorem. We show that currently both high inclination triples exhibit observable unusual retrograde apsidal motion. Moreover, the eclipses will disappear in all but one of the four systems within a few decades. Short-term numerical integrations of the dynamical evolution reveal that both high inclination triples are currently subject to ongoing, large amplitude ($\Delta e\sim0.3$) inner eccentricity variations on centuries-long timescales, in accord with the ZKL theorem. Longer-term integrations predict that two of the four systems may become dynamically unstable on $\sim$ Gyr timescales, while in the other two triples common envelope phases and stellar mergers may occur. Finally we investigate the dynamical properties of a sample of 71 KIC/TIC triples statistically, and find that the mutual inclinations and outer mass ratios are anti-correlated at the 4$\sigma$ level. We discuss the implications for the formation mechanisms of compact triples.

The Crab pulsar and its nebula are among the most studied astrophysical systems, and constitute one of the most promising environments where high energy processes and particle acceleration can be investigated. They are the only objects for which previous X-ray polarisation has been reported. We present here the first Imaging X-ray Polarimetry Explorer (IXPE) observation of the Crab pulsar and nebula. The total pulsar pulsed emission in the [2-8] keV energy range is unpolarised. Significant polarisation up to 15% is detected only in the core of the main peak. The nebula has a total space integrated polarised degree of 20% and polarisation angle of 145deg. The polarised maps show a large level of asymmetry, and regions with polarised degree up to 45-50%. The polarisation pattern suggests a predominantly toroidal magnetic field.

A deep (~70%) and extended (~150 days) eclipse was seen towards the young multiple stellar system V773 Tau in 2010. We interpret it as due to the passage of a circumbinary disk around the B components moving in front of the A components. Our aim is to characterise the orientation and structure of the disk, to refine the orbits of the subcomponents, and to predict when the next eclipse will occur. We combine the photometry from several ground based surveys, construct a model for the light curve of the eclipse, and use high angular resolution imaging to refine the orbits of the three components of the system, A, B and C. Frequency analysis of the light curves, including from the TESS satellite, enables characterisation of the rotational periods of the Aa and Ab stars. A toy model of the circumbinary disk shows that it extends out to approximately 5 au around the B binary and has an inclination of 73 degrees with respect to the orbital plane of AB, where the lower bound of the radius of the disk is constrained by the geometry of the AB orbit and the upper bound is set by the stability of the disk. We identify several frequencies in the photometric data that we attribute to rotational modulation of the Aa and Ab stellar companions. We produce the first determination of the orbit of the more distant C component around the AB system and limit its inclination to 93 degrees. The high inclination and large diameter of the disk, together with the match from theory suggest that B is an almost equal mass, moderately eccentric binary. We identify the rotational periods of the Aa and Ab stars, identify a third frequency in the light curve that we attribute to the orbital period of the stars in the B binary. We predict that the next eclipse will be around 2037, during which both detailed photometric and spectroscopic monitoring will characterise the disk in greater detail.

Type Ia supernovae (SNe Ia) are used as standardisable candles to measure cosmological distances, but differences remain in their corrected luminosities which display a magnitude step as a function of host galaxy properties such as stellar mass and rest-frame $U-R$ colour. Identifying the cause of these steps is key to cosmological analyses and provides insight into SN physics. Here we investigate the effects of SN progenitor ages on their light curve properties using a galaxy-based forward model that we compare to the Dark Energy Survey 5-year SN Ia sample. We trace SN Ia progenitors through time and draw their light-curve width parameters from a bimodal distribution according to their age. We find that an intrinsic luminosity difference between SNe of different ages cannot explain the observed trend between step size and SN colour. The data split by stellar mass are better reproduced by following recent work implementing a step in total-to-selective dust extinction ratio $(R_V)$ between low- and high-mass hosts, although an additional intrinsic luminosity step is still required to explain the data split by host galaxy $U-R$. Modelling the $R_V$ step as a function of galaxy age provides a better match overall. Additional age vs. luminosity steps marginally improve the match to the data, although most of the step is absorbed by the width vs. luminosity coefficient $\alpha$. Furthermore, we find no evidence that $\alpha$ varies with SN age.

This manual is intended to help ALMA and other interferometer users improve images by recognising limitations and how to overcome them and deciding when and how to use self-calibration. The images provided by the ALMA Science Archive are calibrated using standard observing and data processing routines, including a quality assurance process to make sure that the observations meet the proposer's science requirements. This may not represent the full potential of the data, since any interferometry observation can be imaged with a range of resolutions and surface brightness sensitivity. The separation between phase calibration source and target usually limits the target dynamic range to a few hundred (or 50--100 for challenging conditions) but if the noise in the target field has not reached the thermal limit, improvements may be possible using self-calibration. This often requires judgements based on the target properties and is not yet automated for all situations. This manual provides background on the instrumental and atmospheric causes of visibility phase and amplitude errors, their effects on imaging and how to improve the signal to noise ratio and image fidelity by self-calibration. We introduce the conditions for self-calibration to be useful and how to estimate calibration parameter values for a range of observing modes (continuum, spectral line etc.). We also summarise more general error recognition and other techniques to tackle imaging problems. The examples are drawn from ALMA interferometric data processed using CASA, but the principles are generally applicable to most similar cm to sub-mm imaging.

The star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations.

We report on a particular example of noise and data representation interacting to introduce systematic error. Many instruments collect integer digitized values and appy nonlinear coding, in particular square-root coding, to compress the data for transfer or downlink; this can introduce surprising systematic errors when they are decoded for analysis. Square root coding and subsequent decoding typically introduces a variable, $\pm 1$ count value-dependent systematic bias in the data after reconstitution. This is significant when large numbers of measurements (e.g., image pixels) are averaged together. Using direct modeling of the probabiliity distribution of particular coded values in the presence of instrument noise, one may apply Bayes' Theorem to construct a decoding table that reduces this error source to a very small fraction of a digitizer step; in our example, systematic error from square root coding is reduced by a factor of 20 from 0.23 count RMS to 0.013 count RMS. The method is suitable both for new experiments such as the upcoming PUNCH mission, and also for post facto application to existing data sets -- even if the instrument noise properties are only loosely known. Further, the method does not depend on the specifics of the coding formula, and may be applied to other forms of nonlinear coding or representation of data values.

We study the bound orbit conditions for equatorial and eccentric orbits around a Kerr black hole both in the parameter space ($E$, $L$, $a$) representing the energy, angular momentum of the test particle, and spin of the black hole, and also ($e$, $\mu$, $a$) space representing the eccentricity, inverse-latus rectum of the orbit, and spin. We apply these conditions and implement the relativistic precession (RP) model to M82X-1, which is an Intermediate-mass black hole (IMBH) system, where two high-frequency Quasi-Periodic Oscillations (HFQPOs) and a low-frequency QPO were simultaneously observed. Assuming that the QPO frequencies can also be generated by equatorial and eccentric trajectories, we calculate the probability distributions to infer $e$, $a$, and periastron distance, $r_p$, of the orbit giving rise to simultaneous QPOs. We find that an eccentric orbit solution is possible in the region between innermost stable circular orbit (ISCO) and the marginally bound circular orbit (MBCO) for $e=0.2768^{+0.0657}_{-0.0451}$, $a=0.2897 \pm 0.0087$, and $r_p=4.6164^{+0.0694}_{-0.1259}$.

GRB 190829A is the fourth closest gamma-ray burst (GRB) to date ($z=0.0785$). Owing to its wide range of radio, optical, X-ray, and the very-high-energy (VHE) observations by H.E.S.S., it has become an essential new source examined by various models with complementary approaches. We here show in GRB 190829A the double-prompt pulses and the three-multiwavelength afterglows are consistent with the type II binary-driven hypernova (BdHN II) model. The progenitor is a binary composed of a carbon-oxygen (CO) star and a neutron star (NS) companion. The gravitational collapse of the iron core of the CO star produces a supernova (SN) explosion and leaves behind a new neutron star ($\nu$NS) at its center. The accretion of the SN ejecta onto the NS companion and onto the $\nu$NS via matter fallback spins up the NSs and produces the double-peak prompt emission. The synchrotron emission from the expanding SN ejecta with the energy injection from the rapidly spinning $\nu$NS and its subsequently spin-down leads to the afterglow in the radio, optical, and X-ray. We model the sequence of physical and related radiation processes in BdHNe and focus on individuating the binary properties that play the relevant roles.

This document characterizes the actual science performance of the James Webb Space Telescope (JWST), as known on 12 July 2022. Following six months of commissioning to prepare JWST for science operations, the observatory is now fully capable of achieving the discoveries for which it was built. A key part of commissioning activities was characterizing the on-orbit performance of the observatory. Here we summarize how the performance of the spacecraft, telescope, science instruments, and ground system differ from pre-launch expectations. Almost across the board, the science performance of JWST is better than expected. In most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies.

Simulation-Based Inference (SBI) is a promising Bayesian inference framework that alleviates the need for analytic likelihoods to estimate posterior distributions. Recent advances using neural density estimators in SBI algorithms have demonstrated the ability to achieve high-fidelity posteriors, at the expense of a large number of simulations ; which makes their application potentially very time-consuming when using complex physical simulations. In this work we focus on boosting the sample-efficiency of posterior density estimation using the gradients of the simulator. We present a new method to perform Neural Posterior Estimation (NPE) with a differentiable simulator. We demonstrate how gradient information helps constrain the shape of the posterior and improves sample-efficiency.

We present a new probabilistic method for detecting, deblending, and cataloging astronomical sources called the Bayesian Light Source Separator (BLISS). BLISS is based on deep generative models, which embed neural networks within a Bayesian model. For posterior inference, BLISS uses a new form of variational inference known as Forward Amortized Variational Inference. The BLISS inference routine is fast, requiring a single forward pass of the encoder networks on a GPU once the encoder networks are trained. BLISS can perform fully Bayesian inference on megapixel images in seconds, and produces highly accurate catalogs. BLISS is highly extensible, and has the potential to directly answer downstream scientific questions in addition to producing probabilistic catalogs.

We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster respectively than Nested Sampling (NS). We found that in higher dimensional applications the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

(abridged) Spectroscopic observations of high-redshift galaxies slowly reveal the same complexity of the interstellar medium (ISM) as expected from resolved observations in nearby galaxies. While providing a wealth of diagnostics, such high-z spectra are often unresolved, and inferring reliable diagnostics represents a major obstacle. Bright, nearby, unresolved galaxies provide many constraints to design methods to infer ISM properties, but they have so far been limited to deterministic methods and/or with simple topological assumptions. Multicomponent ISM models associate a complex enough distribution of matter and phases with a simple enough topological description to be constrained with probabilistic frameworks. Our goal is to provide a flexible statistical framework that is agnostic to the model grid and that considers either a few discrete components defined by their parameter values and/or statistical distributions of parameters. In this study we present a the MULTIGRIS framework and use it to model a suite of spectral lines in the star-forming regions of the metal-poor dwarf galaxy IZw18. We are able to calculate accurate values for the metallicity, number of ionizing photons, masses of ionized and neutral hydrogen, as well as the dust mass and the dust-to-gas mass ratio. We find a relatively modest amount of H$_2$ which is predominantly CO-dark and traced by C$^+$ rather than C$^0$. Nevertheless, more than $90$% of the [CII] emission is associated with the neutral atomic gas. Our models confirm the necessity to consider an X-ray source. Finally, we investigate the escape fraction of ionizing photons for different energy ranges. While the escape fraction for the main HII region lies around $50-65$%, we show that most of the soft X-ray photons are able to escape and may play a role in the ionization and heating of the circumgalactic or intergalactic medium.

pocoMC is a Python package for accelerated Bayesian inference in astronomy and cosmology. The code is designed to sample efficiently from posterior distributions with non-trivial geometry, including strong multimodality and non-linearity. To this end, pocoMC relies on the Preconditioned Monte Carlo algorithm which utilises a Normalising Flow in order to decorrelate the parameters of the posterior. It facilitates both tasks of parameter estimation and model comparison, focusing especially on computationally expensive applications. It allows fitting arbitrary models defined as a log-likelihood function and a log-prior probability density function in Python. Compared to popular alternatives (e.g. nested sampling) pocoMC can speed up the sampling procedure by orders of magnitude, cutting down the computational cost substantially. Finally, parallelisation to computing clusters manifests linear scaling.

Microwave telescopes require an ever-increasing control of experimental systematics in their quest to measure the Cosmic Microwave Background (CMB) to exquisite levels of precision. One important systematic for ground and balloon-borne experiments is ground pickup, where beam sidelobes detect the thermal emission of the much warmer ground while the main beam is scanning the sky. This generates scan-synchronous noise in experiment timestreams, which is difficult to filter out without also deleting some of the signal from the sky. Therefore, efficient modelling of pickup can help guide the design of experiments and of analysis pipelines. In this work, we present an extension to the beamconv algorithm that enables us to generate time-ordered data (TOD) from beam-convolved sky and ground maps simultaneously. We simulate ground pickup for both a ground-based experiment and a telescope attached to a stratospheric balloon. Ground templates for the balloon experiment are obtained by re-projecting satellite maps of the Earth's microwave emission.

We present a multi-wavelength analysis of two flare-related jets on November 13, 2014, using data from SDO/AIA, RHESSI, Hinode/XRT, and IRIS. Unlike most coronal jets where hard X-ray (HXR) emissions are usually observed near the jet base, in these events HXR emissions are found at several locations, including in the corona. We carry out the first differential emission measure (DEM) analysis that combines both AIA (and XRT when available) bandpass filter data and RHESSI HXR measurements for coronal jets, and obtain self-consistent results across a wide temperature range and into non-thermal energies. In both events, hot plasma first appeared at the jet base, but as the base plasma gradually cooled, hot plasma also appeared near the jet top. Moreover, non-thermal electrons, while only mildly energetic, are found in multiple HXR locations and contain a large amount of total energy. Particularly, the energetic electrons that produced the HXR sources at the jet top were accelerated near the top location, rather than traveling from a reconnection site at the jet base. This means that there was more than one particle acceleration site in each event. Jet velocities are consistent with previous studies, including upward and downward velocities around ~200 km/s and ~100 km/s respectively, and fast outflows of 400-700 km/s. We also examine the energy partition in the later event, and find that the non-thermal energy in accelerated electrons is most significant compared to other energy forms considered. We discuss the interpretations and provide constraints on mechanisms for coronal jet formation.

If the physics behind dark energy and/or dark matter violates the parity symmetry assumed in the standard cosmological paradigm, the linear polarization of the cosmic microwave background (CMB) photons can rotate due to their coupling to the dark sector. Recent 3$\sigma$ hints of this 'cosmic birefringence' in the EB spectrum of the CMB polarization motivates us to pursue new directions to independently validate and characterize the signal. Here, we explore the prospects to probe cosmic birefringence from small-scale fluctuations in the CMB using polarized Sunyaev Zel'dovich (pSZ) tomography. We find that pSZ can be used to infer the redshift dependence of cosmic birefringence and also help calibrate the instrumental polarization orientation. To illustrate the prospects, we show that pSZ tomography may probe an axion-like dark energy model with masses $m_\phi \lesssim 10^{-32}$eV with $\mathcal{O}(0.1)$ degrees of rotation between reionization and recombination.

A long-term (1995-2016) survey for novae in the nearby Andromeda galaxy (M31) was conducted as part of the Research-based Science Education initiative. During the course of the survey 180 nights of observation were completed at Kitt Peak, Arizona. A total of 262 novae were either discovered or confirmed, 40 of which have not been previously reported. Of these, 203 novae formed a spatially-complete sample detected by the KPNO/WIYN~0.9m telescope within in a $20'\times20'$ field centered on the nucleus of M31. An additional 50 novae are part of a spatially-complete sample detected by the KPNO~4m telescope within a larger $36'\times36'$ field. Consistent with previous studies, it was found that the spatial distribution of novae in both surveys follow the galaxy's bulge light somewhat more closely than the overall background light of the galaxy. After correcting for the limiting magnitude and the spatial and temporal coverage of the surveys, a final nova rate in M31 was found to be $R=40^{+7}_{-6}$~yr$^{-1}$, which is considerably lower than recent estimates. When normalized to the $K$-band luminosity of M31, this rate yields a luminosity-specific nova rate, $\nu_K$, of $3.3\pm0.4$~yr$^{-1}~[10^{10}~L_{\odot,K}]^{-1}$. By scaling the M31 nova rate using the relative infrared luminosities of M31 and our Galaxy, a nova rate of $R_\mathrm{G}=29\pm5$ is predicted for the Milky Way.

NGC 362 is a non-common Type II Galactic globular cluster, showing a complex pseudo two-color diagram or 'chromosome map'. The clear separation of its stellar populations in the color-magnitude diagram and the distribution of the giant stars in the chromosome map strongly suggests that NGC 362 could host stars with both cluster-nominal as well as enhanced heavy-element abundances, and one of them could be iron. However, despite previous spectroscopic observations of NGC 362, no such iron variation has been detected. Our main goal is to confirm or disprove this result by searching for any internal variation of [Fe/H] which would give us insight into the formation and evolution of this interesting globular cluster. In this paper, we present the abundance analysis for a sample of 11 red giant branch members based on high-resolution and high S/N spectra obtained with the MIKE echelle spectrograph mounted at the Magellan-Clay telescope. HST and GAIA photometry and astrometry has been used to determine atmospheric parameters and membership. We obtained T$_{\text{eff}}$, log(g) and v$_{\text{t}}$ for our target stars and measured the mean iron content of the sample and its dispersion with three different methods, which lead to [Fe/H]$_1$=-1.10$\pm0.02$, [Fe/H]$_2$=-1.09$\pm0.01$ and [Fe/H]$_3$=-1.10$\pm0.01$, while the internal dispersion turned out to be $\sigma_{[\text{Fe/H}]_1}$=0.06$\pm0.01$, $\sigma_{[\text{Fe/H}]_2}$=0.03$\pm0.01$ and $\sigma_{[\text{Fe/H}]_3}$=0.05$\pm0.01$ respectively. The error analysis gives an internal dispersion due to observational error of 0.05 dex. Comparing the observed dispersion with the internal errors, we conclude that NGC 362 does not show any trace of an internal iron spread.

In modern analysis pipelines, Einstein-Boltzmann Solvers (EBSs) are an invaluable tool for obtaining CMB and matter power spectra. To accelerate the computation of these observables, the CosmicNet strategy is to replace the bottleneck of an EBS, which is the integration of a system of differential equations for linear cosmological perturbations, by neural networks. This strategy offers advantages compared to the direct emulation of the final observables, including small networks that are easy to train in high-dimensional parameter spaces, and which do not depend by on primordial spectrum parameters nor observation-related quantities such as selection functions. In this second CosmicNet paper, we present a more efficient set of networks that are already trained for extended cosmologies beyond LCDM, with massive neutrinos, extra relativistic degrees of freedom, spatial curvature, and dynamical dark energy. We release a new branch of the CLASS code, called CLASSNET, which automatically uses networks within a region of trusted accuracy. We demonstrate the accuracy and performance of CLASSNET by presenting parameter inference runs from Planck, BAO and supernovae data, performed with CLASSNET and the COBAYA inference package. We have eliminated the perturbation module as a bottleneck of the EBS, with a speedup that is even more remarkable in extended cosmologies, where the usual approach would have been more expensive while the network's performance remains the same. We obtain a speedup factor of order 150 for the emulated perturbation module of CLASS. For the whole code, this translates into an overall speedup factor of order 3 when computing CMB harmonic spectra (now dominated by the highly parallelizable and further optimizable line-of-sight integration), and of order 50 when computing matter power spectra (less than 0.1 seconds even in extended cosmologies).

SDSS J1226+2152 is the first extragalactic science target to be observed by JWST. The strongly lensed z=2.9233 Lyman break galaxy, SGAS J122651.3+215220, lensed by a cluster at z=0.4358, was targeted as part of the JWST-ERS program TEMPLATES: Targeting Extremely Magnified Panchromatic Lensed Arcs and their Extended Star formation. Aiming to combine the exquisite capabilities of JWST with the extreme magnification provided by strong gravitational lensing, these observations will peer into galaxies at cosmic noon and probe the building blocks of star formation. Here, we present an Hubble Space Telescope strong lensing analysis, lens model, source-plane interpretation, and the lensing outputs needed to analyze the JWST observations in the context of the source's intrinsic properties. The lens model outputs are made publicly available to the community through the Mikulski Archive for Space Telescopes (MAST) Portal.

Aiming at investigating the collisional dynamics of massive stars in dense stellar systems hosting a black hole, two relativistic generalizations of the Chandrasekhar dynamical friction formula are derived and discussed. We extend the original formalism to the case of relativistic velocities and distributions, and we account for strong encounters with massive objects by including the first post-Newtonian correction (1PN), so that the effects of general relativity, at the lowest order, are accounted.

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass $M_{\rm GC} \sim 10^5\,M_\odot$ forms within $\approx 3\,{\rm Myr}$ when two GMCs with mass $M_{\rm GMC}=3.6\times10^5\,M_\odot$ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of $\sim 2$ when the radiation feedback is neglected, and the difference is slightly more pronounced (16\%) when extreme $\rm Ly\alpha$ feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $\beta=v/c>0.6$. We consider the radio to X-ray synchrotron emission produced by collisionless shocks driven by such fast ejecta into the interstellar medium. Analytic expressions are given for spherical ejecta with broken power-law mass (or energy) distributions, $M(>\gamma\beta)\propto(\gamma\beta)^{-s}$ with $s=s_{\rm KN}$ at $\gamma\beta<\gamma_0\beta_0$ and $s=s_{\rm ft}$ at $\gamma\beta>\gamma_0\beta_0$ (where $\gamma$ is the Lorentz factor). For parameter values characteristic of merger calculation results -- a "shallow" mass distribution, $1<s_{\rm KN}<3$, for the bulk of the ejecta (at $\gamma\beta\approx 0.2$), and a steep, $s_{\rm ft}>5$, "fast tail" mass distribution -- our model provides an accurate (to 10's of percent) description of the evolution of the flux, including at the phase of deceleration to sub-relativistic expansion. This is a significant improvement over earlier results, based on extrapolations of results valid for $\gamma\beta\gg1$ or $\ll1$ to $\gamma\beta\approx1$, which overestimate the flux by an order of magnitude for typical parameter values. It will enable a more reliable inference of ejecta parameters from future measurements of the non-thermal emission. For the merger event GW170817, the existence of a "fast tail" is expected to produce detectable radio and X-ray fluxes over a time scale of $\sim10^4$days.

The curvature scalar invariants of the Riemann tensor are important in General Relativity because they allow a manifestly coordinate invariant characterisation of certain geometrical properties of spacetimes such as, among others, curvature singularities, gravitomagnetism. We calculate explicit analytic expressions for the set of Zakhary-McIntosh curvature invariants for accelerating Kerr-Newman black holes in (anti-)de Sitter spacetime as well as for the Kerr-Newman-(anti-)de Sitter black hole. These black hole metrics belong to the most general type D solution of the Einstein-Maxwell equations with a cosmological constant. Explicit analytic expressions for the Euler-Poincare density invariant, which is relevant for the computation of the Euler-Poincare characteristic $\chi(M)$, and the Kretschmann scalar are also provided for both cases. We perform a detailed plotting of the curvature invariants that reveal a rich structure of the spacetime geometry surrounding the singularity of a rotating, electrically charged and accelerating black hole . These graphs also help us in an exact mathematical way to explore the interior of these black holes. Our explicit closed form expressions show that the above gravitational backgrounds possess a non-trivial Hirzebruch signature density. Possible physical applications of this property for the electromagnetic duality anomaly in curved spacetimes that can spoil helicity conservation are briefly discussed.

We point out that dark matter which is produced non-adiabatically in a phase transition (PT) with fast bubble walls receives a boost in velocity which leads to long free-streaming lengths. We find that this could be observed via the suppressed matter power spectrum for dark matter masses around $10^8 - 10^9$ GeV and energy scales of the PT around $10^{2} - 10^3~{\rm GeV}$. The PT should take place at the border of the supercooled regime, i.e. approximately when the Universe becomes vacuum dominated. This work offers novel physics goals for galaxy surveys, Lyman-$\alpha$, lensing, and 21-cm observations, and connects these to the gravitational waves from such phase transitions, and more speculatively to possible telescope signals of heavy dark matter decays.

At all Universe scales, there is a detectable amount of magnetic field. There are several probable origins for this observed magnetic field, including the possibility of its origin in the early Universe. There are several models for primordial magnetogenesis, and if the inflationary background is taken into account, breaking conformal symmetry is required to generate a sufficient amount of magnetic field. The conformal symmetry breaking is introduced either by new couplings between electromagnetic field and inflaton field or including higher derivative terms to the theory. To unify these different approaches in the literature, we propose an Effective Field Theory (EFT) approach, where EFT parameters describe the magnetogenesis scenario in the early Universe, and different choices of parameters correspond to different models. We show that the vector perturbations do not have temporal evolution; hence, only the gauge field is the relevant gauge-invariant variable for the EFT. We explicitly show that the generation of primordial magnetic fields requires two necessary conditions -- conformal invariance breaking and causal propagation. Hence, conformal invariance breaking is only a necessary condition, \emph{not} a sufficient condition. We confirm this by considering a specific model of primordial magnetogenesis.

Long-range interacting systems irreversibly relax as a result of their finite number of particles, $N$. At order $1/N$, this process is described by the inhomogeneous Balescu--Lenard equation. Yet, this equation exactly vanishes in one-dimensional inhomogeneous systems with a monotonic frequency profile and sustaining only 1:1 resonances. In the limit where collective effects can be neglected, we derive a closed and explicit $1/N^{2}$ collision operator for such systems. We detail its properties highlighting in particular how it satisfies an $H$-theorem for Boltzmann entropy. We also compare its predictions with direct $N$-body simulations. Finally, we exhibit a generic class of long-range interaction potentials for which this $1/N^{2}$ collision operator exactly vanishes.

Leptonic mixing patterns are usually extracted on the basis of groups or algebraic structures. In this paper, we introduce an alternative geometric method to study the correlations between the leptonic mixing parameters. At the 3{\sigma} level of the recent global fit data of neutrino oscillations, the distribution of the scatter points of the angles between the vectors of the magnitude of the leptonic mixing matrix is analysed. We find that the scatter points are concentrated on several special regions. Using the data in these regions, correlations of the leptonic mixing angles and the Dirac CP violating phase are obtained. The implications of the correlations are shown through the predicted flavor ratio of high-energy astrophysical neutrinos (HANs) at Earth.

We systematically investigate the microlensing effect of a charged spherically symmetric wormhole (Ellis wormhole as an example), where the light source is remote from the throat. Remarkably, there will be at most three images by considering the charge part. We study all situations including three images, two images, and one image, respectively. The numerical result shows that the range of total magnification is from $10^5$ to $10^{-2}$ depending on various metrics. However, we cannot distinguish the case that forms three images or only one image as the total magnification is of order $10^5$. Finally, our theoretical investigation could shed new light on exploring the wormhole with the microlensing effect.

Kinetic Inductance Detectors (KIDs) are superconducting energy-resolving detectors, sensitive to single photons from the near-infrared to ultraviolet. We study a hybrid KID design consisting of a beta phase tantalum ($\beta$-Ta) inductor and a NbTiN interdigitated capacitor (IDC). The devices show an average intrinsic quality factor $Q_i$ of 4.3$\times10^5$ $\pm$ 1.3 $\times10^5$. To increase the power captured by the light sensitive inductor, we 3D-print an array of 150$\times$150 $\mu$m resin micro lenses on the backside of the sapphire substrate. The shape deviation between design and printed lenses is smaller than 1$\mu$m, and the alignment accuracy of this process is $\delta_x = +5.8 \pm 0.5$ $\mu$m and $\delta_y = +8.3 \pm 3.3$ $\mu$m. We measure a resolving power for 1545-402 nm that is limited to 4.9 by saturation in the KID's phase response. We can model the saturation in the phase response with the evolution of the number of quasiparticles generated by a photon event. An alternative coordinate system that has a linear response raises the resolving power to 5.9 at 402 nm. We verify the measured resolving power with a two-line measurement using a laser source and a monochromator. We discuss several improvements that can be made to the devices on a route towards KID arrays with high resolving powers.

In this paper, a practical approach to the trajectory design for asteroid exploration missions with CubeSats is presented. When applied trajectories are sought, operative concerns and uncertainties affecting the spacecraft dynamics must be considered during the design process. Otherwise, trajectories that are possible on paper might become unfeasible when real-world constraints are considered. The risk of such eventualities leads to the urge of extending the trajectory design focus on the uncertainties affecting the dynamics and on the operative constraints derived by ground operations. This is especially true when targeting highly perturbed environments such as small bodies with low-cost solutions as CubeSats, whose capabilities in deep-space are still unknown. The case study presented is the Milani CubeSat which will be launched in 2024 with Hera in the frame of the AIDA mission.

We construct a model of cooperative superradiant emission from a highly relativistic multi-particle source. We revise an existing model of the literature for a relativistic two-level particle, and construct from it a Hamiltonian describing relativistic velocity dependent multi-particle superradiance. We adapt the standard diagrammatic framework to compute time evolution and density operators from our Hamiltonian, and demonstrate during the process a departure from standard results and calculation methods. In particular, we demonstrate that the so-called vertical photon result of the literature is modified by the relativistic Lorentz factor of the sample; we also introduce a set of coupled differential equations describing certain propagators in the velocity-dependent small sample framework, which we solve numerically via a hybrid fourth order Runge-Kutta and convolution approach. We demonstrate our methods for the simple case of two highly relativistic particles travelling with slightly differing velocities simulated at varying relativistic mean sample $\beta$ factors, and evaluate velocity coherence requirements for a sample to demonstrate enhanced superradiant emission in the observer frame. We find these coherence requirements to become increasingly restrictive at higher $\beta$ factors, even in the context of standard results of relativistic velocity differential transformations.

In this work we study the long-lived post-merger gravitational wave signature of a boson-star binary coalescence. We use full numerical relativity to simulate the post-merger and track the gravitational afterglow over an extended period of time. We implement recent innovations for the binary initial data, which significantly reduce spurious initial excitations of the scalar field profiles, as well as a measure for the angular momentum that allows us to track the total momentum of the spatial volume, including the curvature contribution. Crucially, we find the afterglow to last much longer than the spin-down timescale. This prolonged gravitational wave afterglow provides a characteristic signal that may distinguish it from other astrophysical sources.