Papers for Friday, Jul 28 2023

We present highly sampled photometry of the supernova (SN) 2023ixf, a Type II SN in M101, beginning 2 days before its first known detection. To gather these data, we enlisted the global Unistellar Network of citizen scientists. These 252 observations from 115 telescopes show the SN's rising brightness associated with shock emergence followed by gradual decay. We measure a peak $M_{V}$ = -18.18 $\pm$ 0.09 mag at 2023-05-25 21:37 UTC in agreement with previously published analyses.

Convolutional neural networks (CNNs) have been shown to both extract more information than the traditional two-point statistics from cosmological fields, and marginalise over astrophysical effects extremely well. However, CNNs require large amounts of training data, which is potentially problematic in the domain of expensive cosmological simulations, and it is difficult to interpret the network. In this work we apply the learnable scattering transform, a kind of convolutional neural network that uses trainable wavelets as filters, to the problem of cosmological inference and marginalisation over astrophysical effects. We present two models based on the scattering transform, one constructed for performance, and one constructed for interpretability, and perform a comparison with a CNN. We find that scattering architectures are able to outperform a CNN, significantly in the case of small training data samples. Additionally we present a lightweight scattering network that is highly interpretable.

This chapter will go through the important nuclear reactions in stellar evolution and explosions, passing through the individual stellar burning stages and also explosive burning conditions. To follow the changes in the composition of nuclear abundances requires the knowledge of the relevant nuclear reaction rates. For light nuclei (entering in early stellar burning stages) the resonance density is generally quite low and the reactions are determined by individual resonances, which are best obtained from experiments. For intermediate mass and heavy nuclei the level density is typically sufficient to apply statistical model approaches. For this reason, while we discuss all burning stages and explosive burning, focusing on the reactions of importance, we will for light nuclei refer to the chapters by M. Wiescher, deBoer & Reifarth (Experimental Nuclear Astrophysics) and P. Descouvement (Theoretical Studies of Low-Energy Nuclear Reactions), which display many examples, experimental methods utilized, and theoretical approaches how to predict nuclear reaction rates for light nuclei. For nuclei with sufficiently high level densities we discuss statistical model methods used in present predictions of nuclear reaction cross sections and thermonuclear rates across the nuclear chart, including also the application to nuclei far from stability and fission modes.

The presence of aerosols is intimately linked to the global energy budget and composition of planet atmospheres. Their ability to reflect incoming light prevents energy from being deposited into the atmosphere, and they shape spectra of exoplanets. We observed one near-infrared secondary eclipse of WASP-80b with the Wide Field Camera 3 aboard the Hubble Space Telescope to provide constraints on the presence and properties of atmospheric aerosols. We detect a broadband eclipse depth of $34\pm10$ ppm for WASP-80b, making this the lowest equilibrium temperature planet for which a secondary eclipse has been detected so far with WFC3. We detect a higher planetary flux than expected from thermal emission alone at $1.6\sigma$ that hints toward the presence of reflecting aerosols on this planet's dayside. We paired the WFC3 data with Spitzer data and explored multiple atmospheric models with and without aerosols to interpret this spectrum. Albeit consistent with a clear dayside atmosphere, we found a slight preference for near-solar metallicities and for dayside clouds over hazes. We exclude soot haze formation rates higher than $10^{-10.7}$ g cm$^{-2}$s$^{-1}$ and tholin formation rates higher than $10^{-12.0}$ g cm$^{-2}$s$^{-1}$ at $3\sigma$. We applied the same atmospheric models to a previously published WFC3/Spitzer transmission spectrum for this planet and find weak haze formation. A single soot haze formation rate best fits both the dayside and the transmission spectra simultaneously. However, we emphasize that no models provide satisfactory fits in terms of chi-square of both spectra simultaneously, indicating longitudinal dissimilarity in the atmosphere's aerosol composition.

Kerr black holes radiate neutrinos in an asymmetric pattern, preferentially in the lower hemisphere relative to the black hole's rotation axis, while antineutrinos are predominantly produced in the upper hemisphere. Leveraging this asymmetric emission, we explore the potential of high-energy, $E_\nu \gtrsim 1$ TeV, neutrino and antineutrino detection to reveal crucial characteristics of an evaporating primordial black hole at the time of its burst when observed near Earth. We improve upon previous calculations by carefully accounting for the non-isotropic particle emission, as Earth occupies a privileged angle relative to the black hole's rotation axis. Additionally, we investigate the angular dependence of primary and secondary photon spectra and assess the evaporating black hole's time evolution during the final explosive stages of its lifetime. Although photon events outnumber neutrinos by about three orders of magnitude, we find that simultaneous measurements of these particles are indispensable for identifying the initial angular momentum and the black hole hemisphere facing Earth. This is particularly important for evaporating black holes within our solar system, at distances $\lesssim 5 \times 10^{-4}$ pc, and observed during the final 100 s of their lifetimes.

The Universe has evolved from an initial diffuse, uniform gas to a complex structure that includes both voids and high-density galaxy clusters connected by gaseous filaments, known as the Cosmic Web, and traced by 3D surveys of galaxies. The filamentary structure contains a significant fraction of the baryonic matter and is predicted to be mostly in the form of a moderately high temperature plasma, the Warm Hot Intergalactic Medium. Plasma at this temperature and ionization level emits mostly in soft X-rays. The filamentary structure, however, is hard to detect because the other sources contributing to the Diffuse X-ray Background are much brighter and, currently, there are very few reported detections of emission from the filaments. We report the first high-confidence level indirect detection of X-ray emission from the Warm Hot Intergalactic Medium. Applying the Power Spectrum Analysis to XMM-Newton and eROSITA data, we separated its contribution from other sources modeled in previous studies. Our result is in good agreement with numerical simulations and fills a critical gap in the picture of the large-scale structure of the Universe, in which filamentary gas, galaxies and dark matter interact and co-evolve.

JWST imaging of the rich galaxy cluster Abell 2744 at $z=0.308$ has been used by the UNCOVER team (Bezanson et al. 2022) to construct mosaic images in the NIRCAM filters. The exceptionally deep images in the ($F115W$, $F150W$, $F200W$) bands reveal a large population of unresolved pointlike sources across the field, the vast majority of which are globular clusters (GCs). To the limits of our photometry, more than 10,000 such objects were measured, most of which are in the halos of the five largest A2744 galaxies but which also include GCs around some satellite galaxies and throughout the IntraCluster Medium. Their luminosity function follows a lognormal shape, with the data reaching to within one magnitude of the classic GCLF turnover point. The colour index ($F115W-F200W$) in particular covers a range of $0.5$ mag, clearly resolving the expected internal spread of GC metallicities. The estimated GC masses are systematically higher than in present-day galaxies, consistent with a large, normal GC population seen at a $3.5~$Gyr earlier stage of dynamical evolution. Lastly, the spatial distribution of the bluer (more metal-poor) GCs resembles the gravitational lensing map of the cluster, consistent with recent theoretical suggestions.

We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<\lambda<3.4\mu$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$\alpha$ emission line, and the Paschen-series broad emission lines in a quasar at z>7. Hot torus dust is clearly detected as an upturn in the continuum emission at $\lambda_{\text{rest}}\simeq1.3\mu$m, leading to a black-body temperature of $T=1413.5^{+5.7}_{-7.4}$K. Compared to similarly-luminous quasars at 0<z<6, the hot dust in J1120+0641 is somewhat elevated in temperature (top 1%). The temperature is more typical among 6<z<6.5 quasars (top 25%), leading us to postulate a weak evolution in the hot dust temperature at z>6 ($2\sigma$ significance). We measure the black hole mass of J1120+0641 based on the H$\alpha$ Balmer line, $M_{\text{BH}}=1.52\pm0.17\cdot 10^9 M_\odot$, which is in good agreement with the previous rest-UV MgII black hole mass measurement. The black hole mass based on the Paschen-series lines is also consistent, indicating no significant extinction in the rest-frame UV measurement. The broad H$\alpha$, Pa-$\alpha$ and Pa-$\beta$ emission lines are consistent with an origin in a common broad-line region (BLR) with density log$N_H/\text{cm}^{-3}\geq 12$, ionisation parameter $-7<$log$U<-4$, and extinction E(B-V)$\lesssim 0.1$mag. These BLR parameters are consistent with similarly-bright quasars at 0<z<4. Overall, we find that both J1120+0641's hot dust torus and hydrogen BLR properties show no significant peculiarity when compared to luminous quasars down to z=0. The quasar accretion structures must have therefore assembled very quickly, as they appear fully "mature" less than 760 million years after the Big Bang.

We study the properties of associations of dwarf galaxies and their dependence on the environment. Associations of dwarf galaxies are extended systems composed exclusively of dwarf galaxies, considering as dwarf galaxies those galaxies less massive than $M_{\star, \rm max} = 10^{9.0}$ ${\rm M}_{\odot}\,h^{-1}$. We identify these particular systems using a semi-analytical model of galaxy formation coupled to a dark matter only simulation in the $\Lambda$ Cold Dark Matter cosmological model. To classify the environment, we estimate eigenvalues from the tidal field of the dark matter particle distribution of the simulation. We find that the majority, two thirds, of associations are located in filaments ($ \sim 67$ per cent), followed by walls ($ \sim 26 $ per cent), while only a small fraction of them are in knots ($ \sim 6 $ per cent) and voids ($ \sim 1 $ per cent). Associations located in more dense environments present significantly higher velocity dispersion than those located in less dense environments, evidencing that the environment plays a fundamental role in their dynamical properties. However, this connection between velocity dispersion and the environment depends exclusively on whether the systems are gravitational bound or unbound, given that it disappears when we consider associations of dwarf galaxies that are gravitationally bound. Although less than a dozen observationally detected associations of dwarf galaxies are currently known, our results are predictions on the eve of forthcoming large surveys of galaxies, which will enable these very particular systems to be identified and studied.

Recent observational studies have uncovered a small number of very metal-poor stars with cold kinematics in the Galactic disc and bulge. However, their origins remain enigmatic. We select a total of 138 Milky Way (MW) analogs from the TNG50 cosmological simulation based on their $z=0$ properties: disky morphology, stellar mass, and local environment. In order to make more predictive statements for the MW, we further limit the spatial volume coverage of stellar populations in galaxies to that targeted by the upcoming 4MOST high-resolution survey of the Galactic disc and bulge. We find that across all galaxies, $\sim$20 per cent of very metal-poor (${\rm [Fe/H]} < -2$) stars belong to the disk, with some analogs reaching 30 per cent. About 50$\pm$10 per cent of the VMP disc stars are, on average, older than 12.5 Gyr and $\sim$70$\pm$10 per cent come from accreted satellites. A large fraction of the VMP stars belong to the halo ($\sim$70) and have a median age of 12 Gyr. Our results with the TNG50 cosmological simulation confirm earlier findings with simulations of fewer individual galaxies, and suggest that the stellar disc of the Milky Way is very likely to host significant amounts of very- and extremely-metal-poor stars that, although mostly of ex situ origin, can also form in situ, reinforcing the idea of the existence of a primordial Galactic disc.

The radial metallicity distribution of the Milky Way's disc is an important observational constraint for models of the formation and evolution of our Galaxy. It informs our understanding of the chemical enrichment of the Galactic disc and the dynamical processes therein, particularly radial migration. We investigate how the metallicity changes with guiding radius in the thin disc using a sample of red-giant stars with robust astrometric, spectroscopic and asteroseismic parameters. Our sample contains $668$ stars with guiding radii $4$ kpc < $R_\mathrm{g}$ < $11$ kpc and asteroseismic ages covering the whole history of the thin disc with precision $\approx 25\%$. We use MCMC analysis to measure the gradient and its intrinsic spread in bins of age and construct a hierarchical Bayesian model to investigate the evolution of these parameters independently of the bins. We find a smooth evolution of the gradient from $\approx -0.07$ dex/kpc in the youngest stars to $\approx -0.04$ dex/kpc in stars older than $10$ Gyr, with no break at intermediate ages. Our results are consistent with those based on asteroseismic ages from CoRoT, with that found in Cepheid variables for stars younger than $1$ Gyr, and with open clusters for stars younger than $6$ Gyr. For older stars we find a significantly lower metallicity in our sample than in the clusters, suggesting a survival bias favouring more metal-rich clusters. We also find that the chemical evolution model of Chiappini (2009) is too metal-poor in the early stages of disc formation. Our results provide strong new constraints for the growth and enrichment of the thin disc and radial migration, which will facilitate new tests of model conditions and physics.

The presence or absence of star clusters in galaxies, and the properties of star cluster populations compared to their host galaxy properties, are important observables for validating models of cluster formation, galaxy formation, and galaxy assembly. In this work, we apply a Bayesian approach to fit two models to data from surveys of young clusters in star forming galaxies. The first model is a logistic regression, which allows us to include galaxies which do not have any young clusters. The second model is a hurdle model, which includes galaxies with zero clusters and also incorporates information about the total mass in the cluster system. We investigate two predictors (star formation rate and total stellar mass in the galaxy) and look at clusters younger than 10 or 100 Myr. We find that in all cases, star formation rate is the better predictor for both the probability of hosting clusters and the total mass in the cluster system. We compare our results to similar models for old globular clusters, and conclude that star cluster formation was more abundant and more efficient at higher redshifts, likely because of the high gas content of galaxies at that time.

Aims. Traditional star-galaxy classification techniques often rely on feature estimation from catalogues, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification's reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. By leveraging convolutional neural networks (CNN) and additional information about the source position, we aim to accurately classify all stars and galaxies within a survey, particularly those with a signal-to-noise ratio (S/N) near the detection limit. Methods. The AutoSourceID-Classifier (ASID-C) algorithm developed here uses 32x32 pixel single filter band source cutouts generated by the previously developed ASID-L code. ASID-C utilizes CNNs to distinguish these cutouts into stars or galaxies, leveraging their strong feature-learning capabilities. Subsequently, we employ a modified Platt Scaling calibration for the output of the CNN. This technique ensures that the derived probabilities are effectively calibrated, delivering precise and reliable results. Results. We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, outperforms similar codes like SourceExtractor. ASID-C opens up new possibilities for accurate celestial object classification, especially for sources with a S/N near the detection limit. Potential applications of ASID-C, like real-time star-galaxy classification and transient's host identification, promise significant contributions to astronomical research.

The light element lithium is formed by nucleosynthesis during the Big Bang. Its abundance can help to define the parameters of the early universe. To find this primordial value, it is necessary to determine Li abundances in the oldest stars because it is readily destroyed by nuclear reactions in stellar interiors. We have made high-resolution (~45,000) spectroscopic observations of five identical unevolved turn-off stars in the 13 Gyr old globular cluster M5. In our analysis we find a range in Li abundance of a factor of two; the spread is five times the individual error. The comparison of these results with those for turn-off stars from five other globular clusters reveals similarly large range in Li. Lithium in M5 and the other clusters all have stars above the field star Li plateau, but none are as high as the predictions for primordial Li. The maximum values for Li are the same in all six clusters. Multiple generations of stars are found in many globular clusters; those later generations are expected to have formed from Li-depleted gas. Such second- and later-generation stars would have no Li. However, only one of the six clusters has a few unevolved stars with upper limits on the Li abundance.

The Cosmic Evolution Early Release Science Survey (CEERS), targeting the Extended Groth Strip extragalactic field, is one of the JWST Director's Discretionary Early Release Science programs. To date, all observations have been executed and include NIRCam/MIRI imaging and NIRSpec/NIRCam spectroscopic exposures. Here, we discuss the MIRI imaging, which includes eight pointings, four of which provide deep imaging with the bluer bands (F560W, F770W) and four with contiguous wavelength coverage in F1000W, F1280W, F1500W, and F1800W, where two of these also include coverage in F770W and F2100W. We present a summary of the data, the data quality, and data reduction. The data reduction is based on the JWST Calibration Pipeline combined with custom modifications and additional steps designed to enhance the output quality, including improvements in astrometry and the removal of detector artifacts. We estimate the image depth of the reduced mosaics, and show that these generally agree with expectations from the Exposure Time Calculator. We compare the MIRI F560W and F770W flux densities for bright sources to measurements from Spitzer/IRAC Ch3 (5.8 $\mu$m) and Ch4 (8.0 $\mu$m), and we find that they agree with systematic differences of $<0.1$ mag. For the redder MIRI bands, we assess their quality by studying the spectral energy distributions (SEDs) of Galactic stars. The SEDs are consistent with the expected Rayleigh-Jeans law with a deviation $\sim 0.03$ mag, indicating that the MIRI colors are reliable. We also discuss all publicly released data products (images and source catalogs), which are available on the CEERS website (https://ceers.github.io/).

Characterizing the physical properties of dust grains in a protoplanetary disk is critical to comprehending the planet formation process. Our study presents ALMA high-resolution observations of the young protoplanetary disk around DG Tau at a 1.3 mm dust continuum. The observations, with a spatial resolution of $\approx 0.04''$, or $\approx5$ au, revealed a geometrically thin and smooth disk without substantial substructures, suggesting that the disk retains the initial conditions of the planet formation. To further analyze the distributions of dust surface density, temperature, and grain size, we conducted a multi-band analysis with several dust models, incorporating ALMA archival data of the 0.87 mm and 3.1 mm dust polarization. The results showed that the Toomre $Q$ parameter is $\lesssim2$ at a 20 au radius, assuming a dust-to-gas mass ratio of 0.01. This implies that a higher dust-to-gas mass ratio is necessary to stabilize the disk. The grain sizes depend on the dust models, and for the DSHARP compact dust, they were found to be smaller than $\sim400$ $\mu$m in the inner region ($r\lesssim20$ au), while exceeding larger than 3 mm in the outer part. Radiative transfer calculations show that the dust scale height is lower than at least one-third of the gas scale height. These distributions of dust enrichment, grain sizes, and weak turbulence strength may have significant implications for the formation of planetesimals through mechanisms such as streaming instability. We also discuss the CO snowline effect and collisional fragmentation in dust coagulation for the origin of the dust size distribution.

We spectroscopically analyzed a sample of 85 Sloan Digital Sky Survey compact, oxygen-dominated galaxies located at redshift $z\sim0.001-0.350$, selected because of their large equivalent width of [O III]$\lambda$5007 (larger than 200\AA). These galaxies might be considered as extreme emission-line galaxies due to their strong [O III]$\lambda$5007 emission line. We detected high-ionization lines, even those related to the presence of Wolf-Rayet stars, in almost all the galaxies studied. We obtained physical properties (electron density and temperature) as well as chemical abundances by using the direct method based on the electron temperature. In this analysis, we obtained three different measurements of T(high): the usual one, T([O III]), but also that of T([S III]) and T([Ar III]) for five and three of the galaxies, respectively. Further, we established a new calibration for T([S III]). We determined oxygen, nitrogen, sulfur, neon, argon, iron, and chlorine abundances when possible, and we compared to the results of other late-type, low-metallicity galaxies, such as blue compact dwarfs, Ims, and green peas. From such a detailed study, we can conclude that the majority of the galaxies in this sample have similar metallicities to the SMC (about $12+\log(O/H)\approx8$dex), and that only 12% of the galaxies are extremely metal-poor, with abundances lower than $7.7$dex. Also, a comparison with some chemical evolution models as well as a brief discussion on the chemical evolution with time is considered.

Past results from the IceCube Collaboration have suggested that the blazar TXS 0506+056 is a potential source of astrophysical neutrinos. However, in the years since there have been numerous updates to event processing and reconstruction, as well as improvements to the statistical methods used to search for astrophysical neutrino sources. These improvements in combination with additional years of data have resulted in the identification of NGC 1068 as a second neutrino source candidate. This talk will re-examine time-dependent neutrino emission from TXS 0506+056 using the most recent northern-sky data sample that was used in the analysis of NGC 1068. The results of using this updated data sample to obtain a significance and flux fit for the 2014 TXS 0506+056 "untriggered" neutrino flare are reported

Outward migration of massive binary stars or black holes in their circumbinary disc is often observed in simulations and it is key to the formation of wide black hole binaries. Using numerical simulations of Population III (Pop III) star formation, we study the angular momentum of Pop III binaries and the torques between stars and gas discs to understand the origin of outward migration and high ellipticity. The outward migration of protostars is produced by gravitational torques exerted on them by their circumstellar minidiscs. The minidiscs, on the other hand, migrate outward mainly by gaining angular momentum by accreting gas from the circumbinary disc. The angular momentum transfer is most efficient for rapidly accreting equal-mass binaries, and weaker when the secondary mass is small or the massive companion evaporates the gas disc via radiative feedback. We conclude that outward migration and the formation of wide equal-mass massive binaries is common in metal-free/metal-poor star formation, mainly driven by their large accretion rates. We expect that the lower gas temperature and accretion rates in metal-enriched circumstellar discs would lead more often to inward migration and closer binary separations. We also observe inward migration for smaller mass Pop III protostars/fragments, leading to the rapid merging of sink particles and likely the formation of close binary black holes that, however, reach separations below the resolution of our simulations. We discuss the implications that Pop III separations and ellipticity may have on the interpretation that gravitational wave signals from merging intermediate-mass black holes come from Pop III remnants.

In the upcoming IceCube-Gen2 extension, the newly developed optical modules will include 4-inch PMTs. For this purpose, the manufacturers Hamamatsu and North Night Vision Technology have developed new PMT models to meet the requirements of the IceCube-Gen2 science case. The specifications include strict requirements on temporal resolution, detection efficiency, and dark noise. We summarize the efforts to measure these performance characteristics and show that both PMT models meet the performance specifications set by IceCube-Gen2. Prototype optical modules based on both PMT models will be deployed with the IceCube Upgrade in order to test them in situ and help decide on a vendor for the Gen2 extension.

Stellar white-light flares are believed to play an essential role on the physical and chemical properties of the atmosphere of the surrounding exoplanets. Here we report an optical monitoring campaign on the nearby flaring system EI~Cnc carried out by the Ground-based Wide Angle Cameras (GWAC) and its dedicated follow-up telescope. A superflare, coming from the brighter component EI~CncA, was detected and observed, in which four components are required to properly model the complex decay light curve. The lower limit of flare energy in the $R-$band is estimated to be $3.3\times10^{32}$ ergs. 27 flares are additionally detected from the GWAC archive data with a total duration of 290 hours. The inferred cumulative flare frequency distribution follows a quite shallow power-law function with a slope of $\beta=-0.50\pm 0.03$ over the energy range between $10^{30}$ and $10^{33}$ erg, which reinforces the trend that stars cooler than M4 show enhanced superflare activity. The flares identified in EI~Cnc enable us to extend the $\tau-E$ relationship previously established in the white-light superflares of solar-type stars down to an energy as low as $\sim10^{30}$erg (i.e., by three orders): $\tau\propto E^{0.42\pm0.02}$, which suggests a common flare mechanism for stars with a type from M to solar-like, and implies an invariant of $B^{1/3}\upsilon_{\rm A}$ in the white-light flares.

We studied the polarization properties of Cygnus X--1 in both hard and soft spectral states with Imaging X-ray Polarimetry Explorer (IXPE) observations. The polarization degree is $\sim 4\%$, and $\sim 2.6-2.8\%$ in the hard and soft states, respectively. The polarization angle is observed along the jet axis and remains the same in both states. Energy-dependent analysis revealed an increasing polarization with the energy. The observed polarization rules out the lamp-post corona and prefers a conical corona if the coronal geometry remains the same in the hard and soft states. The polarization degree is found to be higher for the disk emission, compare to the coronal emission. The scattering of the seed photons inside the corona is likely to reduce the polarization. This could explain the different polarization in different spectral states.

The energy of the transition from Galactic to extra-galactic origin of cosmic rays is one of the major unresolved issues of cosmic-ray physics. However, strong constraints can be obtained from studying the anisotropy in the arrival directions of cosmic rays. The sensitivity to cosmic-ray anisotropy is, in particular, a matter of statistics. Recently, the cosmic ray anisotropy measurements in the TeV to PeV energy range were updated from IceCube using 11 years of data. The IceCube-Gen2 surface array will cover an area about 8 times larger than the existing IceTop surface array with a corresponding increase in statistics and capability to investigate cosmic-ray anisotropy with higher sensitivity. In this contribution, we present details on the performed simulation studies and sensitivity to the cosmic-ray anisotropy signal for the IceCube-Gen2 surface array.

Sub-Neptune planets formed in the protoplanetary disk accreted hydrogen-helium (H,He) envelopes. Planet formation models of sub-Neptunes formed by pebble accretion result in small rocky cores surrounded by polluted H,He envelopes where most of the rock (silicate) is in vapor form at the end of the formation phase. This vapor is expected to condense and rain-out as the planet cools. In this Letter we examine the timescale for the rainout and its effect on the thermal evolution. We calculate the thermal and structural evolution of a 10 Earth masses planet formed by pebble accretion, taking into account material redistribution from silicate rainout (condensation and settling) and from convective mixing. We find that the duration of the rainout in sub-Neptunes is on Gyr timescale and varies with envelope mass: planets with envelopes below 0.75 Earth mass rainout into a core-envelope structure in less than 1 Gyr, while planets in excess of 0.75 Earth mass of H,He preserve some of their envelope pollution for billions of years. The energy released by the rainout inflates the radius with respect to planets that start out from a plain core-envelope structure. This inflation would result in estimates of the H,He contents of observed exoplanets based on the standard core-envelope structure to be too high.We identify a number of planets in the exoplanet census where rainout may operate, which would result in their H,He contents to be overestimated by up to a factor two. Future accurate age measurements by the PLATO mission may allow the identification of planets formed with polluted envelopes.

We report on an analysis of the high-energy muon component in near-vertical extensive air showers detected by the surface array IceTop in coincidence with the in-ice array of the IceCube Neutrino Observatory. In the coincidence measurement, the predominantly electromagnetic signal measured by IceTop is used to estimate the cosmic-ray primary energy, and the energy loss of the muon bundle in the deep in-ice array is used to estimate the number of muons in the shower with energies above 500 GeV ("TeV muons"). The average multiplicity of these TeV muons is determined for cosmic-ray energies between 2.5 PeV and 100 PeV assuming three different hadronic interaction models: Sibyll 2.1, QGSJet-II.04, and EPOS-LHC. For all models considered, the results are found to be in good agreement with the expectations from simulations. A tension exists, however, between the high-energy muon multiplicity and other observables; most importantly the density of GeV muons measured by IceTop using QGSJet-II.04 and EPOS-LHC.

FRB 20210912A is a fast radio burst (FRB), detected and localised to sub-arcsecond precision by the Australian Square Kilometre Array Pathfinder. No host galaxy has been identified for this burst despite the high precision of its localisation and deep optical and infrared follow-up, to 5-$\sigma$ limits of $R=26.7$ mag and $K_\mathrm{s}=24.9$ mag with the Very Large Telescope. The combination of precise radio localisation and deep optical imaging has almost always resulted in the secure identification of a host galaxy, and this is the first case in which the line-of-sight is not obscured by the Galactic disk. The dispersion measure of this burst, $\mathrm{DM_{FRB}}=1233.696\pm0.006~\mathrm{pc}\ \mathrm{cm}^{-3}$, allows for a large source redshift of $z>1$ according to the Macquart relation. It could thus be that the host galaxy is consistent with the known population of FRB hosts, but is too distant to detect in our observations ($z>0.7$ for a host like that of the first repeating FRB source, FRB 20121102A); that it is more nearby with a significant excess in $\mathrm{DM_{host}}$, and thus dimmer than any known FRB host; or, least likely, that the FRB is truly hostless. We consider each possibility, making use of the population of known FRB hosts to frame each scenario. The fact of the missing host has ramifications for the FRB field: even with high-precision localisation and deep follow-up, some FRB hosts may be difficult to detect, with more distant hosts being the less likely to be found. This has implications for FRB cosmology, in which high-redshift detections are valuable.

We reconstructed dark spots in the images of supermassive black holes SgrA* and M87* provided by the Event Horizon Telescope (EHT) collaboration by using the geometrically thin accretion disk model. In this model, the black hole is highlighted by the hot accretion matter up to the very vicinity of the black hole event horizon. The existence of hot accretion matter in the vicinity of black hole event horizons is predicted by the Blandford-Znajek mechanism, which is confirmed by recent general relativistic MHD simulations in supercomputers. A dark spot in the black hole image in the described model is a gravitationally lensed image of an event horizon globe. The lensed images of event horizons are always projected at the celestial sphere inside the awaited positions of the classical black hole shadows, which are invisible in both cases of M87* and SgrA*. We used the sizes of dark spots in the images of SgrA* and M87* for inferring their spins, 0.65<a<0.9 and a>0.75, accordingly.

The IceCube Neutrino Observatory measures high energy atmospheric neutrinos with high statistics. These atmospheric neutrinos are produced in cosmic ray interactions in the atmosphere, mainly by the decay of pions and kaons. The rate of the measured neutrinos is affected by seasonal temperature variations in the stratosphere, which are expected to increase with the energy of the particle. In this contribution, seasonal energy spectra are obtained using a novel spectrum unfolding approach, the Dortmund Spectrum Estimation Algorithm (DSEA+), in which the energy distribution from 125 GeV to 10 TeV is estimated from measured quantities with machine learning algorithms. The seasonal spectral difference to the annual average flux will be discussed based on preliminary results from IceCube's atmospheric muon neutrino data.

IceCube Neutrino Observatory, the cubic kilometer detector embedded in ice of the geographic South Pole, is capable of detecting particles from several GeV up to PeV energies enabling precise neutrino spectrum measurement. The diffuse neutrino flux can be subdivided into three components: astrophysical, from extraterrestrial sources; conventional, from pion and kaon decays in atmospheric Cosmic Ray cascades; and the yet undetected prompt component from the decay of charmed hadrons. A particular focus of this work is to test the predicted angular dependence of the atmospheric neutrino flux using an unfolding method. Unfolding is a set of methods aimed at determining a value from related quantities in a model-independent way, eliminating the influence of several assumptions made in the process. In this work, we unfold the muon neutrino energy spectrum and employ a novel technique for rebinning the observable space to ensure sufficient event numbers within the low statistic region at the highest energies. We present the unfolded energy and zenith angle spectrum reconstructed from IceCube data and compare the result with model expectations and previous measurements.

We have classified the point-like X-ray sources detected by the SRG/eROSITA telescope in the deep Lockman Hole survey. The goal was to separate the sources into Galactic and extragalactic objects. In this work have used the results of our previous cross-match of X-ray sources with optical catalogs. To classify SRG/eROSITA sources we have used the flux ratio $F_{x}/F_{o}$ and information about the source optical extent. As a result, of the 6885 X-ray sources in the eROSITA catalog 357 sources have been classified as Galactic and 5929 and as extragalactic. 539 out of 6885 have been treated as hostless, i.e., having no optical counterparts in the optical catalogs under consideration. 60 have remained unclassified due to the insufficient reliability of optical photometry. Precision and recall for the extragalactic sources are 99.9 and 98.9% (respectively) and 91.6 and 99.7% for the Galactic sources. Using this classification, we have constructed the curves of cumulative number counts for the Galactic and extragalactic sources in the Lockman Hole field. The code that accompanies this paper is available at https://github.com/mbelveder/ero-lh-class.git.

We present a methodological framework for studying galaxy evolution by utilizing Graph Theory and network analysis tools. We study the evolutionary processes of local ultraluminous infrared galaxies (ULIRGs) and quasars and the underlying physical processes, such as star formation and active galactic nucleus (AGN) activity, through the application of Graph Theoretical analysis tools. We extract, process and analyse mid-infrared spectra of local (z < 0.4) ULIRGs and quasars between 5-38 microns through internally developed Python routines, in order to generate similarity graphs, with the nodes representing ULIRGs being grouped together based on the similarity of their spectra. Additionally, we extract and compare physical features from the mid-IR spectra, such as the polycyclic aromatic hydrocarbons (PAHs) emission and silicate depth absorption features, as indicators of the presence of star-forming regions and obscuring dust, in order to understand the underlying physical mechanisms of each evolutionary stage of ULIRGs. Our analysis identifies five groups of local ULIRGs based on their mid-IR spectra, which is quite consistent with the well established fork classification diagram by providing a higher level classification. We demonstrate how graph clustering algorithms and network analysis tools can be utilized as unsupervised learning techniques for revealing direct or indirect relations between various galaxy properties and evolutionary stages, which provides an alternative methodology to previous works for classification in galaxy evolution. Additionally, our methodology compares the output of several graph clustering algorithms in order to demonstrate the best-performing Graph Theoretical tools for studying galaxy evolution.

We critically examine the state of current constraints on the dark energy (DE) equation of state (EoS) $w$. Our study is partially motivated by the observation that, while broadly consistent with the cosmological constant value $w=-1$, several independent probes appear to point towards a slightly phantom EoS ($w \sim -1.03$). We pay attention to the apparent preference for phantom DE from Planck Cosmic Microwave Background (CMB) data alone, whose origin we study in detail and attribute to a wide range of (physical and geometrical) effects. We deem the combination of Planck CMB, Baryon Acoustic Oscillations, Type Ia Supernovae, and Cosmic Chronometers data to be particularly trustworthy, inferring from this final consensus dataset $w=-1.013^{+0.038}_{-0.043}$, in excellent agreement with the cosmological constant value. Overall, despite a few scattered hints, we find no compelling evidence forcing us away from the cosmological constant (yet).

In this work, we measured the polarization properties of the X-rays emitted from the X-ray tubes, which were used during the calibration of the instrument onboard Imaging X-ray Polarimetry Explorer (IXPE). X-ray tubes are used as a source of unpolarized X-rays to calibrate the response of the gas pixel detectors to unpolarized radiation. However, even though the characteristic fluorescent emission lines are unpolarized, continuum bremsstrahlung emission can be polarized based on the geometry of the accelerated electrons and emitted photons. Hence, characterizing the contribution of polarized X-rays from bremsstrahlung emission is of interest, also for future measurements. We find that when accelerated electrons are parallel to the emitted photons, the bremsstrahlung emission is unpolarized, and when they are perpendicular, the polarization increases with energy, as expected from the theoretical predictions. A comparison with the theoretical predictions is also shown.

We investigate the discrepancy between the predicted size of accretion disks (ADs) in quasars and the observed sizes as deduced from gravitational microlensing studies. Specifically, we aim to understand whether the discrepancy is due to an inadequacy of current AD models or whether it can be accounted for by the contribution of diffuse broad-line region (BLR) emission to the observed continuum signal. We employed state-of-the-art emission models for quasars and high-resolution microlensing magnification maps and compared the attributes of their magnification-distribution functions to those obtained for pure Shakura-Sunyaev disk models. We tested the validity of our detailed model predictions by examining their agreement with published microlensing estimates of the half-light radius of the continuum-emitting region in a sample of lensed quasars. Our findings suggest that the steep disk temperature profiles found by microlensing studies are erroneous as the data are largely affected by the BLR, which does not obey a temperature-wavelength relation. We show with a sample of 12 lenses that the mere contribution of the BLR to the continuum signal is able to account for the deduced overestimation factors as well as the implied size-wavelength relation. Our study points to a likely solution to the AD size conundrum in lensed quasars, which is related to the interpretation of the observed signals rather than to disk physics. Our findings significantly weaken the tension between AD theory and observations, and suggest that microlensing can provide a new means to probe the hitherto poorly constrained diffuse BLR emission around accreting black holes.

Recently another long period radio pulsar GPM J1839$-$10 is reported, similar to GLEAM-X J162759.5$-$523504.3. Previously, the energy budget and rotational evolution of long period radio pulsars had been considered. This time, the death line and pulse width for neutron star and white dwarf pulsars are investigated. The pulse width is included as the second criterion for neutron star and white dwarfs pulsars. It is found that: (1) PSR J0250+5854 and PSR J0901$-$4046 etc should be normal radio pulsars. They have narrow pulse width and they lie near the radio emission death line. (2) The two long period radio pulsars GLEAM-X J162759.5$-$523504.3 and GPM J1839$-$10 is unlikely to be normal radio pulsars. Their possible pulse width is relatively large. And they lie far below the fiducial death line on the $P-\dot{P}$ diagram. (3) GLEAM-X J162759.5$-$523504.3 and GPM J1839$-$10 may be magnetars or white dwarf radio pulsars. At present, there are many parameters and uncertainties in both of these two possibilities.

As a fundamental parameter for modern cosmology, the Hubble constant $H_0$ is experiencing a serious crisis. In this paper, we explore an independent approach to measure $H_0$ based on the time-delay cosmography with strong gravitational lensing of a quasar by a galaxy cluster. Specifically we focus on the strong lensing system SDSS J1004+4112 with the maximum image separation of 14.62$''$, the first system of a quasar lensed by a galaxy cluster with five multiple images. Incorporating the latest time-delay measurements, we investigate the lens model dependence from the combination of 16 different lens mass models. We find that the lens model dependence is indeed large, with the combined measurement of the Hubble constant of $H_0=67.5^{+14.5}_{-8.9}km/s/Mpc$ that is obtained by summing posteriors of the Hubble constant from the 16 models with equal weighting. Interestingly, our results show that the value of Hubble constant decreases as the complexity of the perturbation around the lens increases, although weighting based on positional errors of quasar images does not significantly improve the $H_0$ constraint. We find that the 16 different mass models predict largely different shapes of the lensed quasar host galaxy as well as other lensed galaxies behind the cluster. By selecting two mass models that best reproduces those shapes, the constraint on the Hubble constant is significantly tightened to $H_0=59.1^{+3.6}_{-3.5}km/s/Mpc$. While we caution that our analysis still does not fully explore all the possible mass model uncertainty, our results highlight the importance of including as many constraints as possible such as extended shapes of lensed galaxies for obtaining tight constraints on the Hubble constant from cluster-lensed quasar lens systems.

IceCube has discovered a flux of astrophysical neutrinos and presented evidence for the first neutrino sources, a flaring blazar known as TXS 0506+056 and the active galaxy NGC 1068. However, the sources responsible for the majority of the astrophysical neutrino flux remain elusive. In addition to hypothetical sources within our Galaxy, high energy neutrinos are produced when cosmic rays interact at their acceleration sites and during propagation through the interstellar medium. The Galactic plane has therefore long been hypothesized as a neutrino source. In this contribution, new results are presented for searches of neutrino sources utilizing a dataset that builds upon recent advances in deep-learning-based reconstruction methods for neutrino-induced cascades. This work presents the first observation of high-energy neutrinos from the Milky Way Galaxy, rejecting the background-only hypothesis at 4.5$\sigma$. The neutrino signal is consistent with diffuse emission from the Galactic plane, potentially in combination with emission by a population of sources.

In the 2030s, a new era of gravitational-wave (GW) observations will dawn as multiple space-based GW detectors, such as the Laser Interferometer Space Antenna, Taiji and TianQin, open the millihertz window for GW astronomy. These detectors are poised to detect a multitude of GW signals emitted by different sources. It is a challenging task for GW data analysis to recover the parameters of these sources at a low computational cost. Generally, the matched filtering approach entails exploring an extensive parameter space for all resolvable sources, incurring a substantial cost owing to the generation of GW waveform templates. To alleviate the challenge, we make an attempt to perform parameter inference for coalescing massive black hole binaries (MBHBs) using deep learning. The model trained in this work has the capability to produce 50,000 posterior samples for redshifted total mass, mass ratio, coalescence time and luminosity distance of a MBHB in about twenty seconds. Our model can serve as a potent data pre-processing tool, reducing the volume of parameter space by more than four orders of magnitude for MBHB signals with a signal-to-noise ratio larger than 100. Moreover, the model exhibits robustness when handling input data that contains multiple MBHB signals.

During a screening in ASAS-SN database searching candidates of Delta Scuti stars with short period, our attention was drawn to the variable star ASASSN-V J104912.47+274312.7, we considered interesting to follow. It is an ASAS-SN discovery, classified by them as a RRc. We observed it for 16 nights in 2021, obtaining several maxima that allow us to refine its period. After a frequency analysis, we conclude that this star belongs to a small subgroup of RRc stars having the peculiarity to show two close frequencies with ratio greater than 0.96, without 0.61 or 0.68 period ratio mode. This behaviour is very similar to two stars found by Netzel & Smolec (2019) in OGLE data, and one star in NGC 6362, studied by Smolec et al. (2017). This large period ratios are also found as additional mode to standard double mode RRc stars (Px/Pl= 0.61-0.68) by Netzel et al. (2023).

We study the cross-correlation between maps of the unresolved $\gamma$-ray background constructed from the 12-year data release of the Fermi Large-Area Telescope, and the overdensity of galaxies in the redshift range $z\lesssim0.4$ as measured by the 2MASS Photometric Redshift survey and the WISE-SuperCOSMOS photometric survey. A signal is detected at the $8-9\sigma$ level, which we interpret in terms of both astrophysical $\gamma$-ray sources, and WIMP dark matter decay and annihilation. The sensitivity achieved allows us to characterise the energy and redshift dependence of the signal, and we show that the latter is incompatible with a pure dark matter origin. We thus use our measurement to place an upper bound on the WIMP decay rate and the annihilation cross-section, finding constraints that are competitive with those found in other analyses. Our analysis is based on the extraction of clean model-independent observables that can then be used to constrain arbitrary astrophysical and particle physics models. In this sense we produce measurements of the $\gamma$-ray emissivity as a function of redshift and rest-frame energy $\epsilon$, and of a quantity $F(\epsilon)$ encapsulating all WIMP parameters relevant for dark matter decay or annihilation. We make these measurements, together with a full account of their statistical uncertainties, publicly available.

The study of quasi-periodic oscillations (QPOs) and power density spectra (PDS) continuum properties can help shed light on the still illusive emission physics of magnetars and as a window into the interiors of neutron stars using asteroseismology. In this work, we employ a Bayesian method to search for the QPOs in the hundreds of X-ray bursts from SGR J1935+2154 observed by {\it Insight}-HXMT, GECAM and Fermi/GBM from July 2014 to January 2022. Although no definitive QPO signal (significance $>3\sigma$) is detected in individual bursts or the averaged periodogram of the bursts grouped by duration, we identify several bursts exhibiting possible QPO at $\sim$ 40 Hz, which is consistent with that reported in the X-ray burst associated with FRB 200428. We investigate the PDS continuum properties and find that the distribution of the PDS slope in the simple power-law model peaks $\sim$ 2.5, which is consistent with other magnetars but higher than 5/3 commonly seen in gamma-ray bursts. Besides, the distribution of the break frequency in the broken power-law model peaks at $\sim$ 60 Hz. Finally, we report that the power-law index of PDS has an anti-correlation and power-law dependence on the burst duration as well as the minimum variation timescale.

The millisecond pulsar(MSP) is believed to be an old neutron star(NS) having undergone spin-up by the accreting material from the donor. Whereas, the discovery of eccentric millisecond pulsars (eMSPs) in the Galactic field challenges such a scenario producing MSP-white dwarf (WD) only in the circular orbit. As orbital periods and companion mass of these eMSPs are located in a narrow range, a reasonable postulation is that they have the same origin. Although many models have been proposed to interpret their origin, however, the origin of the narrow range of the orbital period is still an open question. The accretion-induced collapse(AIC) of the ONe WD is considered to be an important pathway to form MSP, which was expected to result in the formation of MSP in the circular orbit due to tidal circularization. Here we revisited this scenario by the binary population synthesis including the specific circularization calculation. Our results indicate that binaries with insufficient circularization in this scenario can evolve into the eMSPs. The narrow initial binary parameters required by insufficient circularization can naturally account for the narrow range of the orbital period. Although the evolution of WD's AIC process has not been well understood, the characteristic of a narrow range in the orbital period of eMSPs can still set constraints on the physics of their evolution.

The existence of giant planets on wide orbits ($\stackrel{>}{_\sim}100$AU) challenge planet formation theories; the core accretion scenario has difficulty in forming them, whereas the disc instability model forms an overabundance of them that is not seen observations. We perform $N$-body simulations investigating the effect of close stellar encounters ($\leq 1200$AU) on systems hosting wide-orbit giant planets and the extent at which such interactions may disrupt the initial wide-orbit planet population. We find that the effect of an interaction on the orbit of a planet is stronger for high-mass, low-velocity perturbers, as expected. We find that due to just a single encounter there is a $\sim 17%$ chance that the wide-orbit giant planet is liberated in the field, a $\sim 10$% chance it is scattered significantly outwards, and a $\sim 6$% chance it is significantly scattered inwards. Moreover, there is a $\sim 21\%$ chance that its eccentricity is excited to e>0.1, making it more prone to disruption in subsequent encounters. The results strongly suggest that the effect of even a single stellar encounter is significant in disrupting the primordial wide-orbit giant planet population; in reality the effect will be even more prominent, as in a young star-forming region more such interactions are expected to occur. We conclude that the low occurrence rate of wide-orbit planets revealed by observational surveys does not exclude the possibility that such planetary systems are initially abundant, and therefore the disc-instability model may be a plausible scenario for their formation.

The recently released data by pulsar timing array (PTA) collaborations have amassed substantial evidence corroborating the existence of a stochastic signal consistent with gravitational-wave background at frequencies around the nanoHertz regime. We investigate the situation that the PTA signal originates from scalar-induced gravitational waves~(SIGWs), which serves as a valuable tool to probe the equation of state parameter $w$ during the early stages of the Universe. The joint consideration of the PTA data from the NANOGrav 15-yr data set, PPTA DR3, and EPTA DR2 yields that $w=0.60^{+0.32}_{-0.39}$ (median + $90\%$ credible interval), indicating a period of condensate domination at the production of SIGWs and the preferred $w$ reveals the dynamics of the homogeneous field. We also impose an upper bound on the reheating temperature and the constraint on $w$ reveals valuable information of the inflationary potential and the dynamics at the end of inflation.

Protostellar binaries harbour complex environment morphologies. Observations represent a snapshot in time, and projection and optical depth effects impair our ability to interpret them. Careful comparison with high-resolution models that include the larger star-forming region can help isolate the driving physical processes and give observations context in the time domain. We carry out zoom-in simulations with AU-scale resolution, and for the first time ever we follow the evolution until a circumbinary disk is formed. We investigate the gas dynamics around the young stars and extract disk sizes. Using radiative transfer, we obtain evolutionary tracers of the binary systems. We find that the centrifugal radius in prestellar cores is a poor estimator of the resulting disk size due to angular momentum transport at all scales. For binaries, the disk sizes are regulated periodically by the binary orbit, having larger radii close to the apastron. The bolometric temperature differs systematically between edge-on and face-on views and shows a high frequency time dependence correlated with the binary orbit and a low frequency time dependence with larger episodic accretion events. These oscillations can bring the system appearance to change rapidly from class 0 to class I and for short time periods even bring it to class II. The highly complex structure in early stages, as well as the binary orbit itself, affects the classical interpretation of protostellar classes and direct translation to evolutionary stages has to be done with caution and include other evolutionary indicators such as the extent of envelope material.

Star formation is a multi-scale problem, and only global simulations that account for the connection from the molecular cloud scale gas flow to the accreting protostar can reflect the observed complexity of protostellar systems. Star-forming regions are characterised by supersonic turbulence and as a result, it is not possible to simultaneously design models that account for the larger environment and in detail reproduce observed stellar systems. Instead, the stellar inventories can be matched statistically, and best matches found that approximate specific observations. Observationally, a combination of single-dish telescopes and interferometers are now able to resolve the nearest protostellar objects on all scales from the protostellar core to the inner 10 AU. We present a new non-parametric methodology which uses high-resolution simulations and post-processing methods to match simulations and observations using deep learning. Our goal is to perform a down-selection from large data sets of synthetic images to a ranked list of best-matching candidates with respect to the observation. This is particularly useful for binary and multiple stellar systems that form in turbulent environments. The objective is to accelerate the rate at which we can do such comparisons, remove biases from hand-picking matches, and contribute to identifying the underlying physical processes that drive the creation and evolution of observed protostellar systems.

We present a 2.5-dimensional magnetohydrodynamic simulation of a systematically rotating prominence inside its coronal cavity using the open-source \texttt{MPI-AMRVAC} code. Our simulation starts from a non-adiabatic, gravitationally stratified corona, permeated with a sheared arcade magnetic structure. The flux rope (FR) is formed through converging and shearing footpoints driving, simultaneously applying randomized heating at the bottom. The latter induces a left-right asymmetry of temperature and density distributions with respect to the polarity-inversion line. This asymmetry drives flows along the loops before the FR formation, which gets converted to net rotational motions upon reconnection of the field lines. As the thermal instability within the FR develops, angular momentum conservation about its axis leads to a systematic rotation of both hot coronal and cold condensed plasma. The initial rotational velocity exceeds $60\ \mathrm{km\ s^{-1}}$. The synthesized images confirm the simultaneous rotations of the coronal plasma seen in 211 and 193 \AA\ and condensations seen in 304 \AA. Furthermore, the formation of the dark cavity is evident in 211 and 193 \AA\ images. Our numerical experiment is inspired by observations of so-called giant solar prominence tornadoes, and reveals that asymmetric FR formation can be crucial in triggering rotational motions. We reproduce observed spinning motions inside the coronal cavity, augmenting our understanding of the complex dynamics of rotating prominences.

Trajectory design in highly-perturbed environments like binary asteroids is challenging. It typically requires using realistic, non-autonomous dynamical models in which periodic solutions vanish. In this work, a novel technique to find regions of bounded motion in the perturbed planar bi-elliptic restricted four-body problem is proposed. Lagrangian descriptors are employed to find regions of bounded motion about Dimorphos, the secondary body of the (65803) Didymos binary system. Results show that Lagrangian descriptors successfully reveal phase space organizing structures both in the unperturbed and perturbed planar bi-elliptic restricted four-body problem. With no solar radiation pressure, regions of bounded motion are visually identified, so granting access to a vast selection of bounded orbits about Dimorphos. Conversely, the presence of solar radiation pressure breaks down the majority of structures, leading to a large region of unstable motion with rare exceptions. Compared to other chaos indicators applied to the astrodynamics, Lagrangian descriptors are more convenient since they avoid propagating variational equations.

Solar flares and coronal mass ejections are interrelated phenomena that together are known as solar eruptive events. These are the main drivers of space weather and understanding their origins is a primary goal of Heliophysics. In this white paper, we advocate for the allocation of sufficient resources to bring together experts in observations and modeling to construct and test next generation data-driven models of solar eruptive events. We identify the key components necessary for constructing comprehensive end-to-end models including global scale 3D MHD resolving magnetic field evolution and reconnection, small scale simulations of particle acceleration in reconnection exhausts, kinetic scale transport of flare-accelerated particles into the lower solar atmosphere, and the radiative and hydrodynamics responses of the solar atmosphere to flare heating. Using this modeling framework, long-standing questions regarding how solar eruptive events release energy, accelerate particles, and heat plasma can be explored. To address open questions in solar flare physics, we recommend that NASA and NSF provide sufficient research and analysis funds to bring together a large body of researchers and numerical tools to tackle the end-to-end modeling framework that we outline. Current dedicated theory and modeling funding programs are relatively small scale and infrequent; funding agencies must recognize that modern space physics demands the use of both observations and modeling to make rapid progress.

The upcoming Extremely Large Telescopes (ELTs) are expected to have the collecting area required to detect potential biosignature gases such as molecular oxygen, $\mathrm{O_2}$, in the atmosphere of terrestrial planets around nearby stars. One of the most promising detection methods is transmission spectroscopy. To maximize our capability to detect $\mathrm{O_2}$ using this method, spectral resolutions $\mathrm{R}\geq 300,000$ are required to fully resolve the absorption lines in an Earth-like exoplanet atmosphere and disentangle the signal from telluric lines. Current high-resolution spectrographs typically achieve a spectral resolution of $\mathrm{R}\sim100,000$. Increasing the resolution in seeing limited observations/instruments requires drastically larger optical components, making these instruments even more expensive and hard to fabricate and assemble. Instead, we demonstrate a new approach to high-resolution spectroscopy. We implemented an ultra-high spectral resolution booster to be coupled in front of a high-resolution spectrograph. The instrument is based on a chained Fabry Perot array which generates a hyperfine spectral profile. We present on-sky telluric observations with a lab demonstrator. Depending on the configuration, this two-arm prototype reaches a resolution of R=250,000-350,000. After carefully modeling the prototype's behavior, we propose a Fabry Perot Interferometer (FPI) design for an eight-arm array configuration aimed at ELTs capable of exceeding R=300,000. The novel FPI resolution booster can be plugged in at the front end of an existing R=100,000 spectrograph to overwrite the spectral profile with a higher resolution for exoplanet atmosphere studies.

Photoevaporation is a potential explanation for several features within exoplanet demographics. Atmospheric escape observed in young Neptune-sized exoplanets can provide insight into and characterize which mechanisms drive this evolution and at what times they dominate. AU Mic b is one such exoplanet, slightly larger than Neptune (4.19 Earth radii). It closely orbits a 23 Myr pre-Main Sequence M dwarf with a period of 8.46 days. We obtained two visits of AU Mic b at Lyman-alpha with HST/STIS. One flare within the first HST visit is characterized and removed from our search for a planetary transit. We present a non-detection in our first visit followed by the detection of escaping neutral hydrogen ahead of the planet in our second visit. The outflow absorbed about 30% of the star's Lyman-alpha blue-wing 2.5 hours before the planet's white-light transit. We estimate the highest velocity escaping material has a column density of 10^13.96 cm^-2 and is moving 61.26 km/s away from the host star. AU Mic b's large high energy irradiation could photoionize its escaping neutral hydrogen in 44 minutes, rendering it temporarily unobservable. Our time-variable Lyman-alpha transit ahead of AU Mic b could also be explained by an intermediate stellar wind strength from AU Mic that shapes the escaping material into a leading tail. Future Lyman-alpha observations of this system will confirm and characterize the unique variable nature of its Lyman-alpha transit, which combined with modeling will tune the importance of stellar wind and photoionization.

The effective field theory of large-scale structure allows for a consistent perturbative bias expansion of the rest-frame galaxy density field. In this work, we present a systematic approach to renormalize galaxy bias and stochastic parameters using a finite cutoff scale $\Lambda$. We derive the differential equations of the Wilson-Polchinski renormalization group that describe the evolution of the finite-scale bias parameters with $\Lambda$, analogous to the $\beta$-function running in QFT. We further provide the connection between the finite-cutoff scheme and the renormalization procedure for $n$-point functions that has been used as standard in the literature so far; some inconsistencies in the treatment of renormalized bias in current EFT analyses are pointed out as well. The fixed-cutoff scheme allows us to predict, in a principled way, the finite part of loop contributions which is due to perturbative modes and which, in the standard renormalization approach, is absorbed into counterterms. We expect that this will allow for the robust extraction of (a yet-to-be-determined amount of) additional cosmological information from galaxy clustering, both when using field-level techniques and $n$-point functions.

Understanding the evolution of galaxies provides crucial insights into a broad range of aspects in astrophysics, including structure formation and growth, the nature of dark energy and dark matter, baryonic physics, and more. It is, however, infeasible to track the evolutionary processes of individual galaxies in real time given their long timescales. As a result, galaxy evolution analyses have been mostly based on ensembles of galaxies that are supposed to be from the same population according to usually basic and crude observational criteria. We propose a new strategy of evaluating the evolution of an individual galaxy by identifying its descendant galaxies as guided by cosmological simulations. As a proof of concept, we examined the evolution of the total mass distribution of a target strong lensing galaxy at $z=0.884$ using the proposed strategy. We selected 158 galaxies from the IllustrisTNG300 simulation that we identified as analogs of the target galaxy. We followed their descendants and found 11 observed strong lensing galaxies that match in stellar mass and size with the descendants at their redshifts. The observed and simulated results are discussed, although no conclusive assessment is made given the low statistical significance due to the small sample size. Nevertheless, the test confirms that our proposed strategy is already feasible with existing data and simulations. We expect it to play an even more important role in studying galaxy evolution as more strong lens systems and larger simulations become available with the advent of next-generation survey programs and cosmological simulations.

The 21-cm brightness temperature fluctuation from the Dark Ages ($z \simeq 30-100$) will allow us to probe the inflationary epoch on very small scales ($>0.1 \, \mbox{Mpc}^{-1}$), inaccessible to cosmic microwave background experiments. Combined with the possibility to collect information from different redshift slices, the 21-cm bispectrum has the potential to significantly improve constraints on primordial non-Gaussianity. However, recent work has shown secondary effects source off-diagonal terms in the covariance matrix which can significantly affect forecasted constraints, especially in signals that peak in the squeezed configuration, such as the local bispectrum. In this paper we propose the three-point $\langle 21-21-\rm CMB \rangle$ bispectrum cross-correlation as a new independent observational channel sensitive to local primordial non-Gaussianity. We find that, contrary to the 21-cm bispectrum, secondary contributions are subdominant to the primordial signal for values $f_{\rm NL}^{\rm loc} \sim 1$, resulting in negligible effects from off-diagonal terms in the covariance matrix. We forecast that CMB $T$ and $E$ modes cross-correlated with an ideal cosmic variance-limited 21-cm experiment with a $0.1$ MHz frequency and $0.1$ arc-minute angular resolution could reach $f_{\rm NL}^{\rm loc} \sim 6 \times 10^{-3}$. This forecast suggests cross-correlation between CMB and 21-cm experiments could provide a viable alternative to 21-cm auto-spectra in reaching unprecedented constraints on primordial local non-Gaussianities.

MWC 349A is one of the rare stars known to have hydrogen radio recombination line (RRL) masers. The bright maser emission makes it possible to study the dynamics of the system at milli-arcsecond (mas) precision. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 1.4 mm and 0.8 mm continuum emission of MWC 349A, as well as the H30$\alpha$ and H26$\alpha$ RRLs. Using the most extended array configuration of C43-10 with a maximum baseline of 16.2km, we spatially resolved the H30$\alpha$ line and 1.4mm continuum emission for the first time. In addition to the known H30$\alpha$ and H26$\alpha$ maser emission from a Keplerian disk at LSR velocities from -12 to 28 km s$^{-1}$ and from an ionized wind for velocities between -12 to -40 km s$^{-1}$ and 28 to 60 km s$^{-1}$, we found evidence of a jet along the polar axis at $V_{\mathrm{LSR}}$ from -85 to -40 km s$^{-1}$ and +60 to +100 km s$^{-1}$. These masers are found in a linear structure nearly aligned with the polar axis of the disk. If these masers lie close to the polar axis, their velocities could be as high as 575 km s$^{-1}$, which cannot be explained solely by a single expanding wind as proposed in B\'aez Rubio et al (2013). We suggest that they originate from a high-velocity jet, likely launched by a magnetohydrodynamic wind. The jet appears to rotate in the same direction as the rotation of the disk. A detailed radiative transfer modeling of these emissions will further elucidate the origin of these masers in the wind.

Our current best guess for a unified theory of gravitation and quantum field theory (string theory) generically predicts a set of requirements for a consistently quantized theory, the Swampland criteria. Refined versions of these criteria have recently been shown to be in mild tension with cosmological observations. We summarize the status of the current impact of and constraints on the Swampland conjectures from cosmology, and subject a variety of dark energy quintessence models to recently released cosmological datasets. We find that instead of tightening the tension, the new data allows for slightly more freedom in the Swampland criteria. We further demonstrate that if there is no theoretical argument made to prevent interactions of the moduli fields with the electromagnetic sector, a novel fine-tuning argument arises from the extremely tight current constraints on such interactions. Finally, we conclude with a cautionary tale on model-independent reconstructions of the Swampland criteria from expansion rate data.

This paper explores the detection capability of space-borne detectors to the anisotropic stochastic gravitational-wave background (SGWB) without relying on the low-frequency approximation. To assess the detection performance, we calculate the power-law integrated sensitivity (PLIS) curve. Our results demonstrate that a single detector has limited capabilities in detecting multipole moments beyond the monopole ($l=0$), quadrupole ($l=2$), and hexadecapole ($l=4$). However, when multiple detectors are combined, the presence of multiple pointing directions and the separation between detectors significantly enhance the detection capabilities for the other multipole moments. For instance, when considering the dipole ($l=1$), combining TianQin with TianQin II and LISA with TianQin significantly improves the detection sensitivity by 2-3 orders of magnitude, compared with using a single TianQin and a single LISA, respectively.

We study how light scalar fields can change the stellar landscape by triggering a new phase of nuclear matter. Scalars coupled to nucleons can develop a non-trivial expectation value at finite baryon density. This sourcing of a scalar reduces the nucleon mass and provides an additional energy density and pressure source. Under generic conditions, a new ground state of nuclear matter emerges, with striking implications for the configuration of stellar remnants. Notably, neutron stars in the new ground state can be significantly heavier than QCD equations of state currently predict. We also find hybrid stellar compositions and stable self-bound objects with sizes as small as the Compton wavelength of the scalar. We discuss several specific realizations of this scenario: the QCD axion and lighter generalizations thereof and linearly or quadratically coupled scalar fields effectively equivalent to a class of scalar-tensor modification of gravity. Lastly, we explore phenomenological signatures relevant to electromagnetic and gravitational wave observations of neutron stars, such as atypical compactness and instability gaps in radii.

White dwarfs and neutron stars are far-reaching and multi-faceted laboratories in the hunt for dark matter. We review detection prospects of wave-like, particulate, macroscopic and black hole dark matter that make use of several exceptional properties of compact stars, such as ultra-high densities, deep fermion degeneracies, low temperatures, nucleon superfluidity, strong magnetic fields, high rotational regularity, and significant gravitational wave emissivity. Foundational topics first made explicit in this document include the effect of the ``propellor phase" on neutron star baryonic accretion, and the contribution of Auger and Cooper pair breaking effects to neutron star heating by dark matter capture.

In curved spacetime, Maxwell's equations can be expressed in forms valid in Minkowski background, with the effect of the metric (gravity) appearing as effective polarizations and magnetizations. The electric and magnetic (EM) fields depend on the observer's frame four-vector. We derive Maxwell's equations valid in general curved spacetime using the fields defined in the normal frame, the coordinate frame, and two other non-covariant methods used in the literature. By analyzing the case in the generic frame we show that the EM fields, as well as the charge and current densities, defined in non-covariant ways do not correspond to physical ones measured by an observer. We show that modification of the homogeneous part is inevitable to any observer, and such a modification is difficult to interpret as the effective medium property. The normal frame is the relevant one to use as it gives the EM fields measured by an Eulerian observer.

We focus on a simple, natural and predictive T model of inflation in Supergravity employing as inflaton the Higgs field which leads to the spontaneous breaking of a U(1)_(B-L) symmetry at the SUSY GUT scale. We use a renormalizable superpotential, fixed by a U(1) R symmetry, and a Kahler potential which parameterizes the Kahler manifold SU(2,1)/(SU(2)xU(1))x(SU(2)/U(1)) with scalar curvature R_K=-6/N+2/N_0 where 0<N_0<6. The spectral index ns turns out to be close to its present central observational value and the tensor-to-scalar ratio r increases with N<36. The model can be nicely linked to MSSM offering an explanation of the magnitude of the mu parameter consistently with phenomenological data. It also allows for baryogenesis via non-thermal leptogenesis with gravitino as light as 1 TeV.

We examine the teleparallel formulation of non-minimally coupled scalar Einstein-Gauss-Bonnet gravity. In the teleparallel formulation, gravity is described by torsion instead of curvature, causing the usual Gauss-Bonnet invariant expressed through curvature to decay into two separate invariants built from torsion. Consequently, the teleparallel formulation permits broader possibilities for non-minimal couplings between spacetime geometry and the scalar field. In our teleparallel theory, there are two different branches of equations in spherical symmetry depending on how one solves the antisymmetric part of the field equations, leading to a real and a complex tetrad. We first show that the real tetrad seems to be incompatible with the regularity of the equations at the event horizon, which is a symptom that scalarized black hole solutions beyond the Riemannian Einstein-Gauss-Bonnet theory might not exist. Therefore, we concentrate our study on the complex tetrad. This leads to the emergence of scalarized black hole solutions, where the torsion acts as the scalar field source. Extending our previous work, we study monomial non-minimal couplings of degree one and two, which are intensively studied in conventional, curvature-based, scalar Einstein-Gauss-Bonnet gravity. We discover that the inclusion of torsion can potentially alter the stability of the resulting scalarized black holes. Specifically, our findings indicate that for a quadratic coupling, which is entirely unstable in the pure curvature formulation, the solutions induced by torsion may exhibit stability within certain regions of the parameter space. In a limiting case, we were also able to find black holes with a strong scalar field close to the horizon but with a vanishing scalar charge.

We study the kinetic mixing between the cosmic microwave background (CMB) photon and the birefringent dark photon as a source of cosmic birefringence. We show that indeed the birefringence of the dark photon propagates to the CMB photon, but the resulting birefringence may not be uniform over the sky. Moreover, our investigation sheds light on the essential role played by kinetic mixing in the generation of two fundamental characteristics of the CMB: circular polarization and spectral distortion.

Dark Matter, if represented by a $Z_2$-symmetric scalar field, can manifest as both particles and condensates. In this paper, we study the evolution of an oscillating homogeneous condensate of a $Z_2$-symmetric scalar field in a thermal plasma in an FLRW universe. We focus on the perturbative regime where the oscillation amplitude is sufficiently small so that parametric resonance is inefficient. This perturbative regime necessarily comprises the late stage of the condensate decay and determines its fate. The coupled coarse-grained equations of motion for the condensate, radiation, and spacetime are derived from first principles using nonequilibrium quantum field theory. We obtain analytical expressions for the relevant microscopic quantities that enter the equations of motion and solve the latter numerically. We find that there is always a nonvanishing relic abundance for a $Z_2$-symmetric condensate because its decay rate decreases faster than the Hubble parameter at late times due to either the amplitude-dependence or the temperature-dependence in the condensate decay rate. Consequently, accounting for the condensate contribution to the overall Dark Matter relic density is essential for $Z_2$ scalar singlet Dark Matter. Unlike normal thermal freeze-out for particles, the condensate relic density depends on the initial condition which we take as arbitrary in the present work provided that it falls within the perturbative regime.

We reflect on the possibility of having a matter action that is invariant only under transverse diffeomorphisms. This possibility is particularly interesting for the dark sector, where no restrictions arise based on the weak equivalence principle. In order to implement this idea we consider a scalar field which couples to gravity minimally but via arbitrary functions of the metric determinant. We show that the energy-momentum tensor of the scalar field takes the perfect fluid form when its velocity vector is time-like. We analyze the conservation of this tensor in detail, obtaining a seminal novel result for the energy density of this field in the kinetic dominated regime. Indeed, in this regime the fluid is always adiabatic and we obtain an explicit expression for the speed of sound. Furthermore, to get insight in the gravitational properties of these theories, we consider the fulfillment of the energy conditions, concluding that nontrivial physically reasonable matter violates the strong energy condition in the potential domination regime. On the other hand, we present some shift-symmetric models of particular interest. These are: constant equation of state models (which may provide us with a successful description of dark matter or dark radiation) and models presenting different gravitational domains (characterized by the focusing or possible defocusing of time-like geodesics), as it happens in unified dark matter-energy models.

We consider thin accretion disks in the field of a class of rotating regular black holes. For this purpose, we obtain the radius of the innermost stable circular orbit, $r_{ISCO}$ and efficiency of accretion disk in converting matter into radiation $\eta$ with the aim of modeling the disk's emission spectrum. We consider a simple model for the disk's radiative flux, differential and spectral luminosity and compare the results with those expected from accretion disks around Kerr black holes. As a remarkable result, from our computations we find that both the luminosity of the accretion disk and the efficiency are larger in the geometry of rotating regular black holes for fixed and small values of the spin parameter $j$ with respect to those predicted with the Kerr metric for a black hole of the same mass. These results may have interesting implications for astrophysical black holes.

Ignorance of the initial condition for the axion dynamics in the early Universe has led us to consider an $O(1)$ valued initial amplitude, and that prefers the decay constant, $F_a$, of the QCD axion to be an intermediate scale such as $10^{12}$ GeV in order to explain the dark matter abundance. We explore a cosmological scenario of $F_a$ being much larger than $10^{12}$ GeV by considering the axion and moduli dynamics during inflation to set the initial amplitude. We show that if the volume moduli (radion) of the extra-dimension is stabilized mainly by the QCD contribution to the moduli potential during inflation, the QCD axion with the string-scale decay constant obtains a mass around the inflationary Hubble parameter. This means that the axion rolls down to the $\theta = 0$ minimum during the inflation realizing almost vanishing initial amplitude, and the inflationary quantum fluctuation can be the dominant source of the current number density of axions. We find natural parameter regions where the axion explains the cold dark matter of the Universe, while the constraint on the isocurvature perturbation is avoided. The presence of the axion miniclusters or axion stars are predicted in a wide range of parameters, including the one explains the Subaru-HCS microlensing event.