Papers for Monday, Oct 23 2023

We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emission from NIRCam and Dragonfly imaging at smaller distances, we trace the surface brightness and color profiles of this galaxy over the entire extent of its predicted dark matter halo, from 0.1-100 kpc. Controlled N-body simulations have predicted that minor mergers create "accretion shelves" in the surface brightness profile at large radius. We observe such a feature in Ark 227 at 10-20 kpc, which, according to models, could be caused by a merger with total mass ratio 1:10. The metallicity declines over this radial range, further supporting the minor merger scenario. There is tentative evidence of a second shelf at $\mu_V\approx 35$ mag arcsec$^{-2}$ extending from 50-100 kpc, along with a corresponding drop in metallicity. The stellar mass in this outermost envelope is $\approx10^7M_\odot$. These results suggest that Ark 227 experienced multiple mergers with a spectrum of lower-mass galaxies -- a scenario that is broadly consistent with the hierarchical growth of structure in a cold dark matter-dominated universe. Finally, we identify an ultra-faint dwarf associated with Ark 227 with $M_\ast\approx10^5 M_\odot$ and $\mu_{V,e}=28.1$ mag arcsec$^{-2}$, demonstrating that JWST is capable of detecting very low-mass dwarfs to distances of at least ~30 Mpc.

We illustrate the formation and evolution of the Milky Way over cosmic time, utilizing a sample of 10 million red giant stars with full chemodynamical information, including metallicities and $\alpha$-abundances from low-resolution Gaia XP spectra. The evolution of angular momentum as a function of metallicity - a rough proxy for stellar age, particularly for high-[$\alpha$/Fe] stars - displays three distinct phases: the disordered and chaotic protogalaxy, the kinematically-hot old disk, and the kinematically-cold young disk. The old high-$\alpha$ disk starts at [Fe/H] $\approx -1.0$, 'spinning up' from the nascent protogalaxy, and then exhibits a smooth 'cooldown' toward more ordered and circular orbits at higher metallicities. The young low-$\alpha$ disk is kinematically cold throughout its metallicity range, with its observed properties modulated by a strong radial gradient. We interpret these trends using Milky Way analogs from the TNG50 cosmological simulation, identifying one that closely matches the kinematic evolution of our Galaxy. This halo's protogalaxy spins up into a relatively thin and misaligned high-$\alpha$ disk at early times, which is subsequently heated and torqued by a major gas-rich merger. The merger contributes a large amount of low-metallicity gas and angular momentum, from which the kinematically cold low-$\alpha$ stellar disk is subsequently born. This simulated history parallels several observed features of the Milky Way, particularly the decisive 'GSE' merger that likely occurred at $z \approx 2$. Our results provide an all-sky perspective on the emerging picture of our Galaxy's three-phase formation, impelled by the three physical mechanisms of spinup, merger, and cooldown.

Absorption of stellar X-ray and Extreme Ultraviolet radiation in the upper atmosphere of close-in exoplanets can give rise to hydrodynamic outflows, which may lead to the gradual shedding of their primordial, light element envelopes. Excess absorption by neutral helium atoms in the metastable state has recently emerged as a viable diagnostic of atmospheric escape. Here we present a public module to the 1D photo-ionization hydrodynamic code ATES, designed to calculate the HeI triplet transmission probability for a broad range of planetary parameters. By relaxing the isothermal outflow assumption, the code enables a self-consistent assessment of the HeI triplet absorption depth along with the atmospheric mass loss rate and the outflow temperature profile, which strongly affects the recombination rate of HeII into HeI triplet. We investigate how the transit signal can be expected to depend upon known system parameters, including host spectral type, orbital distance, as well as planet gravity. At variance with previous studies, which identified K-type stars as favorable hosts, we conclude that late M-dwarfs with Neptune-sized planets orbiting at ~0.05-0.1 AU can be expected to yield the strongest transit signal well in excess of 30% for near-cosmological He/H abundances. More generally, we show that the physics which regulates the population and depletion of the metastable state, combined with geometrical effects, can yield somewhat counter-intuitive results, such as a non-monotonic dependence of the transit depth on orbital distance. These are compounded by a strong degeneracy between the stellar EUV flux intensity and the atmospheric He/H abundance, both of which are highly uncertain. Compared against spectroscopy data our modelling suggests that either a large fraction of the targets have helium depleted envelopes, or, that the input stellar EUV spectra are systematically overestimated.

We apply image moment invariant analysis to total intensity and polarimetric images calculated from general relativistic magnetohydrodynamic simulations of accreting black holes. We characterize different properties of the models in our library by their invariant distributions and their evolution in time. We show that they are highly sensitive to different physical effects present in the system which allow for model discrimination. We propose a new model scoring method based on image moment invariants that is uniformly applicable to total intensity and polarimetric images simultaneously. The method does not depend on the type of images considered and its application to other non-ring like images (e.g., jets) is straight forward.

This study presents the results concerning six red giant stars members of the globular cluster NGC 6558. Our analysis utilized high-resolution near-infrared spectra obtained through the CAPOS initiative (the APOgee Survey of Clusters in the Galactic Bulge), which focuses on surveying clusters within the Galactic Bulge, as a component of the Apache Point Observatory Galactic Evolution Experiment II survey (APOGEE-2). We employ the BACCHUS (Brussels Automatic Code for Characterizing High accUracy Spectra) code to provide line-by-line elemental-abundances for Fe-peak (Fe, Ni), $\alpha$-(O, Mg, Si, Ca, Ti), light-(C, N), odd-Z (Al), and the $s$-process element (Ce) for the 4 stars with high signal-to-noise ratios. This is the first reliable measure of the CNO abundances for NGC 6558. Our analysis yields a mean metallicity for NGC 6558 of $\langle$[Fe/H]$\rangle$ = $-$1.15 $\pm$ 0.08, with no evidence for a metallicity spread. We find a Solar Ni abundance, $\langle$[Ni/Fe]$\rangle$ $\sim$ $+$0.01, and a moderate enhancement of $\alpha$-elements, ranging between $+$0.16 to $<+$0.42, and a slight enhancement of the $s$-process element $\langle$[Ce/Fe]$\rangle$ $\sim$ $+$0.19. We also found low levels of $\langle$[Al/Fe]$\rangle \sim $+$0.09$, but with a strong enrichment of nitrogen, [N/Fe]$>+$0.99, along with a low level of carbon, [C/Fe]$<-$0.12. This behaviour of Nitrogen-Carbon is a typical chemical signature for the presence of multiple stellar populations in virtually all GCs; this is the first time that it is reported in NGC 6558. We also observed a remarkable consistency in the behaviour of all the chemical species compared to the other CAPOS bulge GCs of the same metallicity.

The `core-cusp' problem has long been considered a key small scale challenge to the standard LCDM paradigm. Halos in dark matter only simulations exhibit `cuspy' (NFW-like) density profiles, where density continuously increases towards the centre. However, many observations of galaxies (particularly in the dwarf regime) indicate that their dark matter density profiles deviate strongly from this prediction, with much flatter central regions (`cores'). Here, we use NewHorizon (NH), a hydrodynamical cosmological simulation which adopts the standard LCDM cosmology, to investigate core formation and the core-cusp problem, using a statistically significant number of galaxies in a cosmological volume. In NH, halos that contain galaxies in the upper M* > 10^10.2 MSun) and lower (M* < 10^8 MSun) ends of the stellar mass distribution contain cusps. However, halos that contain galaxies with intermediate (10^8 MSun < M* < 10^10.2 MSun) stellar masses are typically cored. These halos typically have DM masses between 10^10.2 MSun and 10^11.5 MSun. Cores are created via supernova-driven gas removal from the central regions of halos, which alters the central gravitational potential, inducing dark matter to migrate to larger radii. While all massive (M* > 10^9.5 MSun) galaxies undergo a cored-phase, in some cases cores can be removed and cusps reformed. This happens if a galaxy undergoes sustained star formation at high redshift, which results in stars (which, unlike the gas, cannot be removed by baryonic feedback) dominating the central gravitational potential. After cosmic star formation peaks, the number of cores and the mass of the halos they are formed in remain constant throughout the simulation, indicating that cores are being routinely formed over cosmic time after a threshold halo mass is reached. The existence of cores is, therefore, not in tension with the standard paradigm.

We used ultra-deep observations obtained with the NIRCam aboard the James Webb Space Telescope to explore the stellar population of NGC 6440: a typical massive, obscured and contaminated globular cluster formed and orbiting within the Galactic bulge. Leveraging the exceptional capabilities of this camera, we sampled the cluster down to ~5 magnitudes below the main-sequence turn-off in the (mF115W , mF115W - mF200W ) colour-magnitude diagram. After carefully accounting for differential extinction and contamination by field interlopers, we find that the main sequence splits into two branches both above and below the characteristic knee. By comparing the morphology of the colour-magnitude diagram with a suitable set of isochrones, we argue that the upper main-sequence bi-modality is likely due to the presence of a He-enriched stellar population with a helium spread of DeltaY = 0.04. The lower main-sequence bi-modality can be attributed to variations in the abundance of water (i.e., oxygen) with Delta[O/Fe] ~ -0.4. This is the first evidence of both helium and oxygen abundance variations in a globular cluster purely based on JWST observations. These results open the window for future in-depth investigations of the multiple population phenomenon in clusters located in the Galactic bulge, which were previously unfeasible with near-UV observations, due to prohibitive reddening and crowding conditions.

The Large Array Survey Telescope (LAST) is a wide-field telescope designed to explore the variable and transient sky with a high cadence and to be a test-bed for cost-effective telescope design. A LAST node is composed of 48 (32 already deployed), 28-cm f/2.2 telescopes. A single telescope has a 7.4 deg^2 field of view and reaches a 5-sigma limiting magnitude of 19.6 (21.0) in 20s (20x20s) (filter-less), while the entire system provides a 355 deg^2 field of view. The basic strategy of LAST is to obtain multiple 20-s consecutive exposures of each field (a visit). Each telescope carries a 61 Mpix camera, and the system produces, on average, about 2.2 Gbit/s. This high data rate is analyzed in near real-time at the observatory site, using limited computing resources (about 700 cores). Given this high data rate, we have developed a new, efficient data reduction and analysis pipeline. The data pipeline includes two major parts: (i) Processing and calibration of single images, followed by a coaddition of the visit's exposures. (ii) Building the reference images and performing image subtraction and transient detection. Here we describe in detail the first part of the pipeline. Among the products of this pipeline are photometrically and astrometrically calibrated single and coadded images, 32-bit mask images marking a wide variety of problems and states of each pixel, source catalogs built from individual and coadded images, Point Spread Function (PSF) photometry, merged source catalogs, proper motion and variability indicators, minor planets detection, calibrated light curves, and matching with external catalogs. The entire pipeline code is made public. Finally, we demonstrate the pipeline performance on real data taken by LAST.

We present the analysis of microlensing event OGLE-2014-BLG-0221, a planetary candidate event discovered in 2014. The photometric light curve is best described by a binary-lens single-source model. Our light curve modeling finds two degenerate models, with event timescales of $t_\mathrm{E}\sim70$ days and $\sim110$ days. These timescales are relatively long, indicating that the discovered system would possess a substantial mass. The two models are similar in their planetary parameters with a Jupiter mass ratio of $q \sim 10^{-3}$ and a separation of $s \sim 1.1$. While the shorter timescale model shows marginal detection of a microlensing parallax signal, the longer timescale model requires a higher order effect of microlensing parallax, lens orbital motion or xallarap to explain the deviation in the light curve. However, the modeling shows significant correlation between the higher order effects and suffers the ecliptic degeneracy that results in a failure to determine the parallax parameters. Bayesian inference is used to estimate the physical parameters of the lens, revealing the lens to be either a late-type star supported by the shorter timescale model or a stellar remnant supported by the longer timescale model. If the lens is a remnant, this would be the second planet found by microlensing around a stellar remnant. Since the models predict different values for relative proper motion and source flux, future high angular resolution follow-up observations (e.g. Keck adaptive optics) are required to rule out either of the models.

We utilize the public GIZMO code to simulate twelve disc galaxies from the NIHAO suite simulated with the GASOLINE code, then compare the corresponding galaxies in the two simulations. We find that while both codes with the same initial conditions and large-scale environments can successfully produce similar disc galaxies, significant differences are still seen in many properties of the galaxies, particularly in the circumgalactic medium (CGM) environment they reside. Specifically, the thermal feedback recipe used in GASOLINE results in ubiquitous long-lasting collimated outflows, primarily driven by high-density hot interstellar medium (ISM) from the galaxy center, and inflows of gas not aligned with the outflow cools rapidly and flows towards the galactic center. In contrast, galaxies from GIZMO code do not exhibit large-scale outflows at low redshifts, but instead display quasi-virialized hot gaseous halos that arise from the strong interaction between inflow of gas and feedback driven outflow. Therefore, the origins of mass and angular momentum of the cold disc in the two simulations are quite different, even though the final morphologies of corresponding galaxies are similar at $z\sim0$. The differences in the distribution of CGM gas are mainly due to different feedback models implemented in the two codes, thus future observations of CGM provide valuable insight into the physics governing the baryon cycle in disc galaxies.

Forthcoming missions probing the absolute intensity of the CMB are expected to be able to measure spectral distortions, which are deviations from its blackbody distribution. As cosmic inflation can induce spectral distortions, these experiments offer a possibility to further test the various promising inflationary proposals, whose predictions need to be carefully determined. After numerically fitting all inflationary observables to match current observations, we compute the predicted spectral distortions of various promising single and multifield inflationary models. The predictions of single-field inflationary models display deviations between 1% and 20% with respect to the standard cosmological model in the observable window, where multi-natural and axion-monodromy inflation stand out in this respect. In the case of multifield inflation, we observe a richer structure of the power spectrum, which, in the case of so-called hybrid attractors, yields spectral distortions about 100 times more intense than the standard signal. These observations open up questions about the relation among our results and other cosmological observables that are also to be probed soon.

We investigate a variety of short warning time, terminal mitigation scenarios via fragmentation for a hypothetical impact of asteroid 2023 NT1, a Near-Earth Object (NEO) that was discovered on July 15, 2023, two days after its closest approach to Earth on July 13. The asteroid passed by Earth within ~0.25 lunar distances with a closest approach distance of ~10$^{5}$ km and speed of 11.27 km/s. Its size remains largely uncertain, with an estimated diameter range of 26 - 58 m and probable diameter estimate (weighted by the NEO size frequency distribution) of 34 m (JPL Sentry, September 12, 2023). The asteroid approached Earth from the direction of the Sun, as did both the Chelyabinsk asteroid in 2013 and comet NEOWISE in 2021. As a result, 2023 NT1 remained undetected until after its closest approach. If it had been on a collision course, it would have had an impact energy of ~1.5 Mt (assuming a spherical asteroid with the probable diameter estimate of 34 m, 2.6 g/cm$^{3}$ uniform density, and impact speed of 15.59 km/s). 2023 NT1 represents a threat that could have caused significant local damage (~3x Chelyabinsk airburst energy). We utilize the PI ("Pulverize It") method for planetary defense to model potential mitigation scenarios of an object like 2023 NT1 through simulations of hypervelocity asteroid disruption and atmospheric ground effects for the case of a terminal defense mode. Simulations suggest that PI is an effective multimodal approach for planetary defense that can operate in extremely short interdiction modes (with intercepts as short as hours prior to impact), in addition to long interdiction time scales with months to years of warning. Our simulations support the proposition that threats like 2023 NT1 can be effectively mitigated with intercepts of one day (or less) prior to impact, yielding minimal to no ground damage, using modest resources and existing technologies.

We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,000 years, yielding median mass, momentum, and energy rates of $\dot{M}$ = 10$^{-4.4}$ M$_{\odot}$ yr$^{-1}$, $\dot{P}$ = 10$^{-3.2}$ M$_{\odot}$ km/s yr$^{-1}$, and $\dot{E}$ = 1 L$_{\odot}$. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (~3$-$5$\sigma$), and no correlations above 3$\sigma$ with clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass - which has been previously described in the literature - may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15-30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force-bolometric luminosity relationship in the literature, and our data extend this relationship up to L $\geq$ 10$^6$ L$_{\odot}$ and $\dot{P}$ $\geq$ 1 M$_{\odot}$ km/s yr$^{-1}$. Our lack of correlation with clump L/M is inconsistent with models of protocluster formation in which all protostars start forming at the same time.

Observational indications of wind-mass transfer from an evolved giant to its distant white dwarf (WD) companion in symbiotic binaries are rare. Here, we present a way to examine the neutral wind from the giant in symbiotic binaries, which is temporarily observable throughout the orbital plane during outbursts. We find that the mass-loss rate from giants in the orbital plane of S-type symbiotic binaries is high, indicating a high wind-mass-transfer efficiency in these systems. We modeled hydrogen column densities in the orbital plane between the observer and the WD for all suitable eclipsing S-type symbiotic binaries during outbursts in any orbital phase. The mass-loss rate from the giant in the orbital plane is on the order of 10$^{-6}\,M_{\odot}\,{\rm yr}^{-1}$, which is a factor of $\sim$10 higher than rates derived from nebular emission produced by the ionized wind from normal giants in symbiotic stars. This finding suggests a substantial focusing of the giant's wind toward the orbital plane and, thus, its effective transfer onto the WD companion. Our finding suggests that wind focusing on the orbital plane may be a common property of winds from giants in S-type symbiotic stars. Such wind-focusing resolves a long-standing problem of the large energetic output from their burning WDs and deficient fueling by the giant via a standard Bondi-Hoyle accretion. It also allows the WD to grow faster in mass, which lends support to the possibility that S-type symbiotic binaries are progenitors of Type Ia supernovae.

We explore neutrino emission from nonrotating, single star models across six initial metallicities and seventy initial masses from the zero-age main sequence to the final fate. Overall, across the mass spectrum, we find metal-poor stellar models tend to have denser, hotter and more massive cores with lower envelope opacities, larger surface luminosities, and larger effective temperatures than their metal-rich counterparts. Across the mass-metallicity plane we identify the sequence (initial CNO $\rightarrow$ $^{14}$N $\rightarrow$ $^{22}$Ne $\rightarrow$ $^{25}$Mg $\rightarrow$ $^{26}$Al $\rightarrow$ $^{26}$Mg $\rightarrow$ $^{30}$P $\rightarrow$ $^{30}$Si) as making primary contributions to the neutrino luminosity at different phases of evolution. For the low-mass models we find neutrino emission from the nitrogen flash and thermal pulse phases of evolution depend strongly on the initial metallicity. For the high-mass models, neutrino emission at He-core ignition and He-shell burning depends strongly on the initial metallicity. Anti-neutrino emission during C, Ne, and O burning shows a strong metallicity dependence with $^{22}$Ne($\alpha$,$n$)$^{25}$Mg providing much of the neutron excess available for inverse-$\beta$ decays. We integrate the stellar tracks over an initial mass function and time to investigate the neutrino emission from a simple stellar population. We find average neutrino emission from simple stellar populations to be 0.5--1.2 MeV electron neutrinos. Lower metallicity stellar populations produce slightly larger neutrino luminosities and average $\beta$ decay energies. This study can provide targets for neutrino detectors from individual stars and stellar populations. We provide convenient fitting formulae and open access to the photon and neutrino tracks for more sophisticated population synthesis models.

This dissertation is written as an annotated collection of selected articles. The dissertation focuses on the modelling of accretion disks in X-ray binaries with a black hole or neutron star. The main objective is to use advanced numerical methods to reveal the fundamental processes that influence the observed spectral and temporal features. Among the significant results of this modelling is the recognition of the puffy accretion disk, a novel type of accretion disk based purely on the results of numerical simulations. Furthermore, the work focuses on modelling the X-ray variability and quasi-periodic oscillations using analytical models of accretion disks, but an advanced description of the oscillations, and a discussion of the implications of this modelling for observational data.

Context. The location of the Solar System constrains the detection of extragalactic sources beyond the Milky Way(MW) plane. The optical observations are hampered in the so--called Zone of Avoidance (ZOA) where stellar crowding and Galactic absorption are severe. Observations at longer wavelengths are needed to discover new background galaxies and complete the census in the ZOA. Aims. The goal of this work is to identify galaxies behind the MW Bulge using near-infrared (NIR) data from the VISTA Variables in V\'ia L\'actea (VVV) survey. Methods. To this end we made use of different VISTA Science Archive (VSA) tools in order to extract relevant information from more than 32 billion catalogued sources in the VVV Bulge region. We find that initial photometric restriction on sources from VSA \texttt{vvvSource} table combined with restrictions on star--galaxy separation parameters obtained from Source Extractor is a successful strategy to achieve acceptable levels of contamination (60\%) and a high completeness (75\%) in the construction of a galaxy target sample. To decontaminate the initial target sample from Galactic sources our methodology also incorporates the visual inspection of false colour RGB images, a crucial quality control which was carried out following a specifically defined procedure. Results. Under this methodology we find 14480 galaxy candidates in the VVV Bulge region making this sample the largest catalogue to date in the ZOA. Moreover these new sources provide a \textbf{fresh} picture of the Universe hidden behind the curtain of stars, dust and gas in the unexplored MW Bulge region. Conclusions. The results from this work further demonstrate the potential of the VVV/VVVX survey to find and study a large number of galaxies and extragalactic structures obscured by the MW, enlightening our knowledge of the Universe in this challenging and impressive region of the sky.

The Neptunian desert is a region in period-radius parameter space with very few Neptune-sized planets at short orbital periods. Amongst these, LTT 9779 b is the only known Neptune with a period shorter than one day to retain a significant H-He atmosphere. If the Neptune desert is the result of X-ray/EUV-driven photoevaporation, it is surprising that the atmosphere of LTT 9779 b survived the intense bombardment of high energy photons from its young host star. However, the star has low measured rotational broadening, which points to the possibility of an anomalously slow spin period and hence a faint X-ray emission history that may have failed to evaporate the planet's atmosphere. We observed LTT 9779 with XMM-Newton and measured an upper limit for its X-ray luminosity that is a factor of fifteen lower than expected for its age. We also simulated the evaporation past of LTT 9779 b and found that the survival of its atmosphere to the present day is consistent with an unusually faint XUV irradiation history that matches both the X-ray and rotation velocity measurements. We conclude that the anomalously low X-ray irradiation of the one Neptune seen to survive in Neptunian desert supports the interpretation of the desert as primarily a result of photoevaporation.

We present an analysis of the JWST/MIRI-MRS spectrum of Sz114, an accreting M5 star surrounded by a large dust disk with a gap at ~39au. The spectrum is molecular-rich: we report the detection of water, CO, CO2, HCN, C2H2, and H2. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water lines with the 17.22 micron emission surpassing that of most late M-star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, flux rations of C2H2/H2O and HCN/H2O in Sz114 also resemble those of earlier-type disks, with a slightly elevated CO2/H2O ratio. While accretion heating can boost all infrared lines, the unusual properties of Sz114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and/or the presence of dust substructures. The latter delay the inward drift of icy pebbles and thus preserve a lower C/O ratio over an extended period. In contrast, late M-star disks-which are typically faint, small in size, and likely lack significant substructures-may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models.

We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses of $10^{9-12}\ {\rm M_{\odot}}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes for the stellar $\overline{\gamma}_*$, dark $\overline{\gamma}_{\rm DM}$, and total $\overline{\gamma}_{\rm T}$ mass from Zhu et al. (2023b). The $\overline{\gamma}_{\rm T}$ approaches a constant value of 2.2 for high $\sigma_{\rm e}$ galaxies, and flattens for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\lesssim2.3$, reaching 1.5 for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\approx1.8$. The total and stellar slopes track each other tightly, with $\overline{\gamma}_{\rm T}\approx\overline{\gamma}_*-0.174$ over the full $\sigma_{\rm e}$ range. This confirms the dominance of stellar matter within $R_{\rm e}$. We also show that there is no perfect conspiracy between baryonic and dark matter, as $\overline{\gamma}_*$ and $\overline{\gamma}_{\rm DM}$ do not vary inversely within the $\sigma_{\rm e}$ range. We find that the central galaxies from TNG50 and TNG100 simulations do not reproduce the observed galaxy mass distribution, which we attribute to the overestimated dark matter fraction, possibly due to a constant IMF and excessive adiabatic contraction effects in the simulations. Finally, we present the stacked dark matter density profiles and show that they are slightly steeper than the pure dark matter simulation prediction of $\overline{\gamma}_{\rm DM}\approx1$, suggesting moderate adiabatic contraction in the central region of galaxies. Our work demonstrate the power of stellar dynamics modelling for probing the interaction between stellar and dark matter and testing galaxy formation theories.

Observations close to the diffraction limit, with high Strehl ratios from Adaptive Optics (AO)-assisted instruments mounted on ground-based telescopes are a reality and will become even more widespread with the next generation instruments that equip 30 meter-class telescopes. This results in a growing interest in tools and methods to accurately reconstruct the observed Point Spread Function (PSF) of AO systems. We will discuss the performance of the PSF reconstruction (PSF-R) software developed in the context of the MICADO instrument of the Extremely Large Telescope. In particular, we have recently implemented a novel algorithm for reconstructing off-axis PSFs. In every case, the PSF is reconstructed from AO telemetry, without making use of science exposures. We will present the results coming from end-to-end simulations and real AO observations, covering a wide range of observing conditions. Specifically, the spatial variation of the PSF has been studied with different AO-reference star magnitudes. The reconstructed PSFs are observed to match the reference ones with a relative error in Strehl ratio and full-width at half maximum below 10% over a field of view of the order of one arcmin, making the proposed PSF-R method an appealing tool to assist observation analysis, and interpretation.

The growth of large-scale structure, together with the geometrical information of cosmic expansion history and cosmological distances, can be used to obtain constraints on the spatial curvature of the universe that probes the early universe physics, whereas modeling the nonlinear growth in a nonflat universe is still challenging due to computational expense of simulations in a high-dimensional cosmological parameter space. In this paper, we develop an approximate method to compute the halo-matter and halo-auto power spectra for nonflat $\Lambda$CDM model, from quantities representing the nonlinear evolution of the corresponding flat $\Lambda$CDM model, based on the separate universe (SU) method. By utilizing the fact that the growth response to long-wavelength fluctuations (equivalently the curvature), $T_{\delta_{\rm b}}(k)$, is approximated by the response to the Hubble parameter, $T_h(k)$, our method allows one to estimate the nonlinear power spectra in a nonflat universe efficiently from the power spectra of the flat universe. We use $N$-body simulations to show that the estimator can provide the halo-matter (halo-auto) power spectrum at $\sim 1\%$ ($\sim 2\%$ ) accuracy up to $k \simeq 3 (1) \, h {\rm Mpc}^{-1}$ even for a model with large curvature $\Omega_K = \pm 0.1$. Using the estimator we can extend the prediction of the existing emulators such as Dark Emulator to nonflat models without degrading their accuracy. Since the response to long-wavelength fluctuations is also a key quantity for estimating the super sample covariance (SSC), we discuss that the approximate identity $T_{\delta_{\rm b}}(k) \approx T_h(k)$ can be used to calculate the SSC terms analytically.

We study a modification of the primordial black hole (PBH) formation model from axion bubbles. We assume that the Peccei-Quinn scalar rolls down in the radial direction from a large field value to the potential minimum during inflation, which suppresses the axion fluctuations and weakens the clustering of PBHs on large scales. We find that the modified model can produce a sufficient number of PBHs that seed the supermassive black holes while avoiding the observational constraints from isocurvature perturbations and angular correlation of the high-redshift quasars.

The 21-cm line emitted by neutral hydrogen is the most promising probe of the Epoch of Reionization (EoR). Multiple radio interferometric instruments are on the cusp of detecting its power spectrum. It is therefore essential to deliver robust theoretical predictions, enabling sound inference of the coeval Universe properties. The nature of this signal traditionally required the modeling of $\mathcal{O}(10^{7-8} \, {\rm Mpc}^3)$ volumes in order to suppress the impact of cosmic variance. However, the recently-proposed Fixed & Paired (F&P) approach (Pontzen et al. 2016) uses carefully-crafted simulation pairs to achieve equal results in smaller volumes. In this work, we thoroughly test the applicability of and improvement granted by this technique to different observables of the 21-cm signal from the EoR. We employ radiation-magneto-hydrodynamics simulations to ensure the most realistic physical description of this epoch, greatly improving over previous studies using a semi-numerical approach without accurate galaxy formation physics and radiative transfer. We estimate the statistical improvement granted by the F&P technique on predictions of the skewness, the power spectrum, the bispectrum and ionized regions size distribution of the 21-cm signal at redshift $7 \leq z \leq 10$. We find that the effective volume of F&P simulations is at least 3.5 times larger than in traditional simulations. This directly translates into an equal improvement in the computational cost (both in terms of time and memory). Finally, we confirm that a combination of different observables like skewness, power spectrum and bispectrum across different redshifts can be utilised to maximise the improvement brought by the F&P technique.

The Committee on Radio Astronomy Frequencies (CRAF) is an Expert Committee of the European Science Foundation. It aims to provide a cost-effective single voice on frequency protection issues for European radio astronomy observatories and research institutes, achieving a significantly greater impact than that achievable by individual national institutions. By working together, European observatories and institutes can profit from synergy effects, cover many more topics, and learn from each other. CRAF was founded in 1988 and has since then been engaged with the International Telecommunication Union (ITU), in particular its Radiocommunication Sector (ITU-R), and the European Conference of Postal and Telecommunications Administrations (CEPT) and its European Communications Committee (ECC). This is the self-evaluation report prepared by CRAF for its periodic review of the years 2011-2021.

The multi planet system HR8799 is the first target observed with MIRI's coronagraphs as part of the MIRI-EC Guaranteed Time Observations exoplanets programme in Nov. 2022. We obtained deep observations in three coronagraphic filters from 10 to 15mic (F1065C, F1140C, F1550C), and one standard imaging filter at 20 mic (F2100W), with the goal to extract the photometry of the four planets, as well as to detect and investigate the distribution of circumstellar dust. Using dedicated observations of a reference star, we tested several algorithms to subtract the stellar diffraction pattern while preserving the fluxes of planets, which can be significantly affected by over-subtraction. Measuring correctly the planet's flux values requires accounting for the attenuation by the coronagraphs as a function of their position, and to estimate the normalisation with respect to the central star. We tested several procedures to derive averaged photometric values and error bars. These observations have enabled us to obtain two main results. First of all, the four planets in the system are well recovered, and their mid-IR fluxes, combined with near-IR flux values from the literature, are compared to two exoplanet atmosphere models, ATMO and Exo-REM. As a main outcome, the MIRI photometric data points imply larger radii (0.86 or 1.07 RJ for planet b) and cooler temperatures (950 or 1100 K for planet b), especially for planet b, in better agreement with evolutionary models. Second of all, these JWST/MIRI coronagraphic data also deliver the first spatially resolved detection of the inner warm debris disk, the radius of which is constrained to about 15 au, with flux densities comparable, but lower than former unresolved spectroscopic measurements with Spitzer. abridged...

We present the analyses of intense X-ray flares detected on the active fast rotator AB Dor using observations from the XMM-Newton. A total of 21 flares are detected, and 13 flares are analysed in detail. The total X-ray energy of these flares is found to be in the range of 10$^{34-36}$ erg, in which the peak flare flux increased up to 34 times from the pre-/post-flaring states for the strongest observed flare. The duration of these flaring events is found to be 0.7 to 5.8 hrs. The quiescent state X-ray spectra are found to be explained by a three-temperature plasma with average temperatures of 0.29, 0.95, and 1.9 keV, respectively. The temperatures, emission measures, and abundances are found to be varying during the flares. The peak flare temperature was found in the 31-89 MK range, whereas the peak emission measure was 10$^{52.5-54.7}$ $cm^{-3}$. The abundances vary during the flares and increase by a factor of $\sim$3 from the quiescent value for the strongest detected flare. The variation in individual abundances follows the inverse-FIP effect in quiescent and flare phases. The X-ray light curves of AB Dor are found to exhibit rotational modulation. The semi-loop lengths of the flaring events are derived in the range of 10$^{9.9-10.7}$ cm, whereas the minimum magnetic field to confine the plasma in the flaring loop is estimated between 200 and 700 G.

Strong gravitational lenses provide unique laboratories for cosmological and astrophysical investigations, but they must first be discovered - a task that can be met with significant contamination by other astrophysical objects and asterisms. Here we review strong lens searches, covering various sources (quasars, galaxies, supernovae, FRBs, GRBs, and GWs), lenses (early- and late-type galaxies, groups, and clusters), datasets (imaging, spectra, and lightcurves), and wavelengths. We first present the physical characteristics of the lens and source populations, highlighting relevant details for constructing targeted searches. Search techniques are described based on the main lensing feature that is required for the technique to work, namely one of: (i) an associated magnification, (ii) multiple spatially-resolved images, (iii) multiple redshifts, or (iv) a non-zero time delay between images. To use the current lens samples for science, and for the design of future searches, we list several selection biases that exist due to these discovery techniques. We conclude by discussing the future of lens searches in upcoming surveys and the new population of lenses that will be discovered.

Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 planetary system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a large sub-Neptune ($P_b=10.9$ days, $R_{p,b}=2.62\pm 0.11\, \mathrm{R}_{\oplus}$) whose orbit is interior to that of the system's smaller planet ($P_c=18.8$ days, $R_{p,c}=2.13\pm 0.12\, \mathrm{R}_{\oplus}$). Here we present revised planetary radii based on new TESS and diffuser-assisted ground-based transit observations and characterize both planetary masses using a set of 145 radial velocity (RV) measurements from HARPS-N ($M_{p,b}=4.23\pm 0.69\, \mathrm{M}_{\oplus}, M_{p,c}=2.88\pm 0.80\, \mathrm{M}_{\oplus}$). Our RV analysis also reveals a third planet candidate ($P_d=32.3$ days, $M_{p,d}\sin{i} = 4.59^{+0.96}_{-0.94}\, \mathrm{M}_{\oplus}$), which if real, would form a chain of 5:3 period ratios, although we show that the system is likely not in a mean motion resonance. Our results indicate that TOI-1266 b and c are the lowest density sub-Neptunes known around M dwarfs and may exhibit distinct bulk compositions of a gas-enveloped terrestrial ($X_{\mathrm{env},b}=5.5\pm 0.7$%) and a water-rich world (WMF$_c=59\pm 14$%), which is supported by hydrodynamic escape models. If distinct bulk compositions for the transiting planets are confirmed by atmospheric characterization, the system's unique architecture would represent an interesting test case for sub-Neptune formation models such as inside-out planet formation at pebble traps.

In 2015, the Planck Collaboration released an all-sky map of the thermal Sunyaev-Zeldovich (SZ) effect, obtained by implementing the Needlet Internal Linear Combination (NILC) method on the Planck PR2 data. The quality of the Planck data has significantly improved since then. The Planck PR4 data release offers upgraded full-sky maps in the LFI and HFI frequency bands with improved systematics and sensitivity. We present a new all-sky thermal SZ Compton $y$-parameter map derived from the Planck PR4 data using NILC and highlight improvements, particularly in noise reduction and handling residual foreground contamination. The PR4 NILC Compton $y$-parameter map has been made publicly available to support further analyses.

It has been observed that several Milky Way (MW) satellite dwarf galaxies are distributed along a coherent planar distribution known as the Vast Polar Structure (VPOS). Here we investigate whether MW satellites located on the VPOS have different physical and orbital properties from those not associated with it. Using the proper motion measurements of the MW satellites from the \textit{Gaia} mission and literature values for their observational parameters, we first discriminate between systems that may or may not be associated with the VPOS, and then compare their chemical and dynamical properties. Comparing the luminosity distributions of the on-plane and off-plane samples, we find an excess of bright satellites observed on the VPOS. Despite this luminosity gap, we do not observe a significant preference for on-plane and off-plane systems to follow different scaling relations. The on-plane systems also show a striking pattern in their radial velocities and orbital phases: co-orbiting satellites are almost all approaching their pericenters, while both counter-orbiting satellites are leaving their last pericenters. This contrasts with the more random distribution of the off-plane sample. The on-plane systems also tend to have the lowest orbital energies for a given value of angular momentum. These results are robust to the assumed MW potential, even in the case of a potential perturbed by the arrival of a massive LMC. Considering them a significant property of the VPOS, we explore several scenarios, all related to the late accretion of satellite systems, which interpret the VPOS as a young structure. We hypothesise that the VPOS formed as a result of the accretion of a group of dwarf galaxies. More accurate proper motions and dedicated studies in the context of cosmological simulations are needed to confirm this scenario.

Spectral observations of Neptune made in 2019 with the MUSE instrument at the Very Large Telescope in Chile have been analysed to determine the spatial variation of aerosol scattering properties and methane abundance in Neptune's atmosphere. The darkening of the South Polar Wave (SPW) at $\sim$ 60$^\circ$S, and dark spots such as the Voyager 2 Great Dark Spot is concluded to be due to a spectrally-dependent darkening ($\lambda < 650$nm) of particles in a deep aerosol layer at $\sim$ 5 bar and presumed to be composed of a mixture of photochemically-generated haze and H$_2$S ice. We also note a regular latitudinal variation of reflectivity at wavelengths of very low methane absorption longer than $\sim$ 650 nm, with bright zones latitudinally separated by $\sim$ 25$^\circ$. This feature, similar to the spectral characteristics of a discrete deep bright spot DBS-2019 found in our data, is found to be consistent with a brightening of the particles in the same $\sim$5-bar aerosol layer at $\lambda > 650 $ nm. We find the properties of an overlying methane/haze aerosol layer at $\sim$ 2 bar are, to first-order, invariant with latitude, while variations in the opacity of an upper tropospheric haze layer reproduce the observed reflectivity at methane-absorbing wavelengths, with higher abundances found at the equator and also in a narrow `zone' at $80^\circ$S. Finally, we find the mean abundance of methane below its condensation level to be 6--7% at the equator reducing to $\sim$3% south of $\sim$25$^\circ$S, although the absolute abundances are model dependent.

The properties of pre-eruptive structures and coronal mass ejections (CMEs) are characterized by those of their footpoints, the latter of which thus attract great interest. However, how to identify the footpoints of pre-eruptive structures and how to identify the footpoints with ground-based instruments, still remain elusive. In this work, we study an arc-shaped structure intruding in the sunspot umbra. It is located close to the (pre-)eruptive flux rope footpoint and is thus expected to help identify the footpoint. We analyse this arc-shaped structure, which we name as "sunspot scar", in a CME event on 2012 July 12 and in two CME events in observationally-inspired MHD simulations performed by OHM and MPI-AMRVAC. The sunspot scar has a more inclined magnetic field with a weaker vertical component and a stronger horizontal component relative to that in the surrounding umbra and manifests as a light bridge in the white light passband. The hot field lines anchored in the sunspot scar are spatially at the transition between the flux rope and the background coronal loops, and temporally in the process of the slipping reconnection which builds up the flux rope. The sunspot scar and its related light bridge mark the edge of the CME flux rope footpoint, and especially, the edge of the pre-eruptive flux rope footpoint in the framework of "pre-eruptive structures being flux ropes". Therefore, they provide a new perspective for the identification of pre-eruptive and CME flux rope footpoints, and also new methods for studying the properties and evolution of pre-eruptive structures and CMEs with photospheric observations only.

Low surface brightness galaxies (LSBGs) which are defined as galaxies that are fainter than the night sky, play a crucial role in understanding galaxy evolution and cosmological models. Upcoming large-scale surveys like Rubin Observatory Legacy Survey of Space and Time (LSST) and Euclid are expected to observe billions of astronomical objects. In this context, using semi-automatic methods to identify LSBGs would be a highly challenging and time-consuming process and demand automated or machine learning-based methods to overcome this challenge. We study the use of transformer models in separating LSBGs from artefacts in the data from the Dark Energy Survey (DES) data release 1. Using the transformer models, we then search for new LSBGs from the DES that the previous searches may have missed. Properties of the newly found LSBGs are investigated, along with an analysis of the properties of the total LSBG sample in DES. We identified 4,083 new LSBGs in DES, adding an additional $\sim17\% $ to the LSBGs already known in DES. This also increased the number density of LSBGs in DES to 5.5 deg$^{-2}$. We performed a clustering analysis of the LSBGs in DES using an angular two-point auto-correlation function and found that LSBGs cluster more strongly than their high surface brightness counterparts. We associated 1310 LSBGs with galaxy clusters and identified 317 among them as ultra-diffuse galaxies (UDGs). We found that these cluster LSBGs are getting bluer and larger in size towards the edge of the clusters when compared with those in the centre. Transformer models have the potential to be on par with convolutional neural networks as state-of-the-art algorithms in analysing astronomical data.

In this study, we have used the {\it Gaia} Third Data Release (Gaia DR3) to investigate an intermediate-age open cluster Collinder 74. Taking into account the stars with membership probabilities over 0.5 and inside the limiting radius of the cluster, we identified 102 most likely cluster members. The mean proper-motion components of Collinder 74 are estimated as ($\mu_{\alpha}\cos \delta, \mu_{\delta})=(0.960 \pm 0.005, -1.526 \pm 0.004$) mas yr$^{-1}$. We detected previously confirmed four blue straggler stars which show flat radial distribution. Colour excess, distance, and age of the cluster were estimated simultaneously by fitting {\sc PARSEC} isochrones to the observational data on {\it Gaia} based colour magnitude diagram. These values were derived as $E(G_{\rm BP}-G_{\rm RP})=0.425\pm 0.046$ mag, $d=2831 \pm118$ pc and $t=1800 \pm 200$ Myr, respectively. The mass function slope was estimated as $\Gamma = 1.34 \pm 0.21$ within the mass range $0.65\leq M/ M_{\odot}\leq 1.58$ which is well matched with that of Salpeter. Stellar mass distribution indicated that the massive and most likely stars are concentrated around the cluster center. The total mass of the cluster was found to be 365 $M/M_{\odot}$ for the stars with probabilities $P>0$. Galactic orbit integration shows that the Collinder 74 follows a boxy pattern outside the solar circle and is a member of the thin-disc component of the Galaxy.

Mature neutron stars are expected exhibit gravity g-modes due to stratification caused by varying composition. These modes are affected by nuclear reactions, leading to complex (damped) mode frequencies and the suppression of high order g-modes, in contrast with the common non-dissipative analysis which leads to an infinite g-mode spectrum. Focusing on the transition between the fast and slow reaction regimes, we examine the impact of nuclear reactions on the g-mode spectrum. The general framework for the analysis is presented along with sample numerical results for the BSk21 equation of state and the standard Urca reactions.

Elementary particles such as quarks and gluons are expected to be fundamental degrees of freedom at ultra high temperatures or densities, while natural phenomena in our daily lives are described in terms of hadronic degrees of freedom. Massive neutron stars and remnants of binary neutron star mergers may contain quark matter, but it is not known how the transition from hadron matter to quark matter occurs. Different transition scenarios predict different gravitational waveforms emitted from binary neutron star mergers. If the difference between the equations of state occurs at sufficiently high density, it is expected that the difference between waveforms mainly appears in the merger or the post-merger phase rather than in the inspiral phase. The typical frequency of gravitational waves after the coalescence is higher than 2 kHz, which is difficult to observe using current detectors. In this study, we performed Bayesian model selection for two representative scenarios and investigated whether observations with future detectors will allow us to identify the correct model. We assume that the relatively low density equation of state around the nuclear saturation density is completely known from accumulated observations. Under this assumption, we find that it is reasonable to expect to be able to identify the correct transition scenario with third-generation detectors or specialized detectors with high sensitivity at high frequencies designed for post-merger signal observation, e.g., NEMO.

We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.27$\pm$0.42 R$_{e}$ for the outer one. We acquired 111 radial-velocity spectra with HARPS-N and 108 radial-velocity spectra with SOPHIE. After careful periodogram analysis, we derive masses for both planets via Gaussian Process regression: 9.13$^{+0.78}_{-0.76}$ M$_{e}$ for TOI-2134b and 41.86$^{+7.69}_{-7.83}$ M$_{e}$ for TOI-2134c. We analysed the photometric and radial-velocity data first separately, then jointly. The inner planet is a mini-Neptune with density consistent with either a water-world or a rocky core planet with a low-mass H/He envelope. The outer planet has a bulk density similar to Saturn's. The outer planet is derived to have a significant eccentricity of 0.67$^{+0.05}_{-0.06}$ from a combination of photometry and RVs. We compute the irradiation of TOI-2134c as 1.45$\pm$0.10 times the bolometric flux received by Earth, positioning it for part of its orbit in the habitable sone of its system. We recommend further RV observations to fully constrain the orbit of TOI-2134c. With an expected Rossiter-McLaughlin (RM) effect amplitude of 7.2$\pm$1.3 m/s, we recommend TOI-2134c for follow-up RM analysis to study the spin-orbit architecture of the system. We calculate the Transmission Spectroscopy Metric, and both planets are suitable for bright-mode NIRCam atmospheric characterisation.

Physical constraints have been suggested to make neural network models more generalizable, act scientifically plausible, and be more data-efficient over unconstrained baselines. In this report, we present preliminary work on evaluating the effects of adding soft physical constraints to computer vision neural networks trained to estimate the conditional density of redshift on input galaxy images for the Sloan Digital Sky Survey. We introduce physically motivated soft constraint terms that are not implemented with differential or integral operators. We frame this work as a simple ablation study where the effect of including soft physical constraints is compared to an unconstrained baseline. We compare networks using standard point estimate metrics for photometric redshift estimation, as well as metrics to evaluate how faithful our conditional density estimate represents the probability over the ensemble of our test dataset. We find no evidence that the implemented soft physical constraints are more effective regularizers than augmentation.

Recent discovery of 20 TeV radiation from the Vela pulsar confirms (tentatively, at the level of crude estimates) the Aristotelian Electrodynamics picture of pulsar radiation: pulsars shine, mostly in GeV, by annihilating colliding Poynting fluxes into curvature radiation near the light cylinder. The observed GeV photons are the curvature radiation of electrons/positrons with Lorentz factors ~10^8. These "super-ultra-relativistic" electrons/positrons must also produce TeV radiation by inverse Compton if low-energy target photons are available, as they are in the Vela pulsar.

The plasma composition of the solar corona is different from that of the solar photosphere. Elements that have a low first ionisation potential (FIP) are preferentially transported to the corona and, therefore, show enhanced abundances in the corona compared to the photosphere. The level of enhancement is measured using the FIP bias parameter. In this work, we use data from the EUV Imaging Spectrometer (EIS) on Hinode to study the plasma composition in an active region following an episode of significant new flux emergence into the pre-existing magnetic environment of the active region. We use two FIP bias diagnostics: Si X 258.375 A/S X 264.233 A (temperature of approximately 1.5 MK) and Ca XIV 193.874 A/Ar XIV 194.396 A (temperature of approximately 4 MK). We observe slightly higher FIP bias values with the Ca/Ar diagnostic than Si/S in the newly emerging loops, and this pattern is much stronger in the preexisting loops (those that had been formed before the flux emergence). This result can be interpreted in the context of the ponderomotive force model, which proposes that the plasma fractionation is generally driven by Alfv\'en waves. Model simulations predict this difference between diagnostics using simple assumptions about the wave properties, particularly that the fractionation is driven by resonant/non-resonant waves in the emerging/preexisting loops. We propose that this results in the different fractionation patterns observed in these two sets of loops.

(Exo-)planets inherit their budget of chemical elements from a protoplanetary disk. The disk temperature determines the phase of each chemical species, which sets the composition of solids and gas available for planet formation. We investigate how gap structures, which are widely seen by recent disk observations, alter the thermal and chemical structure of a disk. Planet-disk interaction is a leading hypothesis of gap formation and so such changes could present a feedback that planets have on planet-forming material. Both the planet gap-opening process and the disk thermal structure are well studied individually, but how the gap-opening process affects disk thermal structure evolution remains an open question. We develop a new modelling method by iterating hydrodynamical and radiative transfer simulations to explore the gap-opening feedback on disk thermal structure. We carry out parameter studies by considering different planet locations rp and planet masses Mp. We find that for the same rp and Mp, our iteration method predicts a wider and deeper gap than the non-iteration method. We also find that the inner disk and gap temperature from the iteration method can vary strongly from the non-iteration or disk without planets, which can further influence dust-trap conditions, iceline locations, and distribution of various ices, such as H2O, CO2, and CO on large dust grains ("pebbles"). Through that, a gap-opening planet can complicate the canonical picture of the non-planet disk C/O ratio and influence the composition of the next generation of planetesimals and planets.

The current and next observation seasons will detect hundreds of gravitational waves (GWs) from compact binary systems coalescence at cosmological distances. When combined with independent electromagnetic measurements, the source redshift will be known, and we will be able to obtain precise measurements of the Hubble constant $H_0$ via the distance-redshift relation. However, most observed mergers are not expected to have electromagnetic counterparts, which prevents a direct redshift measurement. In this scenario, one of the possibilities is to use the dark sirens method that statistically marginalizes over all the potential host galaxies within the GW location volume to provide a probabilistic redshift to the source. Here we presented $H_{0}$ measurements using two new dark sirens compared to previous analyses using DECam data, GW190924$\_$021846 and GW200202$\_$154313. The photometric redshifts of the possible host galaxies of these two events are acquired from the DECam Local Volume Exploration Survey (DELVE) carried out on the Blanco telescope at Cerro Tololo in Chile. The combination of the $H_0$ posterior from GW190924$\_$021846 and GW200202$\_$154313 together with the bright siren GW170817 leads to $H_{0} = 68.84^{+15.51}_{-7.74}\, \rm{km/s/Mpc}$. Including these two dark sirens improves the 68% confidence interval (CI) by 7% over GW170817 alone. This demonstrates that the inclusion of well-localized dark sirens in such analysis improves the precision with which cosmological measurements can be made. Using a sample containing 10 well-localized dark sirens observed during the third LIGO/Virgo observation run, we determine a measurement of $H_{0} = 76.00^{+17.64}_{-13.45}\, \rm{km /s/Mpc}$.

Dark matter and cosmic inflation represent two of the major puzzles in cosmology. They are typically addressed by introducing separate fields with independent dynamics. On the other hand, extra dimensions might play an important role for observable physics. We introduce a five-dimensional model called dark inflaxion that includes an axion-like particle, whose Kaluza-Klein modes can describe inflation and dark matter. We show that a simple yet natural choice of the mass scale of the effective four-dimensional fields of our model can accommodate simultaneously the observable values of inflationary parameters and dark-matter abundance. It will be interesting to explore the consequences and predictions of more generic scenarios within the scope of our model, which include the possibility of multifield inflation and dark matter.

We investigate a generic source of stochastic gravitational wave background (SGWB) due to the parametric resonance of oscillating scalar fields in the early Universe. By systematically analyzing simple benchmark models using lattice simulation and considering a wide range of parameter space, we demonstrate that such a scenario can lead to detectable signals in GW detectors over a broad frequency range and potentially address the recent finding by NANOGrav. Furthermore, these models are found to naturally yield ultra-light dark matter candidates or dark radiation detectable by CMB observatories.

The immiscibility of hydrogen-helium mixture under the temperature and pressure conditions of planetary interiors is crucial for understanding the structures of gas giant planets (e.g., Jupiter and Saturn). While the experimental probe at such extreme conditions is challenging, theoretical simulation is heavily relied in an effort to unravel the mixing behavior of hydrogen and helium. Here we develop a method via a machine learning accelerated molecular dynamics simulation to quantify the physical separation of hydrogen and helium under the conditions of planetary interiors. The immiscibility line achieved with the developed method yields substantially higher demixing temperatures at pressure above 1.5 Mbar than earlier theoretical data, but matches better to the experimental estimate. Our results revise the structures of Jupiter and Saturn where H-He demixing takes place in a large fraction of the interior radii, i.e., 27.5% in Jupiter and 48.3% in Saturn. This direct evidence of an H-He immiscible layer supports the formation of helium rain and explains the helium reduction in atmosphere of Jupiter and Saturn.

In a recent Comment on the paper "Dark matter as a Weyl geometric effect", by Burikham et al., Phys. Rev. D 107, 064008 (2023), posted on arxiv. org as eprint arXiv:2306.11926, it was claimed that the exact solution found in the above mentioned paper by Burikham et al. "is wrong". In this Reply to the Comment we present, in a clear and comprehensive way, a step by step derivation of the exact solution of the vacuum static spherically symmetric field equations of the Weyl geometric gravity theory, and we show that, contrary to the claims in arXiv:2306.11926, the obtained solution is correct, and it satisfies all the equations of motion of the basic theory. Hence, it can be considered as a viable alternative model for the explanation of the behavior of the galactic rotation curves, without invoking the presence of dark matter.

We introduce a new output amplifier for fully-depleted thick p-channel CCDs based on double-gate MOSFETs. The charge amplifier is an n-type MOSFET specifically designed and operated to couple the fully-depleted CCD with high charge-transfer efficiency. The junction coupling between the CCD and MOSFET channels has enabled high sensitivity, demonstrating sub-electron readout noise in one pixel charge measurement. We have also demonstrated the non-destructive readout capability of the device. Achieving single-electron and single-photon per pixel counting in the entire CCD pixel array has been made possible through the averaging of a small number of samples. We have demonstrated fully-depleted CCD readout with better performance than the floating diffusion and floating gate amplifiers available today, in both single and multisampling regimes, boasting at least six times the speed of floating gate amplifiers.

We study phenomenology of a light scalar dark matter (DM). In the model, there are an inert doublet scalar and a singlet Dirac fermion $\psi$, both charged under a global $Z_2$ symmetry. The mass of the lightest inert scalar $H$ can be lighter than 10 GeV by imposing appropriate relations between three scalar quartic couplings. The lightest $Z_2$ odd particle is stable and DM. In this paper, focusing on the parameter space where $H$ is lighter than $\psi$ and is DM, we discuss DM physics related to relic density, direct detection, indirect detection, collider searches and other cosmological observations. We clarify differences from the case where $\psi$ is instead DM, which has been focused on in the previous works.