Papers for Thursday, Feb 08 2024

A number of scenarios have been proposed to explain the low velocity dispersion (and hence possible absence of dark matter) of the low surface brightness galaxies NGC1052-DF2 and NGC1052-DF4. Most of the proposed mechanisms are based on the removal of dark matter via the interaction of these galaxies with other objects. A common feature of these processes is the prediction of very faint tidal tails, which should be revealed by deep imaging ({\mu}g > 30 mag/arcsec2). Using ultra-deep images obtained with the Gemini telescopes, about 1 mag deeper than previously published data, we analyzed the possible presence of tidal tails in both galaxies. We confirm the presence of tidal tails in NGC1052-DF4, but see no evidence for tidal effects in NGC1052-DF2, down to surface brightnesses of {\mu}g=30.9 mag/arcsec2. We therefore conclude that while the absence of dark matter in NGC1052-DF4 could be attributed to the removal of dark matter by gravitational interactions, in the case of NGC1052-DF2 this explanation seems less plausible, and therefore other possibilities such as an incorrect distance measurement or that the system may be rotating could alleviate the dark matter problem.

Upcoming surveys will measure the cosmic microwave background (CMB) weak lensing power spectrum in exquisite detail, allowing for strong constraints on the sum of neutrino masses among other cosmological parameters. Standard CMB lensing power spectrum estimators aim to extract the connected non-Gaussian trispectrum of CMB temperature maps. However, they are generically dominated by a large Gaussian noise bias which thus needs to be subtracted at high accuracy. This is currently done with realistic map simulations of the CMB and noise, whose finite accuracy currently limits our ability to recover the CMB lensing on small-scale. In this paper, we propose a novel estimator which instead avoids this large Gaussian bias. This estimator relies only on the data and avoids the need for bias subtraction with simulations. Thus our bias avoidance method is (1) insensitive to misestimates in simulated CMB and noise models and (2) avoids the large computational cost of standard simulation-based methods like "realization-dependent $N^{(0)}$" (${\rm RDN}^{(0)}$). We show that our estimator is as robust as standard methods in the presence realistic inhomogeneous noise (e.g. from scan strategy) and masking. Moreover, our method can be combined with split-based methods, making it completely insensitive to mode coupling from inhomogeneous atmospheric and detector noise. We derive the corresponding expressions for our estimator when estimating lensing from CMB temperature and polarization. Although in this paper we specifically consider CMB weak lensing power spectrum estimation, we illuminate the relation between our new estimator, ${\rm RDN}^{(0)}$ subtraction, and general optimal trispectrum estimation. Through this discussion we conclude that our estimator is applicable to analogous problems in other fields which rely on estimating connected trispectra/four-point functions like large-scale structure.

The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines observed with the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption lines over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide, warm, ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($\Sigma_{\rm SFR}$), and the outflows in star-forming regions with higher $\Sigma_{\rm SFR}$ tend to be more ionized. Considering the ionized outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $\eta\gtrsim 0.2$. We compare the LMC's outflows with those detected in starburst galaxies, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto \Sigma_{\rm SFR}^{0.22}$ over a range of star-forming conditions ($\Sigma_{\rm SFR}\sim10^{-2}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way.

The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $\mu$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $\mu$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $\mu$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $\mu$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample.

We introduce a new class of hilltop and plateau potentials which can successfully unify inflation and dark energy resulting in Quintessential Inflation (QI). Interestingly these new potentials are related through an inverse transformation. Namely, if $V(\phi) = V_0 \, v(\phi)$ is a plateau potential then the inverse potential $V(\phi) = V_0 \, \left[v(\phi)\right]^{-1}$ describes hilltop QI. A simple example is provided by the KKLT-inspired potential $v(\phi) = \left\lbrack \frac{M^{2n} + \phi^{2n}}{N^{2n} + \phi^{2n}}\right\rbrack \,$. When $M/N \ll 1$ this potential describes plateau QI, while its inverse, $\left[v(\phi)\right]^{-1}$ describes hilltop QI. Other simple models of QI arise for the class of potentials $V(\phi) \sim \exp\left\lbrack\mp f(\phi)\right\rbrack$, where the $-$ ($+$) sign is associated with a plateau (hilltop). A key feature of this new class of QI models is the near absence of small parameters which are usually associated with the presence of dark energy. A forecast for the gravitational wave background generated in these models is provided.

Radiative transfer (RT) models are critical in the interpretation of exoplanetary spectra, in simulating exoplanet climates and when designing the specifications of future flagship observatories. However, most models differ in methodologies and input data, which can lead to significantly different spectra. In this paper, we present the experimental protocol of the MALBEC (Modeling Atmospheric Lines By the Exoplanet Community) project. MALBEC is an exoplanet model intercomparison project (exoMIP) that belongs to the CUISINES (Climates Using Interactive Suites of Intercomparisons Nested for Exoplanet Studies) framework which aims to provide the exoplanet community with a large and diverse set of comparison and validation of models. The proposed protocol tests include a large set of initial participating RT models, a broad range of atmospheres (from Hot Jupiters to temperate terrestrials) and several observation geometries, which would allow us to quantify and compare the differences between different RT models used by the exoplanetary community. Two types of tests are proposed: transit spectroscopy and direct imaging modeling, with results from the proposed tests to be published in dedicated follow-up papers. To encourage the community to join this comparison effort and as an example, we present simulation results for one specific transit case (GJ-1214 b), in which we find notable differences in how the various codes handle the discretization of the atmospheres (e.g., sub-layering), the treatment of molecular opacities (e.g., correlated-k, line-by-line) and the default spectroscopic repositories generally used by each model (e.g., HITRAN, HITEMP, ExoMol).

The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorithms to mine this information-rich dataset to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and Polycyclic Aromatic Hydrocarbons (PAHs), as measured over 150 pc-scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously-unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the ISM, (ii) the PAH 11.3/7.7 \mic~ band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes.

In this paper we review the basics of magneto-rotational properties of neutron stars focusing on spin-up/spin-down behavior at different evolutionary stages. The main goal is to provide equations for the spin frequency changes in various regimes (radio pulsar, propeller, accretor, etc.). Since presently spin behavior of neutron stars at all stages remains a subject of many uncertainties, we review different suggestions made over the years in the literature.

In a previous paper, we reported evidence for quasi-periodicities in the \textsl{TESS} light curves of BL Lacerate and two other blazars found serendipitously in the SDSS AGN catalog. In this work, we find tentative evidence for quasi-periodic features in the \textsl{TESS} observations of five sources in the fourth catalog of the Fermi--LAT (4FGL) sources: J090453.4$-$573503, J2345$-$1555, B0422+004, J002159.2$-$514028, and B0537$-$441. We analysed the \textsl{TESS} light curves of these blazars that we extracted using a customized approach. The quasi-periodic oscillations (QPOs) are searched for using two timing analysis techniques: generalized Lomb-Scargle periodogram and weighted wavelet Z-transform. Their apparent periods lie in the range of 2.8--6.5 days and have at least 3$\sigma$ significance in both of these methods. QPOs at such timescales can originate from the kink instability model which relates the quasi-periodic feature with the growth of kinks in the magnetized relativistic jets. We performed MCMC simulations to obtain the posterior distribution of parameters associated with this model and found the kink period consistent with previous studies.

Ultra-short period compact binaries are important sources of gravitational waves, which include e.g. the progenitors of type Ia supernovae or the progenitors of merger episodes that may lead to massive and magnetic single white dwarfs. J0526+5934 is an unresolved compact binary star with an orbital period of 20.5 minutes that belongs to this category. The visible component of J0526+5934 has been recently claimed to be a hot sub-dwarf star with a CO white dwarf companion. Our aim is to provide strong observational plus theoretical evidence that the primary star is rather an extremely-low mass white dwarf, although the hot subdwarf nature cannot be completely ruled out. We analyse optical spectra together with time-series photometry of the visible component of J0526+5934 to constrain its orbital and stellar parameters. We also employ evolutionary sequences for low-mass white dwarfs to derive independent values of the primary mass. From the analysis of our observational data, we find a stellar mass for the primary star in J0526+5934 of 0.26+-0.05 Msun, which perfectly matches the 0.237+-0.035 Msun independent measurement we derived from the theoretical evolutionary models. This value is considerably lower than the theoretically expected and generally observed mass range of hot subdwarf stars, but falls well within the mass limit values of extremely low-mass white dwarfs. We conclude J0526+5934 is the fifth ultra-short period detached double white dwarf currently known.

We report striking Doppler velocity gradients observed during the well-observed September 10th 2017 solar flare, and argue that they are consistent with the presence of an above-the-looptop termination shock beneath the flare current sheet. Observations from the Hinode Extreme-ultraviolet Imaging Spectrometer (EIS) measure plasma sheet Doppler shifts up to 35 km/s during the late-phase of the event. By comparing these line-of-sight flows with plane-of-sky measurements, we calculate total velocity downflows of 200+ km/s, orientated 6-10{\deg} out of the plane of sky. The observed velocities drop rapidly at the base of the hot plasma sheet seen in extreme ultraviolet, consistent with simulated velocity profiles predicted by our 2.5D magnetohydrodynamics model that features a termination shock at the same location. Finally, the striking velocity deceleration aligns spatially with the suppression of Fe XXIV non-thermal velocities, and a 35--50 keV hard X-ray looptop source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Together, these observations are consistent with the presence of a possible termination shock within the X8.2-class solar flare.

We report here on a project that has developed a practical approach to processing all-sky image collections on cloud platforms, using as an exemplar application the creation of three-color Hierarchical Progressive Survey (HiPS) maps of the 2MASS data set with the Montage Image Mosaic Engine on Amazon Web Services. We will emphasize issues that must be considered by scientists wishing to use cloud platforms to perform such parallel processing, so providing a guide for scientists wishing to exploit cloud platforms for similar large-scale processing. A HiPS map is based on the HEALPix sky-tiling scheme. Progressive zooming of a HiPS map reveals an image sampled at ever smaller or larger spatial scales that are defined by the HEALPix standard. Briefly, the approach used by Montage involves creating a base mosaic at the lowest required HEALPix level, usually chosen to match as closely as possible the spatial sampling of the input images, then cutting out the HiPS cells in PNG format from this mosaic. The process is repeated at successive HEALPix levels to create a nested collection of FITS files, from which PNG files are created that are shown in HiPS viewers. Stretching FITS files to produce PNGs is based on an image histogram. For composite regions (up and including the whole sky), the histograms for each tile can be combined to create a composite histogram for the region. Using this single histogram for each of the individual FITS files means all the PNGs are on the same brightness scale and displaying them side by side in a HiPS viewer produces a continuous uniform map across the entire sky.

For all active instruments, the Keck Observatory Archive (KOA) now ingests raw data from the Keck Telescopes within 1 minute of acquisition, quick-look reduced data within 5 minutes of creation, and science ready reduced data for four instruments as they are created by their automated pipelines. On August 1, 2023, KOA released the Observers Data Access Portal (ODAP), which enables observers at the telescope and their collaborators anywhere in the world to securely monitor and download science, calibration, and quick-look data as they are ingested into the archive. The portal is built using Python Socket IO.WebSockets that ensure metadata appear in the portal as the data themselves are ingested. The portal itself is a dynamic web interface built with React. It enables users to view and customize metadata fields, filter metadata according to data type, and download data as they are ingested or in bulk through wget scripts. Observers have used the ODAP since its release and have provided feedback that will guide future releases.

The strongly coupled hydrodynamic, magnetic, and radiation properties of the plasma in the solar chromosphere makes it a region of the Sun's atmosphere that is poorly understood. We use data obtained with the high-resolution Visible Broadband Imager (VBI) equipped with an H$\beta$ filter and the Visible Spectro-Polarimeter (ViSP) at the Daniel K. Inouye Solar Telescope to investigate the fine-scale structure of the plage chromosphere. To aid the interpretation of the VBI imaging data, we also analyze spectra from the CHROMospheric Imaging Spectrometer on the Swedish Solar Telescope. The analysis of spectral properties, such as enhanced line widths and line depths explains the high contrast of the fibrils relative to the background atmosphere demonstrating that H$\beta$ is an excellent diagnostic for the enigmatic fine-scale structure of the chromosphere. A correlation between the parameters of the H$\beta$ line indicates that opacity broadening created by overdense fibrils could be the main reason for the spectral line broadening observed frequently in chromospheric fine-scale structures. Spectropolarimetric inversions of the ViSP data in the Ca II 8542 {\AA} and Fe I 6301/6302 {\AA} lines are used to construct semiempirical models of the plage atmosphere. Inversion outputs indicate the existence of dense fibrils in the Ca II 8542 {\AA} line. The analyses of the ViSP data show that the morphological characteristics, such as orientation, inclination and length of fibrils are defined by the topology of the magnetic field in the photosphere. Chromospheric maps reveal a prominent magnetic canopy in the area where fibrils are directed towards the observer.

The two most favored explanations of the Fermi Galactic Center gamma-ray excess (GCE) are millisecond pulsars and self annihilation of the smooth dark matter halo of the galaxy. In order to distinguish between these possibilities, we would like to optimally use all information in the available data, including photon direction and energy information. To date, analyses of the GCE have generally treated directional and energy information separately, or have ignored one or the other completely. Here, we develop a method for analyzing the GCE that relies on simulation-based inference with neural posterior models to jointly analyze photon directional and spectral information while correctly accounting for the spatial and energy resolution of the telescope, here assumed to be the Fermi Large Area Telescope (LAT). Our results also have implications for analyses of the diffuse gamma-ray background, which we discuss.

Gamma-ray bursts (GRBs) can be probes of the early universe, but currently, only 26% of GRBs observed by the Neil Gehrels Swift Observatory GRBs have known redshifts ($z$) due to observational limitations. To address this, we estimated the GRB redshift (distance) via a supervised machine learning model that uses optical afterglow observed by Swift and ground-based telescopes. The inferred redshifts are strongly correlated (a Pearson coefficient of 0.93) with the observed redshifts, thus proving the reliability of this method. The inferred and observed redshifts allow us to estimate the number of GRBs occurring at a given redshift (GRB rate) to be 7.6-8 $yr^{-1} Gpc^{-1}$ for $1.9<z<2.3$. Since GRBs come from the collapse of massive stars, we compared this rate with the star formation rate highlighting a discrepancy of a factor of 3 at $z<1$.

Most, if not all, sun-like stars host one or more planets, making multi-planetary systems commonplace in our galaxy. We utilize hundreds of multi-planet simulations to explore the origin of such systems, focusing on their orbital architecture. The first set of simulations assumes in-situ assembly of planetary embryos, while the second explores planetary migration. After applying observational biases to the simulations, we compare them to 250+ observed multi-planetary systems, including 13 systems with planets in the habitable zone. For all of the systems, we calculate two of the so-called statistical measures: the mass concentration ($S_{c}$) and orbital spacing ($S_{s}$). After analytic and empirical analyses, we find that the measures are related to first-order with a power law: $S_{c} \sim S_{s}^\beta$. The in-situ systems exhibit steeper power-law relations relative to the migration systems. We show that different formation scenarios cover different regions in the $S_{s} - S_{c}$ diagram with some overlap. Furthermore, we discover that observed systems with $S_{s} < 30$ are likely dominated by the migration scenario, while those with $S_{s} \geq 30$ are likely dominated by the in-situ scenario. We apply these criteria to determine that a majority (62%) of observed multi-planetary systems formed via migration, whereas most systems with currently observed habitable planets formed via in-situ assembly. This work provides methods of leveraging the statistical measures ($S_{s}$ and $S_{c}$) to disentangle the formation history of observed multi-planetary systems based on their present-day architectures.

We report on the detection of extended $\gamma$-ray emission from lobes in the radio galaxy NGC 6251 using observation data of Fermi Large Area Telescope (Fermi-LAT). The maximum likelihood analysis results show that a radio morphology template provides a better fit than a point-like source description for the observational data at a confidence level of 8.1$\sigma$, and the contribution of lobes accounts for more than 50\% of the total $\gamma$-ray flux. Furthermore, the $\gamma$-ray energy spectra show a significant disparity in shape between the core and lobe regions, with a curved log-parabola shape observed in core region and a power-law form observed in lobes. Neither the core region nor the northwest lobe displays the significant flux variations in the long-term $\gamma$-ray light curves. The broadband spectral energy distributions of both core region and northwest lobe can be will explained with a single-zone leptonic model. The $\gamma$-rays of core region are due to the synchrotron-self-Compton process while the $\gamma$-rays from northwest lobe are interpreted as inverse Compton emission of the cosmic microwave background.

Core-collapse supernovae (CCSNe) are the terminal explosions of massive stars. While most massive stars explode as iron-core-collapse supernovae (FeCCSNe), slightly less massive stars explode as electron-capture supernovae (ECSNe), shaping the low-mass end of CCSNe. ECSNe was proposed $\sim 40$ years ago and first-principles simulations also predict their successful explosions. Observational identification and investigation of ECSNe are important for the completion of stellar evolution theory. To date, only one promising candidate has been proposed, SN 2018zd, other than the historical progenitor of the Crab Nebula, SN 1054. We present representative synthetic light curves of low-mass FeCCSNe and ECSNe exploding with energies in circumstellar media (CSM) estimated with theoretically or observationally plausible methods. The plateaus of the ECSNe are shorter, brighter, and bluer than those of the FeCCSNe. To investigate the robustness of their intrinsic differences, we adopted various explosion energies and CSM. Although they may have similar bolometric light-curve plateaus, ECSNe are bluer than FeCCSNe in the absence of strong CSM interaction, illustrating that multicolor observations are essential to identify ECSNe. This provides a robust indicator of ECSNe because the bluer plateaus stem from the low-density envelopes of their super-asymptotic-giant-branch progenitors. Furthermore, we propose a distance-independent method to identify ECSNe: $(g-r)_{t_{\rm PT}/2} < 0.008 \times t_{\rm PT} - 0.4$, i.e., blue $g-r$ at the middle of the plateau $(g-r)_{t_{\rm PT}/2}$, where $t_{\rm PT}$ is the transition epoch from plateau to tail. Using this method, we identified SN 2018zd as an ECSN, which we believe to be the first ECSN identified with modern observing techniques.

The detection of gravitational wave events has stimulated theoretical modeling of the formation and evolution of double compact objects (DCOs). However, even for the most studied isolated binary evolution channel, there exist large uncertainties in the input parameters and treatments of the binary evolution process. So far, double neutron stars (DNSs) are the only DCOs for which direct observations are available through traditional electromagnetic astronomy. In this work, we adopt a population synthesis method to investigate the formation and evolution of Galactic DNSs. We construct 324 models for the formation of Galactic DNSs, taking into account various possible combinations of critical input parameters and processes such as mass transfer efficiency, supernova type, common envelope efficiency, neutron star kick velocity, and pulsar selection effect. We employ Bayesian analysis to evaluate the adopted models by comparing with observations. We also compare the expected DNS merger rate in the Galaxy with that inferred from the known Galactic population of Pulsar-NS systems. Based on these analyses we derive favorable range of the aforementioned key parameters.

The Be isotopic measurements preliminarily reported by the AMS-02 Collaboration have reached an unprecedented energy of 12 GeV/$n$. As secondary cosmic rays (CRs), the Be isotopes include both stable and unstable species, which are crucial for constraining the propagation parameters of Galactic CRs. However, uncertainties in their production cross sections can skew the interpretation of the CR data, especially when cross-section measurements are of significantly lower quality than CR measurements. In this work, we consider the uncertainties of the cross sections to interpret the Be isotopic data by adopting a cross-section parametrization that fully utilizes the available experimental data. Owing to the high-quality measurements of the $^7$Be production cross section, we innovatively employ $^7$Be instead of $^9$Be to constrain propagation parameters. Notably, the diffusion halo thickness is constrained to $5.67\pm0.76$~kpc, representing a moderate value compared to previous analogous works. Combining the well-constrained CR propagation model and the precise CR measurements of $^9$Be, we conversely constrain the major production cross section of $^9$Be and find that it ought to be remarkably lower than previously thought. Our analysis also questions the reliability of certain cross sections measured by some experiments, potentially marking the first time CR data has been used to identify dubious nucleon production cross sections. The method presented in this work holds promise for analyzing upcoming isotopic data from other nuclei.

High-energy radiation from stars impacts planetary atmospheres deeply affecting their chemistry, providing departures from chemical equilibrium. While the upper atmospheric layers are dominated by ionizations induced by extreme ultraviolet radiation, deeper into the atmosphere molecular abundances are controlled by a characteristic X-ray dominated chemistry, mainly driven by an energetic secondary electron cascade. In this work, we aim at identifying molecular photochemically induced fingerprints in the transmission spectra of a giant planet atmosphere. We have developed a numerical code capable of synthesizing transmission spectra with arbitrary spectral resolution, exploiting updated infrared photoabsorption cross sections. Chemical mixing ratios are computed using a photochemical model, tailored to investigate high energy ionization processes. We find that in case of high levels of stellar activity, synthetic spectra in both low and high resolutions show significant, potentially observable out-of-equilibrium signatures arising mainly from CO, CH$_4$, C$_2$H$_2$, and HCN.

In this short paper we outline the case for a small radio telescope array in the southern hemisphere with operations dedicated to rapid follow-up and monitoring of astrophysical transients. We argue that the science harvest from such a facility would be very large, using AMI-LA as an outstanding example of how such a programme is already being operated in the north with an enormous track record of success. A southern radio transients facility would in turn take pressure off the Square Kilometre Array and the other world class larger arrays with 10-100 times more collecting area, which will never have the programme time available to comprehensively pursue this science. We discuss comparisons with the development of transient surveys and follow up in optical astronomy, and also how single millimetre dishes can contribute to radio transients science in the south. This paper is not a funding proposal aimed at any particular body, but rather a concept and discussion piece, and the authors welcome comments and feedback.

We perform three-dimensional numerical simulations of magnetized relativistic jets propagating in a uniform density environment in order to study the effect of the entrainment and the consequent deceleration, extending a previous work in which magnetic effects were not present. As in previous papers, our aim is to understand the connection between the jet properties and the resulting Fanaroff-Riley classification. We consider jets with different low densities, and therefore low power, and different magnetizations. We find that lower magnetization jets effectively decelerate to sub-relativistic velocities and may then result in an FR~I morphology on larger scales. At the opposite, in the higher magnetization cases, the entrainment and consequent deceleration are substantially reduced. }

The nature of the GeV gamma-ray Galactic center excess (GCE) in the data of Fermi-Large Area Telescope (LAT) is still to be unveiled. We present a new analysis of the inner Galaxy Fermi-LAT data at energies above 10 GeV, based on an innovative method which combines the skyFACT adaptive template fitting with and the 1pPDF pixel-count statistics. We find a strong evidence for the GCE also at high energies, $\sigma > 5$ regardless of the GCE spatial template. Remarkably, our fits prefer the bulge morphological model over the dark matter one at high significance, and show no evidence for an additional dark matter template on top of the bulge component. Through the 1pPDF analysis, we find that the model best describing the gamma-ray data requires a smooth, diffuse GCE following a bulge morphology, together with sub-threshold point sources. The 1pPDF fit reconstructs a consistent population of faint point sources down at least to $10^{-12}$ ph cm$^{-2}$ s$^{-1}$. Between $10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and $10^{-11}$ ph cm$^{-2}$ s$^{-1}$ the 1pPDF measures a number of point sources significantly higher than the ones in the Fermi 4FGL catalog. The robustness of our results brings further support to the attempt of explaining, at least partially, the high-energy tail of the GCE in terms of a population of point sources, likely corresponding to millisecond pulsars.

Solar magnetic fields alter scattering polarization in spectral lines like Sr I at 4607 {\AA} via the Hanle effect, making it a potential diagnostic for small-scale mixed-polarity photospheric magnetic fields. Recently, observational evidence for scattering polarization in the Sr I 4607 {\AA} at the solar disk center was found. Here, we investigate the reliability of the reconstruction method making possible this detection. To this end, we apply it to linear polarization profiles of the Sr I 4607 {\AA} line radiation emerging at the disk center obtained from a detailed 3D radiative transfer calculation in a magneto-hydrodynamic simulation snapshot with a small-scale dynamo contribution. The reconstruction method systematically reduces the scattering amplitudes by up to a factor of two, depending on the noise level. We demonstrate that the decrease can be attributed to two systematic errors: first, the physical constraint that underlies our assumptions regarding the dependence of scattering polarization on the quadrupolar moment of the radiation field, and second, the limitations of our method in accurately determining the sign of the radiation field tensor from the observed intensity image. However, consistently applying the reconstruction process and after taking into account image degradation effects due to the temporally variable image quality, such as imposed by seeing, observed and synthesized polarization signals show remarkable agreement. We thus conclude that the observed scattering polarization at solar disk center is consistent with that emerging from magneto-hydrodynamic model of the solar photosphere with an average magnetic field of 170 G at the visible surface.

In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths.

The search for neutrinos with energies greater than $10^{17}~$eV is being actively pursued. Although normalization of the dominant neutrino flux is highly uncertain, a floor level is guaranteed by the interactions of extragalactic cosmic rays with Milky Way gas. We estimate that this floor level gives an energy flux of $E^2\phi_\nu\simeq 10^{-13^{+0.5}_{-0.5}}~$GeV~cm$^{-2}$~sr$^{-1}$~s$^{-1}$ at $10^{18}~$eV, where uncertainties arise from the modeling of the gas distribution and the experimental determination of the mass composition of ultra-high-energy cosmic rays on Earth. Based on a minimal model of cosmic-ray production to explain the mass-discriminated energy spectra observed on Earth above $5{\times}10^{18}$~eV, we also present generic estimates of the neutrino fluxes expected from extragalactic production that generally exceed the aforementioned guaranteed floor. The prospects for detecting neutrinos above $10^{18}$~eV remain however challenging, unless proton acceleration to the highest energies is at play in a sub-dominant population of cosmic-ray sources or new physical phenomena are at work.

$\Lambda$CDM tensions are by definition model dependent; one sees anomalies through the prism of $\Lambda$CDM. Thus, progress towards tension resolution necessitates checking the consistency of the $\Lambda$CDM model to localise missing physics either in redshift or scale. Since the Universe is dynamical and redshift is a proxy for time, it is imperative to first perform consistency checks involving redshift, then consistency checks involving scale, as the next steps to settle the ``systematics versus new physics" debate and foster informed model building. We present a review of the hierarchy of assumptions underlying the $\Lambda$CDM cosmological model and comment on whether relaxing them can address the tensions. We focus on the lowest lying fruit of identifying missing physics through the identification of redshift dependent $\Lambda$CDM model fitting parameters. We highlight recent progress made on ${S_8:= \sigma_8 \sqrt{\Omega_{\rm m}/0.3}}$ tension and elucidate how similar progress can be made on $H_0$ tension. Our discussions indicate that $H_0$ tension, equivalently a redshift dependent $H_0$, and a redshift dependent $S_8$ imply a problem with background $\Lambda$CDM cosmology.

Dusty, distant, massive ($M_*\gtrsim 10^{11}\,\rm M_\odot$) galaxies are usually found to show a remarkable star-formation activity, contributing on the order of $25\%$ of the cosmic star-formation rate density at $z\approx3$--$5$, and up to $30\%$ at $z\sim7$ from ALMA observations. Nonetheless, they are elusive in classical optical surveys, and current near-infrared surveys are able to detect them only in very small sky areas. Since these objects have low space densities, deep and wide surveys are necessary to obtain statistically relevant results about them. Euclid will be potentially capable of delivering the required information, but, given the lack of spectroscopic features at these distances within its bands, it is still unclear if it will be possible to identify and characterize these objects. The goal of this work is to assess the capability of Euclid, together with ancillary optical and near-infrared data, to identify these distant, dusty and massive galaxies, based on broadband photometry. We used a gradient-boosting algorithm to predict both the redshift and spectral type of objects at high $z$. To perform such an analysis we make use of simulated photometric observations derived using the SPRITZ software. The gradient-boosting algorithm was found to be accurate in predicting both the redshift and spectral type of objects within the Euclid Deep Survey simulated catalog at $z>2$. In particular, we study the analog of HIEROs (i.e. sources with $H-[4.5]>2.25$), combining Euclid and Spitzer data at the depth of the Deep Fields. We found that the dusty population at $3\lesssim z\lesssim 7$ is well identified, with a redshift RMS and OLF of only $0.55$ and $8.5\%$ ($H_E\leq26$), respectively. Our findings suggest that with Euclid we will obtain meaningful insights into the role of massive and dusty galaxies in the cosmic star-formation rate over time.

In 2022 the DART mission spacecraft impacted the asteroid Dimorphos, the secondary body of the binary Didymos system, ejecting a large number of dust particles, rocks and boulders. The ESA Hera mission will reach the system in 2026 for post--impact studies and possible detection of orbiting fragments. We investigate the long term dynamics of the large boulders ejected by DART to test if any of these objects survive in orbit until the arrival of the Hera mission. To model the dynamics of the boulders we use a numerical model which includes the gravity of non-spherical Didymos and Dimorphos, the solar gravity and the radiation pressure. The SPICE kernels are used to define the correct reference frame for the integration. The dynamics of the boulders is highly chaotic and 1% of the initial boulders survive at least for 4 years on quasi--stable orbits. These orbits are characterised by wide oscillations in eccentricity in antiphase with those in inclination (including spin flips), a mechanism similar to the Kozai one. This behaviour may protect these bodies from close encounters with both asteroids. We also compute the distribution on the surfaces of the asteroids of sesquinary impacts which may influence the dust emission, after the initial DART impact, and the surface composition of the asteroids. The probability of observing boulders by the mission Hera is small but not negligible and an almost constant flux of escaping boulders is expected in the coming years since their lifetime after the DART impact covers a large time interval. Most of re--impacts on Dimorphos occur in the hemisphere opposite to the impact site, preferentially close to the equatorial plane.

We consider the phase shift in the gravitational wave signal induced by fast oscillations of scalar dark matter surrounding binary systems, which could be probed by the future experiments LISA and DECIGO. This effect depends on the local matter density and the mass of the dark matter particle. We compare it to the phase shift due to a standard dynamical friction term, which should generically be present. We find that the effect associated with the oscillations only dominates over the dynamical friction for dark matter masses below $10^{-21}$ eV, with masses below $10^{-23}$ eV implying cloud sizes that are too large to be realistic. Moreover, for masses of the order of $10^{-21}$ eV, LISA and DECIGO would only detect this effect for dark matter densities greater than that in the solar system by a factor $10^5$ or $10^4$ respectively. We conclude that this signal can be ignored for most dark matter scenarios unless very dense clouds of very light dark matter are created early in the Universe at a redshift $z\sim 10^4$.

Solar radio U-bursts are generated by electron beams traveling along closed magnetic loops in the solar corona. Low-frequency ($<$ 100 MHz) U-bursts serve as powerful diagnostic tools for studying large-sized coronal loops that extend into the middle corona. However, the positive frequency drift component (descending leg) of U-bursts has received less attention in previous studies, as the descending radio flux is weak. In this study, we utilized LOFAR interferometric solar imaging data from a U-burst that has a significant descending leg component, observed between 10 to 90 MHz on June 5th, 2020. By analyzing the radio source centroid positions, we determined the beam velocities and physical parameters of a large coronal magnetic loop that reached just about 1.3 $\rm{R_{\odot}}$ in altitude. At this altitude, we found the plasma temperature to be around 1.1 MK, the plasma pressure around 0.20 $\rm{mdyn,cm^{-2}}$, and the minimum magnetic field strength around 0.07 G. The similarity in physical properties determined from the image suggests a symmetric loop. The average electron beam velocity on the ascending leg was found to be 0.21 c, while it was 0.14 c on the descending leg. This apparent deceleration is attributed to a decrease in the range of electron energies that resonate with Langmuir waves, likely due to the positive background plasma density gradient along the downward loop leg.

In the context of future large surveys like the Euclid mission, extracting the cosmic web from galaxies at higher redshifts with more statistical power will become feasible, particularly within the group-cluster mass regime. Therefore, it is imperative to enlarge the number of metrics that can used to constrain our cosmological models at these large scales. The number of cosmic filaments surrounding galaxies, groups and clusters, namely the connectivity, has recently emerged as a compelling probe of the large-scale structures, and has been investigated in various observational and numerical analyses. In this first paper, we examine dark matter-only cosmological simulations using the widely used DisPerSE filament finder code under two theories of gravity: the Poisson ($\Lambda$CDM) and the Monge-Amp\`ere models, in order to quantify how alternative models of gravity alter the properties of the cosmic skeleton. We specifically focused on this alternative gravity theory due to its propensity to enhance the formation of anisotropic structures such as filaments, but it also makes them more resistant to collapse, which consequently reduces the formation of halos. Indeed, our findings reveal that replacing the Poisson equation has a significant impact on the hierarchical formation scenario. This is evidenced by examining the redshift evolution of both the slope and the offset of the connectivity. Additionally, we demonstrated that current observations are generally in better agreement with our well-established gravity model. Finally, our study suggests that filament connectivity in the group-cluster regime could serve as a probe of our gravity model at cosmological scales. We also emphasize that our approach could be extended to alternative theories of dark matter, such as warm or fuzzy dark matter, given the extraordinary datasets provided by next-generation surveys.

In this letter, motivated by double-peaked broad Balmer emission lines probably related to tidal disruption events (TDEs), a potential TDE candidate is reported in SDSS J160536+134838 (=SDSS J1605) at $z\sim0.44$ having quasar-like spectrum but with double-peaked broad H$\beta$. The long-term CSS light curve can be naturally described by a main-sequence star of $2.82_{-0.19}^{+0.20}{\rm M_\odot}$ disrupted by the central black hole (BH) of $144_{-21}^{+26}\times10^6{\rm M_\odot}$ in SDSS J1605. Meanwhile, the ASAS-SN light curves afterwards show none apparent trend variability, indicating the bright CSS flare in SDSS J1605 unique and different enough from variability of normal AGN. Furthermore, there is a consistency between the TDE model determined sizes of debris with the sizes of emission regions for the double-peaked broad H$\beta$ described by the accretion disk model, supporting the disk-like BLRs probably related to a central TDE in SDSS J1605. And the virial BH mass $\sim$7 times higher than the TDE model determined value can be naturally explained by R-L relation determined BLRs sizes very larger than the actual distance of emission regions related to TDEs debris in SDSS J1605. Although no clear conclusion on double-peaked broad lines absolutely related to TDEs, the results here provide clues to detect potential TDEs in AGN with double-peaked broad lines.

Ultra-steep spectrum (USS) radio sources with complex filamentary morphologies are a poorly understood subclass of diffuse radio source found in galaxy clusters. They are characterised by power law spectra with spectral indices less than -1.5, and are typically located in merging clusters. We present X-ray and radio observations of the galaxy cluster A272, containing a USS diffuse radio source. The system is an ongoing major cluster merger with an extended region of bright X-ray emission south of the core. Surface brightness analysis yields a $3\sigma$ detection of a merger shock front in this region. We obtain shock Mach numbers $M_\rho = 1.20 \pm 0.09$ and $M_T = 1.7 \pm 0.3$ from the density and temperature jumps, respectively. Optical data reveals that the system is a merger between a northern cool core cluster and a southern non-cool core cluster. We find that the USS source, with spectral index $\alpha^{\text{74 MHz}}_{\text{1.4 GHz}} = -1.9 \pm 0.1$, is located in the bright southern region. Radio observations show that the source has a double-lobed structure with complex filaments, and is centred on the brightest cluster galaxy of the southern subcluster. We provide two suggestions for the origin of this source; the first posits the source as an AGN relic that has been re-energised by the passing of a merger shock front, while the second interprets the complex structure as the result of two overlapping AGN radio outbursts. We also present constraints on the inverse Compton emission at the location of the source.

The interstellar object population of the Milky Way is a product of its stars. However, what is in fact a complex structure in the Solar neighbourhood has traditionally in ISO studies been described as smoothly distributed. Using a debiased stellar population derived from the Gaia DR3 stellar sample, we infer that the velocity distribution of ISOs is far more textured than a smooth Gaussian. The moving groups caused by Galactic resonances dominate the distribution. 1I/`Oumuamua and 2I/Borisov have entirely normal places within these distributions; 1I is within the non-coeval moving group that includes the Matariki (Pleiades) cluster, and 2I within the Coma Berenices moving group. We show that for the composition of planetesimals formed beyond the ice line, these velocity structures also have a chemodynamic component. This variation will be visible on the sky. We predict that this richly textured distribution will be differentiable from smooth Gaussians in samples that are within the expected discovery capacity of the Vera C. Rubin Observatory. Solar neighbourhood ISOs will be of all ages and come from a dynamic mix of many different populations of stars, reflecting their origins from all around the Galactic disk.

Hidden within the Gaia satellite's multiple data releases lies a valuable cache of dark companions. To facilitate the efficient and reliable detection of these companions via combined analyses involving Gaia, Hipparcos, and Tycho-2 catalogs, we introduce an astrometric modeling framework. This method incorporates analytical least square minimization and nonlinear parameter optimization techniques to a set of common calibration sources across the different space-based astrometric catalogues. This enables us to discern the error inflation, astrometric jitter, differential parallax zero-point, and frame rotation of various catalogues relative to Gaia DR3. Our findings yield the most precise Gaia DR2 calibration parameters to date, revealing notable dependencies on magnitude and color. Intriguingly, we identify sub-mas frame rotation between Gaia DR1 and DR3, along with an estimated astrometric jitter of 2.16 mas for the revised Hipparcos catalog. In a thorough comparative analysis with previous studies, we offer recommendations on calibrating and utilizing different catalogs for companion detection. Furthermore, we provide a user-friendly pipeline (https://github.com/ruiyicheng/Download_HIP_Gaia_GOST) for catalog download and bias correction, enhancing accessibility and usability within the scientific community.

Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing GEMS both through confirmed planets, as well as protoplanetary disk observations, and a combination of tests to reconcile these with theoretical predictions. We then introduce the \textit{Searching for GEMS} survey, where we utilize multi-dimensional nonparameteric statistics to simulate hypothetical survey scenarios to predict the required sample size of transiting GEMS with mass measurements to robustly compare their bulk-density with canonical hot-Jupiters orbiting FGK stars. Our Monte-Carlo simulations predict that a robust comparison requires about 40 transiting GEMS (compared to the existing sample of $\sim$ 15) with 5-$\sigma$ mass measurements. Furthermore, we discuss the limitations of existing occurrence estimates for GEMS, and provide a brief description of our planned systematic search to improve the occurrence rate estimates for GEMS.

By obtaining the assumption that planetary dust particles can escape from the gravitational attraction of a planet, we consider the possibility for the dust grains to leave the star's system by means of the radiation pressure. By taking the typical dust parameters into account, we consider their dynamics and show that they can reach the deep cosmos, taking part in panspermia. It has been shown that, during $5$ billion years, the dust grains will reach $10^5$ stellar systems, and by taking the Drake equation into account, it has been shown that the whole galaxy will be full of planetary dust particles.

Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation.

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a $1.4 \, \mathrm{M}_{\odot}$ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state to an accuracy of 10m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10^5)$ years to do so for a soft equation of state.

We present a study of the nearby low-metallicity dwarf galaxy IC 1613, focusing on the search for massive stars and related feedback processes, as well as for faint supernova remnants (SNR) in late stages of evolution. We obtained the deepest images of IC 1613 in the narrow-band H{\alpha}, He II and [S II] emission lines and new long-slit spectroscopy observations using several facilities (6-m BTA, 2.5m SAI MSU, and 150RTT telescopes), in combination with the multi-wavelength archival data from MUSE/VLT, VLA, XMM-Newton, and Swift/XRT. Our deep narrow-band photometry identifies several faint shells in the galaxy, and we further investigate their physical characteristics with the new long-slit spectroscopy observations and the archival multi-wavelength data. Based on energy balance calculations and assumptions about their possible nature, we propose that one of the shells is a possible remnant of a supernova explosion. We study five out of eight Wolf-Rayet (WR) star candidates previously published for this galaxy using the He ii emission line mapping, MUSE/VLT archival spectra, and new long-slit spectra. Our analysis discards the considered WR candidates and finds no new ones. We found P Cyg profiles in H{\alpha} line in two stars, which we classify as Luminous Blue Variable (LBV) star candidates. Overall, the galaxy IC 1613 may have a lower rate of WR star formation than previously suggested.

We study models where a scalar field has derivative and non-derivative couplings to the Ricci tensor and the co-Ricci tensor with a view to inflation. We consider both the metric formulation and the Palatini formulation. In the Palatini case, the couplings to the Ricci tensor and the Ricci scalar give the same result regardless of whether the connection is unconstrained or the non-metricity or the torsion is assumed to vanish. When the co-Ricci tensor is included, the unconstrained case and the zero torsion case are physically different. We reduce all the actions to the Einstein frame with minimally coupled matter, and find the leading order differences between the metric case and the Palatini cases.

We calculate the effect of dark matter on the ringdown waveform and shadow of supermassive black holes at the core of galaxies. Our main focus is on the supermassive black hole at the core of M87, which is large enough to allow for viable observational data. We compare the effects of a dark matter spike to those expected from a galactic halo of the same mass. The radial pressure is shown to be negligible for both the spike and the halo, implying that there is no difference between the isotropic case and the anisotropic case. Our calculation for the halo starts from the Hernquist density function for which the corresponding metric can be obtained analytically in closed form. The effect of the spike is orders of magnitude more significant than the halo as long as the distribution scale of the latter is within a few orders of magnitude of the value expected from observations. Our results indicate that the impact of the spike surrounding M87* on the ringdown waveform may in principle be detectable. Finally, we point out the somewhat surprising fact that existing Event Horizon Telescope observations of black hole shadows are within an order of magnitude from being able to detect, or rule out, the presence of a spike.

Symbolic regression (SR) searches for analytical expressions representing the relationship between a set of explanatory and response variables. Current SR methods assume a single dataset extracted from a single experiment. Nevertheless, frequently, the researcher is confronted with multiple sets of results obtained from experiments conducted with different setups. Traditional SR methods may fail to find the underlying expression since the parameters of each experiment can be different. In this work we present Multi-View Symbolic Regression (MvSR), which takes into account multiple datasets simultaneously, mimicking experimental environments, and outputs a general parametric solution. This approach fits the evaluated expression to each independent dataset and returns a parametric family of functions f(x; \theta) simultaneously capable of accurately fitting all datasets. We demonstrate the effectiveness of MvSR using data generated from known expressions, as well as real-world data from astronomy, chemistry and economy, for which an a priori analytical expression is not available. Results show that MvSR obtains the correct expression more frequently and is robust to hyperparameters change. In real-world data, it is able to grasp the group behaviour, recovering known expressions from the literature as well as promising alternatives, thus enabling the use SR to a large range of experimental scenarios.

The inflationary reheating phase begins when accelerated expansion ends. As all Standard Model particles are coupled to gravity, gravitational interactions will lead to particle production. This includes the thermal bath, dark matter and gravitational radiation. Here, we compute the spectrum of gravitational waves from the inflatoncondensate during the initial phase of reheating. As particular examples of inflation, we consider the Starobinsky model and T-models, all of which are in good phenomenological agreement with CMB anisotropy measurements. The T-models are distinguished by the shape of the potential about its minimum and can be approximated by $V \sim \phi^k$, where $\phi$ is the inflaton. Interestingly, the shape of the gravitational wave spectrum (when observed) can be used to distinguish among the models considered. As we show, the Starobinsky model and T-models with $k=2$, provide very different spectra when compared to models with $k=4$ or $k>4$. Observation of multiple harmonics in the spectrum can be interpreted as a direct measurement of the inflaton mass. Furthermore, the cutoff in frequency can be used to determine the reheating temperature.

K-essence theories are usually studied in the framework of one scalar field $\phi$. Namely, the Lagrangian of K-essence is the function of scalar field $\phi$ and its covariant derivative. However, in this paper, we explore a two-field pure K-essence, i.e. the corresponding Lagrangian is the function of covariant derivatives of two scalar fields without the dependency of scalar fields themselves. That is why we call it pure K-essence. The novelty of this K-essence is that its Lagrangian contains the quotient term of the kinetic energies from the two scalar fields. This results in the presence of many interesting features, for example, the equation of state can be arbitrarily small and arbitrarily large. As a comparison, the range for equation of state of quintessence is from $-1$ to $+1$. Interestingly, this novel K-essence can play the role of inflation field, dark matter and dark energy. Finally, the absence of the scalar fields themselves in the equations of motion makes the study considerable simple such that even the exact black hole solutions can be found.

The pre-merger detection of gravitational waves from the early inspiral phase of compact binary coalescence events would allow the observation of the earlier stages of the merger in the electromagnetic band. This would significantly impact multi-messenger astronomy, giving astronomers potential access to rich new information. Here, we introduce a proof-of-concept deep-learning-based approach to produce pre-merger early-warning alerts for binary black hole systems. We show the possibility of using a Long Short-Term Memory network trained on the whitened detector strain in the time domain to detect and classify compact binary events. In this work, we consider a single advanced Laser Interferometer Gravitational-Wave Observatory detector at design sensitivity and make approximate sensitivity and early warning capability comparisons with approximations to traditional matched filtering approaches. We find that our model is competitive in both aspects, and when applied to a simulated test dataset was able to produce an early alert up to four seconds before the merger.

We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) inhomogeneous and the homogeneous models as post processes by employing snapshots obtained by our self-consistent, realistic Boltzmann simulations in two spatial dimensions (2D). We show that the FFC growth rate is considerably larger in the inhomogeneous model than in the homogeneous model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability under inhomogeneous conditions. On the other hand, in the homogeneous model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collision-induced enhancement occurs on a long time scale in the inhomogeneous model. Despite that collision-induced enhancement does not appear on the short time scale, that should be taken into account to predict the final conversion probability.

Dark matter freeze-in at stronger coupling is operative when the Standard Model (SM) bath temperature never exceeds the dark matter mass. An attractive feature of this scenario is that it can be probed by direct detection experiments as well as at the LHC. In this work, we show how the mechanism can be realized in a simple UV complete framework, emphasizing the role of the maximal temperature of the SM thermal bath. We demonstrate that the maximal temperature can coincide with the reheating temperature or be close to it such that dark matter production is always Boltzmann-suppressed. This possibility is realized, for example, if the inflaton decays primarily into feebly interacting right-handed neutrinos, which subsequently generate the SM thermal bath. In this case, the SM sector temperature remains constant over cosmological times prior to reheating.

The Simplified General Perturbations 4 (SGP4) orbital propagation method is widely used for predicting the positions and velocities of Earth-orbiting objects rapidly and reliably. Despite continuous refinement, SGP models still lack the precision of numerical propagators, which offer significantly smaller errors. This study presents dSGP4, a novel differentiable version of SGP4 implemented using PyTorch. By making SGP4 differentiable, dSGP4 facilitates various space-related applications, including spacecraft orbit determination, state conversion, covariance transformation, state transition matrix computation, and covariance propagation. Additionally, dSGP4's PyTorch implementation allows for embarrassingly parallel orbital propagation across batches of Two-Line Element Sets (TLEs), leveraging the computational power of CPUs, GPUs, and advanced hardware for distributed prediction of satellite positions at future times. Furthermore, dSGP4's differentiability enables integration with modern machine learning techniques. Thus, we propose a novel orbital propagation paradigm, ML-dSGP4, where neural networks are integrated into the orbital propagator. Through stochastic gradient descent, this combined model's inputs, outputs, and parameters can be iteratively refined, surpassing SGP4's precision. Neural networks act as identity operators by default, adhering to SGP4's behavior. However, dSGP4's differentiability allows fine-tuning with ephemeris data, enhancing precision while maintaining computational speed. This empowers satellite operators and researchers to train the model using specific ephemeris or high-precision numerical propagation data, significantly advancing orbital prediction capabilities.

Using an idealised climate model incorporating seasonal forcing, we investigate the impact of rotation rate on the abundance of clouds on an Earth-like aquaplanet, and the resulting impacts upon albedo and seasonality. We show that the cloud distribution varies significantly with season, depending strongly on the rotation rate, and is well explained by the large-scale circulation and atmospheric state. Planetary albedo displays non-monotonic behaviour with rotation rate, peaking at around 1/2$\Omega_E$. Clouds reduce the surface temperature and total precipitation relative to simulations without clouds at all rotation rates, and reduce the dependence of total precipitation on rotation rate, causing non-monotonic behaviour and a local maximum around 1/8$\Omega_E$ ; these effects are related to the impacts of clouds on the net atmospheric and surface radiative energy budgets. Clouds also affect the seasonality. The influence of clouds on the extent of the winter Hadley cell and the intertropical convergence zone is relatively minor at slow rotation rates ($<$1/8$\Omega_E$ ), but becomes more pronounced at intermediate rotation rates, where clouds decrease their maximum latitudes. The timing of seasonal transitions varies with rotation rate, and the addition of clouds reduces the seasonal phase lag.

Collective neutrino oscillations are typically studied using the lowest-order quantum kinetic equation, also known as the mean-field approximation. However, some recent quantum many-body simulations suggest that quantum entanglement among neutrinos may be important and may result in flavor equilibration of the neutrino gas. In this work, we develop new quantum many-body models for neutrino gases in which any pair of neutrinos can interact at most once in their lifetimes. A key parameter of our models is $\gamma=\mu \Delta z$, where $\mu$ is the neutrino coupling strength, which is proportional to the neutrino density, and $\Delta z$ is the duration over which a pair of neutrinos can interact each time. Our models reduce to the mean-field approach in the limit $\gamma\to0$ and achieve flavor equilibration in time $t \gg (\gamma\mu)^{-1}$. These models demonstrate the emergence of coherent flavor oscillations from the particle perspective and may help elucidate the role of quantum entanglement in collective neutrino oscillations.

Mini-magnetospheres are small ion-scale structures that are well-suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device (LAPD) at the University of California - Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the micro-physics of the interaction. By analyzing the generalized Ohm's law terms, we find that the electron-only reconnection is driven by kinetic effects, through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface.