Papers for Tuesday, Aug 08 2023

A two-body gas-phase reaction rate coefficient can be given by the usual Arrhenius-type formula which depends on temperature. The UMIST Database for Astrochemistry is a widely used database for reaction rate coefficients. They provide fittings for coefficients valid over a particular range of temperatures. The permissible upper-temperature limits vary over a wide range: from 100 K to 41000K. A wide range of temperatures occurs in nature; thus, it requires evaluating the rate coefficients at temperatures outside the range of validity. As a result, a simple extrapolation of the rate coefficients can lead to unphysically large values at high temperatures. These result in unrealistic predictions. Here we present a solution to prevent the gas-phase reaction coefficients from diverging at a very high temperature. We implement this into the spectral synthesis code CLOUDY which operates over a wide range of temperatures from CMB to 10$^{10}$ K subject to different astrophysical environments.

We present a study of the long term variability of Jupiter's mid-infrared auroral CH4 emissions. 7.7 - 7.9 micron images of Jupiter recorded by Earth-based telescopes over the last three decades were collated in order to quantify the magnitude and timescales over which the northern auroral hotspot's CH4 emissions varies. We find that the ratio of the radiance of the poleward northern auroral emissions to a lower-latitude zonal mean, henceforth 'Relative Poleward Radiance' or RPR, exhibits a 37% variability over a range of timescales. We searched for patterns of variability in order to test whether seasonally-varying solar insolation, the 11-year solar cycle, or short-term solar wind variability at Jupiter's magnetopause could explain the observed evolution. The variability of the RPR exhibits a weak (r < 0.2) correlation with the solar insolation received at Jupiter's high-northern latitudes. This rules out the hypothesis suggested in previous work (e.g. Sinclair et al., 2017a) that shortwave solar heating of aurorally-produced haze particles is the dominant heating mechanism in the lower stratosphere. We also find the variability exhibits negligible (r < 0.18) correlation with the monthly-mean sunspot number, which rules out variability associated with the solar cycle. On shorter timescales, we find moderate correlations of the RPR with solar wind conditions at Jupiter in the preceding days before images were recorded. For example, we find correlations of r = 0.45 and r = 0.51 of the RPR with the mean and standard deviation on the solar wind dynamical pressure in the preceding 7 days. The moderate correlation suggests that either: 1) only a subset of solar wind compressions lead to brighter, poleward, CH4 emissions and/or 2) a subset of CH4 emission brightening events are driven by internal magnetospheric and independent of the solar wind.

TeV halos have been suggested as a common phenomenon associated with middle-aged pulsars. Based on our recent work on PSR~J0631+1036, we select 3 candidate TeV halos from the first Large High Altitude Air Shower Observatory (LHAASO) catalog of gamma-ray sources. The corresponding pulsars, given by the positional coincidences and property similarities, are PSR J1958+2846, PSR J2028+3332, and PSR J1849$-$0001. We analyze the GeV $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard {\it the Fermi Gamma-ray Space Telescope} for the first two pulsars, as the last is gamma-ray quiet. We remove the pulsed emissions of the pulsars from the source regions from timing analysis, and determine that there are no residual GeV emissions in the regions as any possible counterparts to the TeV sources. Considering the previous observational results for the source regions and comparing the two pulsars to Geminga (and Monogem), the LHAASO-detected TeV sources are likely the pulsars' respective TeV halos. We find that the candidate and identified TeV halos, including that of PSR~J1849$-$0001, have luminosites (at 50\,TeV) approximately proportional to the spin-down energy $\dot{E}$ of the pulsars, and the ratios of the former to the latter are $\sim 6\times 10^{-4}$.

The upcoming generation of optical spectrographs on four meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide high-quality spectra for thousands of galaxies. These data will allow us to examine of the stellar population properties at intermediate redshift, an epoch that remains unexplored by large and deep surveys. We assess our capability to retrieve the mean stellar metallicity in galaxies at different redshifts and S/N, while simultaneously exploiting the UV and optical rest-frame wavelength coverage. The work is based on a comprehensive library of spectral templates of stellar populations, covering a wide range of age and metallicity values and built assuming various SFHs. We simulated realistic observations of a large sample of galaxies carried out with WEAVE at the WHT at different redshifts and S/N values. We measured all the reliable indices on the simulated spectra and on the comparison library. We then adopted a Bayesian approach to obtain the probability distribution of stellar metallicity. The analysis of the spectral indices has shown how some mid-UV indices can provide reliable constraints on stellar metallicity, along with optical indicators. The analysis of the mock observations has shown that even at S/N=10, the metallicity can be derived within 0.3 dex, in particular, for stellar populations older than 2 Gyr. Our results are in good agreement with other theoretical and observational works in the literature and show how the UV indicators can be advantageous in constraining metallicities. This is very promising for the upcoming surveys carried out with new, highly multiplexed, large-field spectrographs, such as StePS at the WEAVE and 4MOST, which will provide spectra of thousands of galaxies covering large spectral ranges at relatively high S/N.

The topology of the large-scale structure of the universe contains valuable information on the underlying cosmological parameters. While persistent homology can extract this topological information, the optimal method for parameter estimation from the tool remains an open question. To address this, we propose a neural network model to map persistence images to cosmological parameters. Through a parameter recovery test, we demonstrate that our model makes accurate and precise estimates, considerably outperforming conventional Bayesian inference approaches.

Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to resolve density fluctuation in the IGM puts a stringent requirement on the resolution of such simulations which in turn limits the volumes which can be modelled, even on most powerful supercomputers. In this work, we present a novel modeling method which combines physics-driven simulations with data-driven generative neural networks to produce outputs that are qualitatively and statistically close to the outputs of hydrodynamical simulations employing 8 times higher resolution. We show that the Ly-$\alpha$ flux field, as well as the underlying hydrodynamic fields, have greatly improved statistical fidelity over a low-resolution simulation. Importantly, the design of our neural network allows for sampling multiple realizations from a given input, enabling us to quantify the model uncertainty. Using test data, we demonstrate that this model uncertainty correlates well with the true error of the Ly-$\alpha$ flux prediction. Ultimately, our approach allows for training on small simulation volumes and applying it to much larger ones, opening the door to producing accurate Ly-$\alpha$ mock skies in volumes of Hubble size, as will be probed with DESI and future spectroscopic sky surveys.

Binary systems in which a neutron star or black hole accretes material from a high-mass star are known as high-mass X-ray binaries (HMXBs). This chapter provides a brief introduction to the physics of wind accretion and an observational view of HMXBs, including their classification, X-ray spectra, X-ray variability, orbital and compact object properties, as well as studies of Galactic and Magellanic HMXB populations. Two classes of X-ray sources whose possible connections to HMXBs have been debated, ultraluminous X-ray sources and gamma-ray binaries, are also discussed. Approximately 300 HMXBs residing either in the Milky Way or the Magellanic Clouds have been discovered. The majority of these HMXBs host wind-accreting neutron stars. Their X-ray properties depend both on the interaction of the accreting material with the neutron star's strong magnetic field and the properties of the donor star's wind. Most HMXBs are classified as either supergiant XBs or Be XBs based on the spectral type of the donor star; these classes exhibit different patterns of X-ray variability and occupy different phase space in diagrams of neutron star spin versus orbital period. While studies of HMXBs in the Milky Way and Magellanic Clouds find that their luminosity functions have similar shapes, an overabundance of Be XBs in the Small Magellanic Cloud points to important variations of the HMXB population with metallicity and age.

The Cameras for Allsky Meteor Surveillance (CAMS) project, funded by NASA starting in 2010, aims to map our meteor showers by triangulating meteor trajectories detected in low-light video cameras from multiple locations across 16 countries in both the northern and southern hemispheres. Its mission is to validate, discover, and predict the upcoming returns of meteor showers. Our research aimed to streamline the data processing by implementing an automated cloud-based AI-enabled pipeline and improve the data visualization to improve the rate of discoveries by involving the public in monitoring the meteor detections. This article describes the process of automating the data ingestion, processing, and insight generation using an interpretable Active Learning and AI pipeline. This work also describes the development of an interactive web portal (the NASA Meteor Shower portal) to facilitate the visualization of meteor radiant maps. To date, CAMS has discovered over 200 new meteor showers and has validated dozens of previously reported showers.

IceCube-Gen2 is a planned extension to the existing IceCube Neutrino Observatory and will provide an order of magnitude increase in the detection rate of cosmic neutrinos by deploying ~10,000 sensors in a volume of ~8 cubic kilometers. As part of the upcoming IceCube Upgrade, we are developing prototype IceCube-Gen2 sensors to test all components in-situ in preparation for mass production required for IceCube-Gen2. The novel IceCube-Gen2 module will contain up to eighteen 4-inch photomultiplier tubes (PMTs). The signals for each PMT are digitized with a 2-channel, 12-bit ADC (low- and high-gain) at a rate of 60 MSps. In addition, each module contains LED flashers for in-ice calibration, an FPGA for performing in-module local coincidence of PMT signals, and onboard $\mu$SD flash memory for buffering data before it is sent to the surface. In this contribution, we discuss the electronics and data acquisition system design.

We present a conceptual design of a high-performance camera system with applications to neutrino detectors, deep sea exploration, and glaciology. The design combines ultra-sensitive cameras with a number of well-calibrated light sources enclosed in a pressure vessel. The instrument will be capable of withstanding extreme environments such as those encountered in Antarctica or the deep ocean, and be deployable as a standalone system that can be retrieved for deep-sea exploration or glaciology. The camera system is designed to be replicated and deployed in multiple detectors, requiring only modest modifications from one detector to another. The instrument combines a number of capabilities essential for neutrino detector calibrations, including characterization of the scattering and absorption properties of the optical medium, measurement of geometries via photogrammetry, and detector surveillance. The ability to deploy the instrument at different detector sites also offers opportunities for cross-calibration efforts. We present the conceptual design of the instrument and describe plans to produce a prototype.

The recent discovery and evidence of neutrino signals from distant sources, TXS 0506+056 and NGC 1068 respectively, provide opportunities to search for rare interactions of neutrinos that they might encounter on their paths. One potential scenario of interest is the interaction between neutrinos and dark matter, which is invisible and expected to be abundantly spread over the Universe. Various astrophysical observations have implied the existence of dark matter. When high-energy neutrinos from extragalactic sources interact with dark matter during their propagation, their spectra might show suppressions at specific energy ranges, where such interactions occur. These attenuation signatures from the interaction might be measurable on Earth with large neutrino telescopes such as the IceCube Neutrino Observatory. This analysis is focused on the search for rare interactions of high-energy neutrinos from the IceCube-identified astrophysical neutrino sources with dark matter in sub-GeV masses and several benchmark mediator cases using the upgoing track-like events. In this poster, sensitivity studies about the interaction of neutrinos and dark matter are presented.

EX Lupi is the prototype by which EXor-type outbursts were defined. It has experienced multiple accretion-related bursts and outbursts throughout the last decades, whose study have greatly extended our knowledge about the effects of these types of events. This star experienced a new burst in 2022. We used multi-band photometry to create color-color and color-magnitude diagrams to exclude the possibility that the brightening could be explained by a decrease in extinction. We obtained VLT/X-shooter spectra to determine the Lacc and Macc during the peak of the burst and after its return to quiescence using 2 methods: empirical relationships between line luminosity and Lacc, and a slab model of the whole spectrum. We examined the 130 year light curve of EX Lupi to provide statistics on the number of outbursts experienced during this period of time. Our analysis of the data taken during the 2022 burst confirmed that a change in extinction is not responsible for the brightening. Our two approaches in calculating the Macc were in agreement, and resulted in values that are 2 orders of magnitude above what had previously been estimated, thus suggesting that EX Lupi is a strong accretor even when in quiescence. We determined that in 2022 March the Macc increased by a factor of 7 with respect to the quiescent level. We also found hints that even though the Macc had returned to almost its pre-outburst levels, certain physical properties of the gas had not returned to the quiescent values. We found that the mass accreted during this three month event was 0.8 lunar masses, which is approximately half of what is accreted during a year of quiescence. We calculated that if EX Lupi remains as active as it has been for the past 130 years, during which it has experienced at least 3 outbursts and 10 bursts, then it will deplete the mass of its circumstellar material in less than 160000 yr.

Metallicity gradients refer to the sloped radial profile of metallicities of gas and stars and are commonly seen in disk galaxies. A well-defined metallicity gradient of the Galactic disk is observed particularly well with classical Cepheids, which are good stellar tracers thanks to their period-luminosity relation allowing precise distance estimation and other advantages. However, the measurement of the inner-disk gradient has been impeded by the incompleteness of previous samples of Cepheids and limitations of optical spectroscopy in observing highly reddened objects. Here we report the metallicities of 16 Cepheids measured with high-resolution spectra in the near-infrared YJ bands. These Cepheids are located at 3-5.6 kpc in the Galactocentric distance, R(GC), and reveal the metallicity gradient in this range for the first time. Their metallicities are mostly between 0.1 and 0.3 dex in [Fe/H] and more or less follow the extrapolation of the metallicity gradient found in the outer part, R(GC) larger than 6.5 kpc. The gradient in the inner disk may be shallower or even flat, but the small sample does not allow to determine the slope precisely. More extensive spectroscopic observations would also be necessary for studying minor populations, if any, with higher or lower metallicities that were reported in previous literature. In addition, three-dimensional velocities of our inner-disk Cepheids show the kinematic pattern that indicates non-circular orbits caused by the Galactic bar, which is consistent with the patterns reported in recent studies on high-mass star-forming regions and red giant branch stars.

We report the detection of carbon monoxide (CO) and dust, formed under hostile conditions, in recurrent nova V745 Sco about 8.7 days after its 2014 outburst. The formation of molecules or dust has not been recorded previously in the ejecta of a recurrent nova. The mass and temperature of the CO and dust are estimated to be T(CO) = 2250 +/- 250 K, M(CO) = (1 to 5) E-8 solar masses, and T(dust) = 1000 +/- 50 K, M(dust) approximately E-8 to E-9 solar masses respectively. At the time of their detection, the shocked gas was at a high temperature of approximately E+7 K as evidenced by the presence of coronal lines. The ejecta were simultaneously irradiated by a large flux of soft X-ray radiation from the central white dwarf. Molecules and dust are not expected to form and survive in such harsh conditions; they are like snowflakes in a furnace. However, it has been posited in other studies that, as the nova ejecta plow through the red giant's wind, a region exists between the forward and reverse shocks that is cool, dense and clumpy where the dust and CO could likely form. We speculate that this site may also be a region of particle acceleration, thereby contributing to the generation of gamma-rays.

We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC(4-3) lines are not detected, the HCO+(4-3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO+/HCN intensity ratio of $\geq$2.2. Based on the photodissociation region modelling, the most likely explanation for the elevated HCO+/HCN ratio is that SDP.81 has low mechanical heating - less than 10% of the total energy budget - and a sub-solar metallicity, Z=0.5 Z$_\odot$. While such conditions might not be representative of the general population of high-redshift dusty galaxies, lower-than-solar metallicity might have a significant impact on gas masses inferred from CO observations. In addition, we report the detection of CO(0-1) absorption from the foreground lensing galaxy and CO(1-0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the southeast.

The nature of Dark Matter remains one of the most important unresolved questions of fundamental physics. Many models, including the Weakly Interacting Massive Particles (WIMPs), assume Dark Matter to be a particle and predict a weak coupling with Standard Model matter. If Dark Matter particles can scatter off nuclei in the vicinity of a massive object, such as a star or a planet, they may lose kinetic energy and become gravitationally trapped in the center of such objects, including Earth. As Dark Matter accumulates in the center of the Earth, self-annihilation of WIMPs into Standard Model particles can result in an excess of neutrinos coming from the center of the Earth and detectable at the IceCube Neutrino Observatory, situated at the geographic South Pole. A search for excess neutrinos from these annihilations has been performed on 10 years of IceCube data, and results have been interpreted in the context of a number of WIMP annihilation channels ($\chi\chi\rightarrow\tau^+\tau^-/W^+W^-/b\bar{b}$) and masses ranging from 10 GeV to 10 TeV. We present the results from this analysis and compare the outcome with previous searches by other experiments. This analysis yields competitive and world-leading results for masses $m_\chi$ > 100 GeV.

We review the non-linear statistics of Primordial Black Holes that form from the collapse of over-densities in a radiation dominated Universe. We focus on the scenario in which large over-densities are generated by rare and Gaussian curvature perturbations during inflation. As new results, we show that the mass spectrum follows a power law determined by the critical exponent of the self-similar collapse up to a power spectrum dependent cut-off and, that the abundance related to very narrow power spectrums is exponentially suppressed. Related to this, we discuss and explicitly show that the Press-Schechter approximation, as well as the statistics of mean profiles, lead to wrong conclusions for the abundance and mass spectrum. Finally, we clarify that the transfer function in the statistics of initial conditions for primordial black holes formation (the abundance) does not play a significant role.

The horizontal irradiance at the sea surface is an informative light pollution indicator to study Artificial Light at Night (ALAN) effects on marine biodiversity (e.g.: zooplankton diel vertical migration). The Posch ratio (PR) for the horizontal irradiance (that is, the ratio of the horizontal irradiance to the zenith radiance) is a useful tool for estimating the irradiance from easily available measurements of the zenith night sky brightness. The PR definition has already been generalized for any pair of linear radiance indicators in any pair of arbitrarily chosen photometric bands, and can also be applied to estimate e.g. the average sky radiance or the radiance at some elevation above the horizon as a function of the radiance in any other direction of the sky. The PR for a single light source depends on the distance from the source, its angular and spectral emission pattern, and the state of the atmosphere. The PR for any set of sources is a linear combination of the individual PRs that each one would produce separately, with weights that can be easily derived from the relative contribution of each source to the zenith radiance. Whereas in populated lands the ALAN PR varies relatively fast from one location to another, due to the particular spatial distribution of lights, in coastal waters and the high seas the light pollution PR is a smooth function of the distance to the shoreline, due to the progressive lack of neighboring sources and the absence of obstacles. In this work we present the fundamental equations of the model and an example of application for the waters surrounding the Iberian Peninsula, North Africa and the West Mediterranean islands.

The science goals of IceCube-Gen2 include multi-messenger astronomy, astroparticle and particle physics. To this end, the observatory will include several detection methods, including a surface array and in-ice optical sensors. The array will have an approximately 8 km$^2$ surface coverage consisting of elevated scintillator panels and radio antennas to detect air showers in the energy range of several 100 TeV to a few EeV. The observatory's design is unique in that the measurements using the surface array can be combined with the observations of $\geq$300 GeV muons, produced in the hadronic cascades, using the optical detectors in the ice. This allows for an enhanced ability to study cosmic-ray and hadronic physics as well as to boost the sensitivity for astrophysical neutrinos from the southern sky by reducing the primary background, atmospheric muons. We will present the baseline design of the surface array and highlight the expected scientific capabilities.

The highly variable blazar OJ~287 is commonly discussed as an example of a binary black hole system. The 130 year long optical light curve is well explained by a model where the central body is a massive black hole of 18.35$\times$10$^9$ solar mass that supports a thin accretion disc. The secondary black hole of 0.15$\times$10$^9$ solar mass impacts the disc twice during its 12 year orbit, and causes observable flares. Recently, it has been argued that an accretion disc with a typical AGN accretion rate and above mentioned central body mass should be at least six magnitudes brighter than OJ~287's host galaxy and would therefore be observationally excluded. Based on the observations of OJ~287's radio jet, detailed in Marscher and Jorstad (2011), and up-to-date accretion disc models of Azadi et al. (2022), we show that the V-band magnitude of the accretion disc is unlikely to exceed the host galaxy brightness by more than one magnitude, and could well be fainter than the host. This is because accretion power is necessary to launch the jet as well as to create electromagnetic radiation, distributed across many wavelengths, and not concentrated especially on the optical V-band. Further, we note that the claimed V-band concentration of accretion power leads to serious problems while interpreting observations of other Active Galactic Nuclei. Therefore, we infer that the mass of the primary black hole and its accretion rate do not need to be smaller than what is determined in the standard model for OJ~287.

It has been long conjectured that a signature of Quantum Gravity will be Lorentz Invariance Violation (LIV) that could be observed at energies much lower than the Planck scale. One possible signature of LIV is an energy-dependent speed of photons. This can be tested with a distant transient source of very high-energy photons. We explore time-of-flight limits on LIV derived from LHAASO's observations of tens of thousands of TeV photons from GRB 221009A, the brightest gamma-ray burst of all time. For a linear (n=1) dependence of the photon velocity on energy, we find a lower limit on the LIV scale of 5.9 (6.2) E_pl for subluminal (superluminal) modes. These are comparable to the stringent limits obtained so far. For a quadratic model (n=2), the limits, which are currently the best available, are much lower, 5.8 (4.6) x 10^{-8} E_pl. Our analysis uses the publicly available LHAASO data which is only in the 0.2-7 TeV range. Higher energy data would enable us to improve these limits by a factor of 3 for $n=1$ and by an order of magnitude for $n=2$.

The error propagation among estimated parameters reflects the correlation among the parameters. We study the capability of machine learning of "learning" the correlation of estimated parameters. We show that machine learning can recover the relation between the uncertainties of different parameters, especially, as predicted by the error propagation formula. Gravitational lensing can be used to probe both astrophysics and cosmology. As a practical application, we show that the machine learning is able to intelligently find the error propagation among the gravitational lens parameters (effective lens mass $M_{L}$ and Einstein radius $\theta_{E}$) in accordance with the theoretical formula for the singular isothermal ellipse (SIE) lens model. The relation of errors of lens mass and Einstein radius, (e.g. the ratio of standard deviations $\mathcal{F}=\sigma_{\hat{ M_{L}}}/ \sigma_{\hat{\theta_{E}}}$) predicted by the deep convolution neural network are consistent with the error propagation formula of SIE lens model. As a proof-of-principle test, a toy model of linear relation with Gaussian noise is presented. We found that the predictions obtained by machine learning indeed indicate the information about the law of error propagation and the distribution of noise. Error propagation plays a crucial role in identifying the physical relation among parameters, rather than a coincidence relation, therefore we anticipate our case study on the error propagation of machine learning predictions could extend to other physical systems on searching the correlation among parameters.

According to various estimates, brown dwarfs (BD) should account for up to 25 percent of all objects in the Galaxy. However, few of them are discovered and well-studied, both individually and as a population. Homogeneous and complete samples of brown dwarfs are needed for these kinds of studies. Due to their weakness, spectral studies of brown dwarfs are rather laborious. For this reason, creating a significant reliable sample of brown dwarfs, confirmed by spectroscopic observations, seems unattainable at the moment. Numerous attempts have been made to search for and create a set of brown dwarfs using their colours as a decision rule applied to a vast amount of survey data. In this work, we use machine learning methods such as Random Forest Classifier, XGBoost, SVM Classifier and TabNet on PanStarrs DR1, 2MASS and WISE data to distinguish L and T brown dwarfs from objects of other spectral and luminosity classes. The explanation of the models is discussed. We also compare our models with classical decision rules, proving their efficiency and relevance.

Focusing on the redshift space observations with plane-parallel approximation and relying on the rotational dependency of the general definition of excursion sets, we introduce the so-called conditional moments of the first derivative ($cmd$) measures for the smoothed matter density field in $3$-Dimension. We derive the perturbative expansion of $cmd$ for real space and redshift space where peculiar velocity disturbs the galaxies' observed locations. Our criteria can successfully recognize the contribution of linear Kaiser and Finger of God effects. Our results demonstrate that the $cmd$ measure has significant sensitivity for pristine constraining the redshift space distortion parameter $\beta=f/b$ and interestingly, the associated normalized quantity in the Gaussian linear Kaiser limit has only $\beta$-dependency. Implementation of the synthetic anisotropic Gaussian field approves the consistency between the theoretical and numerical results. Including the first-order contribution of non-Gaussianity perturbatively in the $cmd$ criterion implies that the N-body simulations for the Quijote suite in redshift space has been mildly skewed with a higher value for the threshold greater than zero. The non-Gaussianity for the perpendicular direction to the line of sight in the redshift space for smoothing scales $R\gtrsim 20$ Mpc/h is almost the same as the real space, while the non-Gaussianity for along to the line of sight direction in redshift space is magnified.

Supernovae (SN) explosions are thought to be an important source of dust in galaxies. At the same time strong shocks from SNe are known as an efficient mechanism of dust destruction via thermal and kinetic sputtering. A critically important question of how these two hypotheses of SNe activity control the dust budget in galaxies is still not quite clearly understood. In this paper we address this question within 3D multi-fluid hydrodynamical simulations, treating separately the SNe injected dust and the dust pre-exicted in ambient interstellar gas. We focus primarily on how the injected and pre-existed dust are destroyed by shock waves and the hot gas in the SN bubble depending of the density of ambient gas. Within our model we estimate an upper limit of the SN-produced dust mass which can be supplied into interstellar medium. For a SN progenitor mass of 30 $M_\odot$ and the ejected dust mass $M_d=1~M_\odot$ we constrain the dust mass that can delivered into the ISM as $\geq 0.13~M_\odot$ provided the injected dust particles are large $a\geq 0.1~\mu$m.

Infrared observations of stellar orbits about Sgr A* probe the mass distribution in the inner parsec of the Galaxy and provide definitive evidence for the existence of a massive black hole. However, the infrared astrometry is relative and is tied to the radio emission from Sgr A* using stellar SiO masers that coincide with infrared-bright stars. To support and improve this two-step astrometry, we present new astrometric observations of 15 stellar SiO masers within 2 pc of Sgr A*. Combined with legacy observations spanning 25.8 years, we re-analyze the relative offsets of these masers from Sgr A* and measure positions and proper motions that are significantly improved compared to the previously published reference frame. Maser positions are corrected for epoch-specific differential aberration, precession, nutation, and solar gravitational deflection. Omitting the supergiant IRS 7, the mean position uncertainties are 0.46 mas and 0.84 mas in RA and Dec., and the mean proper motion uncertainties are 0.07 mas yr$^{-1}$ and 0.12 mas yr$^{-1}$, respectively. At a distance of 8.2 kpc, these correspond to position uncertainties of 3.7 AU and 6.9 AU and proper motion uncertainties of 2.7 km s$^{-1}$ and 4.6 km s$^{-1}$. The reference frame stability, the uncertainty in the variance-weighted mean proper motion of the maser ensemble, is 8 $\mu$as yr$^{-1}$ (0.30 km s$^{-1}$) in RA and 11 $\mu$as yr$^{-1}$ (0.44 km s$^{-1}$) in Dec., which represents a 2.3-fold improvement over previous work and a new benchmark for the maser-based reference frame.

We report non-destructive 3-dimensional imaging and analysis techniques for material returned by the Stardust cometary collector. Our technique utilizes 3-dimensional laser scanning confocal microscopy (3D LSCM) to image whole Stardust tracks, in situ, with attainable resolutions <90 nm/pixel edge. LSCM images illustrate track morphology and fragmentation history; image segmentation techniques provide quantifiable volumetric and dynamic measurements. We present a method for multipart image acquisition and registration in 3-D. Additionally, we present a 3D deconvolution method for aerogel, using a theoretically calculated point spread function for first-order corrections of optical aberrations induced by light diffraction and refractive index mismatches. LSCM is a benchtop technique and is an excellent alternative to synchrotron x-ray computed microtomography for optically transparent media. Our technique, developed over the past 2 years, is a non-invasive, rapid technique for fine-scale imaging of high value returned samples from the Stardust mission, as well as various other samples from the geosciences.

Djerfisherite is an important carrier of potassium in highly reduced enstatite chondrites, where it occurs in sub-round metal-sulfide nodules. These nodules were once free-floating objects in the protoplanetary nebula. Here, we analyze existing and new data to derive an equation of state (EOS) for djerfisherites of K_{6}(Cu,Fe,Ni)^{B} (Fe,Ni,Cu)^{C}_{24} S_{26}Cl structural formula. We use this EOS to calculate the thermal stability of djerfisherite coexisting in equilibrium with a cooling vapor of solar composition enriched in a dust analogous to anhydrous, chondritic interplanetary dust (C-IDP). We find that condensed mineral assemblages closely match those found in enstatite chondrites, with djerfisherite condensing above 1000 K in C-IDP dust enriched systems. Results may have implications for the volatile budgets of terrestrial planets, and the incorporation of K into early-formed, highly reduced, planetary cores. Previous work links enstatite chondrites to the planet Mercury, where the surface has a terrestrial K/Th ratio, high S/Si ratio, and very low FeO content. Mercury's accretion history may yield insights into Earth's.

With simple estimates based on recent observational data and the assumption that among known FRB sources most repeaters with extremely high rates of repetition (super-repeaters) are already identified, we demonstrate that the hypothesis that super-repeaters and one-off events come from the same population of magnetars is not self-contradictory. In this toy model, the super-repeater stage has a duration of about a few years and the period when one-off events are mostly emitted corresponds to the active life of a magnetar $\sim $few$\times 10^3$~years. Intervals between strong events (observed as one-off FRBs) from the same source are $\sim10$~years, corresponding to the expected time interval between giant flares.

Over the history of planetary exploration, atmospheric entry vehicles have been used to deliver probes and landers to Venus, Mars, Jupiter, and Titan. While the entry vehicles are tools for furthering scientific exploration, by delivering probes and landers which perform in-situ exploration, the entry vehicle and trajectory design in itself is of significant interest. Entering an atmosphere subjects the vehicle to deceleration and aerodynamic heating loads which the vehicle must withstand to deliver the probe or the lander inside the atmosphere. The conditions encountered depend on the destination, the atmosphere-relative entry speed, the vehicle type, ballistic coefficient, the vehicle geometry, and the entry-flight path angle. The driving constraints are the peak aerodynamic deceleration, peak heat rate, and the total heat load encountered during the critical phase of the entry. Jupiter presents the most extreme entry conditions, while Titan presents the most benign entry environment. This study presents a survey of the design trade space for atmospheric entry missions across the Solar System using carpet plots, along with benchmarks from historical missions, which can serve as a useful reference for designing future missions.

More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new halo density profile model for galaxies, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro-Frenk-White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterised density profile model is fit and tested using the DarkSage semi-analytic model of galaxy formation. We find that our density profile model can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the DarkSage simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this density profile. We use the SDSS DR10 galaxy group catalogue to validate this halo density profile model. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function.

There are expected to be $\sim 10^8$ isolated black holes (BHs) in the Milky Way. OGLE-2011-BLG-0462/MOA-2011-BLG-191 (OB110462) is the only such BH with a mass measurement to date. However, its mass is disputed: Lam et al. (2022a,b) measured a lower mass of $1.6 - 4.4 M_\odot$, while Sahu et al. (2022); Mr\'{o}z et al. (2022) measured a higher mass of $5.8 - 8.7 M_\odot$. We re-analyze OB110462, including new data from the Hubble Space Telescope (HST) and re-reduced Optical Gravitational Lensing Experiment (OGLE) photometry. We also re-reduce and re-analyze the HST dataset with newly available software. We find significantly different ($\sim 1$ mas) HST astrometry than Lam et al. (2022a,b) in the de-magnified epochs due to the amount of positional bias induced by a bright star $\sim$0.4 arcsec from OB110462. After modeling the updated photometric and astrometric datasets, we find the lens of OB110462 is a $6.0^{+1.2}_{-1.0} M_\odot$ BH. Future observations with the Nancy Grace Roman Space Telescope, which will have an astrometric precision comparable or better to HST but a field of view $100\times$ larger, will be able to measure hundreds of isolated BH masses via microlensing. This will enable the measurement of the BH mass distribution and improve understanding of massive stellar evolution and BH formation channels.

RAFGL2591 is a massive star-forming complex in the Cygnus-X region comprising of a cluster of embedded protostars and young stellar objects located at a distance of 3.33 kpc. We investigate low-frequency radio emission from the protostellar jet associated with RAFGL2591 using the Giant Metrewave Radio Telescope (GMRT) at 325, 610 and 1280 MHz. For the first time, we have detected radio jet lobes in the E-W direction, labelled as GMRT-1 and GMRT-2. While GMRT-1 displays a flat radio spectral index of $\alpha$ = -0.10 , GMRT-2 shows a steeply negative value $\alpha$ = -0.62 suggestive of non-thermal emission. H$_2$ emission maps show the presence of numerous knots, arcs and extended emission towards the East-West jet, excited by the protostar VLA 3. In addition, we report a few H$_2$ knots in the North-East and South-West for the first time. The radio lobes (GMRT-1, GMRT-2) and H$_2$ emission towards this region are understood in the context of the prominent East-West jet as well as its lesser-known sibling jet in the North-East and South-West direction. To model the radio emission from the lobes, we have employed a numerical model including both thermal and non-thermal emission and found number densities towards these lobes in the range 100 - 1000 cm$^{-3}$ . The misalignment of the East-West jet lobes exhibits a reflection symmetry with a bending of $\sim$ 20$\circ$ . We attempt to understand this misalignment through precession caused by a binary partner and/or a supersonic side wind from source(s) in the vicinity.

The forecast on a time scale of 1 or 2 hours is crucial for all kind of new generation facilities (ELTs) instrumentation supported by the adaptive optics that will be mainly operated in Service Mode. In a recent study (Masciadri et al. 2020) we have showed that we can forecast the seeing and atmospheric parameters at such short time scales using an autoregressive method achieving unprecedented model accuracies with a substantial gain with respect to forecasts performed the day before (i.e. on longer time scales) obtained with an atmospheric mesoscale model. Equally we showed a gain with respect to the method by persistence using simply real-time measurements in situ on the same short time scale (1-2 hours). The auto-regressive method makes use of the forecasts done with mesoscale atmospheric models and real-time measurements and since 2019 has been implemented in the operational forecast system ALTA Center supporting LBT observations. In this contribution we apply the same approach to the VLT site extending the method to all the main astroclimatic parameters i.e. the seeing, the wavefront coherence time, the isoplanatic angle and the ground layer fraction. We prove that such a method offers unprecedented forecast accuracies for all the astroclimatic parameters with clear gains with respect to the prediction by persistence. Preliminary calculations indicate also better accuracies than those obtained with the machine learning based approach using in situ measurements. We will apply soon such a method in an operational forecast system that we conceived for the VLT.

We introduce a method for extracting spectral information from energy-resolved light curves folded at the neutron star spin period (known as pulse profiles) in accreting X-ray binaries. Spectra of these sources are sometimes characterized by features superimposed on a smooth continuum, such as iron emission lines and cyclotron resonant scattering features. We address here the question on how to derive quantitative constraints on such features from energy-dependent changes in the pulse profiles. We developed a robust method for determining in each energy-selected bin the value of the pulsed fraction using the fast Fourier transform opportunely truncated at the number of harmonics needed to satisfactorily describe the actual profile. We determined the uncertainty on this value by sampling through Monte Carlo simulations a total of 1000 faked profiles. We rebinned the energy-resolved pulse profiles to have a constant minimum signal-to-noise ratio throughout the whole energy band. Finally we characterize the dependence of the energy-resolved pulsed fraction using a phenomenological polynomial model and search for features corresponding to spectral signatures of iron emission or cyclotron lines using Gaussian line profiles. We apply our method to a representative sample of NuSTAR observations of well-known accreting X-ray pulsars. We show that, with this method, it is possible to characterize the pulsed fraction spectra, and to constrain the position and widths of such features with a precision comparable with the spectral results. We also explore how harmonic decomposition, correlation, and lag spectra might be used as additional probes for detection and characterization of such features.

Mitigation of the impact of foreground contributions to measurements of Cosmic Microwave Background (CMB) polarization is a crucial step in modern CMB data analysis and is of particular importance for a detection of large-scale CMB $B$ modes. A large variety of techniques, based on different assumptions and aiming at either a full component separation or merely cleaning the foreground signals from the CMB maps, have been described in the literature. In this work, we consider this problem within a unified framework based on the maximum likelihood principle, under the assumption that the signal at each frequency can be represented as a linear mixture of sky templates. We discuss the impact of various additional assumptions on the final outcome of the procedure. We find that the component separation problem can be fully solved in two specific situations: when we either know the frequency scaling of all the components or can correctly model them with a limited number of unknown parameters, as is the case in the parametric component separation techniques; or when we either know the statistical properties of all the components, the foregrounds and CMB, or can correctly model them with a limited number of parameters, as for instance in SMICA-like approaches. However, we also show that much less stringent assumptions are sufficient if we only aim at recovering the cleaned CMB signal. In particular, we discuss a ``minimally informed'' non-parametric method based on maximum likelihood. The method only assumes that the component properties are independent on the sky direction, at least over some region of the sky, and that the CMB covariance is known up to some limited number of parameters. We apply this method to recover the CMB $B$ modes polarization signal in the context of forthcoming CMB experiments and compare its performance with that of the standard parametric... (abridged)

We introduce a new chemical enrichment and stellar feedback model into GIZMO, using the SIMBA sub-grid models as a base. Based on the state-of-the-art chemical evolution model of Kobayashi et al., SIMBA-C tracks 34 elements from H$\rightarrow$Ge and removes SIMBA's instantaneous recycling approximation. Furthermore, we make some minor improvements to SIMBA's base feedback models. SIMBA-C provides significant improvements on key diagnostics such as the knee of the $z=0$ galaxy stellar mass function, the faint end of the main sequence, and the ability to track black holes in dwarf galaxies. SIMBA-C also matches better with recent observations of the mass-metallicity relation at $z=0,2$. By not assuming instantaneous recycling, SIMBA-C provides a much better match to galactic abundance ratio measures such as [O/Fe] and [N/O]. SIMBA-C thus opens up new avenues to constrain feedback models using detailed chemical abundance measures across cosmic time.

We measure the [${\alpha}$/Fe] abundances for 183 quiescent galaxies at z = 0.60 - 0.75 with stellar masses ranging 10.4 \leq log10 10.4 $\leq$ log10 (M$_*$ /M$_\odot$) $\leq$ 11.6 selected from the LEGA-C survey. We estimate [${\alpha}$/Fe] from the ratio of the spectral indices Mgb (${\lambda} \sim 5177$ {\AA}) and Fe4383, compared to predictions of simple stellar population models. We find that 91% of quiescents in our sample have supersolar [${\alpha}$/Fe], with an average value of [${\alpha}$/Fe] = +0.24 $\pm$ 0.01. We find no significant correlation between [${\alpha}$/Fe] and stellar metallicity, mass, velocity dispersion, and average formation time. Galaxies that formed the bulk of their stellar mass on time scales shorter than 1 Gyr follow the same [${\alpha}$/Fe] distribution as those which formed on longer time scales. In comparison to local early-type galaxies and to stacked spectra of quiescent galaxies at z = 0.38 and z = 0.07, we find that the average [${\alpha}$/Fe] has not changed between z = 0.75 and the present time. Our work shows that the vast majority of massive quiescent galaxies at z $\sim$ 0.7 are ${\alpha}$-enhanced, and that no detectable evolution of the average [${\alpha}$/Fe] has taken place over the last $\sim$ 6.5 Gyr.

Aims : We aimed at constraining the accretion-ejection phenomena around the strongly-accreting Northern component of the S CrA young binary system (S CrA N) by deriving its magnetic field topology and its magnetospheric properties, and by detecting ejection signatures, if any. Methods : We led a two-week observing campaign on S CrA N with the ESPaDOnS optical spectropolarimeter at the Canada-France-Hawaii Telescope. We recorded 12 Stokes I and V spectra over 14 nights. We computed the corresponding Least-Square Deconvolution (LSD) profiles of the photospheric lines and performed Zeeman-Doppler Imaging (ZDI). We analysed the kinematics of noticeable emission lines, namely He I $\lambda 5876$ and the four first lines of the Balmer series, known to trace the accretion process. Conclusions : The findings from spectropolarimetry are complementary to those provided by optical long-baseline interferometry, allowing us to construct a coherent view of the innermost regions of a young, strongly accreting star. Yet, the strong and complex magnetic field reconstructed for S CrA N is inconsistent with the observed magnetic signatures of the emission lines associated to the post-shock region. We recommend a multi-technique, synchronized campaign of several days to put more constrains on a system that varies on a $\sim$ 1 day timescale.

Core-collapse supernovae (CCSNe) are prime candidates for gravitational-wave detectors. The analysis of their complex waveforms can potentially provide information on the physical processes operating during the collapse of the iron cores of massive stars. In this work we analyze the early-bounce rapidly rotating CCSN signals reported in the waveform catalog of Richers et al 2017, which comprises over 1800 axisymmetric simulations extending up to about $10$~ms of post-bounce evolution. It was previously established that for a large range of progenitors, the amplitude of the bounce signal, $\Delta h$, is proportional to the ratio of rotational-kinetic energy to potential energy, T/|W|, and the peak frequency, $f_{\rm peak}$, is proportional to the square root of the central rest-mass density. In this work, we exploit these relations to suggest that it could be possible to use such waveforms to infer protoneutron star properties from a future gravitational wave observation, if the distance and inclination are well known. Our approach relies on the ability to describe a subset of the waveforms in the early post-bounce phase in a simple template form depending only on two parameters, $\Delta h$ and $f_{\rm peak}$. We use this template to perform a Bayesian inference analysis of waveform injections in Gaussian colored noise for a network of three gravitational wave detectors formed by Advanced LIGO and Advanced Virgo. We show that it is possible to recover the injected parameters, peak frequency and amplitude, with an accuracy better than 10% for more than 50% of the detectable signals (given known distance and inclination angle). However, inference on waveforms from outside the Richers catalog is not reliable, indicating a need for carefully verified waveforms of the first 10 ms after bounce of rapidly rotating supernovae of different progenitors with agreement between different codes.

We could accurately predict the shadow path and successfully observe an occultation of a bright star by Chiron on 2022 December 15. The Kottamia Astronomical Observatory in Egypt did not detect the occultation by the solid body, but we detected three extinction features in the light curve that had symmetrical counterparts with respect to the central time of the occultation. One of the features is broad and shallow, whereas the other two features are sharper with a maximum extinction of $\sim$25$\%$ at the achieved spatial resolution of 19 km per data point. From the Wise observatory in Israel, we detected the occultation caused by the main body and several extinction features surrounding the body. When all the secondary features are plotted in the sky plane we find that they can be caused by a broad $\sim$580 km disk with concentrations at radii of 325 \pm 16 km and 423 \pm 11 km surrounding Chiron. At least one of these structures appears to be outside the Roche limit. The ecliptic coordinates of the pole of the disk are $\lambda$ = 151$^\circ~\pm$ 8$^\circ$ and $\beta$ = 18$^\circ~\pm$ 11$^\circ$, in agreement with previous results. We also show our long-term photometry indicating that Chiron had suffered a brightness outburst of at least 0.6 mag between March and September 2021 and that Chiron was still somewhat brighter at the occultation date than at its nominal pre-outburst phase. The outermost extinction features might be consistent with a bound or temporarily bound structure associated with the brightness increase. However, the nature of the brightness outburst is unclear, and it is also unclear whether the dust or ice released in the outburst could be feeding a putative ring structure or if it emanated from it.

In this Letter, we give a new geometrical interpretation of HI intensity mapping foreground filters in harmonic space, for both single-dish and interferometer mode surveys. We derive the foreground-filtered HI auto power spectrum and then extend this to the cross-power spectrum of HI with CMB lensing. Foreground filtering leads to a loss of isotropy in Fourier space, resulting in harmonic space non-diagonal correlations, which we show are small compared to the diagonal ones. On large scales, foreground filters lead to a major loss of power in the HI$\,\times\,$CMB lensing correlations.

We present ALMA [C II] 158 $\mu \rm{m}$ and dust continuum observations of the $z=6.79$ quasar J0109--3047 at a resolution of $0."045$ ($\sim$300 pc). The dust and [C II] emission are enclosed within a $\sim 500\, \rm{pc}$ radius, with the central beam ($r<144\ \rm{pc}$) accounting for $\sim$25\% (8\%) of the total continuum ([C II]) emission. The far--infrared luminosity density increases radially from $\sim$5 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ to a central value of $\sim$70 $\times 10^{11} L_\odot\ \rm{kpc}^{-2}$ (SFRD $\sim$50-700 $M_\odot\ \rm{yr}^{-1}\ \rm{kpc}^{-2}$). The [C II] kinematics are dispersion-dominated with a constant velocity dispersion of $137 \pm 6 \,\rm{km\ s}^{-1}$. The constant dispersion implies that the underlying mass distribution is not centrally peaked, consistent with the expectations of a flat gas mass profile. The lack of an upturn in velocity dispersion within the central beam is inconsistent with a black hole mass greater than $M_{\rm{BH}}<6.5\times 10^{8}\ M_\odot\ (2\sigma$ level), unless highly fine-tuned changes in the ISM properties conspire to produce a decrease of the gas mass in the central beam comparable to the black hole mass. Our observations therefore imply either that a) the black hole is less massive than previously measured or b) the central peak of the far-infrared and [C II] emission are not tracing the location of the black hole, as suggested by the tentative offset between the near-infrared position of the quasar and the ALMA continuum emission.

We universally search for evidence of kinematic and spatial correlation of supernova remnant (SNR) and molecular cloud (MC) associations for nearly all SNRs in the coverage of the MWISP CO survey, i.e. 149 SNRs, 170 SNR candidates, and 18 pure pulsar wind nebulae (PWNe) in 1 deg < l < 230 deg and -5.5 deg < b < 5.5 deg. Based on high quality and unbiased 12CO/13CO/C18O (J = 1--0) survey data, we apply automatic algorithms to identify broad lines and spatial correlations for molecular gas in each SNR region. The 91 percent of SNR-MC associations detected previously are identified in this paper by CO line emission. Overall, there could be as high as 80 percent of SNRs associated with MCs. The proportion of SNRs associated with MCs is high within the Galactic longitude less than ~50 deg. Kinematic distances of all SNRs that are associated with MCs are estimated based on systemic velocities of associated MCs. The radius of SNRs associated with MCs follows a lognormal distribution, which peaks at ~8.1 pc. The progenitor initial mass of these SNRs follows a power-law distribution with an index of ~-2.3 that is consistent with the Salpeter index of -2.35. We find that SNR-MC associations are mainly distributed in a thin disk along the Galactic plane, while a small amount distributed in a thick disk. With the height of these SNRs from the Galactic plane below ~45 pc, the distribution of the average radius relative to the height of them is roughly flat, and the average radius increases with the height when above ~45 pc.

Sunspots host various oscillations and wave phenomena like umbral flashes, umbral oscillations, running penumbral waves, and coronal waves. All fan loops rooted in sunspot umbra constantly show a 3-min period propagating slow magnetoacoustic waves in the corona. However, their origin in the lower atmosphere is still unclear. In this work, we studied these oscillations in detail along a clean fan loop system rooted in active region AR12553 for a duration of 4-hour on June 16, 2016 observed by Interface Region Imaging Spectrograph (IRIS) and Solar Dynamics Observatory (SDO). We traced foot-points of several fan loops by identifying their locations at different atmospheric heights from the corona to the photosphere. We found presence of 3-min oscillations at foot-points of all the loops and at all atmospheric heights. We further traced origin of these waves by utilising their amplitude modulation characteristics while propagating in the solar atmosphere. We found several amplitude modulation periods in the range of 9-14 min, 20-24 min, and 30-40 min of these 3-min waves at all heights. Based on our findings, we interpret that 3-min slow magnetoacoustic waves propagating in coronal fan loops are driven by 3-min oscillations observed at the photospheric foot-points of these fan loops in the umbral region. We also explored any connection between 3-min and 5-min oscillations observed at the photospheric foot-points of these loops and found them to be weakly coupled. Results provide clear evidence of magnetic coupling of the solar atmosphere through propagation of 3-min waves along fan loops at different atmospheric heights.

Rare, energetic (long) thermonuclear (Type I) X-ray bursts are classified either as intermediate-duration or superbursts, based on their duration. Intermediate-duration bursts lasting a few to tens of minutes are thought to arise from the thermonuclear runaway of a relatively thick (10^10 g/cm2) helium layer, while superbursts lasting hours are attributed to the detonation of an underlying carbon layer. We present a catalogue of 84 long thermonuclear bursts from 40 low-mass X-ray binaries, and defined from a new set of criteria distinguishing them from the more frequent short bursts. The three criteria are: (1) a total energy release larger than 10^40 erg, (2) a photospheric radius expansion phase longer than 10 s, and (3) a burst time-scale longer than 70 s. This work is based on a comprehensive systematic analysis of 70 bursts found with INTEGRAL, RXTE, Swift, BeppoSAX, MAXI, and NICER, as well as 14 long bursts from the literature that were detected with earlier generations of X-ray instruments. For each burst, we measure its peak flux and fluence, which eventually allows us to confirm the distinction between intermediate-duration bursts and superbursts. Additionally, we list 18 bursts that only partially meet the above inclusion criteria, possibly bridging the gap between normal and intermediate-duration bursts. With this catalogue, we significantly increase the number of long-duration bursts included in the MINBAR and thereby provide a substantial sample of these rare X-ray bursts for further study.

In circumstellar disks around T Tauri stars, visible and near-infrared stellar irradiation is intercepted by dust at the disk's optical surface and reprocessed into thermal infrared; this subsequently undergoes radiative diffusion through the optically thick bulk of the disk. The gas component -- overwhelmingly dominant by mass, but contributing little to the opacity -- is heated primarily by gas-grain collisions. In hydrodynamical simulations, however, typical models for this heating process (local isothermality, $\beta$-cooling, two-temperature radiation hydrodynamics) incorporate simplifying assumptions that limit their ranges of validity. To build on these methods, we develop a ``three-temperature" numerical scheme, which self-consistently models energy exchange between gas, dust, and radiation, as a part of the PLUTO radiation-hydrodynamics code. With a range of test problems in 0D, 1D, 2D, and 3D, we demonstrate the efficacy of our method, and make the case for its applicability to a wide range of problems in disk physics, including hydrodynamic instabilities and disk-planet interaction.

We investigate the low-frequency spectral emission from a network of superconducting cosmic string loops in hopes of explaining the observed radio synchrotron background. After considering constraints from a variety of astrophysical and cosmological measurements, we identify a best-fit solution with string tension $G\mu \simeq 6.5 \times 10^{-12}$ and current $\mathcal{I} \simeq 2.5 \times 10^6$ GeV. This model yields a convincing fit to the data and may be testable in the near future by spectral distortion (TMS, BISOU) and 21 cm experiments (HERA, SKA, REACH). We also find that soft photon heating protects us against current constraints from global $21$ cm experiments.

Machine Learning (ML) algorithms are becoming popular in cosmology for extracting valuable information from cosmological data. In this paper, we evaluate the performance of a Convolutional Neural Network (CNN) trained on matter density snapshots to distinguish clustering Dark Energy (DE) from the cosmological constant scenario and to detect the speed of sound ($c_s$) associated with clustering DE. We compare the CNN results with those from a Random Forest (RF) algorithm trained on power spectra. Varying the dark energy equation of state parameter $w_{\rm{DE}}$ within the range of -0.7 to -0.99, while keeping $c_s^2 = 1$, we find that the CNN approach results in a significant improvement in accuracy over the RF algorithm. The improvement in classification accuracy can be as high as $40\%$ depending on the physical scales involved. We also investigate the ML algorithms' ability to detect the impact of the speed of sound by choosing $c_s^2$ from the set $\{1, 10^{-2}, 10^{-4}, 10^{-7}\}$ while maintaining a constant $w_{\rm DE}$ for three different cases: $w_{\rm DE} \in \{-0.7, -0.8, -0.9\}$. Our results suggest that distinguishing between various values of $c_s^2$ and the case where $c_s^2=1$ is challenging, particularly at small scales and when $w_{\rm{DE}}\approx -1$. However, as we consider larger scales, the accuracy of $c_s^2$ detection improves. Notably, the CNN algorithm consistently outperforms the RF algorithm, leading to an approximate $20\%$ enhancement in $c_s^2$ detection accuracy in some cases.

Thick disks are a prevalent feature observed in numerous disk galaxies including our own Milky Way. Their significance has been reported to vary widely, ranging from a few to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation which has high spatial and stellar mass resolutions to investigate the issue of thick disk mass fraction. We also use the NewHorizon2 simulation that was run on the same initial conditions but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M* > 10^10 M_{\odot} in NewHorizon, the contribution of the thick disk was found to be 34 \pm 15% in r-band luminosity or 48 \pm 13% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick disk stars based on a single [{\alpha}/Fe] cut is severely contaminated by stars of different kinematic properties while missing a bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic mid-plane with high circularity and get heated with time via fluctuation in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties.

1ES1927+654 is a nearby active galactic nucleus that has shown an enigmatic outburst in optical/UV followed by X-rays, exhibiting strange variability patterns at timescales of months-years. Here we report the unusual X-ray, UV, and radio variability of the source in its post-flare state (Jan 2022- May 2023). Firstly, we detect an increase in the soft X-ray (0.3-2 keV) flux from May 2022- May 2023 by almost a factor of five, which we call the bright-soft-state. The hard X-ray 2-10 keV flux increased by a factor of two, while the UV flux density did not show any significant changes ($\le 30\%$) in the same period. The integrated energy pumped into the soft and hard X-ray during this period of eleven months is $\sim 3.57\times 10^{50}$ erg and $5.9\times 10^{49}$ erg, respectively. From the energetics, it is evident that whatever is producing the soft excess (SE) is pumping out more energy than either the UV or hard X-ray source. Since the energy source presumably is ultimately the accretion of matter onto the SMBH, the SE emitting region must be receiving the majority of this energy. In addition, the source does not follow the typical disc-corona relation found in AGNs, neither in the initial flare (in 2017-2019) nor in the current bright soft state (2022-2023). We found that the core (<1 pc) radio emission at 5 GHz gradually increased till March 2022 but showed a dip in August 2022. The G\"udel Benz relation ($L_{\rm radio}/L_{\rm X-ray}\sim 10^{-5}$), however, is still within the expected range for radio-quiet AGN and further follow-up radio observations are currently being undertaken.

Transverse oscillations that do not show significant damping in solar coronal loops are found to be ubiquitous. Recently, the discovery of high-frequency transverse oscillations in small-scale loops has been accelerated by the Extreme Ultraviolet Imager onboard Solar Orbiter. We perform a meta-analysis by considering the oscillation parameters reported in the literature. Motivated by the power law of the velocity power spectrum of propagating transverse waves detected with CoMP, we consider the distribution of energy fluxes as a function of oscillation frequencies and the distribution of the number of oscillations as a function of energy fluxes and energies. These distributions are described as a power law. We propose that the power law slope ($\delta=-1.40$) of energy fluxes depending on frequencies could be used for determining whether high-frequency oscillations dominate the total heating ($\delta < 1$) or not ($\delta > 1$). In addition, we found that the oscillation number distribution depending on energy fluxes has a power law slope of $\alpha=1.00$, being less than 2, which means that oscillations with high energy fluxes provide the dominant contribution to the total heating. It is shown that, on average, higher energy fluxes are generated from higher frequency oscillations. The total energy generated by transverse oscillations ranges from about $10^{20}$ to $10^{25}$ erg, corresponding to the energies for nanoflare ($10^{24}-10^{27}$ erg), picoflare ($10^{21}-10^{24}$ erg), and femtoflare ($10^{18}-10^{21}$ erg). The respective slope results imply that high-frequency oscillations could provide the dominant contribution to total coronal heating generated by decayless transverse oscillations.

The cold dark matter (CDM) model predicts galaxies have 100 times more dark matter mass than stars. Nevertheless, recent observations report the existence of dark-matter-deficient galaxies with less dark matter than expected. To solve this problem, we investigate the physical processes of galaxy formation in head-on collisions between gas-containing dark matter subhaloes (DMSHs). Analytical estimation of the collision frequency between DMSHs associated with a massive host halo indicates that collisions frequently occur within 1/10th of the virial radius of the host halo, with a collision timescale of about 10 Myr, and the most frequent relative velocity increases with increasing radius. Using analytical models and numerical simulations, we show the bifurcation channel of the formation of dark-matter-dominated and dark-matter-deficient galaxies. In the case of low-velocity collisions, a dark-matter-dominated galaxy is formed by the merging of two DMSHs. In the case of moderate-velocity collisions, the two DMSHs penetrate each other. However the gas medium collides, and star formation begins as the gas density increases, forming a dwarf galaxy without dark matter at the collision surface. In the case of high-velocity collisions, shock-breakout occurs due to the shock waves generated at the collision surface reaching the gas surface, and no galaxy forms. For example, the simulation demonstrates that a pair of DMSHs with a mass of 10^9 Msun containing gas of 0.1 solar metallicity forms a dark-matter-deficient galaxy with a stellar mass of 10^7 Msun for a relative velocity of 200 km/s.

Context. The advent of wide-band receiver systems on interferometer arrays enables one to undertake high-sensitivity and high-resolution radio continuum surveys of the Galactic plane in a reasonable amount of telescope time. However, to date, there are only a few such studies of the first quadrant of the Milky Way that have been carried out at frequencies below 1 GHz. The Giant Metrewave Radio Telescope (GMRT) has recently upgraded its receivers with wide-band capabilities (now called the uGMRT) and provides a good opportunity to conduct high resolution surveys, while also being sensitive to the extended structures. Aims. We wish to assess the feasibility of conducting a large-scale snapshot survey, the Metrewave Galactic Plane with the uGMRT Survey (MeGaPluG), to simultaneously map extended sources and compact objects at an angular resolution lower than $10''$ and a point source sensitivity of 0.15 mJy/beam. Methods. We performed an unbiased survey of a small portion of the Galactic plane, covering the W43/W44 regions ($l=29^\circ-35^\circ$ and $|b|<1^\circ$) in two frequency bands: 300$-$500 MHz and 550$-$750 MHz. The 200 MHz wide-band receivers on the uGMRT are employed to observe the target field in several pointings, spending nearly 14 minutes on each pointing in two separate scans. We developed an automated pipeline for the calibration, and a semi-automated self-calibration procedure is used to image each pointing using multi-scale CLEAN and outlier fields. Results. We produced continuum mosaics of the surveyed region at a final common resolution of $25''$ in the two bands that have central frequencies of 400 MHz and 650 MHz, with a point source sensitivity better than 5 mJy/beam. We plan to cover a larger footprint of the Galactic plane in the near future based on the lessons learnt from this study. (Abridged)

In granular systems, thermal cycling causes compaction, creep, penetration of dense objects, and ratcheting of grains against each other. On asteroid surfaces, thermal cycling is high amplitude and can happen billions of times in a few million years. We use a 1-dimensional thermophysical conductivity model to estimate the relative displacement of grains in proximity to one another, caused by variations in thermal conductivity or shadows. We find that grains would experience relative displacements of order a few microns during each thermal cycle. If thermal cycling causes diffusive transport, then the asteroid's few centimeters deep thermal skin depth could flow a few centimeters in a million years. Thermal cycling could cause long-distance flows on sloped surfaces, allowing fine materials to collect in depressions.

The mode profile of a coupled optical cavity often exhibits a resonant doublet, which arises from the strong coupling between its sub-cavities. Traditional readout methods rely on setting fields of different frequencies to be resonant in either sub-cavity, which is challenging in the case of strong coupling. In this regime, the coupled cavity behaves as a single resonator, and a field must be resonant in all its parts. Consequently, specialized sensing schemes are necessary to control strongly coupled cavities. To address this issue, we propose a novel technique for the relative measurement of the degrees of freedom of a strongly coupled cavity. Our approach enables simultaneous frequency stabilization and fine-tuning of frequency splitting in the resonant doublet. Overall, our proposed technique offers a promising solution to control the properties of coupled cavities, facilitating advanced applications in the fields of gravitational-wave detection, quantum cavity optomechanics, and other related areas.

There is a high demand for nuclear data in multidisciplinary subject like nuclear astrophysics. The two areas of nuclear physics which are most clearly related to one another are stellar evolution and nucleosynthesis. The necessity for nuclear data for astrophysical applications puts experimental methods as well as reliability and predicative ability of current nuclear models to the test. Despite recent, considerable advances, there are still significant issues and mysteries. Only a few characteristics of nuclear astrophysics are covered in the current work which include $^{20}$Ne(n,$\gamma$)$^{21}$Ne, $^{52}$Fe(n,$\gamma$)$^{53}$Fe, $^{53}$Fe(n,$\gamma$)$^{54}$Fe, $^{54}$Fe(n,$\gamma$)$^{55}$Fe and $^{55}$Fe(n,$\gamma$)$^{56}$Fe reactions which are important in stellar nucleosynthesis. The reaction rates are calculated using nuclear statistical model. These rates are subsequently fitted to polynomials of temperature T$_9$ in order to facilitate calculations for stellar nucleosynthesis.

The primary goal of the novel photosensor technology called ABALONETM (U.S. Pat. 9,064,678) has been to enable the next-generation of astroparticle physics experiments to open new windows on the universe with large-area detectors of unprecedented performance, radio-purity, robustness, and integration flexibility. This foundational technology provides the means for modern, scalable and cost-effective production. The wide range of possible application-specific detector configurations can also make a difference in homeland security and medical imaging. Thus far, we have reduced to practice two such detector configurations: (i) the composite sandwich Panel assembly that hosts matrices of closely packed ABALONE units (U.S. Pat. 10,823,861), and, quite differently, the Tandem Detector Modules that host two back-to-back oriented ABALONE photosensors. This article for the first time presents the latter configuration, specifically designed to provide $2\times2\pi$ angular acceptance for neutrino astronomy and other large-area radiation detectors.

Atom and matter interferometers are precise quantum sensing experiments that can probe differential forces along separated spacetime paths. Various atom and matter interferometer experiments have been proposed to study dark matter, gravitational waves, and exotic new physics. Increasingly, these experimental concepts have proposed space-based designs to maximize interrogation times and baselines. However, decoherence and phase shifts caused by astrophysical backgrounds could largely undermine or destroy the target sensitivity of the experiments. We calculate the decoherence effects induced by solar photons, the solar wind, cosmic rays, solar neutrinos and zodiacal dust on space-based atom and matter interferometers. We find that, in future space-based atom and matter interferometers, the solar wind generically produces decoherence beyond the quantum noise limit, without proper shielding. In addition, solar photons are also an important background for matter interferometers.

This paper presents a post hoc analysis of a deep learning-based full-disk solar flare prediction model. We used hourly full-disk line-of-sight magnetogram images and selected binary prediction mode to predict the occurrence of $\geq$M1.0-class flares within 24 hours. We leveraged custom data augmentation and sample weighting to counter the inherent class-imbalance problem and used true skill statistic and Heidke skill score as evaluation metrics. Recent advancements in gradient-based attention methods allow us to interpret models by sending gradient signals to assign the burden of the decision on the input features. We interpret our model using three post hoc attention methods: (i) Guided Gradient-weighted Class Activation Mapping, (ii) Deep Shapley Additive Explanations, and (iii) Integrated Gradients. Our analysis shows that full-disk predictions of solar flares align with characteristics related to the active regions. The key findings of this study are: (1) We demonstrate that our full disk model can tangibly locate and predict near-limb solar flares, which is a critical feature for operational flare forecasting, (2) Our candidate model achieves an average TSS=0.51$\pm$0.05 and HSS=0.38$\pm$0.08, and (3) Our evaluation suggests that these models can learn conspicuous features corresponding to active regions from full-disk magnetograms.

The nonlinear wave propagation in large extra spatial dimensions (on and above $d=2$) is investigated in the context of nonlinear electrodynamics theories that depends exclusively on the invariant $\mathcal{F}$. In this vein, we consider propagating waves under the influence of external uniform electric and magnetic fields. Features related to the blackbody radiation in the presence of a background constant electric field such as the generalization of the spectral energy density distribution and the Stefan-Boltzmann law are obtained. Interestingly enough, anisotropic contributions to the frequency spectrum appear in connection to the nonlinearity of the electromagnetic field. In addition, the long wavelength regime and the Wien's displacement law in this situation are studied. The corresponding thermodynamics quantities at thermal equilibrium, such as energy, pressure, entropy and heat capacity densities are contemplated as well.

The problem of electromagnetic emission generation in plasma with electron beams is relevant both for practical applications and for interpretation of radio emission processes in astrophysical systems. In this work, we consider the case of counter-propagating electron beams injection into plasma. Such systems may occur in cosmic plasma in the case of closely spaced particle acceleration regions, and they can also be implemented in laboratory facilities. Using particle-in-cell numerical simulations, we have shown that high beam-to-radiation conversion efficiency can be achieved in the case when beams excite small scale oblique plasma oscillations. In this case, the radiation is generated in the vicinity of the second harmonic of the plasma frequency. For such waves the surrounding plasma is transparent. It has been found that the efficiency and spectrum of the radiation are not dependent on the thickness of the beams. It has been also shown that parameters of the system necessary for efficient radiation generation by the discussed mechanism can be found using the exact linear theory of beam-plasma instability.

The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission has been orbiting Mars since 2014 and now has tens of thousands of orbits which we use to characterize Mars' dynamic space environment. Through global field line tracing with MAVEN magnetic field data we find an altitude dependent draping morphology that differs from expectations of induced magnetospheres in the vertical ($\hat Z$ Mars Sun-state, MSO) direction. We quantify this difference from the classical picture of induced magnetospheres with a Bayesian multiple linear regression model to predict the draped field as a function of the upstream interplanetary magnetic field (IMF), remanent crustal fields, and a previously underestimated induced effect. From our model we conclude that unexpected twists in high altitude dayside draping ($>$800 km) are a result of the IMF angle in the X, Y MSO plane. We propose that this is a natural outcome of current theories of induced magnetospheres but has been underestimated due to approximations of the IMF as solely $\pm \hat Y$ directed. We additionally estimate that distortions in low altitude ($<$800 km) dayside draping along $\hat Z$ are directly related to remanent crustal fields. We show dayside draping propagates down tail and previously reported inner magnetotail twists are likely caused by the crustal field of Mars, while the outer tail morphology is governed by an induced response to the IMF direction. We conclude with an updated understanding of induced magnetospheres which details dayside draping for multiple directions of the incoming IMF and discuss the repercussions of this draping for magnetotail morphology.

The atmospheric neutrino flux at very high energies is dominated by prompt neutrinos, mostly contributed by the decays of charmed hadrons produced in the forward direction by cosmic ray interactions with air nuclei. Theoretical predictions of the prompt atmospheric neutrino flux have large uncertainties mainly related to charm hadron production. Prompt neutrinos can also be studied through high-energy colliders. In particular, two ongoing forward experiments and the proposed Forward Physics Facility at the LHC can detect forward prompt neutrinos. We will present the kinematic regions relevant to the prompt atmospheric neutrino flux in terms of collider kinematic variables, the collision energy $\sqrt{s}$ and the center-of-mass rapidity of charm hadrons $y$, and discuss implications of the forward experiments at the LHC on the theoretical predictions of the prompt atmospheric neutrino flux.

Pixelated CdZnTe detectors are a promising imaging-spectrometer for gamma-ray astrophysics due to their combination of relatively high energy resolution with room temperature operation negating the need for cryogenic cooling. This reduces the size, weight, and power requirements for telescope-based radiation detectors. Nevertheless, operating CdZnTe in orbit will expose it to the harsh radiation environment of space. This work, therefore, studies the effects of $61 \ \mathrm{MeV}$ protons on $2 \times 2 \times 1 \ \mathrm{cm}^3$ pixelated CdZnTe and quantifies proton-induced radiation damage of fluences up to $2.6 \times 10^8 \ \mathrm{p/cm^2}$. In addition, we studied the effects of irradiation on two separate instruments: one was biased and operational during irradiation while the other remained unbiased. Following final irradiation, the $662 \ \mathrm{keV}$ centroid and nominal $1\%$ resolution of the detectors were degraded to $642.7 \ \mathrm{keV}, 4.9 \% \ ( \mathrm{FWHM})$ and $653.8 \ \mathrm{keV}, 1.75 \% \ (\mathrm{FWHM})$ for the biased and unbiased systems respectively. We therefore observe a possible bias dependency on proton-induced radiation damage in CdZnTe. This work also reports on the resulting activation and recovery of the instrument following room temperature and $60^{\circ}\mathrm{C}$ annealing.

Nuclear isomers play a key role in the creation of the elements in the universe and have a number of fascinating potential applications related to the controlled release of nuclear energy on demand. Particularly, $^{93m}$Mo isomer is a good candidate for studying the depletion of nuclear isomer via nuclear excitation by electron capture. For such purposes, efficient approach for $^{93m}$Mo production needs to be explored. In the present work, we demonstrate experimentally an efficient production of $^{93m}$Mo through $^{93}$Nb(p, n) reaction induced by intense laser pulse. When a ps-duration, 100-J laser pulse is employed, the $^{93m}$Mo isomer at 2425 keV (21/2$^+$, $T_{1/2}$ = 6.85 h) are generated with a high yield of $1.8\times10^6$ particles/shot. The resulting peak efficiency is expected to be $10^{17}$ particles/s, which is at least five orders of magnitudes higher than using classical proton accelerator. The effects of production and destruction of $^{93m}$Mo on the controversial astrophysical p-isotope $^{92}$Mo are studied. It is found that the $^{93}$Nb(p, n)-$^{93m}$Mo reaction is an important production path for ^{93m}Mo seed nucleus, and the influence of ^{93m}Mo-^{92}Mo reaction flow on ^{92}Mo production cannot be ignored. In addition, we propose to directly measure the astrophysical rate of (p, n) reaction using laser-induced proton beam since the latter one fits the Maxwell-Boltzmann distribution well. We conclude that laser-induced proton beam opens a new path to produce nuclear isomers with high peak efficiency towards the understanding of p-nuclei nucleosythesis.

Superradiantly unstable ultralight particles around a classical rotating black hole (BH) can form an exponentially growing bosonic cloud, which have been shown to provide an astrophysical probe to detect ultralight particles and constrain their mass. However, the classical BH picture has been questioned, and different theoretical alternatives have been proposed. Exotic compact objects (ECOs) are horizonless alternatives to BHs featuring a reflective surface (with a reflectivity $\mathcal{K}$) in place of the event horizon. In this work, we study superradiant instabilities around ECOs, particularly focusing on the influence of the boundary reflection. We calculate the growth rate of superradiant instabilities around ECOs, and show that the result can be related to the BH case by a correction factor $g_{\mathcal{K}}$, for which we find an explicit analytical expression and a clear physical interpretation. Additionally, we consider the time evolution of superradiant instabilities and find that the boundary reflection can either shorten or prolong the growth timescale. As a result, the boundary reflection alters the superradiance exclusion region on the Regge plane, potentially affecting constraints on the mass of ultralight particles. For a mildly reflective surface ($|\mathcal{K}|\lesssim 0.5$), the exclusion region is not substantially changed, while significant effects from the boundary reflection can occur for an extreme reflectivity ($|\mathcal{K}|\gtrsim0.9$).

Dense environments hosting compact binary mergers can leave an imprint on the gravitational-wave emission which, in turn, can be used to identify the characteristics of the environment. To demonstrate such scenario, we consider a simple setup of binary black holes with an environment consisting of a scalar-field bubble. We use this as a proxy for more realistic environments and as an example of the simplest physics beyond the standard model. We perform Bayesian inference on the numerical relativity waveforms using state-of-the-art waveform templates for black-hole mergers. In particular, we perform parameter estimation and model selection on signals from black-hole mergers with different mass ratios, total mass and loudness, and hosted by scalar-field bubbles of varying field amplitude. We find that sub-dominant gravitational wave modes emitted during the coalescence and ringdown are key to identifying environmental effects. In particular, we find that for face-on signals dominated by the quadrupole mode, the environment is only detectable if both the ringdown and the late inspiral/early merger fall in the detector band, so that inconsistencies can be found between the inferred binary parameters and those of the final black hole. For edge-on mergers we find that the environment can be detected even if only the ringdown is in band, thanks to the information encoded in the quasi-normal mode structure of the final black-hole.

We analyze the evolution of the mass density contrast in spherical perturbations of flat Friedman-Lemaitre-Robertson-Walker cosmologies. Both dark matter and dark energy are included. In the absence of dark energy the evolution equation coincides with that obtained by Bonnor within the ``Newtonian cosmology''.

The IceCube experiment has substantial simulation needs and is in continuous search for the most cost-effective ways to satisfy them. The most CPU-intensive part relies on CORSIKA, a cosmic ray air shower simulation. Historically, IceCube relied exclusively on x86-based CPUs, like Intel Xeon and AMD EPYC, but recently server-class ARM-based CPUs are also becoming available, both on-prem and in the cloud. In this paper we present our experience in running a sample CORSIKA simulation on both ARM and x86 CPUs available through Google Kubernetes Engine (GKE). We used the production binaries for the x86 instances, but had to build the binaries for ARM instances from source code, which turned out to be mostly painless. Our benchmarks show that ARM-based CPUs in GKE were not only the most cost-effective but were also the fastest in absolute terms in all the tested configurations. While the advantage is not drastic, about 20% in cost-effectiveness and less than 10% in absolute terms, it is still large enough to warrant an investment in ARM support for IceCube.

We explore how quantum gravity effects, manifested through the breaking of discrete symmetry responsible for both Dark Matter and Domain Walls, can have observational effects through CMB observations and gravitational waves. To illustrate the idea we consider a simple model with two scalar fields and two $\mathcal{Z}_2$ symmetries, one being responsible for Dark Matter stability, and the other spontaneously broken and responsible for Domain Walls, where both symmetries are assumed to be explicitly broken by quantum gravity effects. We show the recent gravitational wave spectrum observed by several pulsar timing array projects can help constrain such effects.