Papers for Monday, Dec 04 2023

A critically important process affecting the climate evolution and potential habitability of an exoplanet is atmospheric escape, in which high-energy radiation from a star drives the escape of hydrogen atoms and other light elements from a planet's atmosphere. L 98-59 is a benchmark system for studying such atmospheric processes, with three transiting terrestrial-size planets receiving Venus-like instellations (4-25 S$_\oplus$) from their M3 host star. We use the VPLanet model to simulate the evolution of the L 98-59 system and the atmospheric escape of its inner three small planets, given different assumed initial water quantities. We find that, regardless of their initial water content, all three planets accumulate significant quantities of oxygen due to efficient water photolysis and hydrogen loss. All three planets also receive enough XUV flux to drive rapid water loss, which considerably affects their developing climates and atmospheres. Even in scenarios of low initial water content, our results suggest that the James Webb Space Telescope (JWST) will be sensitive to observations of retained oxygen on the L 98-59 planets in its future scheduled observations, with planets b and c being the most likely targets to possess an extended atmosphere. Our results constrain the atmospheric evolution of these small rocky planets, and they provide context for current and future observations of the L 98-59 system to generalize our understanding of multi-terrestrial planet systems.

According to a standard initial mass function, stars in the range 7-12 Msun constitute ~50% (by number) of the stars more massive than ~7 Msun, but, in spite of this, their evolutionary properties, and in particular their final fate, are still scarcely studied. In this paper we present a detailed study of the evolutionary properties of solar metallicity, non rotating stars in the range 7-15 Msun, from the pre main sequence phase up to the presupernova stage or up to an advanced stage of the thermally pulsing phase, depending on the initial mass. We find that (1) the 7.00 Msun develops a degenerate CO core and evolves as a classical AGB star in the sense that it does not ignite the C burning reactions; (2) stars with the initial mass M >= 9.22 Msun end their life as core collapse supernovae; (3) stars in the range 7.50 <= M/Msun <= 9.20 develop a degenerate ONeMg core and evolve through the thermally pulsing SAGB phase; 4) stars in the mass range 7.50 <= M/Msun <= 8.00 end their life as hybrid CO/ONeMg- or ONeMg- WD; (5) stars with the initial mass in the range 8.50 <= M/Msun <= 9.20 may potentially explode as electron capture supernovae.

The intergalactic helium became fully ionized by the end of cosmic noon ($z\sim2$). Similarly to the reionization of hydrogen, helium reionization is expected to be patchy, driven by luminous quasars that ionize the intergalactic gas in their surrounding environment. Probing the morphology of ionized electrons during this epoch can provide crucial information about early structure formation, including the clustering and luminosities of quasars, the accretion rates, variability, and lifetimes of active galactic nuclei, as well as the growth and evolution of supermassive black holes. In this study, we present how measurements of the cosmic microwave background (CMB) can be used to reconstruct the optical-depth fluctuations resulting from patchy helium reionization. As helium reionization occurred at lower redshifts, upcoming probes of large-scale structure surveys will present a significant opportunity to enhance the prospects of probing this epoch by their combined analysis with the CMB. Using a joint information-matrix analysis of hydrogen and helium reionization, we show that near-future galaxy and CMB surveys will have enough statistical power to detect optical-depth fluctuations due to doubly-ionized helium, providing a way of measuring the redshift and duration of helium reionization to high significance. We also show that modeling uncertainties in helium reionization can impact the measurement precision of parameters characterizing hydrogen reionization.

The solar physics community is entering a golden era that is ripe with next-generation ground- and space-based facilities. With ever-increasing resolving power stemming from the newest observational telescopes, it becomes more challenging to obtain (near-)simultaneous measurements at high spatial, temporal and spectral resolutions, while operating at the diffraction limit of these new facilities. Hence, in recent years there has been increased interest in the capabilities integral field units (IFUs) offer towards obtaining the trifecta of high spatial, temporal and spectral resolutions contemporaneously. To date, IFUs developed for solar physics research have focused on mid-optical and infrared measurements. Here, we present an IFU prototype that has been designed for operation within the near-ultraviolet to mid-optical wavelength range, hence providing additional spectral coverage to the instrument suites developed to date. The IFU was constructed as a low-budget proof-of-concept for the upcoming 2m class Indian National Large Solar Telescope and employs circular cross-section fibres to guide light into a Czerny-Turner configuration spectrograph, with the resulting spectra captured using a high quantum efficiency scientific CMOS camera. Mapping of each input fibre allows for the reconstruction of two-dimensional spectral images, with frame rates exceeding 20 per second possible while operating in a non-polarimetric configuration. The science verification data presented here highlights the suitability of fibre-fed IFUs operating at near-ultraviolet wavelengths for solar physics research. Importantly, the successful demonstration of this type of instrument paves the way for further technological developments to make a future variant suitable for upcoming ground-based and space-borne telescope facilities.

Narrow Fe K$\alpha$ fluorescent emission lines arising at $\sim$kpc-scale separations from the nucleus have only been detected in a few AGN. The detections require that the extended line emission be spatially resolved and sufficiently bright. Compared to narrow Fe K$\alpha$ lines arising closer to the nucleus, they have much lower fluxes but show substantially larger equivalent widths, EW$_{\rm Fe K\alpha}$. We show that, in the optically-thin limit, a purely analytical argument naturally predicts large, EW$_{\rm FeK\alpha}\sim$1 keV, values for such lines, regardless of the details of equivalent hydrogen column density, $N_H$, or reprocessor geometry. Monte Carlo simulations corroborate this result and show that the simple analytic EW$_{\rm FeK\alpha}$ prescription holds up to higher $N_H$ approaching the Compton-thick regime. We compare to $Chandra$ observations from the literature and discuss that our results are consistent with the large EW$_{\rm FeK\alpha}$ values reported for local AGN, for which the line is detected in extended, up to $\sim$kpc-scale, regions. We argue that large EW$_{\rm FeK\alpha}$ from kpc-scale regions in AGN should be ubiquitous, because they do not depend on the absolute luminosity of the central X-ray source, and are measured only against the scattered continuum. We predict values to be of the order of $\sim$1 keV or larger, even for covering factors $\ll$1, and for arbitrarily small column densities. We propose that the large-scale molecular material that is now routinely being detected with the Atacama Large Millimeter/Submillimeter Array (ALMA) may act as an extended X-ray scattering reprocessor giving rise to $\sim$kpc-scale Fe K$\alpha$ emission.

Next-generation radio interferometers like the Square Kilometer Array have the potential to unlock scientific discoveries thanks to their unprecedented angular resolution and sensitivity. One key to unlocking their potential resides in handling the deluge and complexity of incoming data. This challenge requires building radio interferometric imaging methods that can cope with the massive data sizes and provide high-quality image reconstructions with uncertainty quantification (UQ). This work proposes a method coined QuantifAI to address UQ in radio-interferometric imaging with data-driven (learned) priors for high-dimensional settings. Our model, rooted in the Bayesian framework, uses a physically motivated model for the likelihood. The model exploits a data-driven convex prior, which can encode complex information learned implicitly from simulations and guarantee the log-concavity of the posterior. We leverage probability concentration phenomena of high-dimensional log-concave posteriors that let us obtain information about the posterior, avoiding MCMC sampling techniques. We rely on convex optimisation methods to compute the MAP estimation, which is known to be faster and better scale with dimension than MCMC sampling strategies. Our method allows us to compute local credible intervals, i.e., Bayesian error bars, and perform hypothesis testing of structure on the reconstructed image. In addition, we propose a novel blazing-fast method to compute pixel-wise uncertainties at different scales. We demonstrate our method by reconstructing radio-interferometric images in a simulated setting and carrying out fast and scalable UQ, which we validate with MCMC sampling. Our method shows an improved image quality and more meaningful uncertainties than the benchmark method based on a sparsity-promoting prior. QuantifAI's source code: https://github.com/astro-informatics/QuantifAI.

We present here broadband radio observations of the CSM interacting SN2019oys. SN2019oys was first detected in the optical and was classified as a Type Ib SN. Then, about $\sim 100$ days after discovery, it showed an optical rebrightening and a spectral transition to a spectrum dominated by strong narrow emission lines, which suggests strong interaction with a distant, dense, CSM shell. We modeled the broadband, multi-epoch, radio spectra, covering 2.2 to 36 GHz and spanning from 22 to 1425 days after optical discovery, as a synchrotron emitting source. Using this modeling we characterized the shockwave and the mass-loss rate of the progenitor. Our broadband radio observations show strong synchrotron emission. This emission, as observed 201 and 221 days after optical discovery, exhibits signs of free-free absorption from the material in front of the shock traveling in the CSM. In addition, the steep power law of the optically thin regime points towards synchrotron cooling of the radiating electrons. Analyzing these spectra in the context of the SN-CSM interaction model gives a shock velocity of 14,000 $\rm km \, s^{-1}$, and an electron number density of $2.6 \times 10^5 \, \rm cm^{-3}$ at a distance of $2.6 \times 10^{16}$ cm. This translates to a high mass-loss rate from the progenitor massive star of $6.7 \times 10^{-4} \, \rm M_{\odot} yr^{-1}$ for an assumed wind of 100 $\rm km s^{-1}$ (assuming constant mass-loss rate in steady winds). The late-time radio spectra, 392 and 557 days after optical discovery, are showing broad spectral peaks. We show that this can be explained by introducing a non-homogeneous CSM structure.

Optical surveys have become increasingly adept at identifying candidate Tidal Disruption Events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. We here present $\texttt{tdescore}$, a simple photometric classifier that is trained using a systematic census of $\sim$3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing $<$2% of the total. $\texttt{tdescore}$ is nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with completeness of 77.0% and a purity of 80.3%. $\texttt{tdescore}$ is thus substantially better than any available TDE photometric classifier scheme in the literature, and performs comparably well to the single-epoch spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multi-wavelength catalogue crossmatching. In a novel extension, we use 'SHapley Additive exPlanations' (SHAP) to provide a human-readable justification for each individual $\texttt{tdescore}$ classification, enabling users to understand and form opinions about the underlying classifier reasoning. $\texttt{tdescore}$ serves as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory.

We present high-resolution $K$-band emission spectra of the quintessential hot Jupiter HD 189733 b from the Keck Planet Imager and Characterizer (KPIC). Using a Bayesian retrieval framework, we fit the dayside pressure-temperature profile, orbital kinematics, mass-mixing ratios of H$_2$O, CO, CH$_4$, NH$_3$, HCN, and H$_2$S, and the $\rm ^{13}CO/^{12}CO$ ratio. We measure mass fractions of $\rm \log H_2O = -2.0^{+0.4}_{-0.4}$ and $\rm \log CO = -2.2^{+0.5}_{-0.5}$, and place upper limits on the remaining species. Notably, we find $\rm \log CH_4 < -4.5$ at 99\% confidence, despite its anticipated presence at the equilibrium temperature of HD 189733 b assuming local thermal equilibrium. We make a tentative ($\sim3\sigma$) detection of $\rm ^{13}CO$, and the retrieved posteriors suggest a $\rm ^{12}C/^{13}C$ ratio similar to or substantially less than the local interstellar value. The possible $\rm ^{13}C$ enrichment would be consistent with accretion of fractionated material in ices or in the protoplanetary disk midplane. The retrieved abundances correspond to a substantially sub-stellar atmospheric $\rm C/O = 0.3\pm0.1$, while the carbon and oxygen abundances are stellar to slightly super-stellar, consistent with core-accretion models which predict an inverse correlation between C/O and metallicity. The specific combination of low C/O and high metallicity suggests significant accretion of solid material may have occurred late in the formation process of HD 189733 b.

The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper reports the results of the ZTF Source Classification Project (SCoPe), which trains neural network and XGBoost machine learning (ML) algorithms to perform dichotomous classification of variable ZTF sources using a manually constructed training set containing 170,632 light curves. We find that several classifiers achieve high precision and recall scores, suggesting the reliability of their predictions for 112,476,749 light curves across 40 ZTF fields. We also identify the most important features for XGB classification and compare the performance of the two ML algorithms, finding a pattern of higher precision among XGB classifiers. The resulting classification catalog is available to the public, and the software developed for SCoPe is open-source and adaptable to future time-domain surveys.

We study the statistics of dwarf galaxy populations as a function of environment in the cold dark matter (CDM) and warm dark matter (WDM) cosmogonies, using hydrodynamical simulations starting from initial conditions with matched phases but differing power spectra, and evolved with the EAGLE galaxy formation model. We measure the abundance of dwarf galaxies within 3~Mpc of DM haloes with a present-day halo mass similar to that of the Milky Way (MW), and find that the radial distribution of galaxies $M_{*}>10^7$\Msun is nearly identical for WDM and CDM. However, the cumulative mass function becomes shallower for WDM at lower masses, yielding 50~per~cent fewer dwarf galaxies of $M_{*}\gtrsim10^{5}$~\Msun than CDM. The suppression of low-mass halo counts in WDM relative to CDM increases significantly from high-density regions to low-density regions for haloes in the region of the half-mode mass, $M_\rm{hm}$. The luminous fraction in the two models also diverges from the overdense to the underdense regions for $M>2M_\rm{hm}$, as the increased collapse delay at small densities pushes the collapse to after the reionization threshold. However, the stellar mass--halo mass relation of WDM haloes relative to CDM increases towards lower-density regions. Finally, we conclude that the suppression of galaxies with $M_{*}\gtrsim10^5$\Msun between WDM and CDM is independent of density: the suppression of halo counts and the luminous fraction is balanced by an enhancement in stellar mass--halo mass relation.

Photonic lanterns (PLs) are tapered waveguides that gradually transition from a multi-mode fiber geometry to a bundle of single-mode fibers. In astronomical applications, PLs can efficiently couple multi-mode telescope light into a multi-mode fiber entrance and convert it into multiple single-mode beams. The output beams are highly stable and suitable for feeding into high-resolution spectrographs or photonic chip beam combiners. For instance, by using relative intensities in the output cores as a function of wavelength, PLs can enable spectroastrometry. In addition, by interfering beams in the output cores with a beam combiner in the backend, PLs can be used for high-throughput interferometric imaging. When used on an Extremely Large Telescope (ELT), with its increased sensitivity and angular resolution, the imaging and spectroastrometric capabilities of PLs will be extended to higher contrast and smaller angular scales. We study the potential spectroastrometry and imaging science cases of PLs on ELTs, including study of exomoons, broad-line regions of quasars, and inner circumstellar disks.

Blue supergiants (BSGs) are key objects for understanding the evolution of massive stars. However, discrepancies between theoretical predictions and empirical observations have opened up important questions yet to be answered. Studying statistically significant and unbiased samples of BSGs can help to improve the situation. We aim to perform a homogeneous and comprehensive quantitative spectroscopic analysis of a large sample of Galactic BSGs from the IACOB spectroscopic database. We derive the projected rotational velocity ($v\sin i$) and macroturbulent broadening ($v_{\rm mac}$) using IACOB-BROAD. We used FASTWIND computations to derive effective temperatures ($T_{\rm eff}$), surface gravities, microturbulences ($\xi$), Si and He surface abundances, and the wind-strength parameter. We provide estimates of the above-mentioned quantities for the largest sample of Galactic BSGs spectroscopically analyzed to date, comprising 538 O9-B5 type stars. We find a drop in the relative number of BSGs at ~21 kK, coinciding with a scarcity of fast rotating stars below that temperature. We speculate that this feature might be roughly delineating the location of the empirical Terminal-Age-Main-Sequence in the 15-85Msol range. By investigating the main characteristics of the $v\sin i$ distribution of O stars and BSGs as a function of $T_{\rm eff}$, we propose that an efficient mechanism transporting angular momentum from the stellar core to the surface might be operating along the main sequence. We find correlations between $\xi$, $v_{\rm mac}$ and the spectroscopic luminosity. We also find that no more than 20% of the BSGs have atmospheres enriched in He, and suggest that the origin of this specific sub-sample of BSGs might be in binary evolution. We do not find clear empirical evidence of an increase in the wind-strength over the wind bi-stability region towards lower $T_{\rm eff}$.

We present photometry and spectroscopy of the slowly evolving superluminous Type IIn SN2015da. SN2015da is extraordinary for its very high peak luminosity, and also for sustaining a high luminosity for several years. Even at 8\,yr after explosion, SN2015da remains as luminous as the peak of a normal SNII-P. The total radiated energy integrated over this time period (with no bolometric correction) is at least 1.6 FOE. Including a mild bolometric correction, adding kinetic energy of the expanding cold dense shell of swept-up circumstellar material (CSM), and accounting for asymmetry, the total explosion kinetic energy was likely 5-10 FOE. Powering the light curve with CSM interaction requires an energetic explosion and 20 Msun of H-rich CSM, which in turn implies a massive progenitor system above 30 Msun. Narrow P Cyg features show steady CSM expansion at 90 km/s, requiring a high average mass-loss rate of roughly 0.1 Msun/yr sustained for 2 centuries before explosion (although ramping up toward explosion time). No current theoretical model for single-star pre-SN mass loss can account for this. The slow CSM, combined with broad wings of H$\alpha$ indicating H-rich material in the unshocked ejecta, disfavor a pulsational pair instability model for the pre-SN mass loss. Instead, violent pre-SN binary interaction is a likely cuprit. Finally, SN2015da exhibits the characteristic asymmetric blueshift in its emission lines from shortly after peak until the present epoch, adding another well-studied superluminous SNeIIn with unambiguous evidence of post-shock dust formation.

The Scalar Field Dark Matter model has been known in various ways throughout its history; Fuzzy, BEC, Wave, Ultralight, Axion-like Dark Matter, etc. All of them consist in proposing that the dark matter of the universe is a spinless field $\Phi$ that follows the Klein-Gordon (KG) equation of motion $\Box\Phi-dV/d\Phi=0$, for a given scalar field potential $V$. The difference between different models is sometimes the choice of the scalar field potential $V$. In the literature we find that people usually work in the nonrelativistic, weak-field limit of the KG equation where it transforms into the Schr\"odinger equation and the Einstein equations into the Poisson equation, reducing the KG-Einstein system, to the Schr\"odinger-Poisson system. In this paper, we review some of the most interesting achievements of this model from the historical point of view and its comparison with observations, showing that this model could be the last answer to the question about the nature of dark matter in the universe.

Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source within the error region of the $\gamma$-ray source {4FGL}~J1646.5$-$4406. The variability is consistent with the eclipse of a compact, steep-spectrum source behind ablated material from a companion in a $\sim 5.3\,$h binary orbit. Based on the eclipse properties and spatial coincidence with {4FGL} J1646.5$-$4406, we argue that the source is likely a recycled pulsar in a "redback" binary system. Using properties of the eclipses from ASKAP and Murchison Widefield Array observations, we provide broad constraints on the properties of the eclipse medium. We identified a potential optical/infra-red counterpart in archival data consistent with a variable low-mass star. Using the Parkes Radio Telescope "Murriyang" and MeerKAT, we searched extensively for radio pulsations but yielded no viable detections of pulsed emission. We suggest that the non-detection of pulses is due to scattering in the intra-binary material, but scattering from the ISM can also plausibly explain the pulse non-detections if the interstellar dispersion measure exceeds $\sim$600$\,$pc$\,$cm$^{-3}$. Orbital constraints derived from optical observations of the counterpart would be highly valuable for future $\gamma$-ray pulsation searches, which may confirm the source nature as a pulsar.

We investigate the effect of radial anisotropy on the rate of orbit decay in parabolic encounters of identical spherical galaxies. Our galaxy models have Hernquist density profiles and Osipkov--Merritt velocity distributions. We find that radially anisotropic models merge in as little as half the time of their isotropic counterparts. Anisotropic models are more susceptible to tidal deformation; this accelerates the transfer of orbital angular momentum to internal degrees of freedom. Even during the initial approach, the anisotropic models become more distorted, and arrive at pericentre already having lost substantial amounts of angular momentum. Our results may have implications for estimates of merger rates and persistence of tidal tails.

This paper presents analytic results on the anatomy of nested sampling, from which a technique is developed to estimate the run-time of the algorithm that works for any nested sampling implementation. We test these methods on both toy models and true cosmological nested sampling runs. The method gives an order-of-magnitude prediction of the end point at all times, forecasting the true endpoint within standard error around the halfway point.

Extremely metal-poor galaxies (XMPGs) at relatively low redshift are excellent laboratories for studying galaxy formation and evolution in the early universe. Much effort has been spent on identifying them from large-scale spectroscopic surveys or spectroscopic follow-up observations. Previous work has identified a few hundred XMPGs. In this work, we obtain a large sample of 223 XMPGs at $z<1$ from the early data of the Dark Energy Spectroscopic Instrument (DESI). The oxygen abundance is determined using the direct $T_{\rm e}$ method based on the detection of the [O III]$\lambda$4363 line. The sample includes 95 confirmed XMPGs based on the oxygen abundance uncertainty; remaining 128 galaxies are regarded as XMPG candidates. These XMPGs are only 0.01% of the total DESI observed galaxies. Their coordinates and other proprieties are provided in the paper. The most XMPG has an oxygen abundance of $\sim 1/34 Z_{\odot}$, stellar mass of about $1.5\times10^7 M_{\odot}$ and star formation rate of 0.22 $M_{\odot}$ yr$^{-1}$. The two most XMPGs present distinct morphologies suggesting different formation mechanisms. The local environmental investigation shows that XMPGs preferentially reside in relatively low-density regions. Many of them fall below the stellar mass-metallicity relations (MZRs) of normal star-forming galaxies. From a comparison of the MZR with theoretical simulations, it appears that XMPGs are good analogs to high-redshift star-forming galaxies. The nature of these XMPG populations will be further investigated in detail with larger and more complete samples from the on-going DESI survey.

Creating radio galaxy catalogues from next-generation deep surveys requires automated identification of associated components of extended sources and their corresponding infrared hosts. In this paper, we introduce RadioGalaxyNET, a multimodal dataset, and a suite of novel computer vision algorithms designed to automate the detection and localization of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4,155 instances of galaxies in 2,800 images with both radio and infrared channels. Each instance provides information about the extended radio galaxy class, its corresponding bounding box encompassing all components, the pixel-level segmentation mask, and the keypoint position of its corresponding infrared host galaxy. RadioGalaxyNET is the first dataset to include images from the highly sensitive Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope, corresponding infrared images, and instance-level annotations for galaxy detection. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to simultaneously detect radio galaxies and the positions of infrared hosts.

LAMOST J024048.51+195226.9 (J0240+1952) was recently identified as the second AE Aquarii (AE Aqr)-type cataclysmic variable, possessing the fastest known rotating white dwarf. We performed a Very Long Baseline Interferometry (VLBI) observation of J0240+1952 utilizing the European VLBI Network at 1.7\,GHz, to obtain the first view of the radio morphology on mas scale. Our high-resolution VLBI image clearly shows that the radio emission is compact on mas scale ($\lesssim2$\,AU), with no evidence for a radio jet or extended emission. The compact radio source has an average flux density of $\sim0.37$\,mJy, and its brightness temperature is given at $\gtrsim2.3\times10^{7}$\,K, confirming a non-thermal origin. The emission exhibits irregular variations on a time-scale of tens of minutes, similar to the radio flares seen in AE Aqr. The measured VLBI position of J0240+1952 is consistent with that derived from \textit{Gaia}. Our results favour the model in which the radio emission is attributed to a superposition of synchrotron radiation from expanding magnetized blobs of this system.

Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $\rho \propto r^{-1}$ and of $\rho \propto r^{-1.5}$. The galaxy shape noise is modeled based on the Year-1 Subaru Hyper Suprime-Cam (HSC) Survey. After reviewing mathematical details of our algorithm and dark matter halo models, we determine an optimal value of the coefficient of the adaptive LASSO regression penalty term for single halo reconstruction. We successfully measure halo masses for massive triaxial halos; the mass determination accuracy is 5 percent for halos with $M = 10^{14.6}~M_\odot$ at $0.0625\leq z \leq 0.2425$, and 5 percent for those with $10^{14.8}~M_\odot$ at $0.0625\leq z \leq 0.4675$, and 20 percent for $M= 10^{15.0} ~M_\odot$ and $M=10^{15.2}~M_\odot$ in the redshift range $0.0625\leq z \leq 0.4675$. The redshift estimate accuracy is consistently below $\Delta z /z \leq 0.05$ for the above halo masses in the range $0.1525\leq z \leq 0.4675$. We also demonstrate that the orientation of triaxial halos and systematic error in our halo model do not affect reconstruction result significantly. Finally, we present results from reconstruction of mass distribution using shear catalogs produced by multiple halos, to show \splinv{}'s capability using realistic shear maps from ongoing and future galaxy surveys.

The Full-disk magnetograph is a main scientific payload onboard the Advanced Space based Solar Observatory (ASO-S/FMG) that through Stokes parameter observation to measures the vector magnetic field. The accuracy of magnetic-field values is an important aspect of checking the quality of the FMG magnetic-field measurement. According to the design of the FMG, the linear calibration method under the weak-field approximation is the preferred scheme for magnetic-field calibration. However, the spacecraft orbital velocity can affect the position of observed spectral lines, then result in a change of the polarization-signal strength. Thus, the magnetic field is modulated by the orbit velocity of the spacecraft. In this article, through cross calibration between FMG and HMI (Helioseismic and Magnetic Imager onboard the Solar Dynamic Observatory), the effects of spacecraft orbital velocity on the coefficient of magnetic-field calibration are investigated. By comparing the magnetic field of FMG and HMI with spacecraft orbital velocity as an auxiliary reference, the revised linear-calibration coefficients that depend on spacecraft orbital velocity are obtained. Magnetic field of FMG corrected by the revised calibration coefficients removing the effect of spacecraft orbital velocity will be more accurate and suitable for scientific research.

The probability that the primordial black hole (PBH) binaries formed in the early Universe can be affected by the Hubble expansion of background, which is non-negligible when the number density of PBHs is very low (it is actually this case for supermassive PBHs). In this paper, taking into account the effect of cosmic expansion on the comoving distance of PBH pairs, we worked out the merger rate of PBHs with any extended mass function. The torques by all PBHs and linear density perturbations are also considered. It is found that the merger rate of PBH, $M\gtrsim 10^6M_\odot$, binaries is significantly lower for $f_\text{pbh}\lesssim 0.01$ than expected.

We present a new total intensity image of M31 at 1.248 GHz, observed with the Five-hundred-meter Aperture Spherical radio telescope (FAST) with an angular resolution of 4 arcmin and a sensitivity of about 16 mK. The new FAST image clearly reveals weak emission outside the ring due to its high sensitivity on large-scale structures. We derive a scale length of 2.7 kpc for the cosmic ray electrons and find that the cosmic ray electrons propagate mainly through diffusion by comparing the scale length at 4.8 GHz. The spectral index of the total intensity varies along the ring, which can be attributed to the variation of the spectra of synchrotron emission. This variation is likely caused by the change of star formation rates along the ring. We find that the azimuthal profile of the non-thermal emission can be interpreted by an axisymmetric large-scale magnetic field with varying pitch angle along the ring, indicating a complicated magnetic field configuration in M31.

In contrast to our understanding of density field tracers, the modelling of direct statistics pertaining to the cosmic velocity field remains open to significant opportunities for improvement. The lack of accurate modelling for the non-linear domain of pairwise velocities restricts our capacity to fully exploit the information encoded in this observable. We present a robust approach for modelling the mean infall velocities, $v_{12}(r,a)$, with broad applicability spanning sub-megaparsec scales and cosmologies extending beyond the standard LCDM paradigm. Our approach involves solving the full pair-conservation equation using accurate non-linear power spectrum descriptions. To assess the robustness of our model, we extend it to cosmologies beyond the standard LCDM, in particular, the Hu-Sawicki $f(R)$-gravity and Dvali-Gabadadze-Porrati (DGP) modified gravity models. Remarkably, our predictions for pairwise velocities of dark matter particles at kilo-parsec scales exhibit excellent agreement with N-body simulations throughout the entire dynamical range ($0.1 \lesssim \xi \lesssim 1000$, or $r\geq0.4$Mpc/h). Furthermore, we show that different gravity models leave distinct signatures in the shape and dynamics of the mean pairwise velocities, providing a potent test of cosmological gravity laws.

Some ideas for onboard processing and photometry with an astrometry satellite are presented, especially designed for GaiaNIR which may be launched about 2045 as a successor of Gaia. - Increased sensitivity, reduced image overlap, and simpler PSF calibration in GaiaNIR will result if the proposed initial processing of data from the detectors is implemented, because the across-scan smearing will become insignificant. - Filter photometry is required for high angular resolution as needed for astrometric and astrophysical reasons. Low-dispersion spectra are questioned because they fail at high star density. This will be a much greater problem with GaiaNIR than it is with Gaia because of the larger number of stars expected. It was the aim to collect in this note all arguments about GaiaNIR photometry which can be stated with words only, in correspondence with readers. The remaining work to be done for the definition of photometric equipment on the satellite requires further quantitative assessments and comparison of various options. Finally, 1) an advantage of filters is that the photometric observations can also be used for astrometry, 2) the XP spectra in Gaia will give very good astrophysical data for about 400 million single stars with G <~ 18.5 mag, but filters would have been better for all fainter and for all multiple stars, and 3) it is presently not clear which advantages for astrophysics low-dispersion spectra in the NIR might have over filters.

In July 2021, faculty from the UAlbany Department of Physics participated in a week-long field expedition with the organization UAPx to collect data on UAPs in Avalon, California, located on Catalina Island, and nearby. This paper reviews both the hardware and software techniques which this collaboration employed, and contains a frank discussion of the successes and failures, with a section about how to apply lessons learned to future expeditions. Both observable-light and infrared cameras were deployed, as well as sensors for other (non-EM) emissions. A pixel-subtraction method was augmented with other similarly simple methods to provide initial identification of objects in the sky and/or the sea crossing the cameras' fields of view. The first results will be presented based upon approximately one hour in total of triggered visible/night-vision-mode video and over 600 hours of untriggered (far) IR video recorded, as well as 55 hours of (background) radiation measurements. Following multiple explanatory resolutions of several ambiguities that were potentially anomalous at first, we focus on the primary remaining ambiguity captured at approximately 4am Pacific Time on Friday, July 16: a dark spot in the visible/near-IR camera possibly coincident with ionizing radiation that has thus far resisted a prosaic explanation. We conclude with quantitative suggestions for serious researchers in this still-nascent field of hard-science-based UAP studies, with an ultimate goal of identifying UAPs without confirmation bias toward either mundane or speculative conclusions.

Thanks to the mid-infrared capability offered by Spitzer, systematic searches of dust in SNe have been carried out over the past decade. Studies have revealed the presence of a substantial amount of dust over a broad range of SN subtypes. How normal SNe present mid-IR excess at later time and turn out to be dusty SNe can be affected by several factors, such as mass-loss history and envelope structure of progenitors and their explosion environment. All these can be combined and related to their environmental properties. A systematic analysis of SNe that exploded under a dusty environment could be of critical importance to measure the properties of the dust-veiled exploding stars, and whether such an intense dust production process is associated with the local environment. In this work, we firstly use the IFS data to study the environmental properties of dusty SNe compared to those of normal ones, and analyze correlations between the environmental properties and their dust parameters. We find that dusty SNe have a larger proportion located at higher SFR regions compared to the normal types. The occurrence of dusty SNe is less dependent on metallicity, with the oxygen abundance spanning from subsolar to oversolar metallicity. We also find the host extinction of dusty SNe scatters a lot, with about 40% of dusty SN located at extremely low extinction environments, and another 30% of them with considerably high host extinction of E(B-V)>0.6 mag.

We investigate the stochastic gravitational waves background arising from the first-order QCD chiral phase transition, considering three distinct sources: bubble collisions, sound waves, and fluid turbulence. Within the framework of the Polyakov-Nambu-Jona-Lasinio (PNJL) model, we calculate the parameters governing the intensity of the gravitational wave phase transition and determine their magnitudes along the adiabatic evolutionary path. We introduce the effective bag constant $B_{\mathrm{eff}}$ related to the dynamical evolution of quarks to evaluate the intensity of the phase transition. By calculating expanded potential at the point of false vacuum, we find that all the bubbles are in the mode of runaway, leading the velocity of the bubble wall to the speed of light. The resulting gravitational wave energy spectrum is estimated, revealing a characteristic amplitude of the generated gravitational waves within the centihertz frequency range. We present the gravitational wave spectrum and compare it with the sensitivity range of detectors, and find that the gravitational wave spectra generated by these sources have the potential to be detected by future detectors such as BBO and $\mu$ARES.

White dwarf (WD) stars are considered cosmic laboratories to study the physics of dense plasma. Furthermore, the use of WD stars as cosmic clocks to date stellar populations and main sequence companions demands an appropriate understanding of the WD physics in order to provide precise ages for these stars. We aim at studying exsolution in the interior of WD stars, a process in which a crystallized ionic binary mixture separates into two solid solutions with different fractions of the constituents. Depending on the parent solid mixture composition, this process can release or absorb heat, thus leading to a delay or a speed-up of WD cooling. Relying on accurate phase diagrams for exsolution, we have modeled this process in hydrogen-rich WDs with both carbon-oxygen and oxygen-neon core composition, with masses ranging from 0.53 to 1.29Msun and from 1.10 to 1.29Msun, respectively. Exsolution is a slow process that takes place at low luminosities (log(L/Lsun)$\lesssim$-2.75) and effective temperatures (Teff$\lesssim$18 000K) in WDs. We find that exsolution begins at brighter luminosities in CO than in ONe WDs of the same mass. Massive WDs undergo exsolution at brighter luminosities than their lower-mass counterparts. The net effect of exsolution on the WD cooling times depends on the stellar mass and the exact chemical profile. For standard core chemical profiles and preferred assumptions regarding miscibility gap microphysics, the cooling delay can be as large as ~0.35 Gyrs at L/Lsun ~ -5. We have neglected a chemical redistribution possibly associated with this process, which could lead to a further cooling delay. Exsolution has a marginal effect on the WD cooling times and, accordingly, we find no WD branches on the Gaia color magnitude diagram associated with it. However, exsolution in massive WDs can alter the faint end of the WD luminosity function, thus impacting WD cosmochronology.

We present a magnetic flux rope (FR) extraction tool for solar coronal magnetic field modelling data, which builds upon the methodology from Wagner et al. (2023). We apply the scheme to magnetic field simulations of active regions AR12473 and AR11176. We compare the method to its predecessor and study the 3D movement of the newly extracted FRs up to heights of 200 and 300 Mm, respectively. The extraction method is based on the twist parameter and a variety of mathematical morphology algorithms, including the opening transform and the morphological gradient. We highlight the differences between the methods by investigating the circularity of the FRs in the plane we extract from. The simulations for the active regions are carried out with a time-dependent data-driven magnetofrictional model (TMFM; Pomoell et al. (2019)). We investigate the FR trajectories by tracking their apex throughout the full simulation time span. We demonstrate that this upgraded methodology provides the user with more tools and less a-priori assumptions about the FR shape that, in turn, leads to a more accurate set of field lines. The propagation analysis yields that the erupting FR from AR12473 showcases stronger dynamics than the AR11176 FR and a significant deflection during its ascent through the domain. The AR11176 FR appears more stable, though there still is a notable deflection. This confirms that at these low coronal heights, FRs do undergo significant changes in the direction of their propagation even for less dynamic cases. The modelling results are also verified with observations, with AR12473 being indeed dynamic and eruptive, while AR11176 only features an eruption outside of our simulation time window.

We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the Flagship I N-body simulation at $z=(0.9,1.2,1.5,1.8)$, which have been populated with H$\alpha$ galaxies leading to catalogues of millions of objects within a volume of about $58\,h^{-3}\,{\rm Gpc}^3$. Our analysis suggests that both models can be used to provide a robust inference of the parameters $(h, \omega_{\rm c})$ in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber $k_{\rm max}=0.45\,h\,{\rm Mpc}^{-1}$, even with a measurement precision well below the percent level. In particular, this is true for a configuration with six free nuisance parameters, including local and non-local bias parameters, a matter counterterm, and a correction to the shot-noise contribution. Fixing either tidal bias parameters to physically-motivated relations still leads to unbiased cosmological constraints. We finally repeat our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, as purity and completeness, showing that we can get consistent cosmological constraints over this range of scales and redshifts.

We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from $LiteBIRD$ satellite simulations, which complement and expand previous studies in the literature. We find that $LiteBIRD$ will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the $TB$ and $EB$ angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of $LiteBIRD$ will reside in $BB$ angular power spectra rather than in $TB$ and $EB$ correlations.

Low Surface Brightness Galaxies (LSBGs) are excellent probes of quenching and other environmental processes near massive galaxies. We study an extensive sample of LSBGs near massive hosts in the local universe that are distributed across a diverse range of environments. The LSBGs with surface-brightness $\mu_{\rm eff,g}> $24.2 mag arcsec$^{-2}$ are drawn from the Dark Energy Survey Year 3 catalog while the hosts with masses $9.0< log(M_{\star}/M_{\odot})< 11.0$ comparable to the Milky Way and the Large Magellanic Cloud are selected from the z0MGS sample. We study the projected radial density profiles of LSBGs as a function of their color and surface brightness around hosts in both the rich Fornax-Eridanus cluster environment and the low-density field. We detect an overdensity with respect to the background density, out to 2.5 times the virial radius for both hosts in the cluster environment and the isolated field galaxies. When the LSBG sample is split by $g-i$ color or surface brightness $\mu_{\rm eff,g}$, we find the LSBGs closer to their hosts are significantly redder and brighter, like their high surface-brightness counterparts. The LSBGs form a clear 'red sequence' in both the cluster and isolated environments that is visible beyond the virial radius of the hosts. This suggests a pre-processing of infalling LSBGs and a quenched backsplash population around both host samples. However, the relative prominence of the 'blue cloud' feature implies that pre-processing is ongoing near the isolated hosts compared to the cluster hosts.

In this study, we explore the dynamical stability of magnetar bursts within the context of the chaos-randomness phase space for the first time, aiming to uncover unique behaviors compared to various astrophysical transients, including fast radio bursts (FRBs). We analyze burst energy time series data from active magnetar sources SGR J1550-5418 and SGR J1935+2154, focusing on burst arrival time and energy differences between consecutive events. We find a distinct separation in the time domain, where magnetar bursts exhibit significantly lower randomness compared to FRBs, solar flares, and earthquakes, with a slightly higher degree of chaos. In the energy domain, magnetar bursts exhibit a broad consistency with other phenomena, primarily due to the wide distribution of chaos-randomness observed across different bursts and sources. Intriguingly, contrary to expectations from the FRB-magnetar connection, the arrival time patterns of magnetar bursts in our analysis do not exhibit significant proximity to repeating FRBs in the chaos-randomness plane. This finding may challenge the hypothesis that FRBs are associated with typical magnetar bursts but indirectly supports the evidence that FRBs may primarily be linked to special magnetar bursts like peculiar X-ray bursts from SGR J1935+2154 observed simultaneously with Galactic FRB 200428.

We present the results of a joint {\it XMM-Newton} and {\it NuSTAR} campaign on the active galactic nucleus (AGN) IC 4329A, consisting of 9 $\times$ 20 ks {\it XMM-Newton} observations, and 5 $\times$ 20 ks {\it NuSTAR} observations within nine days, performed in August 2021. Within each observation, the AGN is not very variable, with the fractional variability never exceeding 5%. Flux variations are observed between the different observations, on timescales of days, with a 30% ratio between the minimum and the maximum 2-10 keV flux. These variations follow the softer-when-brighter behavior typically observed in AGN. In all observations, a soft excess is clearly present. Consistently with previous observations, the X-ray spectra of the source exhibit a cut-off energy between 140 and 250 keV, constant within the error in the different observations. We detected a narrow iron \ka line consistent with being constant during the monitoring, and likely originating in a distant neutral medium. We do not find evidence of a broad component of the iron line, suggesting that the inner disk does not produce strong reflection. We find that the reflection component is weak ($R_{\rm max}=0.009\pm0.002$). We also found the presence of a warm absorber component together with an ultra-fast outflow. Looking at their energetic, these outflows have enough mechanical power to exercise a significant feedback impact on the AGN surrounding environment.

We examine the production of the Hoyle and associated excited states from the viewpoint of pocket resonances in the reaction of an $\alpha$-particle on a ground state prolate $^8$Be nucleus within the optical model coupled-channel framework. The predicted reaction cross sections, as a function of the center-of-mass energy $E_{\rm cm}$, show prominent resonances, including the Hoyle resonance. The positions and widths of these resonances are sensitive to the target deformation ($\beta_2$ parameter) and the parity of the nuclear surface potential $-$ deeper for the even-parity $L$ partial waves relative to those for the odd-parity $L$ partial waves at the surface region because of the Bose-Einstein exchange of the $\alpha$-bosons. Decomposing the reaction cross sections to different partial waves, we find that the resonance energies and widths reasonably agree with the available experimental data and previous hyperspherical calculations for the $0_2^+$ (Hoyle state), $0_3^+$, $1_1^-$ and $3_1^-$ states of $^{12}$C, except for the narrow theoretical width of the $2_2^+$ state. Analyzing the wavefunctions and the resonance widths, we identify the narrow and sharp $0_2^+$, $3_1^-$ and $2_2^+$ resonances as pocket resonances -- resonances which occur below the potential barrier, while the broad $0_3^+$ and $1_1^-$ resonances as above-the-barrier resonances. For astrophysical applications, we also evaluate the astrophysical $S(E_{\rm cm})$-factor for $E_{\rm cm}$ $<$ 1.0 MeV, for the fusion of $\alpha$+$^8$Be into the $^{12}$C$(2^+)$ state based on our estimated $s$-wave $\alpha$+$^8$Be reaction cross section and the associated $\gamma$- and $\alpha$-decay widths for the decay of $^{12}$C excited states in the potential pocket.

Among the binary sources of interest for LISA some are quasimonochromatic, in the sense that the change in the gravitational wave frequency $\Delta f\lesssim 1\;\mbox{yr}^{-1}$ during the observation time. We study these sources in the frequency domain (FD) by taking into account the Doppler shift induced by LISA's motion and the LISA pattern functions, and we compare our results with previous work in the time domain. We also discuss the transition from the quasimonochromatic case to the stationary phase approximation commonly used in Fisher matrix calculations, which applies when $\Delta f\gtrsim 1\;\mbox{yr}^{-1}$.

Temporary matter domination and late entropy dilution, injected by a "long-lived" particle in the early universe, serves as a standard mechanism for yielding the correct dark matter relic density. We recently pointed out the cosmological significance of diluting particle's partial decay into dark matter. When repopulated in such a way, dark matter carries higher momentum than its thermal counterpart, resulting in a suppression of the linear matter power spectrum that is constrained by the large scale structure observations. In this work, we study the impact of such constraints on the minimal left-right symmetric model that accounts for the origin of neutrino mass. We map out a systematic anatomy of possible dilution scenarios with viable parameter spaces, allowed by cosmology and various astrophysical and terrestrial constraints. We show that to accommodate the observed dark matter relic abundance the spontaneous left-right symmetry breaking scale must be above PeV and cosmology will continue to provide the most sensitive probes of it. In case the dilutor is one of the heavier right-handed neutrinos, it can be much lighter and lie near the electroweak scale.

In this work, we made an extensive study about the possible presence of anisotropies in strange stars. To accomplish this task, we use three different configurations for the strange matter: the unpaired matter, a two-flavor super-conducting (2SC) strange matter, and a fully three-flavor super-conducting strange matter (CFL). For each configuration, we calculate the relevant quantities for the strange stars, such as the mass-radius relation, the dimensionless tidal parameter, the moment of inertia, and the surface curvature for different degrees of anisotropies. Whenever possible, we compare our results with constraints found in the literature, especially focusing on the existence of very massive pulsars (PSR J0952-0607), as well as very light compact objects (HESS J1731-347).

The unification of quantum mechanics and general relativity has long been elusive. Only recently have empirical predictions of various possible theories of quantum gravity been put to test. The dawn of multi-messenger high-energy astrophysics has been tremendously beneficial, as it allows us to study particles with much higher energies and travelling much longer distances than possible in terrestrial experiments, but more progress is needed on several fronts. A thorough appraisal of current strategies and experimental frameworks, regarding quantum gravity phenomenology, is provided here. Our aim is twofold: a description of tentative multimessenger explorations, plus a focus on future detection experiments. As the outlook of the network of researchers that formed through the COST Action CA18108 "Quantum gravity phenomenology in the multi-messenger approach (QG-MM)", in this work we give an overview of the desiderata that future theoretical frameworks, observational facilities, and data-sharing policies should satisfy in order to advance the cause of quantum gravity phenomenology.

The background of primordial gravitational waves (PGWs) predicted by the inflationary scenario induce incoherence of the electromagnetic field propagating over cosmological distances. We propose a new schema to constrain the underlying inflationary parameters, in particular the tensor-to-scalar ratio r, based on angular size-redshift \theta-z measurement accomplished by very long baseline interferometry (VLBI) surveys. VLBI observations rely on the van-Citter Zernike theorem, which expresses the coherence length in terms of its wavelength and angular size, i.e., d_coh < \lambda/\theta. In this study, we show that the interaction of the radio signal (involved in VLBI) with PGWs, along the propagation from a source located at redshift z_* to the Earth, leads to the blurring of the visibility. The blurring effect is evaluated for the highly-squeezed PGW, where it turns out that in order not to ruin the visibility, the projected baseline of the interferometer must be smaller than \xi_sq(z), which is inversely proportional to the the tensor-to-scalar ratio through \xi_sq \propto r^{-1/2}. Hence, VLBI observations based on interference pattern lead to a constraint on r imposed by the fact that \xi_sq(z) is greater than d_coh of the emitted radiation. In order to evaluate the constraint, we use a sample of compact radio quasars observed in VLBI and located at redshift range 0.46<z<2.73. We obtain a stringent upper-bound on the tensor-to-scalar ratio, r<2.10^{-6}, far beyond present constraints on this parameter. Further issues and caveats that potentially affect the results are reviewed. In particular, the possible effect of quantum-to-classical transition of PGWs is discussed. Ultimately, the background of primordial tensor perturbations may be more constrained with the help of the high-precision VLBI measurement of angular size-redshift of more distant sources.

The tidal acceleration experienced by an object at the event horizon of a black hole decreases as one over the square of the black hole's mass. As such there is a maximum mass at which a black hole can tidally disrupt an object outside of its event horizon and potentially produce observable emission. This maximum mass is known as the ``Hills mass'', and in full general relativity is a function of both the black hole's spin $a_\bullet$ and the inclination angle of the incoming object's orbit with respect to the black hole's spin axis $\psi$. In this paper we demonstrate that the Hills mass can be represented by a simple analytical function of $a_\bullet$ and $\psi$, the first general solution of this problem. This general solution is found by utilising the symmetries of a class of critical Kerr metric orbits known as the innermost bound spherical orbits. Interestingly, at fixed black hole spin the maximum Hills mass can lie at incoming orbital inclinations outside of the black hole's equatorial plane $\psi \neq \pi/2$. When compared to previous results in the literature this effect can lead to an increase in the maximum Hills mass (at fixed spin) by as much as a factor of $\sqrt{11/5} \simeq 1.48$ for a maximally rotating black hole. We then demonstrate how Bayesian inference, coupled with an estimate of the mass of a black hole in a tidal disruption event, can be used to place conservative constraints on that black hole's spin. We provide a publicly available code tidalspin which computes these spin distributions.

In this contribution to the proceedings of the 182nd Nobel Symposium, I reflect on the concept of "discovery" as it is used by physicists and astronomers. In particular, I comment on how the scientific community distinguishes discoveries from propositions that are supported only by lesser forms of evidence, emphasizing the social nature of this process and remarking on the subjective factors that go into making such judgements. I advocate for an approach that is intentionally Bayesian in nature, in which individuals are encouraged to evaluate and publicly state their priors and to update them systematically. I close by applying these practices to the case example of the Galactic Center Gamma-Ray Excess.

Is the Martian ionosphere symmetric in local time? Does the ionosphere change between aphelion, when there is no dust in the atmosphere, compared to perihelion, when the temperature rises and there is dust in the air even if no storm is blowing? We utilize data from the MAVEN Radio Occultation Science Experiment (Withers et al., 2020) - with unprecedented coverage in solar zenith angle - to answer these questions and isolate the effects that local time and season induce on the ionosphere at solar minimum. 219 out of the 1228 electron density profiles of the Martian ionosphere collected by MAVEN ROSE between July 2016 and December 2022 show, besides the ever-present dayside M2 layer, a distinct M1 layer right below. This allowed us to study the behavior of the M2 and M1 peak density and altitude as a function of solar zenith angle, and, also, for the first time, local time, and Martian season. We find that the M1 layer at low SZA can occur at altitudes lower than 100 km; that the peak altitudes and densities of both the M2 and M1 layers change more with season than at the dusk ionosphere; and that the M2 peak density decreases at a faster rate than the M1 with SZA. This study provides a baseline to accurately characterize the photoproduced Martian ionosphere at solar minimum.

The extreme sensitivity required for direct observation of gravitational waves by the Advanced LIGO detectors means that environmental noise can potentially contaminate gravitational wave signals. Consequently, environmental monitoring efforts have been undertaken and novel noise mitigation techniques have been developed which have helped keep environmental artifacts from influencing gravitational wave detections for the $90$ gravitational wave events detected from 2015--2020 by the aLIGO detectors. The increasing rate of gravitational wave detections due to detector sensitivity improvements requires sophisticated, reliable and automated ways to monitor and assess the degree of environmental coupling between gravitational wave detectors and their surroundings. We introduce a computational tool, PEMcheck, for quantifying the degree of environmental coupling present in gravitational wave signals using data from the network of environmental monitoring sensors. We study its performance when applied to the $79$ gravitational waves confidently detected in LIGO's third observing run and test its performance in the case of extreme environmental contamination of gravitational wave data. We find that PEMcheck's automated analysis identifies only a small number of gravitational waves that merit further study by environmental noise experts due to possible contamination, a substantial improvement over the manual vetting that occurred for every gravitational wave candidate in previous observing runs. Overall, PEMcheck works as intended. Consequently, PEMcheck will play a critical role in event validation during LIGO's fourth observing run.