Papers for Wednesday, Oct 04 2023

We propose that the ultralight dark matter (ULDM) model, in which dark matter particles have a tiny mass of $m=O(10^{-22})eV$, has characteristic scales for physical quantities of observed galaxies such as mass, size, acceleration, mass flux, and angular momentum from quantum mechanics. The typical angular momentum per dark matter particle is $\hbar$ and the typical physical quantities are functions of specific angular momentum $\hbar/m$ and average background density of the particles. If we use the Compton wavelength instead for the length scale, we can obtain bounds for these physical quantities. For example, there is an upper bound for acceleration of ULDM dominated objects, $a_c={c^3 m}/{\hbar}$. We suggest that the physical scales of galaxies depend on the time of their formation and that these characteristic scales are related to some mysteries of observed galaxies. Future observations from the James Webb Space Telescope and NANOgrav can provide evidences for the presence and evolution of these scales.

Due to be launched late 2026, the PLATO mission will bring the study of main-sequence solar-type and low-mass stars into a new era. In particular, PLATO will provide the community with a stellar sample with solar-type oscillations and activity-induced brightness modulation of unequalled size. We present here the main features of the analysis module that will be dedicated to measure stellar surface rotation and activity in the PLATO Stellar Analysis System.

A critical discourse is provided on the current status of the astrophysics of galaxies in view of open fundamental questions on the law of gravitation and the physics-driven variation of the galaxy-wide stellar initial mass function (GWIMF). The Einstein/Newtonian plus cold or warm dark-matter-based models face many significant but unresolved tensions (e.g. planes of satellites, the highly organised and symmetrical structure of the Local Group, the local Gpc-scale void). The accumulating nature of these indicates quite compellingly the need for a different theoretical framework. An example is the prediction, made in 1997, of the existence of satellite galaxies with near-exact properties to the Hercules dwarf spheroidal satellite galaxy discovered in 2007 if it is a tidal dwarf galaxy void of dark matter. Such results indicate that once the evidence is accepted that dark matter particles have no role in galaxies with gravitation being effectively Milgromian, then the resulting theoretical understanding of galaxies becomes much simpler and highly predictive with remarkable successes. Along with these scientific advancements, the recent observational data solidify the existence of one of the most important relations in star-formation theory, namely the mmax-Mecl relation. The data rule out stochastic star formation and confirm that the IMF becomes increasingly top-heavy with increasing star-formation rate density. Galaxies and their embedded-star-cluster building blocks are therefore self-regulated dynamical systems which are computationally accessible. Finding a theory for the formation of galaxies involves the development of a new cosmological model, which may differ substantially from the dark-matter based one. Very significant new opportunities have thus emerged for inquisitive and daring researchers which may appear risky now but almost certainly lead to major breakthroughs.

We report the discovery of Pavo, a faint ($M_V = -10.0$), star-forming, irregular, and extremely isolated dwarf galaxy at $D\approx2$ Mpc. Pavo was identified in Dark Energy Camera Legacy Survey imaging via a novel approach that combines low surface brightness galaxy search algorithms and machine learning candidate classifications. Follow-up imaging with the Inamori-Magellan Areal Camera & Spectrograph on the 6.5 m Magellan Baade telescope revealed a color--magnitude diagram (CMD) with an old stellar population, in addition to the young population that dominates the integrated light, and a tip-of-the-red-giant-branch distance estimate of $1.99^{+0.20}_{-0.22}$ Mpc. The blue population of stars in the CMD is consistent with the youngest stars having formed no later than 150 Myr ago. We also detected no H$\alpha$ emission with SOAR telescope imaging, suggesting we may be witnessing a temporary low in Pavo's star formation. We estimate the total stellar mass of Pavo to be $\log M_\ast/\mathrm{M_\odot} = 5.6 \pm 0.2$ and measure an upper limit on its HI gas mass of $1.0 \times 10^6\,\mathrm{M_\odot}$ based on the HIPASS survey. Given these properties, Pavo's closest analog is Leo P ($D=1.6$ Mpc), previously the only known isolated, star-forming, Local Volume dwarf galaxy in this mass range. However, Pavo appears to be even more isolated, with no other known galaxy residing within over 600 kpc. As surveys and search techniques continue to improve, we anticipate an entire population of analogous objects being detected just outside the Local Group.

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies.

We apply for the first time the Monte Carlo star cluster modeling method to study tidal tail and stellar stream formation from globular clusters, assuming a circular orbit in a smooth Galactic potential. Approximating energetically unbound bodies (potential escapers; PEs) as collisionless enables this fast but spherically symmetric method to capture asymmetric tidal phenomena with unprecedented detail. Beyond reproducing known stream features, including epicyclic overdensities, we show how 'returning tidal tails' may form after the stream fully circumnavigates the Galaxy back to the cluster, enhancing the stream's velocity dispersion. While a realistically clumpy, time-dependent Galactic potential may disrupt such tails, they warrant scrutiny as potentially excellent constraints on the Galactic potential's history and substructure. Re-examining the escape timescale $\Delta t$ of PEs, we find new behavior related to chaotic scattering in the three-body problem; the $\Delta t$ distribution features sharp plateaus corresponding to distinct locally smooth patches of the chaotic saddle separating the phase space basins of escape. We study for the first time $\Delta t$ in an evolving cluster, finding that $\Delta t\sim(E_{\rm J}^{-0.1},E_{\rm J}^{-0.4})$ for PEs with (low, high) Jacobi energy $E_{\rm J}$, flatter than for a static cluster ($E_{\rm J}^{-2}$). Accounting for cluster mass loss and internal evolution -- and (roughly) for ongoing relaxation among PEs -- lowers the median $\Delta t$ from ${\sim}10\,$Gyr to ${\lesssim}100\,$Myr. We finally outline future improvements to escape physics in the Monte Carlo method intended to enable both the first large-parameter-space studies of tidal tail/stellar stream formation from full globular cluster simulations and detailed comparisons to stream observations.

We conduct gas and dust hydrodynamical simulations of protoplanetary discs with one and two embedded planets to determine the impact that a second planet located further out in the disc has on the potential for subsequent planet formation in the region locally exterior to the inner planet. We show how the presence of a second planet has a strong influence on the collection of solid material near the inner planet, particularly when the outer planet is massive enough to generate a maximum in the disc's pressure profile. This effect in general acts to reduce the amount of material that can collect in a pressure bump generated by the inner planet. When viewing the inner pressure bump as a location for potential subsequent planet formation of a third planet, we therefore expect that the mass of such a planet will be smaller than it would be in the case without the outer planet, resulting in a small planet being sandwiched between its neighbours - this is in contrast to the expected trend of increasing planet mass with radial distance from the host star. We show that several planetary systems have been observed that do not show this trend but instead have a smaller planet sandwiched in between two more massive planets. We present the idea that such an architecture could be the result of the subsequent formation of a middle planet after its two neighbours formed at some earlier stage.

In this work, for the first time in literature, we compare the predictions of non-minimally coupled Natural and Coleman-Weinberg potentials in the $n_s-r$ plane against the constraints from the latest cosmological data in an extended $\Lambda$CDM model where we include non-standard self-interactions among massive neutrinos, mediated by a heavy scalar or vector boson. For the inflationary potentials, we consider two different formulations in gravity that are non-minimally coupled to the scalar field of the inflaton: \textit{Metric and Palatini.} We only consider the self-interaction to be present among $\tau$-neutrinos and only at moderate strengths. This is because strong interactions among $\tau$-neutrinos, or any strength self-interaction among electron- and muon-neutrinos, as well as any strength flavor-universal interactions, are strongly disfavoured from particle physics experiments. In terms of cosmological data, we use the latest public CMB datasets from Planck and BICEP/Keck collaborations, along with other data from CMB lensing, BAO, RSD, and SNe Ia luminosity distance measurements. We find that there are some situations where predictions from the inflationary models are ruled out at more than 2$\sigma$ by the minimal $\Lambda$CDM$+r$ model, but they are allowed in the self-interacting neutrino scenario.

We present CARMA, the Cluster Ages to Reconstruct the Milky Way Assembly project, that aims at determining precise and accurate age measurements for the entire system of known Galactic globular clusters and at using them to trace the most significant merger events experienced by the Milky Way. The strength of CARMA relies on the use of homogeneous photometry, theoretical isochrones, and statistical methods, that will enable to define a systematic-free chronological scale for the complete sample of Milky Way globulars. In this paper we describe the CARMA framework in detail, and present a first application on a sample of six metal-rich globular clusters with the aim of putting the final word on the debated origin of NGC6388 and NGC6441. Our results demonstrate that this pair of clusters is coeval with other four systems having a clear in-situ origin. Moreover, their location in the age-metallicity plane matches the one occupied by in-situ field stars. Such an accurate age comparison enabled by the CARMA methodology rules out the possibility that NGC6388 and NGC6441 have been accreted as part of a past merger event.

We use the SIMBA suite of cosmological hydrodynamical simulations to investigate the importance of various stellar and AGN feedback mechanisms in partitioning the cosmic baryons between the intergalactic (IGM) and circumgalactic (CGM) media in the $z\leq 1$ Universe. We identify the AGN jets as the most prominent mechanism for the redistribution of baryons between the CGM and large-scale IGM. In contrast to the full feedback models, deactivating AGN jets results in $\approx15-20$ per cent drop in fraction of baryons that reside in the IGM and a consequent increase of baryon fraction in the CGM by $\approx 50$ per cent. Additionally, we find that stellar feedback mechanisms modify the baryon partition between IGM and CGM on a $10$ per cent level. We further examine the physical properties of simulated haloes in different mass bins, and their response to various feedback models. On average, a sixfold decrease in the CGM mass fraction due to the inclusion of feedback from AGN-driven jets is detected in $10^{12}M_{\odot} \leq M_{\rm 200} \leq 10^{14}M_{\odot}$ haloes. Examination of the average radial gas density profiles of $M_{200} > 10^{12}M_{\odot}$ haloes reveals up to an order of magnitude decrease in gas densities due to the AGN jet feedback. We compare gas density profiles from SIMBA simulations to the predictions of the analytical modified NFW model, and show that the latter provides a reasonable approximation of the gas radial profiles within the virial radii of the full range of halo masses when rescaled by the CGM fraction of the halo. The relative partitioning of cosmic baryons and, subsequently, the feedback models can be constrained observationally with fast radio bursts (FRBs) in upcoming surveys.

We fit the UV/optical lightcurves of the Seyfert 1 galaxy Mrk 817 to produce maps of the accretion disk temperature fluctuations $\delta T$ resolved in time and radius. The $\delta T$ maps are dominated by coherent radial structures that move slowly ($v \ll c$) inwards and outwards, which conflicts with the idea that disk variability is driven only by reverberation. Instead, these slow-moving temperature fluctuations are likely due to variability intrinsic to the disk. We test how modifying the input lightcurves by smoothing and subtracting them changes the resulting $\delta T$ maps and find that most of the temperature fluctuations exist over relatively long timescales ($\sim$100s of days). We show how detrending AGN lightcurves can be used to separate the flux variations driven by the slow-moving temperature fluctuations from those driven by reverberation. We also simulate contamination of the continuum emission from the disk by continuum emission from the broad line region (BLR), which is expected to have spectral features localized in wavelength, such as the Balmer break contaminating the $U$ band. We find that a disk with a smooth temperature profile cannot produce a signal localized in wavelength and that any BLR contamination should appear as residuals in our model lightcurves. Given the observed residuals, we estimate that only $\sim$20% of the variable flux in the $U$ and $u$ lightcurves can be due to BLR contamination. Finally, we discus how these maps not only describe the data, but can make predictions about other aspects of AGN variability.

Current population synthesis modeling suggests that 30-50% of Type II supernovae originate from binary progenitors, however, the identification of a binary progenitor is challenging. One indicator of a binary progenitor is that the surrounding stellar population is too old to contain a massive single star.Measurements of the progenitor mass of SN 2017eaw are starkly divided between observations made temporally close to core-collapse which show a progenitor mass of 13-15 solar masses (final helium core mass of 4.4 to 6.0 solar masses - which is a more informative property than initial mass) and those from the stellar population surrounding the SN which find M<10.8 solar masses (helium core mass <3.4 solar masses). In this paper, we reanalyze the surrounding stellar population with improved astrometry and photometry, finding a median age of 16.8 (+3.2, -1.0) Myr for all stars younger than 50 Myr (helium core mass of 4.7 solar masses) and 85.9 (+3.2, -6.5) Myr for stars younger than 150 Myr. 16.8 Myr is now consistent with the helium core mass range derived from the temporally near explosion observations for single stars. Applying the combined constraints to population synthesis models, we determine that the probability of the progenitor of SN 2017eaw being an initially single-star is 65% compared to 35% for prior binary interaction. 85.9 Myr is inconsistent with any formation scenarios. We demonstrate that combining progenitor age constraints with helium core mass estimates from red supergiant SED modeling, late-time spectra, and indirectly from light curve modeling can help to differentiate single and binary progenitor scenarios and provide a framework for the application of this technique to future observations.

With the goal of untangling the origin of extended main-sequence turnoffs (eMSTOs) and extended red clumps (eRCs) in star clusters, in this work we present the study of the intermediate-age cluster NGC 419, situated along the Bridge of the Small Magellanic Cloud. To this aim, we analyzed multi-epoch, high angular resolution observations acquired with the Hubble Space Telescope for this dynamically young cluster, which enabled the determination of precise proper motions and therefore the assessment of the cluster membership for each individual star in the field of view. With this unprecedented information at hand, we first studied the radial distribution of kinematically selected member stars in different eMSTO subregions. The absence of segregation supports the rotation scenario as the cause for the turnoff color extension and disfavors the presence of a prolonged period of star formation in the cluster. A similar analysis on the eRC of NGC 419 confirms the absence of segregation, providing further evidence against an age spread, which is at odds with previous investigations. Even so, the currently available evolutionary models including stellar rotation fail at reproducing the two photometric features simultaneously. We argue that either shortcomings in these models or a different origin for the red clump feature, such as a nonstandard differential mass loss along the red giant branch phase, are the only way to reconcile our observational findings with theoretical expectations.

Light echoes give us a unique perspective on the nature of supernovae and non-terminal stellar explosions. Spectroscopy of light echoes can reveal details on the kinematics of the ejecta, probe asymmetry, and reveal details on its interaction with circumstellar matter, thus expanding our understanding of these transient events. However, the spectral features arise from a complex interplay between the source photons, the reflecting dust geometry, and the instrumental setup and observing conditions. In this work we present an improved method for modeling these effects in light echo spectra, one that relaxes the simplifying assumption of a light curve weighted sum, and instead estimates the true relative contribution of each phase. We discuss our logic, the gains we obtain over light echo analysis method(s) used in the past, and prospects for further improvements. Lastly, we show how the new method improves our analysis of echoes from Tycho's supernova (SN 1572) as an example.

Current models of galaxy formation require strong feedback from active galactic nuclei (AGN) to explain the observed lack of star formation in massive galaxies since z~2 but direct evidence of this energy input is limited. We use the SIMBA cosmological galaxy formation simulations to assess the ability of thermal Sunyaev-Zel'dovich (tSZ) measurements to provide such evidence, by mapping the pressure structure of the circumgalactic medium around massive z~0.2-1.5 galaxies. We undertake a stacking approach to calculate the total tSZ signal and its radial profile in simulations with varying assumptions of AGN feedback, and we assess its observability with current and future telescopes. By convolving our predictions with the 2.1' beam of the Atacama Cosmology Telescope (ACT), we show that current observations at z~1 are consistent with SIMBA's fiducial treatment of AGN feedback, and inconsistent with SIMBA models without feedback. At z~0.5, observational signals lie between SIMBA run with and without AGN feedback, suggesting AGN in SIMBA may inject too much energy at late times. By convolving our data with a 9.5'' beam corresponding to the TolTEC camera on the Large Millimeter Telescope Alfonso Serrano (LMT), we predict a unique profile for AGN feedback that can be distinguished with future higher-resolution measurements. Finally, we explore a novel approach to quantify the non-spherically symmetric features surrounding our galaxies by plotting radial profiles representing the component of the stack with m-fold symmetry.

The recently reported observation of VFTS 243 is the first example of a massive black-hole binary system with negligible binary interaction following black-hole formation. The black-hole mass ($\approx 10\ M_{\odot}$) and near-circular orbit ($e\approx 0.02$) of VFTS 243 suggest that the progenitor star experienced complete collapse, with energy-momentum being lost predominantly through neutrinos. VFTS 243 enables us to constrain the natal kick and neutrino-emission asymmetry during black-hole formation. At 68% C.L., the natal kick velocity (mass decrement) is $\lesssim 10$ km/s ($\lesssim 1.0\ M_{\odot}$). Most likely $\approx 0.3\ M_{\odot}$ were ejected, presumably in neutrinos, and the black hole experienced a natal kick of $4$ km/s. The neutrino-emission asymmetry is $\lesssim 4$%, with best fit values of $\sim$0-0.2%. Such a small neutrino natal kick accompanying black-hole formation is in agreement with theoretical predictions.

We report on X-ray (NICER/NuSTAR/MAXI/Swift) and radio (MeerKAT) timing and spectroscopic analysis from a three-month monitoring campaign in 2022 of a high-intensity outburst of the dipping neutron star low-mass X-ray binary 1A 1744-361. The 0.5-6.8 keV NICER X-ray hardness-intensity and color-color diagrams of the observations throughout the outburst suggests that 1A 1744-361 spent most of its outburst in an atoll state, but we show that the source exhibited Z-state-like properties at the peak of the outburst, similar to a small sample of other sources in the atoll state. A timing analysis with NICER data revealed several instances of an $\approx8$ Hz quasi-periodic oscillation (QPO; fractional rms amplitudes of ~5%) around the peak of the outburst, the first from this source, which we connect to the normal branch QPOs (NBOs) seen in the Z state. Our observations of 1A 1744-361 are fully consistent with the idea of the mass accretion rate being the main distinguishing parameter between atoll and Z states. Radio monitoring data by MeerKAT suggests that the source was at its radio-brightest during the outburst peak, and that the source transitioned from the 'island' spectral state to the 'banana' state within ~3 days of the outburst onset, launching transient jet ejecta. The observations present the strongest evidence for radio flaring, including jet ejecta, during the island-to-banana spectral state transition at low accretion rates (atoll state). The source also exhibited Fe XXV, Fe XXVI K$\alpha$, and K$\beta$ X-ray absorption lines, whose origins likely lie in an accretion disk atmosphere.

We present JWST MIRI 5.6, 10 and 21 micron observations of the candidate failed supernova N6946-BH1 along with HST WFPC/IR 1.1 and 1.6 micron data and ongoing optical monitoring data with the LBT. There is a very red, dusty source at the location of the candidate which has only ~10-15% of the luminosity of the progenitor star. The source is very faint in the HST near-IR observations (~10^3 Lsun) and is not optically variable to a limit of ~10^3 Lsun at R-band. The dust is likely silicate and probably has to be dominated by very large grains as predicted for dust formed in a failed SN. The required visual optical depths are modest, so we should begin to have a direct view of the source in the near-IR in only a few years.

Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale and necessitates the use of robust numerical simulations. We aim to model the birth and early evolution of a low-mass protostar. We also seek to describe the interior structure of the protostar and the radiative behavior of its accretion shock front. We carry out a high resolution numerical simulation of the collapse of a gravitationally unstable $1$ $\mathrm{M_{\odot}}$ dense molecular cloud core to stellar densities using three-dimensional radiation hydrodynamics under the gray flux-limited diffusion approximation. We follow the initial isothermal phase, the first adiabatic contraction, the second gravitational collapse triggered by the dissociation of $\mathrm{H}_{2}$ molecules, and $\approx 247$ days of the subsequent main accretion phase. We find that the sub-critical radiative behavior of the protostar's shock front causes it to swell as it accretes matter. We also find that the protostar is turbulent from the moment of its inception despite its radiative stability. This turbulence causes significant entropy mixing inside the protostar, which regulates the swelling. Furthermore, we find that the protostar is not fully ionized at birth, but the relative amount of ionized material within it increases as it accretes matter from its surroundings. Finally, we report in the results of the first 3D calculations involving a frequency-dependent treatment of radiative transfer, which has not produced any major differences with its gray counterpart.

The success of the Transiting Exoplanet Survey Satellite (TESS) mission has led to the discovery of an abundance of Venus Zone (VZ) terrestrial planets that orbit relatively bright host stars. Atmospheric observations of these planets play a crucial role in understanding the evolutionary history of terrestrial planets, past habitable states, and the divergence of Venus and Earth climates. The transmission spectrum of a Venus-like exoplanet can be difficult to distinguish from that of an Earth-like exoplanet however, which could severely limit what can be learned from studying exoVenuses. In this work we further investigate differences in transmission between hypothetical exoEarths and exoVenuses, both with varying amounts of atmospheric carbon dioxide (CO$_2$). The exoEarths and exoVenuses were modelled assuming they orbit TRAPPIST-1 on the runaway greenhouse boundary. We simulated James Webb Space Telescope (JWST) Near-Infrared Spectrograph (NIRSpec) PRISM transit observations of both sets of planets between 0.6-5.2 $\mu$m, and quantified the detectability of major absorption features in their transmission spectra. The exoEarth spectra include several large methane (CH$_4$) features that can be detected in as few as 6 transits. The CH$_4$ feature at 3.4 $\mu$m is the optimal for feature for discerning an exoEarth from an exoVenus since it is easily detectable and does not overlap with CO$_2$ features. The sulfur dioxide (SO$_2$) feature at 4.0 $\mu$m is the best indicator of an exoVenus, but it is detectable in atmospheres with reduced CO$_2$ abundance.

The superbubbles (SBs) carved in the interstellar medium by stellar winds and supernovae (SNe) are filled with hot ($T>$10$^{6}$~K) gas that produces soft X-ray emission (0.3-2.0\,keV). Models that assume a constant density medium and central SNe events fail to reproduce the soft X-ray luminosity that is observed in some SBs. We address this problem by generating models that trace the history of SNe in the SB and that produce off-centre SNe, that account for the missing soft X-ray emission. We test the models against archival, radio, optical, and X-ray observations of the SB DEM\,L50 located in the Large Magellanic Cloud. The soft X-ray properties of DEM L50, including its high luminosity, makes it a perfect candidate to test our models. Furthermore, the multiple wave-band observations of this object will help us to assess how well our models can reproduce other SB properties beside its soft X-ray properties. We find that a configuration where DEM L50 forms at the edge of a filament reproduces the observed soft X-ray luminosity, optical morphology, shell velocity, and swept-up mass of neutral gas. This configuration is supported by IR observations of the LMC. In addition, we find that off-centre SNe, which enhance soft X-ray emission, naturally occur for all of the initial ambient conditions we tested in our models. Finally, we show that an off-centre SN can explain the observed soft X-ray luminosity of DEM L50, and that the resulting luminosity is consistent with a plasma in non-equilibrium ionisation.

GW190521 is a remarkable gravitational-wave signal on multiple fronts: its source is the most massive black hole binary identified to date and could have spins misaligned with its orbit, leading to spin-induced precession -- an astrophysically consequential property linked to the binary's origin. However, due to its large mass, GW190521 was only observed during its final 3-4 cycles, making precession constraints puzzling and giving rise to alternative interpretations, such as eccentricity. Motivated by these complications, we trace the observational imprints of precession on GW190521 by dissecting the data with a novel time domain technique, allowing us to explore the morphology and interplay of the few observed cycles. We find that precession inference hinges on a quiet portion of the pre-merger data that is suppressed relative to the merger-ringdown. Neither pre-merger nor post-merger data alone are the sole driver of inference, but rather their combination: in the quasi-circular scenario, precession emerges as a mechanism to accommodate the lack of a stronger pre-merger signal in light of the observed post-merger. In terms of source dynamics, the pre-merger suppression arises from a tilting of the binary with respect to the observer. Establishing such a consistent picture between the source dynamics and the observed data is crucial for characterizing the growing number of massive binary observations and bolstering the robustness of ensuing astrophysical claims.

This study considers the characteristics of planetary systems with giant planets based on a population-level analysis of the California Legacy Survey planet catalog. We identified three characteristics common to hot Jupiters. First, while not all hot Jupiters have a detected outer giant planet companion ($M \sin i$ = 0.3--30 $M_{\textrm{Jup}}$), such companions are ubiquitous when survey completeness corrections are applied for orbital periods out to 40,000 days. Giant harboring systems without a hot Jupiter also host at least one outer giant planet companion per system. Second, the mass distributions of hot Jupiters and other giant planets are indistinguishable. However, within a planetary system that includes a hot Jupiter, the outer giant planet companions are at least $3\times$ more massive than the inner hot Jupiters. Third, the eccentricity distribution of the outer companions in hot Jupiter systems (with an average model eccentricity of $\langle e\rangle=0.34\pm0.05$) is different from the corresponding outer planets in planetary systems without hot Jupiters ($\langle e\rangle=0.19\pm0.02$). We conclude that the existence of two gas giants, where the outermost planet has an eccentricity $\ge0.2$ and is $3\times$ more massive, are key factors in the production of a hot Jupiter. Our simple model based on these factors predicts that $\sim$10\% of warm and cold Jupiter systems will by chance meet these assembly criteria, which is consistent with our measurement of $16\pm6\%$ relative occurrence of hot Jupiter systems to all giant-harboring systems. We find that these three features favor coplanar high-eccentricity migration as the dominant mechanism for hot Jupiter formation.

Gaia DR3 provided a first release of RP spectra and astrophysical parameters for ultracool dwarfs. We used these Gaia RP spectra and astrophysical parameters to select the most outlying ultracool dwarfs. These objects have spectral types of M7 or later and might be young brown dwarfs or low metallicity objects. This work aimed to find ultracool dwarfs which have Gaia RP spectra significantly different to the typical population. However, the intrinsic faintness of these ultracool dwarfs in Gaia means that their spectra were typically rather low signal-to-noise in Gaia DR3. This study is intended as a proof-of-concept for future iterations of the Gaia data releases. Based on well studied subdwarfs and young objects, we created a spectral type-specific color ratio, defined using Gaia RP spectra; this ratio is then used to determine which objects are outliers. We then used the objects kinematics and photometry external to Gaia to cut down the list of outliers into a list of 'prime candidates'. We produce a list of 58 Gaia RP spectra outliers, seven of which we deem as prime candidates. Of these, six are likely subdwarfs and one is a known young stellar object. Four of six subdwarf candidates were known as subdwarfs already. The two other subdwarf candidates: 2MASS J03405673+2633447 (sdM8.5) and 2MASS J01204397+6623543 (sdM9), are new classifications.

During cosmic recombination, charged particles bind into neutral atoms and the mean free path of photons rapidly increases, resulting in the familiar diffusion damping of primordial radiation temperature variations. An additional effect is a small photon spectrum distortion, because photons arriving from a particular sky direction were originally in thermal equilibrium at various spatial locations with different temperatures; the combination of these different blackbody temperature distributions results in a spectrum with a Compton $y$-distortion. Using the approximation that photons had zero mean free path prior to their second-to-last scattering, we derive an expression for the resulting $y$-distortion, and compute the angular correlation function of the diffusion $y$-distortion and its cross-correlation with the square of the photon temperature fluctuation. Detection of the cross-correlation is within reach of existing arcminute-resolution microwave background experiments such as the Atacama Cosmology Telescope and the South Pole Telescope.

We present a study of the mid-IR (MIR) emission of quiescent galaxies (QGs) beyond the local universe. Using deep $JWST$ imaging in the SMACS-0723 cluster field we identify a mass limited ($M_{*} >10^{9}$M$_{\odot}$) sample of intermediate redshift QGs ($0.2<z<0.7$) and perform modeling of their rest-frame UV to MIR photometry. We find that QGs exhibit a range of MIR spectra that are composed of a stellar continuum and a dust component that is 1-2 orders of magnitude fainter to that of star-forming galaxies. The observed scatter in the MIR spectra, especially at $\lambda_{\rm rest} > 5 \mu$m, can be attributed to different dust continuum levels and/or the presence of Polycyclic Aromatic Hydrocarbons (PAHs) features. The latter would indicate enhanced 11.3- and 12.7 $\mu$m PAHs strengths with respect to those at 6.2- and 7.7$ \mu$m, consistent with the observed spectra of local ellipticals and indicative of soft radiation fields. Finally, we augment the average UV-to-MIR spectrum of the population with cold dust and gas emission in the far-IR/mm and construct a panchromatic UV-to-radio SED that can serve as a template for the future exploration of the interstellar medium of $z>0$ QGs with ALMA and $JWST$.

The Green Bank Telescope (GBT) Diffuse Ionized Gas Survey (GDIGS) traces ionized gas in the Galactic midplane by observing radio recombination line (RRL) emission from 4-8 GHz. The nominal survey zone is $32.3^{\circ}> {\ell} > -5^{\circ}$, $|b|<0.5^{\circ}$. Here, we analyze GDIGS Hn${\alpha}$ ionized gas emission toward discrete sources. Using GDIGS data, we identify the velocity of 35 H II regions that have multiple detected RRL velocity components. We identify and characterize RRL emission from 88 H II regions that previously lacked measured ionized gas velocities. We also identify and characterize RRL emission from eight locations that appear to be previously-unidentified H II regions and 30 locations of RRL emission that do not appear to be H II regions based on their lack of mid-infrared emission. This latter group may be a compact component of the Galactic Diffuse Ionized Gas (DIG). There are an additional 10 discrete sources that have anomalously high RRL velocities for their locations in the Galactic plane. We compare these objects' RRL data to 13CO, H I and mid-infrared data, and find that these sources do not have the expected 24 ${\mu}$m emission characteristic of H II regions. Based on this comparison we do not think these objects are H II regions, but we are unable to classify them as a known type of object.

In this letter, we present a new method that quantitatively identifies the occurrence probability of equations of state (EoS) beyond "standard" EoS models that disfavor sharp and strong phase-transitions, based on neutron star mass and radius observations. The radii of two neutron stars with different masses are naturally correlated, in part because both of them are sensitive to the symmetry energy of the EoS. We show the radii of two neutron stars observed by NICER (PSR J0740+6620 and PSR 0030+0451) are correlated if these two neutron stars are built upon EoSs with no sharp and strong first-order phase transitions. We further show that the linear correlation of the neutron star radii can be significantly weakened, when strong and sharp first-order phase transitions occur. We propose a new quantity, ${D}_{\mathrm{L}}$, which measures the extent to which the linear correlation of the radii of two neutron stars is weakened. Our method gives a 48% identification probability (with a 5% false alarm rate) that the NICER observations indicate the necessity for a sharp and strong phase transition. Future observations can confirm or rule out this identification. Our method is generalizable to any pair of neutron star masses and can be employed with other sets of observations in the future.

Advancements in space exploration and sample return technology present a unique opportunity to leverage sample return capsules (SRCs) towards studying atmospheric entry of meteoroids and asteroids. Specifically engineered for the secure transport of valuable extraterrestrial samples from interplanetary space to Earth, SRCs offer unexpected benefits that reach beyond their intended purpose. As SRCs enter the Earth's atmosphere at hypervelocity, they are analogous to naturally occurring meteoroids and thus, for all intents and purposes, can be considered artificial meteors. Furthermore, SRCs are capable of generating shockwaves upon reaching the lower transitional flow regime, and thus can be detected by strategically positioned geophysical instrumentation. NASA's OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC is one of only a handful of artificial objects to re-enter the Earth's atmosphere from interplanetary space since the end of the Apollo era and it will provide an unprecedented observational opportunity. This review summarizes past infrasound and seismic observational studies of SRC re-entries since the end of the Apollo era and presents their utility towards the better characterization of meteoroid flight through the atmosphere.

We present polarization measurements in the $2-8 \thinspace \mathrm{keV}$ band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be $\Pi_{X} = 17.9 \pm 2.8 \%$ at an electric-vector position angle $\psi_X = 25.0 \pm 4.6^{\circ}$ using a spectropolarimetric fit from joint IXPE and XMM-Newton observations. There is no evidence for the polarization degree or angle varying significantly with energy or time on both short time scales (hours) or longer time scales (days). The contemporaneous polarization degree at optical wavelengths was $>$7$\times$ lower, making 1ES 0229+200 the most strongly chromatic blazar yet observed. This high X-ray polarization compared to the optical provides further support that X-ray emission in high-peaked blazars originates in shock-accelerated, energy-stratified electron populations, but is in tension with many recent modeling efforts attempting to reproduce the spectral energy distribution of 1ES 0229+200 which attribute the extremely high energy synchrotron and Compton peaks to Fermi acceleration in the vicinity of strongly turbulent magnetic fields.

It has long been discussed whether stellar feedback in the form of winds and/or radiation can shred the nascent molecular cloud, thereby controlling the star formation rate. However, directly probing and quantifying the impact of stellar feedback on the neutral gas of the nascent clouds is challenging. We present an investigation doing exactly that toward the RCW 79 HII region using the ionized carbon line at 158 $\mu$m ([CII]) from the FEEDBACK Legacy Survey. We combine this data with information on the dozen ionizing O stars responsible for the evolution of the region, and observe in [CII] for the first time both blue- and red-shifted mostly neutral high-velocity gas which reaches velocities up to 25 km s$^{-1}$ relative to the bulk emission of the molecular cloud. This high-velocity gas mostly contains neutral gas and partly forms a fragmented shell, similar to recently found shells in a few Galactic HII regions. However, this shell does not account for all of the observed neutral high-velocity gas. We also find high-velocity gas streaming out of the nascent cloud through holes and obtain a range of dynamical timescales below 1.0 Myr for the high-velocity gas which is well below the 2.3$\pm$0.5 Myr age of the OB cluster. This suggests a different scenario for the evolution of RCW 79, where the high-velocity gas is not solely stemming from a spherical expanding bubble, but also from gas recently ablated at the edge of the turbulent molecular cloud into the surrounding interstellar medium through low-pressure holes or chimneys. The resulting mass ejection rate estimate for the cloud is 0.9-3.5$\times$10$^{-2}$ M$_{\odot}$~yr$^{-1}$, which leads to short erosion timescales, i.e. $<$5 Myr, for the nascent molecular cloud. This finding provides direct observational evidence of rapid molecular cloud dispersal.

Chromospheric jets are plausible agents of energy and mass transport in the solar chromosphere, although their driving mechanisms have not yet been elucidated. Magnetic field measurements are key for distinguishing the driving mechanisms of chromospheric jets. We performed a full Stokes synthesis in the infrared range with a realistic radiative magnetohydrodynamics simulation that generated a chromospheric jet to predict spectro-polarimetric observations from the Sunrise Chromospheric Infrared spectro-Polarimeter (SCIP) onboard the SUNRISE III balloon telescope. The jet was launched by the collision between the transition region and an upflow driven by the ascending motion of the twisted magnetic field at the envelope of the flux tube. This motion is consistent with upwardly propagating non-linear Alfvenic waves. The upflow could be detected as continuous Doppler signals in the CaII 849.8 nm line at the envelope where the dark line core intensity and strong linear polarisation coexist. The axis of the flux tube was bright in both FeI 846.8 nm and CaII 849.8 nm lines with down-flowing plasma inside it. The structure, time evolution, and Stokes signals predicted in our study will improve the physical interpretation of future spectro-polarimetric observations with SUNRISE III/SCIP.

Photoionization models frequently assume constant temperature or density within HII regions. We investigate this assumption by measuring the detailed temperature and density structures of four HII regions in the Large Magellanic Cloud and the Small Magellanic Cloud, using integral-field spectroscopic data from the Wide-Field Spectrograph on the ANU 2.3m telescope. We analyse the distribution of emission-lines of low-ionization species, intermediate-ionization species and high-ionization species. We present the complex electron temperature and density structures within HII regions. All four nebulae present a negative gradient in the electron density profile. Both positive and negative temperature gradients are observed in the nebulae. We create a series of nebula models with a constant ISM pressure and varying temperature and density distributions. Comparison of the line ratios between our HII regions and models suggests that none of the simple nebula models can reproduce the observed temperature and density structures. Comparison between the models and the data suggests that the ISM pressure of nebulae in LMC and SMC is between log(P/k)=6-7.5. Complex internal structures of the nebulae highlight the importance of future Monte-Carlo photoionization codes for accurate nebula modeling, which include a comprehensive consideration of arbitrary geometries of HII regions.

The origin of the Galactic Center gamma-ray excess has not been conclusively determined after over a decade of careful study. The two most widely discussed possibilities are a population of millisecond pulsars (MSPs), and annihilation of dark matter particles. In contrast with annihilating dark matter, MSPs are expected to produce periodic emission. We show that even though the number of photons contributing to the excess is small, there is potentially sufficient information in the data from Fermi to detect a periodic MSP signal. Such a detection would definitively prove that at least some fraction of the excess is due to MSPs. We argue that this conclusion is robust to potential timing perturbations of the gamma-ray photons, such as those due to Earth's orbit, even if the number of parameters that must be used to model the perturbations is $\sim 7$.

Short period exoplanets on circular orbits are thought to be tidally locked into synchronous rotation. If tidally locked, these planets must possess permanent day- and nightsides, with extreme irradiation on the dayside and none on the nightside. However, so far the tidal locking hypothesis for exoplanets is supported by little to no empirical evidence. Previous work showed that the super-Earth LHS 3844b likely has no atmosphere, which makes it ideal for constraining the planet's rotation. Here we revisit the Spitzer phase curve of LHS 3844b with a thermal model of an atmosphere-less planet and analyze the impact of non-synchronous rotation, eccentricity, tidal dissipation, and surface composition. Based on the lack of observed strong tidal heating we rule out rapid non-synchronous rotation (including a Mercury-like 3:2 spin-orbit resonance) and constrain the planet's eccentricity to less than 0.001 (more circular than Io's orbit). In addition, LHS 3844b's phase curve implies that the planet either still experiences weak tidal heating via a small-but-nonzero eccentricity (requiring an undetected orbital companion), or that its surface has been darkened by space weathering; of these two scenarios we consider space weathering more likely. Our results thus support the hypothesis that short period rocky exoplanets are tidally locked, and further show that space weathering can significantly modify the surfaces of atmosphere-less exoplanets.

Ultra-High Energy (UHE, $>$0.1\,PeV) $\gamma$-ray Astronomy is rapidly evolving into an expanding branch of the $\gamma$-ray astronomy with the surprising discovery of 12 PeVatrons and the detection of a handful of photons above 1 PeV. Nearly all known celestial object types that have emissions in the TeV band are found also emitting UHE photons. UHE $\gamma$-rays have a well-defined horizon inside our galaxy due to the absorption of infrared and cosmic microwave backgrounds in the universe. With the last 30 years, traditional cosmic ray (CR) detection techniques allow the detection of UHE $\gamma$-rays, and opened up the last observation window. For leptonic sources, UHE radiation is in the deep Klein-Nishina regime which is largely suppressed. Therefore UHE $\gamma$-ray detection will help to locate and identify hadronic radiation sources, tracing the historic pursuit for the origin of CRs around the knee of the spectrum. The Crab Nebula is again the focus of attention with measured photon emissions above 1\,PeV. In the absence of hadronic processes, this may indicate the existence of an extreme accelerator of e$^+$/e$^-$. Utilization of the CR extensive air shower detection techniques broadens the field of view of the source observations, enabling the measurement of UHE radiation surrounding the sources. These observations can probe the particle propagation inside and outside the accelerators and the subsequent injection/escape into the interstellar medium.

We use semi-analytical models to study the effects of primordial black hole (PBH) accretion on the cosmic radiation background during the epoch of reionization ($z\gtrsim 6$). We consider PBHs floating in the intergalactic medium (IGM), and located inside haloes, where star formation can occur. For stars with a mass $\gtrsim 25 \rm\ M_{\odot}$, formed in suitable host haloes, we assume they quickly burn out and form stellar remnant black holes (SRBHs). Since SRBHs also accrete material from their surroundings, we consider them to have similar radiation feedback as PBHs in the halo environment. To estimate the background radiation level more accurately, we take into account the impact of PBHs on structure formation, allowing an improved modeling of the halo mass function. We consider the radiation feedback from a broad suite of black holes: PBHs, SRBHs, high-mass X-ray binaries (HMXBs), and supermassive black holes (SMBHs). We find that at $z\gtrsim 30$, the radiation background energy density is generated by PBHs accreting in the IGM, whereas at lower redshifts, the accretion feedback power from haloes dominates. We also analyze the total power density by modeling the accretion spectral energy distribution (SED), and break it down into select wavebands. In the UV band, we find that for $f_{\rm PBH} \lesssim 10^{-3}$, the H-ionizing and Lyman-$\alpha$ fluxes from PBH accretion feedback do not violate existing constraints on the timing of reionization, and on the effective Wouthuysen-Field coupling of the 21-cm spin temperature of neutral hydrogen to the kinetic temperature of the IGM. However, in the X-ray band, with the same abundance, PBHs contribute significantly and could account for the unresolved part of the cosmic X-ray background.

Dusty Star-Forming Galaxies (DSFGs) are amongst the most massive and active star-forming galaxies during the cosmic noon. Theoretical studies have proposed various formation mechanisms of DSFGs, including major merger-driven starbursts and secular star-forming disks. Here, we report J0107a, a bright ($\sim8$ mJy at observed-frame 888 $\mu$m) DSFG at $z=2.467$ that appears to be a gas-rich massive disk and might be an extreme case of the secular disk scenario. J0107a has a stellar mass $M_\star\sim10^{11}M_\odot$, molecular gas mass $M_\mathrm{mol}\gtrsim10^{11}M_\odot$, and a star formation rate (SFR) of $\sim500M_\odot$ yr$^{-1}$. J0107a does not have a gas-rich companion. The rest-frame 1.28 $\mu$m JWST NIRCam image of J0107a shows a grand-design spiral with a prominent stellar bar extending $\sim15$ kpc. ALMA band 7 continuum map reveals that the dust emission originates from both the central starburst and the stellar bar. 3D disk modeling of the CO(4-3) emission line indicates a dynamically cold disk with rotation-to-dispersion ratio $V_\mathrm{max}/\sigma\sim8$. The results suggest a bright DSFG may have a non-merger origin, and its vigorous star formation may be triggered by bar and/or rapid gas inflow.

It is an open issue how the surrounding environment of supernova remnant shocks affect nonthermal X-rays from accelerated electrons, with or without interacting dense material. We have conducted spatially resolved X-ray spectroscopy of the shock-cloud interacting region of RCW 86 with XMM-Newton. It is found that bright soft X-ray filaments surround the dense cloud observed with 12CO and HI emission lines. These filaments are brighter in thermal X-ray emission, and fainter and possibly softer in synchrotron X-rays, compared to those without interaction. Our results show that the shock decelerates due to the interaction with clouds, which results in an enhancements of thermal X-ray emission. This could possibly also explain the softer X-ray synchrotron component, because it implies that those shocks that move through a low density environment, and therefore decelerate much less, can be more efficient accelerators. This is similar to SN 1006 and Tycho, and is in contrast to RX J1713.7-3946. This difference among remnants may be due to the clumpiness of dense material interacting with the shock, which should be examined with future observations.

(Aims) Near-infrared imaging polarimetry at high-angular resolutions has revealed intriguing distribution of circumstellar dust towards FU Ori-type objects (FUors). These dust grains are probably associated with either an accretion disk or an infalling envelope. Follow-up observations in the mid-infrared would lead us to a better understanding of the hierarchy of the mass accretion processes onto FUors (i.e., envelope and disk accretion), which hold keys for understanding the mechanism of their accretion outbursts and the growth of low-mass young stellar objects in general. (Method) We have developed a semi-analytic method to estimate the mid-infrared intensity distributions using the observed polarized intensity (PI) distributions in H-band (lambda=1.65 micron). We have derived intensity distributions for two FUors, FU Ori and V1735 Cyg, at three wavelengths (lambda=3.5/4.8/12 micron) for various cases, i.e. with a star or a flat compact self-luminous disk as an illuminating source; an optically thick disk or an optically thin envelope for circumstellar dust grains; and three different dust models. (Results) We have been able to obtain self-consistent results for many cases and regions, in particular when the viewing angle of the disk/envelope is zero (face-on). Our calculations suggest that the mid-infrared extended emission at the above wavelengths is dominated by the single scattering process. The contribution of thermal emission is negligible unless we add an additional heating mechanism such as adiabatic heating in spiral structures and/or fragments. The uncertain nature of the central illuminating source, the distribution of circumstellar dust grains and the optical properties of dust grains yield uncertainties in the intensity levels on orders of magnitude, e.g., 20-800, for the disk/envelope aspect ratios of ~0.2 and lambda=3-13 micron.

Many models of dark matter (DM) are now widely considered and probed intensively with accelerators, underground detectors, and astrophysical experiments. Among the various approaches, high-energy astrophysical observations are extremely useful to complement laboratory searches for some DM candidates. In the near future, the Cherenkov Telescope Array (CTA) should enable us to access much heavier weakly interacting massive particles, as well as a broad range of other DM candidates. In this talk, we describe DM searches with CTA.

In this document we describe a model of an array of water Cherenkov detectors proposed to study ultra-high-energy gamma rays in the southern hemisphere, and a continuous model of secondary particles produced on the ground from gamma and proton showers. We use the model of the detector and the parametrization of showers for the identification of the most promising configuration of detector elements, using a likelihood ratio test statistic to classify showers and a stochastic gradient descent technique to maximize a utility function describing the measurement precision on the gamma-ray flux.

Nuclear Stellar Discs have been observed in the vast majority of barred disc galaxies including the Milky Way. Their intense star formation is sustained by the intense gas inflows driven by their surrounding bars and frequently supports a large-scale galactic fountain. Despite their central role in galaxy evolution, their chemical evolution remains largely unexplored. Here we present the first systematic, multizonal modelling of the chemical evolution of these nuclear stellar discs. We argue that the chemical composition of nuclear stellar discs is best understood relative to the bar tips from which their gas is drawn. We make predictions of the detailed abundance profiles within the nuclear stellar disc under different accretion scenarios from the galactic bar and show by which observable differences they can be decided. We show that with their difference to normal disc star formation, nuclear discs offer a unique laboratory to break parameter degeneracies in chemical evolution models. This allows us to identify the effects of the main gas parameters and disentangle them from the global enrichment history. We show this for the example of the ejection fractions from the outer and nuclear disk.

We characterize the evolution of the rest-frame 1500 $\unicode{xC5}$ UV luminosity Function (UVLF) from AstroSat/UVIT F154W and N242W imaging in the Great Observatories Origins Survey North (GOODS-N) field. With deep FUV observations, we construct the UVLF for galaxies at z$<0.13$ and subsequently characterize it with a Schechter function fit. The fitted parameters are consistent with previous determinations. With deep NUV observations, we are able to construct the UVLF in seven redshift bins in the range z $\sim$ 0.8 - 0.4, with galaxies identified till $\sim$2 mag fainter than previous surveys, owing to the high angular-resolution of UVIT. The fitted Schechter function parameters are obtained for these UVLFs. At z $\sim$ 0.8 - 0.7, we also utilize Hubble Space Telescope (HST) F275W observations in the GOODS-N field to construct the UVLF in 2 redshift bins, whose fitted Schechter function parameters are then found to be consistent with that determined from UVIT at z $\sim$ 0.75. We thus probe the variation of the fitted UVLF parameters over z $\sim$ 0.8 - 0.4, a span of $\sim$2.7 Gyr in age. We find that the slope of the Schechter function, $\alpha$, is at its steepest at z $\sim$ 0.65, implying highest star-formation at this instant with galaxies being relatively more passive before and after this time (also confirmed with FUV star-formation rate determinations). We infer that this is a short-lived instance of increased cosmic star-formation even though cosmic star-formation may be winding-down over longer timespan at this redshift range.

We present the measurement and characterization of the brighter-fatter effect (BFE) on a NASA Roman Space Telescope development Teledyne H4RG-10 near-infrared detector using laboratory measurements with projected point sources. After correcting for other interpixel non-linearity effects such as classical non-linearity and inter-pixel capacitance, we quantify the magnitude of the BFE by calculating the fractional area change per electron of charge contrast. We also introduce a mathematical framework to compare our results with the BFE measured on similar devices using autocorrelations from flat-field images. We find an agreement of 18 +/- 5% between the two methods. We identify potential sources of discrepancy and discuss future investigations to characterize and address them.

Understanding the complex evolution of protoplanetary disks (PPDs) and their dispersal via energetic stellar radiation are prominent challenges in astrophysics. It has recently been established that specifically the X-ray luminosity from the central protostar can significantly heat the surface of the disk, causing powerful photoevaporative winds that eject a considerable fraction of the disc's mass. Recent work in the field has moreover shown the importance of global PPD simulations that simultaneously take into account non-ideal magnetohydrodynamic (MHD) effects and detailed thermochemistry. Our motivation with the current paper lies in combining these two aspects and figure out how they interact. Focus is put on the Hall Effect (HE) and the impact it has on the overall field topology and mass loss/accretion rates. Utilizing a novel X-ray temperature parametrisation, we perform 2D-axisymmetric MHD simulations with the NIRVANA fluid code, covering all non-ideal effects. We find that, in the aligned orientation, the HE causes prominent inward displacement of the poloidal field lines that can increase the accretion rate through a laminar Maxwell stress. We find that outflows are mainly driven by photoevaporation -- unless the magnetic field strength is considerable (i.e., $\beta_p\leq 10^{3}$) or the X-ray luminosity low enough (i.e., $\log{L_X}\leq 29.3$). Inferred mass loss rate are in the range of the expected values $10^{-8}$ to $10^{-7}M_{\odot}yr^{-1}$. For comparison, we have also performed pure hydrodynamic (HD) runs and compared them with the equivalent MHD runs. Here we have found that the magnetic field does indeed contribute to the mass loss rate, albeit only discernibly so for low enough $L_X$ (i.e., $\log{L_X}\leq 30.8$). For values higher than that, the wind mass loss predicted from the MHD set converges to the ones predicted from pure HD.

Bounds are derived on the axion-like particle (ALP) to two-photon coupling in the mass range $2.65-5.27$ eV. The bounds are obtained by searching for the signal from ALP decay in the Multi Unit Spectroscopic Explorer (MUSE) observations of five dwarf spheroidal galaxies, under the assumption that ALPs constitute the dark matter component of the haloes. These bounds are of the same order and improve on the robustness of those of Reference~\cite{Regis}, and currently represent the strongest bounds within the considered mass range.

Approximated forms of the RII and RIII redistribution matrices are frequently applied to simplify the numerical solution of the radiative transfer problem for polarized radiation, taking partial frequency redistribution (PRD) effects into account. A widely used approximation for RIII is to consider its expression under the assumption of complete frequency redistribution (CRD) in the observer frame (RIII CRD). The adequacy of this approximation for modeling the intensity profiles has been firmly established. By contrast, its suitability for modeling scattering polarization signals has only been analyzed in a few studies, considering simplified settings. In this work, we aim at quantitatively assessing the impact and the range of validity of the RIII CRD approximation in the modeling of scattering polarization. Methods. We first present an analytic comparison between RIII and RIII CRD. We then compare the results of radiative transfer calculations, out of local thermodynamic equilibrium, performed with RIII and RIII CRD in realistic 1D atmospheric models. We focus on the chromospheric Ca i line at 4227 A and on the photospheric Sr i line at 4607 A.

Outgassing is a central process during the formation and evolution of terrestrial planets and their atmospheres both within and beyond the solar system. Although terrestrial planets' early atmospheres likely form via outgassing during planetary accretion, the connection between a planet's bulk composition and its initial atmospheric properties is not well understood. One way to inform this connection is to analyze the outgassing compositions of meteorites, and in particular carbonaceous chondrites, because they are some of the most volatile-rich, primitive materials (in terms of their bulk compositions) that are available for direct study. In addition, they may serve as compositional analogs for the building block materials of terrestrial planets in our solar system and around other Sun-like stars. This study builds upon previous outgassing experiments that monitored the abundances of volatile species (e.g., H2O, CO, and CO2) released from the Murchison meteorite. To gain a more complete understanding of Murchison's outgassing composition, we perform a series of heating experiments under atmospheric pressure (1 bar) and vacuum (1E-9 bar) conditions on samples of the Murchison meteorite and subsequent bulk element analysis to inform the outgassing trends of a suite of major elements in Murchison (e.g., Fe, Mg, Zn, and S). Under both pressure conditions, sulfur outgases significantly at the highest temperatures (800C - 1000C). For the samples heated under vacuum conditions, we also detect outgassing of zinc. Combined with prior outgassing experiments, this study provides important insights into the volatile depletion patterns of undifferentiated planetesimals and the early outgassing compositions of terrestrial exoplanets.

We report the detection of outflowing molecular gas at the center of the nearby radio galaxy NGC6328 (z=0.014), which has a gigahertz-peaked spectrum radio core and a compact (2 pc) young double radio lobe tracing jet. Utilizing Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations, as well as a novel code developed to fit the 3D gas distribution and kinematics, to study the molecular gas kinematics, we find that the bulk of the gas is situated within a highly warped disk structure, most likely the result of a past merger event. Our analysis further uncovers, within the inner regions of the gas distribution (R<300 pc) and at a position angle aligning with that of the radio jet (150 degrees), the existence of two anti-diametric molecular gas structures kinematically detached from the main disk. These structures most likely trace a jet-induced cold gas outflow with a total lower limit mass of $2\times 10^6\,\mathrm{M_\odot}$ mass, corresponding to an outflow rate of $2\,\mathrm{M_\odot\,yr^{-1}}$ and a kinetic power of $2.7\times 10^{40}\,\mathrm{erg\,s^{-1}}$. The energy required to maintain such a molecular outflow is aligned with the mechanical power of the jet.

We present a case study of an M4-class flare on 28 March 2022, near Solar Orbiter's first science perihelion (0.33 AU). Solar Orbiter was 83.5{\deg} west of the Sun-Earth line, making the event appear near the eastern limb, while Earth-orbiting spacecraft observed it near the disk center. The timing and location of the STIX X-ray sources were related to the plasma evolution observed in the EUV by the Extreme Ultraviolet Imager (EUI) on Solar Orbiter and the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory, and to the chromospheric response observed in 1600 {\AA} by AIA. We performed differential emission measure (DEM) analysis to further characterize the flaring plasma at different subvolumes. The pre-flare magnetic field configuration was analyzed using a nonlinear force-free (NLFF) extrapolation. In addition to the two classical hard X-ray (HXR) footpoints at the ends of the flaring loops, later in the event we observe a nonthermal HXR source at one of the anchor points of the erupting filament. The evolution of the AIA 1600~{\AA} footpoints indicates that this change in footpoint location represents a discontinuity in an otherwise continuous westward motion of the footpoints throughout the flare. The NLFF extrapolation suggests that strongly sheared field lines close to, or possibly even part of, the erupting filament reconnected with a weakly sheared arcade during the first HXR peak. The remainder of these field lines reconnected later in the event, producing the HXR peak at the southern filament footpoint. Our results show that the reconnection between field lines with very different shear in the early phase of the flare plays a crucial role in understanding the motion of the HXR footpoint during later parts of the flare. This generalizes simpler models, such as whipping reconnection, which only consider reconnection propagating along uniformly sheared arcades.

Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution data from The HI/OH/Recombination line survey of the Milky Way (THOR), we aim to constrain the kinematic and physical properties of the cold atomic hydrogen gas phase toward the inner Galactic plane. HI self-absorption (HISA) has proven to be a viable method to detect cold atomic hydrogen clouds in the Galactic plane. With the help of a newly developed self-absorption extraction routine (astroSABER), we build upon previous case studies to identify HI self-absorption toward a sample of Giant Molecular Filaments (GMFs). We find the cold atomic gas to be spatially correlated with the molecular gas on a global scale. The column densities of the cold atomic gas traced by HISA are usually of the order of $10^{20}\rm\,cm^{-2}$ whereas those of molecular hydrogen traced by $\rm^{13}CO$ are at least an order of magnitude higher. The HISA column densities are attributed to a cold gas component that accounts for a fraction of $\sim$5% of the total atomic gas budget within the clouds. The HISA column density distributions show pronounced log-normal shapes that are broader than those traced by HI emission. The cold atomic gas is found to be moderately supersonic with Mach numbers of a $\sim$few. In contrast, highly supersonic dynamics drive the molecular gas within most filaments.

Upcoming deep optical surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time will scan the sky to unprecedented depths and detect billions of galaxies. This amount of detections will however cause the apparent superposition of galaxies on the images, called blending, and generate a new systematic error due to the confusion of sources. As consequences, the measurements of individual galaxies properties such as their redshifts or shapes will be impacted, and some galaxies will not be detected. However, galaxy shapes are key quantities, used to estimate masses of large scale structures, such as galaxy clusters, through weak gravitational lensing. This work presents a new catalog matching algorithm, called friendly, for the detection and characterization of blends in simulated LSST data for the DESC Data Challenge 2. By identifying a specific type of blends, we show that removing them from the data may partially correct the amplitude of the $\Delta\Sigma$ weak lensing profile that could be biased low by around 20% due to blending. This would result in impacting clusters weak lensing mass estimate and cosmology.

Forward shocks by radio jets, driven into the intracluster medium, are one of the indicators that can be used to evaluate the power of the jet. Meanwhile high-angular-resolution X-ray observations show the Mach numbers of powerful radio jets are smaller compared to that of theoretical and numerical studies, $\mathcal{M_{\rm obs}} < 2$. Our aim is to systematically investigate various factors, such as projection effects and temperature non-equilibration between protons and electrons, that influence the Mach number estimate in a powerful jet. Using two-temperature magnetohydrodynamic simulation data for the Cygnus A radio jets, whose Mach number is approximately 6, we construct mock X-ray maps of simulated jets from various viewing angles. Further, we evaluate the shock Mach number from density/temperature jump using the same method of X-ray observations. Our results demonstrate that measurements from density jump significantly underestimate the Mach numbers, $\mathcal{M} < 2$, around the jet head at a low viewing angle, $\lessapprox 50^{\circ}$. The observed post-shock temperature is strongly reduced by the projection effect, as our jet is in the cluster center where the gas density is high. On the other hand, the temperature jump is almost unity, even if thermal electrons are in instant equilibration with protons. Upon comparison, we find that shock property of our model at viewing angle of $<$ $55^{\circ}$ is in a good agreement with that of Cygnus A observations.

Non-radiating protons in the radio lobes have an essential role to form the jet morphology which is shown by recent radio and X-ray observations. However, since protons and electrons are not always in energy equilibrium due to weak Coulomb coupling, it is difficult to estimate the energy contribution of protons for inflation of radio lobes. The main focus of this study is to examine the effect of the variable model for electron heating by turbulence and shock waves on the thermal energy distribution of electron and proton. We performed two-temperature three-dimensional magnetohydrodynamic simulations of sub-relativistic jets in the galaxy cluster while varying jet magnetization parameters. Because the energy partition rate between electrons and protons in shock and turbulence is determined by plasma kinetic scale physics, our global simulations include electron instantaneous heating sub-grid models for shock waves and turbulence. We find that most of the bulk kinetic energy of the jet is converted into thermal energy of protons through both shocks and turbulence. Thus, protons are energetically dominant. Meanwhile, thermal electrons stored in the lobe evolve toward energy equipartition with magnetic energy through turbulent dissipation. We further estimated the radio power and the mechanical jet power of radio lobes following the same method as used for radio and X-ray observations, and compared these powers with that of the observed radio jets. The two-temperature model quantitatively explains the radiatively inefficient radio cavities, but cannot reproduce the radiatively efficient cavity, even for strong magnetized jets. This implies that a significant population of pair-plasma is needed in the radiatively efficient radio cavities.

Magnetic activity and rotation are related to the age of low-mass main sequence stars. To further constrain these relations, we study a sample of 574 main sequence stars members of common proper motion pairs with white dwarfs, identified thanks to Gaia astrometry. We use the white dwarfs as age indicators, while the activity indexes and rotational velocities are obtained from the main sequence companions using standard procedures. We find that stars older than 5 Gyr do not display Halpha nor Caii H&K emission unless they are fast rotators due to tidal locking from the presence of unseen companions and that the rotational velocities tend to decrease over time, thus supporting the so-called gyrochronology. However, we also find moderately old stars (~2-6 Gyr) that are active presumably because they rotate faster than they should for their given ages. This indicates that they may be suffering from weakened magnetic braking or that they possibly evolved through wind accretion processes in the past. The activity fractions that we measure for all stars younger than 5 Gyr range between ~10-40 per cent. This is line with the expectations, since our sample is composed of F, G, K and early M stars, which are thought to have short (<2 Gyr) activity lifetimes. Finally, we observe that the Halpha fractional luminosities and the R'HK indexes for our sample of (slowly rotating) stars show a spread (-4>log(LHalpha/Lbol); log(R'HK)> -5) typically found in inactive M stars or weakly active/inactive F, G, K stars.

The Cherenkov Telescope Array Observatory (CTAO) is a next-generation ground-based gamma-ray observatory that will study the universe at very high energy using atmospheric Cherenkov light. CTAO will comprise over 67 telescopes of three different sizes, located in the northern and southern hemispheres. Among these, the Medium-Sized Telescope (MST) will play a crucial role in CTAO's observations, providing excellent sensitivity and angular resolution for gamma rays in the energy range of 100 GeV to 5 TeV. The MST is based on a modified single-mirror Davies-Cotton design, featuring a segmented mirror with a diameter of 12 meters, a total reflective surface of 88 square meters, and a focal length of 16 meters. It will cover an approximately 8-degree field of view and be equipped with two different cameras: NectarCAM and FlashCam, at the northern and southern CTAO sites, respectively. The MST's design is optimized for efficient observation of extended sources, including supernova remnants and pulsar wind nebulae, as well as the study of gamma-ray bursts and active galactic nuclei. Currently, the MST is in the midst of production and testing stages with the aim of being installed in 2025 for the CTAO Pathfinder project. In this project, one MST telescope will be deployed at each CTAO site to provide on-site experience with pre-production components. This approach facilitates cost and risk reduction before starting serial production. This contribution provides an overview of the current status and plans of the MST's construction at both the northern and southern CTAO sites, as well as details on the telescope and camera designs and their expected performance.

Contamination by Radio Frequency Interference (RFI) is a ubiquitous challenge for radio astronomy. In particular, transient RFI is difficult to detect and avoid, especially in large data sets with many time bins. In this work, we present a Bayesian methodology for time-dependent, transient anomaly mitigation. In general, the computation time for correcting for transient anomalies in time-separated data sets grows proportionally with the number of time bins. We demonstrate that utilising likelihood reweighting can allow our Bayesian anomaly mitigation method to be performed with a computation time close to independent of the number of time bins. In particular, we identify a factor of 25 improvement in computation time for a test case with 2000 time bins. We also demonstrate how this method enables the flagging threshold to be fit for as a free parameter, fully automating the mitigation process. We find that this threshold fitting also prevents overcorrecting of the data in the case of wide priors. Finally, we investigate the potential of the methodology as a transient detector. We demonstrate that the method is able to reliably flag an individual anomalous data point out of 302,000 provided the SNR > 10.

The Front-End Board (FEB) is a key component of the NectarCAM camera, which has been developed for the Medium-Sized-Telescopes (MST) of the Cherenkov Telescope Array Observatory (CTAO). The FEB is responsible for reading and converting the signals from the camera's photo-multiplier tubes (PMTs) into digital data, as well as generating module level trigger signals. This contribution provides an overview of the design and performances of a new version of the FEB that utilizes an improved version of the NECTAr chip. The NECTAr chip includes a switched capacitor array for sampling signals at 1 GHz, and a 12-bit analog-to-digital converter (ADC) for digitizing each sample when the trigger signal is received. The integration of this advanced NECTAr chip significantly reduces the deadtime of NectarCAM by an order of magnitude as compared to the previous version. This contribution also presents the results of laboratory testing of the new FEB, including measurements of timing performance, linearity, dynamic range, and deadtime.

We release photometric redshifts, reaching $\sim$0.7, for $\sim$14M galaxies at $r\leq 20$ in the 11,500 deg$^2$ of the SDSS north and south galactic caps. These estimates were inferred from a convolution neural network (CNN) trained on $ugriz$ stamp images of galaxies labelled with a spectroscopic redshift from the SDSS, GAMA and BOSS surveys. Representative training sets of $\sim$370k galaxies were constructed from the much larger combined spectroscopic data to limit biases, particularly those arising from the over-representation of Luminous Red Galaxies. The CNN outputs a redshift classification that offers all the benefits of a well-behaved PDF, with a width efficiently signaling unreliable estimates due to poor photometry or stellar sources. The dispersion, mean bias and rate of catastrophic failures of the median point estimate are of order $\sigma_{\rm MAD}=0.014$, <$\Delta z_{\rm norm}$>$=0.0015$, $\eta(|\Delta z_{\rm norm}|>0.05)=4\%$ on a representative test sample at $r<19.8$, out-performing currently published estimates. The distributions in narrow intervals of magnitudes of the redshifts inferred for the photometric sample are in good agreement with the results of tomographic analyses. The inferred redshifts also match the photometric redshifts of the redMaPPer galaxy clusters for the probable cluster members. The CNN input and output are available at: https://deepdip.iap.fr/treyer+2023.

Although spectroscopic observations of RR Lyrae stars have now been carried out for almost a century, it is only recently that it has been identified that the hydrogen line exhibits three successive emissions in each pulsation cycle. The purpose of this study is to clarify the physical origin of these three emissions and their connection with atmospheric dynamics and the influence of Blazhko modulation on their intensity. From the value of the blueshift of the main H$\alpha$ emission, the velocity of the hypersonic shock front was estimated between 100 and 150 $\pm 10$ km$\,$s$^{-1}$ (Mach number between 10 and 15). It has been established that the shock velocity increases from the minimum Blazhko to its maximum, and after that, gradually decreases to the Blazhko minimum to start growing again. This observational result is consistent with the shock model proposed in 2013 to explain the Blazhko effect. The intensity of the H$\alpha$ emission increases with the shock velocity up to a maximum value around 137 km$\,$s$^{-1}$ and then decreases as the shock velocity increases further. This effect would be the consequence of the more and more important ionization of the atoms in the radiative shock wake. The second (blueshifted) H$\alpha$ emission, is the consequence of an approximately constant supersonic compression (Mach number between 2 and 3) of the upper atmosphere falling on the photospheric layers, during 3 to 16% of the pulsation period. Finally, the third H$\alpha$ emission (P-Cygni profile) would be the consequence of the expansion of the high atmosphere induced by the shock wave during its final weakening.

Photometric redshift estimation plays a crucial role in modern cosmological surveys for studying the universe's large-scale structures and the evolution of galaxies. Deep learning has emerged as a powerful method to produce accurate photometric redshift estimates from multi-band images of galaxies. Here, we introduce a multimodal approach consisting of the parallel processing of several subsets of image bands prior, the outputs of which are then merged for further processing through a convolutional neural network (CNN). We evaluate the performance of our method using three surveys: the Sloan Digital Sky Survey (SDSS), The Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) and Hyper Suprime-Cam (HSC). By improving the model's ability to capture information embedded in the correlation between different bands, our technique surpasses the state-of-the-art photometric redshift precision. We find that the positive gain does not depend on the specific architecture of the CNN and that it increases with the number of photometric filters available.

We report measurements of H I content in 11 nearby, actively star-forming, blue compact dwarf galaxies (BCDs) from 21 cm observations with the Arecibo telescope. These BCDs, selected by their red (W2[4.6 $\mu$m]$-$W3[12 $\mu$m]$>$3.8 mag) and bright mid-infrared (MIR) emission (W4[22 $\mu$m]$<$ 7.6 mag), have high specific star formation rates (median sSFR $\sim$10$^{-7.8}$ yr$^{-1}$), similar to high redshift galaxies. H I emission was detected in six sources. We analyze our new detections in the context of previous H I observations of 218 dwarf irregulars (dIs) and BCDs in the literature. The $M_{\rm HI}$-$M_{\ast}$ relation resulting from our observations confirms the dominating fraction of H I gas among baryons in galaxies with lower stellar masses. This Arecibo BCD sample has significantly lower median H I depletion timescales ($\tau_{\rm HI}\sim$ 0.3 Gyr) than other dIs/BCDs ($\sim$ 6.3 Gyr) in the literature. The majority of the sources (10/11) in the Arecibo sample are very red in W1[3.4 $\mu$m]$-$W2[4.6 $\mu$m] colour ($>$ 0.8 mag) implying the presence of warm dust. We investigate the relation of $\tau_{\rm HI}$ with stellar mass ($M_{\ast}$) and sSFR. We find that $\tau_{\rm HI}$ is significantly anti-correlated with $M_{\ast}$ for higher sSFR ($>$10$^{-8.5}$ yr$^{-1}$) and with sSFR for higher stellar mass ($>10^{7.5}\,{\rm M}_{\odot}$) dwarf galaxies. The high sSFR for the BCDs in the Arecibo observed sample is mainly due to their high atomic gas star formation efficiency (SFE) or low $\tau_{\rm HI}$. The low $\tau_{\rm HI}$ or high SFE in these sources is possibly due to runaway star formation in compact and dense super star clusters.

The IceCube neutrino observatory detected two distinct flares of high-energy neutrinos from the direction of the blazar TXS 0506+056: a $\sim 300$ TeV single neutrino on September 22, 2017 and a $3.5\sigma$ signature of a dozen TeV neutrinos in 2014/2015. In a previous work, it was shown that these two episodes of neutrino emission could be due to an inspiral of a supermassive binary black hole (SMBBH) close to its merger at the core of TXS 0506+056. Such an inspiral can lead to quasi-periodic particle emission due to jet precession close to the final coalescence. This model made predictions on when the next neutrino emission episode must occur. On September 18, 2022, IceCube detected an additional, $\sim 170$ TeV neutrino in directional coincidence with the blazar TXS 0506+056, being consistent with the model prediction. Additionally, in April 2021, the Baikal Collaboration reported the detection of a $224\pm 75$ TeV neutrino, with TXS 0506+056 being in the uncertainty range of the event direction. We show that these four distinct flares of neutrino emission from TXS 0506+056 are consistent with a precessing jet scenario, driven by an inspiraling SMBBH. Using improved modeling, we are now able to constrain the total mass together with the mass ratio for the binary. We predict when the next neutrino flares from TXS 0506+056 should be happening. Finally, we estimate the detection potential of the Laser-interferometer Space Antenna (LISA) for the merger in the future.

We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Faraday dispersion (sigma_phi). From THOR, we find median sigma_phi of 15.9 +/- 3.2 rad m^-2 for G46.8-0.3 and 17.6 +/- 1.6 rad m^-2 for G39.2-0.3. When comparing to polarization at 6cm, we find evidence for sigma_phi >~ 30 rad m^-2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium (ISM) using the SF of extra-Galactic radio sources (EGRS) and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales.

Nova Her 2021 was observed with TESS 12.62 days after its most recent outburst in June 12.537 2021. This cataclysmic variable belongs to the intermediate polar class, with an spin period of $\sim$501 s and orbital period of 0.1529 days. During TESS observations of Sector 40, the orbital period of 0.1529(1) days is detected significantly 17 days after the onset of the outburst. A modulation, of unknown origin, with a period of $\sim$0.537 days is present in the data from day 13 until day 17.

We present a comprehensive spectro-polarimetric study of persistent Black hole X-ray binary 4U $1957\!+\!115$ with {\it IXPE} and {\it NICER} observations. The source is observed in disk dominated thermal state with disk temperature, kT$_{in}\approx1.4$ keV. The emission during thermal state from the source is found to be moderately polarized and \textit{IXPE} measures a degree of polarization (PD) $= 1.95\pm0.37\%$ ($>4.79\sigma$) along with a polarization angle (PA) $= -42.78^{\circ}\pm5.41^{\circ}$ in the energy range of $2-8$ keV. PD is found to be an increasing function of energy, whereas PA indicates switching within the energy range which could be due to high inclination and the returning radiation within the system. Simultaneous energy spectra ($0.6-10$ keV) from {\it NICER} are modelled to study the spectral properties. Furthermore, the spin parameter of the black hole is estimated with spectro-polarimetric data as a$_{\ast}=0.988\pm0.001\,(1\sigma)$, which is corroborated by {\it NICER} observations. Finally, we discuss the implications of our findings.

High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of-view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g., LEM, HUBS, and SuperDIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g., O VII/VIII and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to $r_{200}$ and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g., temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology.

We present a direct imaging study of V892 Tau, a young Herbig Ae/Be star with a close-in stellar companion and circumbinary disk. Our observations consist of images acquired via Keck 2/NIRC2 with non-redundant masking and the pyramid wavefront sensor at K$^\prime$ band (2.12$\mu$m) and L$^\prime$ band (3.78$\mu$m). Sensitivity to low-mass accreting companions and cool disk material is high at L$^\prime$ band, while complimentary observations at K$^\prime$ band probe hotter material with higher angular resolution. These multi-wavelength, multi-epoch data allow us to differentiate the secondary stellar emission from disk emission and deeply probe the structure of the circumbinary disk at small angular separations. We constrain architectural properties of the system by fitting geometric disk and companion models to the K$^\prime$ and L$^\prime$ band data. From these models, we constrain the astrometric and photometric properties of the stellar binary and update the orbit, placing the tightest estimates to date on the V892 Tau orbital parameters. We also constrain the geometric structure of the circumbinary disk, and resolve a circumprimary disk for the first time.

We present a method to identify likely visual binaries in Gaia eDR3 that does not rely on parallax or proper motion. This method utilizes the various PSF sizes of 2MASS/Gaia, where at $<2.5$" two stars may be unresolved in 2MASS but resolved by Gaia. Due to this, if close neighbors listed in Gaia are a resolved pair, the associated 2MASS source will have a predictable excess in the J-band that depends on the $\Delta G$ of the pair. We demonstrate that the expected relationship between 2MASS excess and $\Delta G$ differs for chance alignments, as compared to true binary systems, when parameters like magnitude and location on the sky are also considered. Using these multidimensional distributions, we compute the likelihood of a close pair of stars to be a chance alignment, resulting in a total(clean) catalog of 68,725(50,230) likely binaries within 200 pc with a completeness rate of $\sim75\%$($\sim64\%$) and contamination rate of $\sim14\%$($\sim0.4\%$). Within this, we find 590 previously unidentified binaries from Gaia eDR3 with projected physical separations $<30$ AU, where 138 systems were previously identified, and for $s<10$ AU we find that 4 out of 15 new likely binaries have not yet been observed with high-resolution imaging. We also demonstrate the potential of our catalog to determine physical separation distributions and binary fraction estimates, from this increase in low separation binaries. Overall, this catalog provides a good complement for the study of local binary populations by probing smaller physical separations and mass ratios, and provides prime targets for speckle monitoring.

The recent detection of tera-electronvolt (TeV) photons from the record-breaking gamma-ray burst GRB 221009A during its prompt phase poses challenges for constraining its Lorentz factor. We re-evaluate the constraints on the jet Lorentz factor considering a two-zone model, wherein the TeV photons originate from the external shock region while the lower energy MeV photons come from the internal prompt emission region. By properly accounting for the evolution of the MeV photon spectrum and light curve, we calculate the optical depth for TeV photons and derive a minimum Lorentz factor about 300. It is consistent with the afterglow modeling for the TeV emission.

We have found the 3-dimensional structures of stable, long-lived planetary vortices. Detailed observations of the horizontal velocities of Jovian vortices exist at only one height in the atmosphere, making their vertical structures poorly understood. We solve the 3-dimensional anelastic equations with a high-resolution pseudo-spectral method using observed Jovian hydrostatic atmospheric temperatures and zonal flow. We examine several families of vortices. We find that {\it constant-vorticity} vortices, which have nearly-uniform vorticity as a function of height and horizontal areas that go to zero at their tops and bottoms, converge to stable vortices that look like the Great Red Spot (GRS). In contrast, the {\it constant-area} vortices proposed in previous studies, which have nearly-uniform areas as a function of height and vertical vorticities that go to zero at their tops and bottoms, are far from equilibrium, break apart, and converge to {\it constant-vorticity} vortices. Our final vortices show other unexpected properties. Vortices that are initially non-hollow become hollow (i.e., have local minima of vertical vorticity at their center), which is a feature of the GRS that cannot be explained with 2-dimensional simulations. The central axes of the final vortices align with the planetary-spin axis even if they initially align with the local direction of gravity. We show how the small ratio of the magnitude of the vertical-to-horizontal velocity scales with the Rossby number and the vertical aspect ratio of the vortices, and analytically prove that the horizontal mid-plane of a stable vortex must lie at a height above the top of the convective zone.

I present a highly efficient integration method for scattering calculations, which can reduce the evaluation time of particularly challenging analyses by a factor of $10^7$. By projecting each item in the integrand onto vector spaces spanned by simple basis functions, a multidimensional numerical integral is replaced by a much easier exercise in matrix multiplication. This method can be generalized to any integrand that depends linearly on multiple input functions. In this paper I apply the method to dark matter (DM) direct detection with anisotropic target materials, where the DM scattering rate depends not only on the DM particle model, the astrophysical DM velocity distribution, and the properties of the target material, but now also on the $SO(3)$ orientation of the detector. The vector space calculation factorizes the astrophysics, the DM particle model, and the SM material properties from each other, while describing detector rotations by simple matrix multiplication acting on the basis functions. Difficult analyses, which otherwise would take decades or centuries on a computing cluster, can now be performed in a few hours on a laptop.

A new integration method drastically improves the efficiency of the dark matter direct detection calculation. In this work I introduce a complete, orthogonal basis of spherical wavelet-harmonic functions, designed for the new vector space integration method. This factorizes the numeric calculation into a part that depends only on the astrophysical velocity distribution; a second part, depending only on the detector form factor; and a scattering matrix defined on the basis functions, which depends on the details of the dark matter (DM) particle model (e.g. its mass). For common spin-independent DM-Standard Model interactions, this scattering matrix can be evaluated analytically in the wavelet-harmonic basis. This factorization is particularly helpful for the more complicated analyses that have become necessary in recent years, especially those involving anisotropic detector materials or more realistic models of the local DM velocity distribution. With the new method, analyses studying large numbers of detector orientations and DM particle models can be performed about 10 million times faster. This paper derives several analytic results for the spherical wavelets, including an extrapolation in the space of wavelet coefficients, and a generalization of the vector space method to a much broader class of linear functional integrals. Both results are highly relevant outside the field of DM direct detection.

Detection of gravitational waves (GWs) paves the beginning of a new era of gravitational wave astronomy. Black holes (BHs) in their ringdown phase provide the cleanest signal of emitted GWs that imprint the fundamental nature of BHs under low energy perturbation. Apart from GWs, any complementary signature of ringing BHs can be of paramount importance. Motivated by this we analyzed the scattering of electromagnetic waves in such a background and demonstrated that the absorption cross section of a ringing Schwarzschild BH can be superradiant. Moreover, we have found out that such superradiant phenomena are transient in nature with a characteristic time scale equal to the GW oscillation time scale. We further point out that the existing ground-based Low Frequency Array (LOFAR), radio telescopes, may be able to detect such transient signals from BHs of mass range $M\sim 10^{-1} - 10^{-2} M_{\odot}$, which should necessarily be of primordial origin. Our present result opens up an intriguing possibility of observing the black hole merging phenomena through electromagnetic waves.

It is now experimentally established that neutrino has finite mass, and investigating its cosmological implication is of significant importance also for fundamental physics. In particular, as a matter component, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and conversely, observations of the LSS can give constraints on the neutrino mass. Performing large numerical simulations of the LSS formation including massive neutrino along with conventional cold dark matter is thus an important task. For this, calculating the neutrino distribution in the phase space by solving the Vlasov equation is a more suitable approach than conventional $N$-body simulations, but it requires solving the PDE in the $(6+1)$-dimensional space and is thus computationally demanding: configuring $n_\mathrm{gr}$ grid points in each coordinate and $n_t$ time grid points leads to $O(n_tn_\mathrm{gr}^6)$ complexity. We propose an efficient quantum algorithm for this task. Linearizing the Vlasov equation by neglecting the relatively weak self-gravity of the neutrino, we perform the Hamiltonian simulation to produce quantum states that encode the phase space distribution of neutrino. We also propose a way to extract the power spectrum of the neutrino density perturbations as classical data from the quantum state by quantum amplitude estimation with accuracy $\epsilon$ and query complexity of order $\widetilde{O}((n_\mathrm{gr} + n_t)/\epsilon)$. As far as we know, this is the first quantum algorithm for the LSS simulation that outputs the quantity of practical interest with guaranteed accuracy.

Hierarchical Bayesian inference can simultaneously account for both measurement uncertainty and selection effects within astronomical catalogs. In particular, the hierarchy imposed encodes beliefs about the interdependence of the physical processes that generate the observed data. We show that several proposed approximations within the literature actually correspond to inferences that are incompatible with any physical detection process, which can be described by a directed acyclic graph (DAG). This generically leads to biases and is associated with the assumption that detectability is independent of the observed data given the true source parameters. We show several examples of how this error can affect astrophysical inferences based on catalogs of coalescing binaries observed through gravitational waves, including misestimating the redshift evolution of the merger rate as well as incorrectly inferring that General Relativity is the correct theory of gravity when it is not. In general, one cannot directly fit for the ``detected distribution'' and ``divide out'' the selection effects in post-processing. Similarly, when comparing theoretical predictions to observations, it is better to simulate detected data (including both measurement noise and selection effects) rather than comparing estimates of the detected distributions of event parameters (which include only selection effects). While the biases introduced by model misspecification from incorrect assumptions may be smaller than statistical uncertainty for moderate catalog sizes (O(100) events), they will nevertheless pose a significant barrier to precision measurements of astrophysical populations.

The even light QCD axion called the $Z_{\mathcal N}$ axion can both solve the strong CP problem and account for the dark matter (DM). We point out that the double level crossings can naturally take place in the mass mixing between the $Z_{\mathcal N}$ axion and axionlike particle (ALP). The first level crossing occurs much earlier than the QCD phase transition, while the second level crossing occurs exactly during the QCD phase transition if it exists. The $Z_{\mathcal N}$ axion DM relic density can be suppressed or enhanced in the single level crossing, and suppressed in the double level crossings, depending on the onset of axion oscillations. The future axion experiments can detect these level crossing effects.

First-order phase transitions in the early universe might produce a detectable background of gravitational waves. As these phase transitions can be generated by new physics, it is important to quantify these effects. Many pure Yang-Mills gauge theories are known to undergo first-order deconfinement phase transitions, with properties that can be studied with lattice simulations. Despite the recent surge of interest in $Sp(2N)$ gauge theories as a candidate for models of physics beyond the standard model, studies of these theories at finite temperature are still very limited. In this contribution we will present preliminary results of an ongoing numerical investigation of the thermodynamic properties of the deconfinement phase transition in $Sp(4)$ Yang-Mills theory, using the linear logarithmic relaxation algorithm. This method enables us to obtain a highly accurate determination of the density of states, allowing for a precise reconstruction of thermodynamic observables. In particular, it gives access to otherwise difficult to determine quantities such as the free energy of the system, even along metastable and unstable branches, hence providing an additional direct observable to study the dynamics of the phase transition.