Papers for Friday, Jun 16 2023

The accurate modelling of the Point Spread Function (PSF) is of paramount importance in astronomical observations, as it allows for the correction of distortions and blurring caused by the telescope and atmosphere. PSF modelling is crucial for accurately measuring celestial objects' properties. The last decades brought us a steady increase in the power and complexity of astronomical telescopes and instruments. Upcoming galaxy surveys like Euclid and LSST will observe an unprecedented amount and quality of data. Modelling the PSF for these new facilities and surveys requires novel modelling techniques that can cope with the ever-tightening error requirements. The purpose of this review is three-fold. First, we introduce the optical background required for a more physically-motivated PSF modelling and propose an observational model that can be reused for future developments. Second, we provide an overview of the different physical contributors of the PSF, including the optic- and detector-level contributors and the atmosphere. We expect that the overview will help better understand the modelled effects. Third, we discuss the different methods for PSF modelling from the parametric and non-parametric families for ground- and space-based telescopes, with their advantages and limitations. Validation methods for PSF models are then addressed, with several metrics related to weak lensing studies discussed in detail. Finally, we explore current challenges and future directions in PSF modelling for astronomical telescopes.

Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present the first 3D simulations of massive star convection extending from the star's center to near its surface, with realistic stellar luminosities. Using these simulations, we make the first prediction of photometric variability due to convectively-driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, so the red noise must be generated by an alternative process.

The next Milky Way supernova will be an epochal event in multi-messenger astronomy, critical to tests of supernovae, neutrinos, and new physics. Realizing this potential depends on having realistic simulations of core collapse. We investigate the neutrino predictions of nearly all modern models (1-, 2-, and 3-d) over the first $\simeq$1 s, making the first detailed comparisons of these models to each other and to the SN 1987A neutrino data. Even with different methods and inputs, the models generally agree with each other. However, even considering the low neutrino counts, the models generally disagree with data. What can cause this? We show that neither neutrino oscillations nor different progenitor masses appear to be a sufficient solution. We outline urgently needed work.

In this paper we explore how different regularization prescriptions affect the counterterms in the renormalization of the galaxy bias expansion. We work in the context of primordial local non-Gaussianity including non-linear gravitational evolution. We carry out the one-loop renormalization of the field $\delta_\rho^2$ (i.e. the square of the matter overdensity field) up to third order in gravitational evolution. Three regularization schemes are considered and their impact on the values of the counterterms is studied. We explicitly verify that the coefficients of the non-boost invariant operators are regularization scheme independent.

The outer envelopes of massive ($M\gtrsim10\,M_{\odot}$) stars exhibit large increases in opacities from forests of lines and ionization transitions (particularly from iron and helium) that trigger near-surface convection zones. One-dimensional models predict density inversions and supersonic motions that must be resolved with computationally intensive 3D radiation hydrodynamic (RHD) modeling. Only in the last decade have computational tools advanced to the point where ab initio 3D models of these turbulent envelopes can be calculated, enabling us to present five 3D RHD Athena++ models (four previously published and one new 13$M_{\odot}$ model). When convective motions are sub-sonic, we find excellent agreement between 3D and 1D velocity magnitudes, stellar structure, and photospheric quantities. However when convective velocities approach the sound speed, hydrostatic balance fails as the turbulent pressure can account for 80% of the force balance. As predicted by Henyey, we show that this additional pressure support leads to a modified temperature gradient which reduces the superadiabaticity where convection is occurring. In addition, all five models display significant overshooting from the convection in the Fe convection zone. As a result, the turbulent velocities at the surface are indicative of those in the Fe zone. There are no confined convection zones as seen in 1D models. In particular, helium convection zones seen in 1D models are significantly modified. Stochastic low frequency brightness variability is also present in the 13$M_{\odot}$ model with comparable amplitude and characteristic frequency to observed stars.

We have measured, at various wavelengths, the spiral arm pitch angles of a sample of distant spiral galaxies from the Hubble Space Telescope eXtreme Deep Field (XDF). According to density wave theory, we should detect colour jumps from red-to-blue across the spiral arms. Colour jumps are a consequence of large-scale shocks, which also generate the classic blue-to-red age/colour gradients, and have only been detected until now in nearby spiral galaxies. Our results indicate that colour jumps and gradients have been occurring in distant galaxies for at least the last 8 Gyr, in agreement with density wave theory.

We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Galaxy survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity ($V_{\rm{rot}}$) and velocity dispersion within 1$R_{\text{e}}$ ($\sigma_{\text{e}}$) for a sub-set of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI Galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured $V/\sigma$ by nearly a factor two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the LEGA-C survey and, looking forward, E-ELT HARMONI. Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift.

This paper introduces the class of omnipotent dark energy (DE) models characterized by non-monotonic energy densities that are capable of attaining negative values with corresponding equation of state parameters featuring phantom divide line (PDL) crossings and singularities. These non-trivial features are phenomenologically motivated by findings of previous studies that reconstruct cosmological functions from observations, and the success of extensions of $\Lambda$CDM, whose actual or effective DE density is omnipotent, in alleviating the observational discordance within $\Lambda$CDM. As an example, we focus on one embodiment of omnipotent DE, viz., the DE parametrization introduced in arXiv:2005.12587 (DMS20). By updating and extending the data sets used in the original paper where it was introduced, we confirm the effectiveness of DMS20 in alleviating the observational discrepancies. Additionally, we uncover that its negative DE density feature, importance of which was not previously investigated, plays a crucial role in alleviating the tensions, along with the PDL crossing feature that the parametrization presupposes. In particular, we find that there is a positive correlation between the $H_0$ parameter and the scale~($a_{\rm p}$) at which DE density transitions from negative to positive, in agreement with previous studies that incorporate this transition feature. For our full data set, the model yields $H_0=70.05\pm0.64$ (68\% CL) relaxing the $H_0$ tension with a preference of crossing to negative DE densities ($a_{\rm p}>0$ at 99\% CL), along with the constraint $a_m=0.922^{+0.041}_{-0.035}$ on the scale of the presupposed PDL crossing.

We present preliminary results from our spectroscopic survey of low-luminosity early-type galaxies in the Coma cluster conducted with the Binospec spectrograph at the 6.5~m MMT. From spatially-resolved profiles of internal kinematics and stellar population properties complemented with high-resolution images, we placed several low-luminosity dEs on the fundamental plane in the low-luminosity extension of the available literature data. We also discovered unusual kpc-sized kinematically-decoupled cores in several dwarf galaxies, which had been probably formed before these galaxies entered the cluster.

Ultra-diffuse galaxies (UDGs) are spatially extended, low surface brightness stellar systems with regular elliptical-like morphology found in large numbers in galaxy clusters and groups. Studies of the internal dynamics and dark matter content of UDGs have been hampered by their low surface brightnesses. We identified a sample of low-mass early-type post-starburst galaxies, `future UDGs' in the Coma cluster still populated with young stars, which will passively evolve into UDGs in the next 5$-$10 Gyr. We collected deep observations for a sample of low-mass early-type galaxies in the Coma cluster using MMT Binospec, which includes present-day and future UDGs. We derived their dark matter content within a half-light radius (70$-$95%) and total dynamical masses ($M_{200}=5.5\cdot10^9-1.4\cdot10^{11} M_{\odot}$) assuming the Burkert density profile and assess how different proposed evolutionary channels affect dark and visible matter in UDGs. We also discuss observational methodology of present and future UDG studies.

We have studied the pulse-to-pulse variability in PSR J1820--0427 and its frequency dependence using high-quality, wide-band observations made from the upgraded Giant Metrewave Radio Telescope (uGMRT; 300-750 MHz) and the Murchison Widefield Array ($\sim$170-200 MHz). The low-frequency data reveal a previously unreported feature in the average profile (at 185 MHz) after accounting for the effects of temporal broadening arising from multi-path scattering due to the Interstellar Medium (ISM). We advance a new method for flux density calibration of beamformed data from the uGMRT and use it to measure the single pulse flux densities across the uGMRT band. Combined with previously published measurements, these flux densities are best fit with a power-law spectrum with a low-frequency turnover. We also use calibrated flux densities to explore the relationship between pulse-to-pulse variability and the spectral index of individual pulses. Our analysis reveals a large scatter in the single-pulse spectral indices and a general tendency for brighter pulses to show a steepening of the spectral index. We also examine the frequency-dependence of the pulse-fluence distribution and its relation to the Stochastic Growth Theory.

The Kerr metric that describes the spacetime of a spinning supermassive black hole (SMBH) is axisymmetric, implying that the nearly parabolic geodesics on which stars approach the SMBH depend on the inclination angle $\iota$ of the orbital angular momentum with respect to the SMBH spin. This inclination affects both the geodesic deviation that determines whether a star is tidally disrupted and whether the tidal debris survives direct capture by the event horizon to produce an observable tidal disruption event (TDE). The steady-state TDE rate is the rate at which stars are scattered into the loss cone determined by these spin- and inclination-dependent effects. As the anisotropy of this loss-cone refilling is highly uncertain, we consider the two extreme limits in which stellar inclination is preserved (IP) or isotropized (ISO). We calculate the inclination distribution in these two limits and find a prograde bias in the IP limit because of the strong retrograde bias for direct capture. However, we find a retrograde bias in the ISO limit for intermediate SMBH masses when the empty loss cone suppresses capture and allows the weaker retrograde bias of geodesic deviation to dominate. We also calculate the total TDE rates and maximum SMBH mass $M_{\rm \bullet, max}$ for tidal disruption in these two limits. In the IP limit, we find a highly spin-dependent capture cutoff in the TDE rate and $M_{\rm \bullet, max} \approx 10^{8.45} M_\odot$ for maximal SMBH spin. In the ISO limit, we find a strong spin-dependent enhancement in the TDE rate at intermediate SMBH masses, a weakly spin-dependent capture cutoff above $M_\bullet \approx 10^{7.5} M_\odot$, and $M_{\rm \bullet, max} \approx 10^{7.95} M_\odot$ for maximal SMBH spin.

We present a method for identifying radio stellar sources using their proper-motion. We demonstrate this method using the FIRST, VLASS, RACS-low and RACS-mid radio surveys, and astrometric information from Gaia Data Release 3. We find eight stellar radio sources using this method, two of which have not previously been identified in the literature as radio stars. We determine that this method probes distances of ~90pc when we use FIRST and RACS-mid, and ~250pc when we use FIRST and VLASS. We investigate the time baselines required by current and future radio sky surveys to detect the eight sources we found, with the SKA (6.7 GHz) requiring <3 years between observations to find all eight sources. We also identify nine previously known and 43 candidate variable radio stellar sources that are detected in FIRST (1.4 GHz) but are not detected in RACS-mid (1.37 GHz). This shows that many stellar radio sources are variable, and that surveys with multiple epochs can detect a more complete sample of stellar radio sources.

The X-ray pulsar RX J0440.9+4431 went through a giant outburst in 2022 and reached a record-high flux of nearly 2.3 Crab, as observed by Swift/BAT. We study the evolution of different spectral and timing properties of the source using NICER observations during the outburst. The pulse period is found to decrease from 208 s to 205 s and the pulse profile evolved significantly during the outburst with energy and luminosity. The emission mechanism and beaming patterns may change significantly, which is related to the state transition of the source. The hardness ratio shows significant evolution during the outburst, and the hardness intensity diagram also shows two different branches. The HID turns towards the diagonal branch from the horizontal branch above the critical luminosity. The observed photon index shows a negative correlation with X-ray flux below the critical luminosity, which turns into a positive correlation above the critical luminosity. This indicates a spectral transition from the sub-critical to the super-critical regime. The magnetic field is estimated to be nearly 1.8 $\times$ 10$^{12}$ G using critical luminosity of $\sim$2.7 $\times$ 10$^{37}$ erg s$^{-1}$. NICER spectra can be described using a cutoff power-law model with a blackbody component and an additional 6.4 keV iron fluorescence line. The iron emission line evolves from a narrow to a broad feature with the increase in luminosity. The iron line flux is strongly correlated with the X-ray flux. Based on luminosity, the Fe band featured two emission lines at 6.4 and 6.67 keV originating from neutral and highly ionized Fe-atoms.

The Cosmological Recombination Radiation (CRR) is one of the guaranteed $\Lambda$CDM Spectral Distortion (SD) signals. Even if very small in amplitude, it provides a direct probe of the three recombination eras, opening the path for testing one of the key pillars in our cosmological interpretation of the measured CMB anisotropies. Here we develop a new emulator, CRRfast, to quickly and accurately represent the CRR for a wide range of cosmologies, using the state-of-the-art CosmoSpec code as a reference. CRRfast has been made publicly available both as stand-alone code and as part of CLASS, thereby completing the set of $\Lambda$CDM sources of SDs that can be modeled with CLASS. With this newly-developed pipeline we investigate the full constraining power of SDs within $\Lambda$CDM and highlight possible future applications to experimental design optimization. Furthermore, we show that the inhomogeneous evolution of the recombination process imprints second-order contributions to the CRR spectrum, leading to a broadening and shifting of the CRR features. These second-order terms are naturally captured by the emulator and allow us to evaluate the $\Lambda$CDM contributions to the average CRR as well as to illustrate the effect of perturbed recombination due to Primordial Magnetic Fields (PMFs). As it turns out, while the $\Lambda$CDM variance effects can be neglected, they could be significantly enhanced in the beyond-$\Lambda$CDM models. In particular in the case of PMFs we demonstrate that through these non-linear terms the parameter space relevant to the Hubble tension could be tested with future CMB spectrometers.

Post-starburst galaxies are sources that had the last major episode of star formation about 1 Gyr before the epoch of the observations and are on their way to quiescence. It is important to study such galaxies at redshift z > 1, during their main quenching phase, and estimate their molecular gas content to constrain the processes responsible for the cessation of star formation. We present CO(3-2) ALMA observations of two massive (Mstar ~ 5 x 10^10 Msun) post-starburst galaxies at z > 1. We measure their molecular gas fraction to be f_H2 = M_H2/Mstar ~ 8% - 16%, consistent with z < 1 post-starburst galaxies from the literature. The star formation efficiency of our targets is ~ 10x lower than that of star-forming galaxies at similar redshift, and they are outliers of the f_H2 - specific star formation rate (sSFR) relation of star-forming galaxies, as they have larger f_H2 than expected given their sSFR. The gas fraction of post-starbursts from our sample and the literature correlates with the Dn4000 spectral index, a proxy of the stellar population age. This suggests that their gas content decreases after the last major burst of star formation. Finally, one of our targets is undergoing a major merger phase with two highly star-forming companions. This hints at a picture where a perturber event (e.g., major merger) quenches star formation without completely removing the molecular gas.

Magnetars are young neutron stars powered by the strongest magnetic fields in the Universe (10$^{13-15}$ G). Their transient X-ray emission usually manifests as short (a few hundred milliseconds), bright, energetic ($\sim$ 10$^{40-41}$ erg) X-ray bursts. Since its discovery in 2014, magnetar J1935+2154 has become one of the most prolific magnetars, exhibiting very active bursting episodes, and other fascinating events such as pulse timing anti-glitches and Fast Radio Bursts. Here, we present evidence for possible 42 Hz (24 ms) quasi-periodic oscillations in the $\nu F_{\nu}$ spectrum peak energy (Ep) identified in a unique burst detected with the Fermi Gamma-ray Burst Monitor in January 2022. While quasi-periodic oscillations have been previously reported in the intensity of magnetar burst lightcurves, quasi-periodic oscillations in the Ep have not. We also find an additional event from the same outburst that appears to exhibit similar character in Ep, albeit of lower statistical quality. For these two exceptional transients, such Ep oscillations can be explained by magnetospheric density and pressure perturbations. For burst-emitting plasma consisting purely of $e^+e^-$ pairs, these acoustic modes propagate along a highly magnetized flux tube of length up to around $L\sim 130$ neutron star radii, with $L$ being lower if ions are present in the emission zone. Detailed time-resolved analyses of other magnetar bursts are encouraged to evaluate the rarity of these events and their underlying mechanisms.

Type Ia supernovae remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of sub-types of Type Ia supernovae have been identified in recent decades. Due to the massive computational load required, inference of the internal structure of Type Ia supernovae ejecta directly from observations using simulations has previously been computationally intractable. However, deep-learning emulators for radiation transport simulations have alleviated such barriers. We perform abundance tomography on 40 Type Ia supernovae from optical spectra using the radiative transfer code TARDIS accelerated by the probabilistic DALEK deep-learning emulator. We apply a parametric model of potential ejecta structures to comparatively investigate abundance distributions and internal ionization fractions of intermediate-mass elements between normal and 1991T-like Type Ia supernovae. Our inference shows that 1991T-like Type Ia supernovae are under-abundant in the typical intermediate mass elements that heavily contribute to the spectral line formation seen in normal Type Ia supernovae at early times. Additionally, we find that the intermediate-mass elements present in 1991T-like Type Ia supernovae are highly ionized compared to those in the normal Type Ia population. Finally, we conclude that the transition between normal and 1991T-like Type Ia supernovae appears to be continuous observationally and that the observed differences come out of a combination of both abundance and ionization fractions in these supernovae populations.

Recent Gaia observations suggest that some hypervelocity stars (HVSs) might originate from outside the Galaxy. We ask if these HVSs could come from as far as Andromeda. Therefore, we simulate HVSs originating in Andromeda with initial conditions based on attributes of high-velocity stars measured in the Milky Way and a simple model for the gravitational potential of Andromeda and the Milky Way. We evaluate the validity of this scenario based on the simulation results. While we expect the vast majority of HVSs in our Galaxy will originate here, we expect the number of stars present from Andromeda at any one time to be between twelve and 3910, depending upon model assumptions. Further, we analyse the properties of HVSs that are able to reach the Milky Way and discuss whether they could be detected experimentally based on recent constraints set on the ejection rate of HVSs from the Milky Way centre.

Young stellar populations provide a powerful record that traces millions of years of star formation history in the solar neighborhood. Using a revised form of the SPYGLASS young star identification methodology, we produce an expanded census of nearby young stars (Age $<50$ Myr). We then use the HDBSCAN clustering algorithm to produce a new SPYGLASS Catalog of Young Associations (SCYA), which reveals 116 young associations within 1 kpc. More than 25\% of these groups are largely new discoveries, as 20 are substantively different from any previous definition, and 10 have no equivalent in the literature. The new associations reveal a yet undiscovered demographic of small associations with little connection to larger structures. Some of the groups we identify are especially unique for their high transverse velocities, which can differ from the solar velocity by 30-50 km s$^{-1}$, and for their positions, which can reach up to 300 pc above the galactic plane. These features may suggest a unique origin, matching existing evidence of infalling gas parcels interacting with the disk ISM. Our clustering also suggests links between often-separated populations, hinting to direct structural connections between Orion Complex and Perseus OB2, and between the subregions of Vela. The $\sim$30 Myr old Cepheus-Hercules association is another emerging large-scale structure, with a size and population comparable to Sco-Cen. Cep-Her and other similarly-aged structures are also found clustered along extended structures perpendicular to known spiral arm structure, suggesting that arm-aligned star formation patterns have only recently become dominant in the solar neighborhood.

We present results from a sample of 106 high-luminosity IRAS sources observed with the Very Large Array in the B and C configurations. 96 sources were observed in the X-band and 52 in the K-band, with 42 of them observed at both wavelengths. We also used previously published observations in the C-band for 14 of them. The detection rate of sources with 3.6~cm continuum emission was $\sim25\%$, while only 10\% have emission at 1.3~cm. In order to investigate the nature of these sources, their physical parameters were calculated mainly using the 3.6~cm continuum emission, and for sources detected at two wavelengths, we used the best fit of three HII region models with different geometries. As a final result, we present a catalog of the detected sources, which includes their basic physical parameters for further analysis. The catalog contains 17 ultracompact HII regions and 3 compact HII regions.

We analyze the structure and dynamics of the plasma atmospheres and Coulomb-liquid oceans on neutron stars. Salient dynamical parameters are identified and their values estimated for the governing set of magnetohydrodynamics equations. Neutron star atmospheres and oceans are strongly stratified and, depending on the rotation period, contain a multitude of long-lived vortices (spots) and/or narrow zonal jets (free-shear zones) in the large plasma-beta regime-i.e., $\beta_p \gg 1$ (hydrodynamic regime). In contrast, when $\beta_p \lesssim 1$ (magnetohydrodynamic regime), the flow is dominated by a global lattice of effectively fixed magnetic islands (plasmoids)-without any jets. Understanding the spatio-temporal variability of dynamic atmospheres and oceans on neutron stars is crucial for interpreting observations of their X-ray emissions.

It was recently suggested that "cosmologically coupled" black holes with masses that increase in proportion to the volume of the Universe might constitute the physical basis of dark energy. We take this claim at face value and discuss its potential astrophysical implications. We show that the gravitational wave emission in binary systems would be significantly enhanced so that the number of black hole mergers would exceed the observed rate by orders of magnitude, with typical masses much larger than those seen by the LIGO-Virgo-KAGRA network. Separately, if the mass growth happens at fixed angular momentum, the supermassive black holes in matter-deficient elliptical galaxies should be slowly rotating. Finally, cosmological coupling would stabilize small black holes against Hawking radiation-induced evaporation.

We study the radial and non-radial oscillations of Cross stars (CrSs), i.e., stars with a quark matter crust and a hadronic matter core in an inverted order compared to conventional hybrid stars. We draw comparisons of their oscillation modes with those of neutron stars, quark stars, and conventional hybrid stars. We find that the stellar stability analysis from the fundamental mode of radial oscillations, and the $g$, $f$ modes of non-radial oscillations are quite similar to those of conventional hybrid stars. However, due to the inverted stellar structure, the first non-radial $p$ mode of CrSs behaves in an inverted way and sits in a higher frequency domain compared to that of conventional hybrid stars. These results provide a direct way to discriminate CrSs from other types of compact stars via gravitational-wave probes. Specifically, compact stars emitting $g$-mode gravitational waves within the $0.5$--$1$ kHz range should be CrSs or conventional hybrid stars rather than neutron stars or pure quark stars, and a further GW detection of the first $p$ mode above 8 kHz or an identification of a decreasing trend of frequencies versus star masses associated with it will help identify the compact object to be a CrS rather than a conventional hybrid star.

We numerically studied magnetic reconnection in a high $\beta$ hydrogen-helium plasma at different altitudes from the photosphere to the upper chromosphere. The time dependent ionization degrees were included to get more realistic diffusivities and viscosity, and appropriate radiative cooling models were applied. Our numerical results indicate that the plasmoid instability always plays a vital role in speeding up magnetic reconnection at different atmospheric layers. In addition, both the strong radiative cooling and the magnetic diffusion caused by the electron-neutral collision ($\eta_{en}$) can significantly accelerate magnetic reconnection below the middle chromosphere. On the other hand, both the ambipolar diffusion and the viscosity result in higher temperature and plasma pressure in the reconnection region in the upper chromosphere, which then hinder the fast reconnection process from developing. The local compression heating triggered by turbulent reconnection mediated with plasmoids is the dominant heating mechanism in the unstable reconnection stage at different atmospheric layers, but the viscous heating and the ambipolar diffusion heating are equally important in the upper chromosphere. The Joule heating contributed by $\eta_{en}$ dominates during the early quasi-steady reconnection stage below the middle chromosphere, the strong radiative cooling also leads to much stronger compression heating and more generation of thermal energy in this region. Though the plasma $\beta$ is the same in all the simulation cases at different altitudes, the temperature increase is more significant in the upper chromosphere with much lower density and weaker radiative cooling.

The sun is highly complex in nature and its observatory imagery features is one of the most important sources of information about the sun activity, space and Earth's weather conditions. The NASA, solar Dynamics Observatory captures approximately 70,000 images of the sun activity in a day and the continuous visual inspection of this solar observatory images is challenging. In this study, we developed a technique of tracking the sun's activity using 2D circular kernel time series transformation, statistical and entropy measures, with machine learning approaches. The technique involves transforming the solar observatory image section into 1-Dimensional time series (1-DTS) while the statistical and entropy measures (Approach 1) and direct classification (Approach 2) is used to capture the extraction features from the 1-DTS for machine learning classification into 'solar storm' and 'no storm'. We found that the potential accuracy of the model in tracking the activity of the sun is approximately 0.981 for Approach 1 and 0.999 for Approach 2. The stability of the developed approach to rotational transformation of the solar observatory image is evident. When training on the original dataset for Approach 1, the match index (T90) of the distribution of solar storm areas reaches T90 ~ 0.993, and T90 ~ 0.951 for Approach 2. In addition, when using the extended training base, the match indices increased to T90 ~ 0.994 and T90 ~ 1, respectively. This model consistently classifies areas with swirling magnetic lines associated with solar storms and is robust to image rotation, glare, and optical artifacts.

Context. The Si IV lines at 1394 \r{A} and 1403 \r{A} form in the solar atmosphere at a temperature of $\sim10^{4.8}$ K. They are usually considered optically thin, but their opacity can be enhanced during solar flares. Traditionally, the intensity ratio of these lines are used as an indicator of the optical thickness. However, observations have shown a wavelength-dependent intensity ratio profile $r(\Delta\lambda)$ of the the 1394 \r{A} to 1403 \r{A} lines. Aims. We aim to study the variation of the intensity ratio profile in solar flares and the physical reasons behind it. Method. The Si IV lines and their intensity ratio profiles are calculated from the one-dimensional radiative hydrodynamics flare model with non-thermal electron heating. Result. During flares, $r(\Delta\lambda)$ is smaller than 2 at the line core but larger than 2 at the line wings. We attribute the deviation of the ratio from 2 to two effects: the resonance scattering effect and the opacity effect. Resonance scattering increases the population ratio of the upper levels of the two lines, and as a result, increases $r(\Delta\lambda)$, in all wavelengths. The opacity effect decreases $r(\Delta\lambda)$, especially at the line core where the opacity is larger. These two effects compete with each other and cause the U-shape of $r(\Delta\lambda)$.

Ice is ubiquitous in the interstellar medium. We model the formation of the main constituents of interstellar ices, including H2O, CO2 , CO, and CH3 OH. We strive to understand what physical or chemical parameters influence the final composition of the ice and how they benchmark to what has already been observed, with the aim of applying these models to the preparation and analysis of JWST observations. We used the Nautilus gas-grain model, which computes the gas and ice composition as a function of time for a set of physical conditions, starting from an initial gas phase composition. All important processes (gas-phase reactions, gas-grain interactions, and grain surface processes) are included and solved with the rate equation approximation. We first ran an astrochemical code for fixed conditions of temperature and density mapped in the cold core L429-C to benchmark the chemistry. One key parameter was revealed to be the dust temperature. When the dust temperature is higher than 12 K, CO2 will form efficiently at the expense of H2O, while at temperatures below 12 K, it will not form. Whatever hypothesis we assumed for the chemistry (within realistic conditions), the static simulations failed to reproduce the observed trends of interstellar ices in our target core. In a second step, we simulated the chemical evolution of parcels of gas undergoing different physical and chemical situations throughout the molecular cloud evolution and starting a few 1e7 yr prior to the core formation (dynamical simulations). Our dynamical simulations satisfactorily reproduce the main trends already observed for interstellar ices. Moreover, we predict that the apparent constant ratio of CO2/H2O observed to date is probably not true for regions of low AV , and that the history of the evolution of clouds plays an essential role, even prior to their formation.

There are some strange quasars with multiple Gaia detections or observed with abnormal astrometric characteristics, such as with large proper motions or significant astrometric noises. Those strange quasars could be potential candidates of quasar-star pairs, dual quasars (DQs), or lensed quasars (LQs). Searching for both DQs and LQs is of great importance in many fields of astrophysics. Here in this work, we select 143 SDSS spectroscopically confirmed quasars that have multiple Gaia EDR3 detections within 1 arcsec of the SDSS quasar' position. We apply several optical identification methods to classify this sample. We firstly exclude 65 quasar-star pairs via their stellar features including their parallaxes and proper motions, stellar features in the SDSS spectra, or via the colour-colour diagram. Based on the spectral-fitting results, we find 2 DQ candidates, one of which presents a double-peaked [O III] emission line feature and the other shows a broad $H_{\beta}$ velocity offset ($\sim$ 870 $ km s^{-1} $) relative to the [O III] $\lambda$5007 line. Via the colour difference method, we further find 56 LQ candidates with similar colours in their multiple images. We also cross-match 143 objects with the HST archive and find 19 targets with archival HST images. Our classification results of those 19 targets are mainly consistent with previous works.

Outflows and winds launched from young stars play a crucial role in the evolution of protostars and the early stages of planet formation. However, the specific details of the mechanism behind these phenomena, including how they affect the protoplanetary disk structure, are still debated. We present {\it JWST} NIRSpec Integral Field Unit (IFU) observations of atomic and H$_2$ lines from 1 -- 5.1 $\mu$m toward the low-mass protostar TMC1A. For the first time, a collimated atomic jet is detected from TMC1A in the [Fe II] line at 1.644 $\mu$m along with corresponding extended H$_2$ 2.12 $\mu$m emission. Towards the protostar, we detected spectrally broad H I and He I emissions with velocities up to 300 km/s that can be explained by a combination of protostellar accretion and a wide-angle wind. The 2$\mu$m continuum dust emission, H I, He I, and O I all show emission from the illuminated outflow cavity wall and scattered line emission. These observations demonstrate the potential of {\it JWST} to characterize and reveal new information about the hot inner regions of nearby protostars. In this case, a previously undetected atomic wind and ionized jet in a well-known outflow.

We consider a model-independent approach to constrain the equivalence redshift, $z_{eq}$, at which dark energy, baryons and cold dark matter equate their magnitudes. To this aim, in the context of a homogeneous and isotropic universe, we first consider a generic model where the dark energy contribution is provided by an unknown function of barotropic fluids. Afterwards, we compute the deceleration and jerk parameters, evaluating at our epoch, namely $z=0$, and at $z=z_{eq}$. Thus, by Taylor expanding around current time the Hubble, luminosity and angular distances, we substitute the theoretical expressions obtained from the aforementioned generic dark energy model, defining a correspondence between quantities evaluated at $z=0$ and $z=z_{eq}$. In so doing, we directly fit these quantities by means of current data sets, involving the most recent Pantheon type IA supernovae, baryonic acoustic oscillation and differential Hubble points. We consider two hierarchies in our fitting procedures and compare our findings in the spatially-flat universe first and including spatial curvature, later. We assess constraints on the overall equation of state of the universe and its first derivative. We compare our results with those predicted by the standard $\Lambda$CDM paradigm. Specifically, our findings are in agreement at the 2$\sigma$ confidence level, assuming a constant dark energy term. However, our analysis does not rule out the possibility of a slight evolution of dark energy, indicating a small deviation from the scenario of a pure cosmological constant. In particular, the possible departures appear consistent with a phenomenological $\omega$CDM model, rather than more complicated dark energy parameterizations.

With the greater power to infer the state of stellar interiors provided by asteroseismology, it has become possible to study the survival of initially convective cores within stars during their main-sequence evolution. Standard theories of stellar evolution predict that convective cores in sub-solar mass stars have lifetimes below 1 Gyr. However, a recent asteroseismic study of the star Kepler-444 concluded that the initial convective core had survived for nearly 8 Gyr. The goal of this paper is to study the convective-core evolution of Kepler-444 and to investigate its proposed longevity. We modify the input physics of stellar models to induce longer convective-core lifetimes and vary the associated parameter across a dense grid of evolutionary tracks. The observations of metallicity, effective temperature, mean density, and asteroseismic frequency ratios are fitted to the models using the BASTA pipeline. We explore different choices of constraints, from which a long convective-core lifetime is only recovered for a few specific combinations: mainly from the inclusion of potentially unreliable frequencies and/or excluding the covariances between the frequency ratios; while for the classical parameters, the derived luminosity has the largest influence. For all choices of observables, our analysis reliably constrains the convective-core lifetime of Kepler-444 to be short, with a median around 0.6 Gyr and a $1\sigma$ upper bound around 3.5 Gyr.

GW190521 is the most massive merging binary black hole (BBH) system detected so far. At least one of the component BHs was measured to lie within the pair-instability supernova (PISN) mass gap ($\sim 50-135\;{\rm M}_{\odot}$), making its formation a mystery. However, the transient observed signal allows alternative posterior distributions. There was suggestion that GW190521 could be an intermediate-mass ratio inspiral (IMRI), with the component masses $m_1\sim 170\;{\rm M}_{\odot}$ and $m_2\sim 16 \;{\rm M}_{\odot}$, happening to straddle the PISN mass gap. Under this framework, we perform binary population synthesis to explore the formation of GW190521-like systems via isolated binary evolution. We numerically calculate the binding energy parameter for massive stars at different metallicities, and employ them in our calculation for common envelope evolution. Our results prefer that the progenitor binaries formed in metal-poor environment with $\rm Z\leq0.0016$. The predicted merger rate density within redshift $z=1.1$ is $\sim 4\times 10^{-5}-5\times 10^{-2} \,\rm Gpc^{-3}yr^{-1}$. We expect that such events are potentially observable by upcoming both space and ground-based gravitational wave detectors.

We invoke the estimates of the amplitudes of the velocity perturbations $f_R$ and $f_\theta$ caused by the influence of a spiral density wave that have been obtained by us previously from three stellar samples. These include Galactic masers with measured VLBI trigonometric parallaxes and proper motions, OB2 stars, and Cepheids. From these data we have obtained new estimates of the spiral pattern speed in the Galaxy $\Omega_p:$ $24.61\pm2.06$, $24.71\pm1.29$ and $25.98\pm1.37$~km s$^{-1}$ kpc$^{-1}$ from the samples of masers, OB2 stars, and Cepheids, respectively. The corotation radii for these three samples $R_{\rm cor}/R_0$ are $1.16\pm0.09$, $1.15\pm0.06$ and $1.09\pm0.06,$ suggesting that the corotation circle is located between the Sun and the Perseus arm segment.

The origin and evolution of fluorine in the Milky Way galaxy is still in debate. In particular, the increase of the [F/Fe] in metal-rich stars found from near-IR HF-lines is challenging to explain theoretically. We determine the fluorine abundances from 50 M giants in the solar neighbourhood spanning a broad range of metallicities (-0.9<[Fe/H]<0.25 dex). These stars are cool enough to have an array of HF lines in the K band. We observed the stars with the IGRINS and investigate each of ten HF molecular lines in detail. Based on a detailed line-by-line analysis of ten HF lines, we find that the R19, R18 and R16 lines should primarily be used for abundance analysis. The R15, R14 and R13 lines can also be used, but the trends based on these lines show increasing dependencies with the stellar parameters. The strongest HF lines, namely R12, R11, R9 and R7 should be avoided since the abundances from them show significant trends with the stellar parameters, and a high sensitivity to variations in the microturbulence, especially for coolest metal-rich stars. This leads to a huge scatter and high fluorine abundances for supersolar metallicity stars, not seen in the trends from the weaker lines for the same stars. When estimating the final mean fluorine abundance trend versus metallicity, we neglect the fluorine abundances from the four strongest lines (R7, R9, R11 and R12) for all stars and use only those derived from R16, R18, and R19 for the coolest metal-rich stars. We confirm the flat trend of [F/Fe] found in other studies in the metallicity range of -1.0<[Fe/H]<0.0. We also find a slight enhancement at supersolar metallicities (0<[Fe/H]<0.15) but we cannot confirm the upward trend seen at [Fe/H]>0.25. We need more observations of M giants at super solar metallicities with a spectrometer like IGRINS to confirm if the metal-rich fluorine abundance upturn is real or not.

GJ 229B, the first type-T brown dwarf to be discovered, has presented a tension between comparisons with evolutionary models and the larger-than-expected mass and radius values derived from spectroscopic and astrometric observations. We examine the hypothesis that GJ 229B is actually a binary sub-stellar object by using two grid-based fits using evolutionary models to explore the range of mass ratios of the possible binary components. We find that the best-fit component values are most consistent with a roughly 2:1 binary mass ratio and an age range of 2-6 Gyr. The observed temperatures, masses, and apparent radii match expected values from evolutionary models for a binary much better than a single-object model, but more detailed observations and modeling are needed to definitively confirm the binary hypothesis.

In this article, I review past, current, and future advances on the study of radio-loud AGN (RLAGN; radio-loud quasars and radio galaxies) lifecycles exclusively in the remnant and restarting phases. I focus on their dynamics and energetics as inferred from radio observations while discussing their radiative lifetimes, population statistics, and trends in their physical characteristics. I briefly summarize multi-wavelength observations, particularly X-rays, that have enabled studies of the large-scale environments of RLAGN in order to understand their role in feedback. Furthermore, I discuss analytic and numerical simulations that predict key properties of remnant and restarting sources as found in wide-area surveys, and discuss the prospects of future surveys that may shed further light on these elusive subpopulations of RLAGN.

This paper reports the first detection of polarization in the X-rays for atoll-source 4U 1820-303, obtained with the Imaging X-ray Polarimetry Explorer (IXPE) at 99.999% confidence level (CL). Simultaneous polarimetric measurements were also performed in the radio with the Australia Telescope Compact Array (ATCA). The IXPE observations of 4U 1820-303 were coordinated with Swift-XRT, NICER, and NuSTAR aiming to obtain an accurate X-ray spectral model covering a broad energy interval. The source shows a significant polarization above 4 keV, with a polarization degree of 2.0(0.5)% and a polarization angle of -55(7) deg in the 4-7 keV energy range, and a polarization degree of 10(2)% and a polarization angle of -67(7) deg in the 7-8 keV energy bin. This polarization also shows a clear energy trend with polarization degree increasing with energy and a hint for a position-angle change of about 90 deg at 96% CL around 4 keV. The spectro-polarimetric fit indicates that the accretion disk is polarized orthogonally to the hard spectral component, which is presumably produced in the boundary/spreading layer. We do not detect linear polarization from the radio counterpart, with a 99.97% upper limit of 50% at 7.25 GHz.

Aims. A transiting planet candidate with a sub-Neptune radius orbiting the nearby ($d$ = 51.9$\pm$0.07 pc) M1.5 V star TOI-1470 with a period of $\sim$2.5 d was announced by the NASA Transiting Exoplanet Survey Satellite (TESS), which observed the field of TOI-1470 in four different sectors. We aim to validate its planetary nature using precise radial velocities (RVs) taken with the CARMENES spectrograph. Methods. We obtained 44 RV measurements with CARMENES spanning eight months between 3 June 2020 and 17 January 2021. For a better characterization of the parent star activity, we also collected contemporaneous optical photometric observations at the Joan Or\'o and Sierra Nevada Observatories, and we retrieved archival photometry from the literature. We used ground-based photometric observations from MuSCAT and also from MuSCAT2 and MuSCAT3 to confirm the planetary transit signals. We performed a combined photometric and spectroscopic analysis by including Gaussian processes and Keplerian orbits to simultaneously account for the stellar activity and planetary signals. Results. We estimate that TOI-1470 has a rotation period of 29$\pm$3 d based on photometric and spectroscopic data. The combined analysis confirms the discovery of the announced transiting planet, TOI-1470 b, with an orbital period of 2.527093$\pm$0.000003 d, a mass of $7.32^{+1.21}_{-1.24}$ M$_{\oplus}$, and a radius of $2.18^{+0.04}_{-0.04}$ R$_{\oplus}$. We also discover a second transiting planet that was not announced previously by TESS, TOI-1470 c, with an orbital period of 18.08816$\pm$0.00006 d, a mass of $7.24^{+2.87}_{-2.77}$ M$_{\oplus}$, and a radius of $2.47^{+0.02}_{-0.02}$ R$_{\oplus}$. The two planets are placed on the same side of the radius valley of M dwarfs and lie between TOI-1470 and the inner border of its habitable zone.

We report ALMA high-angular resolution (~ 50 au) observations of the binary system SVS13-A. More specifically, we analyse deuterated water (HDO) and sulfur dioxide (SO2) emission. The molecular emission is associated with both the components of the binary system, VLA4A and VLA4B. The spatial distribution is compared to that of formamide (NH2CHO), previously analysed in the system. Deuterated water reveals an additional emitting component spatially coincident with the dust accretion streamer, at a distance larger than 120 au from the protostars, and at blue-shifted velocities (> 3 km/s from the systemic velocities). We investigate the origin of the molecular emission in the streamer, in light of thermal sublimation temperatures calculated using updated binding energies (BE) distributions. We propose that the observed emission is produced by an accretion shock at the interface between the accretion streamer and the disk of VLA4A. Thermal desorption is not completely excluded in case the source is actively experiencing an accretion burst.

In this study, we use Keck/KCWI spectroscopy to measure the age, metallicity and recessional velocity of NGC~1052-DF9 (DF9), a dwarf galaxy in the NGC~1052 group. We compare these properties to those of two other galaxies in the group, NGC~1052-DF2 and NGC~1052-DF4, which have low dark matter content. The three galaxies are proposed constituents of a trail of galaxies recently hypothesised to have formed as part of a ``bullet dwarf'' collision. We show that the ages and total metallicities of the three galaxies are within uncertainties of one another which may be expected if they share a related formation pathway. However, the recessional velocity we recover for DF9 (1680 $\pm$ 10 km s$^{-1}$) is higher than predicted for a linearly projected interpretation of the ``bullet dwarf'' trail. DF9 is then either not part of the trail or the correlation of galaxy velocities along the trail is not linear in 2D projection due to their 3D geometry. After examining other proposed formation pathways for the galaxies, none provide a wholly satisfactory explanation for all of their known properties. We conclude further work is required to understand the formation of this interesting group of galaxies.

The analysis of the available TESS light curves of $\alpha$ Sex (HD 87887) reveals low-frequency pulsations with a period of about 9.1 hours in this spectroscopic A0 III standard star. The IUE observations in December 1992 reveal large flux variations both in the far UV and in the mid UV which are accompanied by variations of the brightness in the V band recorded by the Fine Error Sensor on board IUE. The ultraviolet variability could be due to an eclipse by an hitherto undetected companion of smaller radius, possibly 2.5 R$_{\odot}$ but this needs confirmation by further monitoring possibly with TESS. An abundance determination yields solar abundances for most elements. Only carbon and strontium are underabundant and titanium, vanadium and barium mildly overabundant. Identification is provided for most of the lines absorbing more than 2% in the optical spectrum of $\alpha$ Sex. Stellar evolution modeling shows that $\alpha$ Sex is near the terminal-age main sequence, and its mass, radius and age are estimated to be $M = 2.57 \pm 0.32$ M$_{\odot}$, $R = 3.07 \pm 0.90$ R$_{\odot}$, $A = 385 \pm 77$ Myr, respectively.

During very early age of neutron stars, the core cools down faster compared to the crust creating a large thermal gradient in the interior of the star. During $10-100$ years, a cooling wave propagates from the core to the crust causing the interior of the star to thermalize. During this duration thermal properties of the core material is of great importance to understand the dynamics of the interior of the star. The heat capacity and thermal conductivity of the core depends on the behaviour of matter inside the core. We investigate these two properties in case of magnetars. Due to presence of large magnetic field, the proton superconductivity is quenched partially inside the magnetars depending upon the comparative values of upper critical field and the strength of the magnetic field present. This produces non-uniformity in the behaviour of matter throughout the star. Moreover, such non-uniformity arises from the variation of nature of the pairing and values of the pairing gap energy. We find that the heat capacity is substantially reduced due to the presence of superfluidity. On the other hand, the thermal conductivity of neutron is enhanced due to proton superconductivity and gets reduced due to neutron superfluidity. Hence, the variation of the thermal properties due to superfluidity in presence of magnetic field is different at different radius inside the star. However, in all the cases the %minimum maximum variation is of the order one. This affects the thermal relaxation time of the star and eventually its the thermal evolution.

We study the region around Collinder 121 (Cr 121) using newly available 6-dimensional data from the Gaia DR3 catalogue. Situated in the third quadrant, near the galactic plane, Collinder 121 lies in the region of Canis Major centred around l = 236 degrees, b = -10 degrees. Previous studies have suggested that the stellar associations in this region comprise an OB association (CMa OB2) lying at about 740 pc with a more distant open cluster (Cr 121) at approximately 1170 pc. Despite these studies, the precise nature of Collinder 121 remains uncertain. This study investigates the region bounded by the box l = 225 to 245 degrees, b = 0.00 to -20.00 degrees to a depth of 700 pc from 500 to 1200 pc which fully encompasses the region discussed in the literature. Using Gaia DR3 data, we do not find associations at the distances given in the literature. Instead, using the HDBSCAN machine learning algorithm, we find a major association of OB stars centred around 803 pc. Within this association we find four smaller subgroups that may be indicative of a larger association and which are located at a mean distance of 827 pc. Proper motion studies find coherence between these four subgroups and show a distinctive east to west increase in the size of the velocity vectors which supports contemporary studies that show similar trends in OB populations in Cygnus and within the Carina spiral Arm. Therefore, we hypothesize that Cr 121 and CMa OB2 are the same cluster, consistent with the 1977 study by Hoogerwerf.

We present a new orbit and mass solution for the four small satellites of Pluto: Styx, Nix, Kerberos, and Hydra. We have reanalyzed all available observations of the Pluto system obtained by the Hubble Space Telescope (HST) from 2005 to 2019 with the ACS, WFPC2, and WFC3 instruments, as well as the New Horizons LORRI images taken on approach to Pluto in 2015. We have used this high precision astrometry to produce updated orbits and mass estimates with uncertainties for all four of the small satellites. We find that the masses of Nix and Hydra are smaller than previously published estimates, with a dynamical mass of 1.8$\pm$0.4$\times$10$^{-3}$ km$^3$/s$^2$ (2.7$\pm$0.6$\times$10$^{16}$ kg) for Nix and 2.0$\pm$0.2$\times$10$^{-3}$ km$^3$/s$^2$ (3.0$\pm$0.3$\times$10$^{16}$ kg) for Hydra. These masses are 60% and 63% of the mean estimates by Brozovic et al. (2015), although still consistent with their 1-sigma uncertainties, and correspond to densities of 1.0$\pm$0.2 g/cm$^3$ for Nix and 1.2$\pm$0.2 g/cm$^3$ for Hydra given the moon volume estimates from Porter et al (2021). Although these densities are consistent with a range of ice-rock compositions, depending on the unknown bulk porosity in the moon interiors, the moons' high albedos and predominantly icy surfaces are most easily explained if their interiors are ice-rich. The tiny masses of Kerberos and Sytx remain very poorly constrained; we find 1-$\sigma$ upper limits for the dynamical mass of Styx to be 3$\times$10$^{-5}$ km$^3$/s$^2$ (5$\times$10$^{14}$ kg) and for Kerberos 5$\times$10$^{-5}$ km$^3$/s$^2$ (8$\times$10$^{14}$ kg), consistent with densities of $<$2.1 g/cm$^3$ for both bodies.

Observational data suggest that the ice shell on Enceladus is thicker at the equator than at the pole, indicating an equator-to-pole ice flow. If the ice shell is in an equilibrium state, the mass transport of the ice flow must be balanced by the freezing and melting of the ice shell, which in turn is modulated by the ocean heat transport. Here we use a numerical ocean model to study the ice-ocean interaction and ocean circulation on Enceladus with different salinities. We find that salinity fundamentally determines the ocean stratification. A stratified layer forms in the low salinity ocean, affecting the ocean circulation and heat transport. However, in the absence of tidal heating in the ice shell, the ocean heat transport is found to always be equatorward, resulting in freezing at the pole and melting at lower latitudes, which cannot maintain the ice shell geometry against the equator-to-pole ice flow. The simulation results suggest that either the ice shell on Enceladus is not in an equilibrium state, or tidal dissipation in the ice shell is important in maintaining the ice shell geometry.

Optical rest-frame spectroscopic diagnostics are usually employed to distinguish between star formation and AGN-powered emission. However, this method is biased against dusty sources, hampering a complete census of the AGN population across cosmic epochs. To mitigate this effect, it is crucial to observe at longer wavelengths in the rest-frame near-infrared (near-IR), which is less affected by dust attenuation and can thus provide a better description of the intrinsic properties of galaxies. AGN diagnostics in this regime have not been fully exploited so far, due to the scarcity of near-IR observations of both AGNs and star-forming galaxies, especially at redshifts higher than 0.5. Using Cloudy photoionization models, we identify new AGN - star formation diagnostics based on the ratio of bright near-infrared emission lines, namely [SIII] 9530 Angstrom, [CI] 9850 Angstrom, [PII] 1.188 $\mu m$, [FeII] $1.257 \mu m$, and [FeII] $1.64 \mu m$ to Paschen lines (either Pa$\gamma$ or Pa$\beta$), providing simple, analytical classification criteria. We apply these diagnostics to a sample of 64 star-forming galaxies and AGNs at 0 < z < 1, and 65 sources at 1 < z < 3 recently observed with JWST-NIRSpec in CEERS. We find that the classification inferred from the near-infrared is broadly consistent with the optical one based on the BPT and the [SII]/H$\alpha$ ratio. However, in the near-infrared, we find $\sim 60 \%$ more AGNs than in the optical (13 instead of 8), with 5 sources classified as 'hidden' AGNs, showing a larger AGN contribution at longer wavelengths, possibly due to the presence of optically thick dust. The diagnostics we present provide a promising tool to find and characterize AGNs from z=0 to z=3 with low and medium-resolution near-IR spectrographs in future surveys.

We present the results of our identification of 14 X-ray sources detected in the eastern Galactic sky ($0<l<180 \circ$ ) in the 4-12 keV energy band on the combined map of the first five all-sky surveys (from December 2019 to March 2022) with the Mikhail Pavlinsky ART-XC telescope onboard the SRG observatory. All 14 sources are reliably detected by the SRG/eROSITA telescope in the 0.2-8 keV energy band. Six of them have been detected in X-rays for the first time, while the remaining ones have already been known previously as X-ray sources, but their nature has remained unknown. We have taken optical spectra for 12 sources with the 1.6-m AZT-33IK telescope at the Sayan Observatory (the Institute of Solar-Terrestrial Physics, the Siberian Branch of the Russian Academy of Sciences). For two more objects we have analyzed the archival spectra taken during the 6dF survey. All objects have turned out to be Seyfert galaxies (one NLSy1, three Sy1, four Sy1.9, and six Sy2) at redshifts $z=0.015-0.238$. Based on data from the eROSITA and ART-XC telescopes onboard the SRG observatory, we have obtained X-ray spectra for all objects in the energy range 0.2-12 keV. In four of them the intrinsic absorption exceeds $N_{\rm H}>10^{22}$ cm$^{-2}$ at a 90% confidence level, with one of them being probably heavily obscured ($N_{\rm H}>5\times 10^{22}$ cm$^{-2}$ with 90% confidence). This paper continues our series of publications on the identification of hard X-ray sources detected during the all-sky survey with the SRG orbital X-ray observatory.

Planets orbiting close to hot stars experience intense extreme-ultraviolet radiation, potentially leading to atmosphere evaporation and to thermal dissociation of molecules. However, this extreme regime remains mainly unexplored due to observational challenges. Only a single known ultra-hot giant planet, KELT-9b, receives enough ultraviolet radiation for molecular dissociation, with a day-side temperature of ~4,600K. An alternative approach uses irradiated brown dwarfs as hot-Jupiter analogues. With atmospheres and radii similar to those of giant planets, brown dwarfs orbiting close to hot Earth-sized white-dwarf stars can be directly detected above the glare of the star. Here we report observations revealing an extremely irradiated low-mass companion to the hot white dwarf WD0032-317. Our analysis indicates a day-side temperature of ~8,000K, and a day-to-night temperature difference of ~6,000K. The amount of extreme-ultraviolet radiation (with wavelengths 100-912\r{A}) received by WD0032-317B is equivalent to that received by planets orbiting close to stars as hot as a late B-type stars, and about 5,600 times higher than that of KELT-9b. With a mass of ~75-88 Jupiter masses, this near-hydrogen-burning-limit object is potentially one of the most massive brown dwarfs known.

We use the \texttt{GRUMPY} galaxy formation model based on a suite of zoom-in, high-resolution, dissipationless $\Lambda$ Cold Dark Matter ($\Lambda$CDM) simulations of the Milky Way (MW) sized haloes to examine total matter density within the half-mass radius of stellar distribution, $\rho_{\rm tot}(<r_{1/2})$, of satellite dwarf galaxies around the MW hosts and their mass assembly histories. We compare model results to $\rho_{\rm tot}(<r_{1/2})$ estimates for observed dwarf satellites of the Milky Way spanning their entire luminosity range. We show that observed MW dwarf satellites exhibit a trend of decreasing $\rho_{\rm tot}(<r_{1/2})$ with increasing stellar mass. This trend is in general agreement with the trend predicted by the model. None of the observed satellites are overly dense compared to the results of our $\Lambda$CDM-based model. We also show that although the halo mass of many satellite galaxies is comparable to the halo mass of the MW progenitor at $z> 10$, at these early epochs halos that survive as satellites to $z=0$ are located many virial radii away from the MW progenitors and thus do not have a chance to merge with it. Our results show that neither the densities estimated in observed Milky Way satellites nor their mass assembly histories pose a challenge to the $\Lambda$CDM model. In fact, the broad agreement between density trends with the stellar mass of the observed and model galaxies can be considered as yet another success of the model.

Based on high-quality APOGEE DR17 and Gaia DR3 data for 1,742 red giants stars within 5 kpc of the Sun and not rotating with the Galactic disc ($V_\phi <$ 100 km s$^{-1}$), we use the nonlinear technique of unsupervised analysis t-SNE to detect coherent structures in the space of ten chemical-abundance ratios: [Fe/H], [O/Fe], [Mg/Fe], [Si/Fe], [Ca/Fe], [C/Fe], [N/Fe], [Al/Fe], [Mn/Fe], and [Ni/Fe]. Additionally, we obtain orbital parameters for each star using the non-axisymmetric gravitational potential {\tt GravPot16}. Seven structures are detected, including the Splash, Gaia-Sausage-Enceladus (GSE), the high-$\alpha$ heated-disc population, N-C-O peculiar stars, and inner disk-like stars, plus two other groups that did not match anything previously reported in the literature, here named Galileo 5 and Galileo 6 (G5 and G6). These two groups overlap with Splash in [Fe/H], G5 being lower metallicity than G6, both between GSE and Splash in the [Mg/Mn] versus [Al/Fe] plane, G5 in the $\alpha$-rich in-situ locus, and G6 on the border of the $\alpha$-poor in-situ one; nonetheless their low [Ni/Fe] hints to a possible ex-situ origin. Their orbital energy distributions are between the Splash and GSE, with G5 being slightly more energetic than G6. We verified the robustness of all the obtained groups by exploring a large range of t-SNE parameters, applying it to various subsets of data, and also measuring the effect of abundance errors through Monte Carlo tests.

We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a period of $\approx$1000 days and an amplitude of $\Delta m \approx 0.6$ mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IR $J$ and $K_{s}$ bands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-SN outbursts in this time period. The IR brightness ($M_{K_s} = -10.7$ mag) and color ($J-K_{s} = 1.6$ mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature ($T_{\mathrm{eff}} = 3500_{-1400}^{+800}$ K) and luminosity ($\log L/L_{\odot} = 5.1\pm0.2$). This places the candidate among the most luminous Type II SN progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass of $M_{\mathrm{init}} = 17\pm4$ $M_{\odot}$. We estimate the pre-SN mass-loss rate of the star between 3-19 years before explosion from the SED modeling at $\dot M \approx 3\times10^{-5}$ to $3\times10^{-4}$ $M_{\odot}$ yr$^{-1}$, perhaps pointing to enhanced mass loss in a pulsation-driven wind.

The overcooling of cool core clusters is a persistent puzzle in the astrophysics of galaxy clusters. We propose that it may naturally be resolved via interactions between the baryons of the intracluster medium (ICM) and its dark matter (DM). DM-baryon interactions can inject heat into the ICM to offset bremmstrahlung cooling, but these interactions are also strongly constrained by existing experiments and astrophysical observations. We survey existing constraints and combine these with the energetic needs of an observed sample of cool core clusters. We find that a robust parameter space exists for baryon-DM scattering solutions to the cooling flow problem, provided that only a sub-component of DM interacts strongly with the baryons. Interestingly, baryon-DM scattering is a thermally stable heating source so long as the baryon temperature is greater than $1/3-1/2$ the DM temperature, a condition that seems to be satisfied observationally.

The morphological diversity of galaxies is a relevant probe of galaxy evolution and cosmological structure formation. However, in large sky surveys, even the morphological classification of galaxies into two classes, like late-type (LT) and early-type (ET), still represents a significant challenge. In this work we present a Deep Learning (DL) based morphological catalog built from images obtained by the Southern Photometric Local Universe Survey (S-PLUS) Data Release 3 (DR3). Our DL method achieves an precision rate of 98.5$\%$ in accurately distinguishing between spiral, as part of the larger category of late type (LT) galaxies, and elliptical, belonging to early type (ET) galaxies. Additionally, we have implemented a secondary classifier that evaluates the quality of each galaxy stamp, which allows to select only high-quality images when studying properties of galaxies on the basis of their DL morphology. From our LT/ET catalog of galaxies, we recover the expected color--magnitude diagram in which LT galaxies display bluer colors than ET ones. Furthermore, we also investigate the clustering of galaxies based on their morphology, along with their relationship to the surrounding environment. As a result, we deliver a full morphological catalog with $164314$ objects complete up to $r_{petro}<18$, covering $\sim 1800$ deg$^2$, including a significant area of the Southern hemisphere that was not covered by previous morphology catalogues.

We present a CCD $VI$ time-series analysis of the globular cluster NGC 6139 and its variable star population. Using the astrometric data available in $Gaia$-DR3 we performed a membership analysis that enabled the construction of a clean Colour-Magnitude Diagram (CMD). Variable stars in the field of the cluster reported by $Gaia$-DR3 and newly discovered ones in this paper are classified and their membership is critically evaluated. We report two cluster member RRc (V12, V15) and four SR (V13, V14, V17, V18) not previously detected and assign variable names to V11 and V16 detected by $Gaia$ as they proof to be cluster members. Light curves and periods for non-member $Gaia$ eclipsing binaries, semi regular variables and newly detected RR Lyrae stars are provided. Fourier decomposition of the light curves of the cluster member RRab and RRc stars leads to the values [Fe/H]=$-1.63$ dex, and distance of 9.63$\pm$0.68 kpc. The Oosterhoff type II nature of the cluster is confirmed. We adopted the mean reddening $E(B-V)$=0.786 mag and performed a differential reddening analysis based on the dispersion of the red giant branch. The differential map allowed a mild correction of the CMD.

The X-ray Silicon-On-Insulator (SOI) pixel sensor named XRPIX has been developed for the future X-ray astronomical satellite FORCE. XRPIX is capable of a wide-band X-ray imaging spectroscopy from below 1 keV to a few tens of keV with a good timing resolution of a few tens of $\mu$s. However, it had a major issue with its radiation tolerance to the total ionizing dose (TID) effect because of its thick buried oxide layer due to the SOI structure. Although new device structures introducing pinned depleted diodes dramatically improved radiation tolerance, it remained unknown how radiation effects degrade the sensor performance. Thus, this paper reports the results of a study of the degradation mechanism of XRPIX due to radiation using device simulations. In particular, mechanisms of increases in dark current and readout noise are investigated by simulation, taking into account the positive charge accumulation in the oxide layer and the increase in the surface recombination velocity at the interface between the sensor layer and the oxide layer. As a result, it is found that the depletion of the buried p-well at the interface increases the dark current, and that the increase in the sense-node capacitance increases the readout noise.

The abundance of refractory elements in giant planets can provide key insights into their formation histories. Due to the Solar System giants' low temperatures, refractory elements condense below the cloud deck limiting sensing capabilities to only highly volatile elements. Recently, ultra-hot giant exoplanets have allowed for some refractory elements to be measured showing abundances broadly consistent with the solar nebula with titanium likely condensed out of the photosphere. Here we report precise abundance constraints of 14 major refractory elements on the ultra-hot giant planet WASP-76b that show distinct deviations from proto-solar, and a sharp onset in condensation temperature. In particular, we find nickel to be enriched, a possible sign of the accretion of a differentiated object's core during the planet's evolution. Elements with condensation temperatures below 1,550 K otherwise closely match those of the Sun before sharply transitioning to being strongly depleted above 1,550 K, well explained by nightside cold-trapping. We further unambiguously detect vanadium oxide on WASP-76b, a molecule long hypothesized to drive atmospheric thermal inversions, and also observe a global east-west asymmetry in its absorption signals. Overall, our findings indicate that giant planets have a mostly stellar-like refractory elemental content and suggest that temperature sequences of hot Jupiter spectra can show abrupt transitions wherein a mineral species is either present, or completely absent if a cold-trap exists below its condensation temperature.

We outline a process to design, architect, and implement the next-generation Event Horizon Telescope (ngEHT): a transformative enhancement to the EHT, that will form a networked global array of radio dishes capable of making high-fidelity real-time movies of supermassive black holes (SMBH) and their emanating jets. This builds upon the EHT principally by deploying additional modest-diameter dishes to optimized geographic locations to enhance the current global mm/submm wavelength Very Long Baseline Interferometric (VLBI) array, which has, to date, utilized mostly pre-existing radio telescopes. The ngEHT program further focuses on observing at three frequencies simultaneously for increased sensitivity and Fourier spatial frequency coverage. Here, the concept, science goals, design considerations, station siting and instrument prototyping is discussed, and a preliminary reference array to be implemented in phases is described.

X-ray polarization from the Imaging X-ray Polarimetry Explorer (IXPE) provides an important new probe of the geometry of the pulsar emission zone and of particle acceleration in the surrounding pulsar wind nebula (PWN). However, with IXPE's modest ~20-30" spatial resolution, separation of the pulsar signal from the nebula is a challenge. Conventional analysis defines an "off" phase window as pure nebular emission and subtracts its polarization to isolate the phase-varying pulsar ("on-off fitting"). We present a more sensitive scheme that uses external measurements of the nebula structure and pulsar light curve to isolate their contributions to the phase- and spatially-varying polarization via least-squares regression ("simultaneous fitting"). Tests with simulation data show ~30% improvement in pulse phase polarization uncertainties, decreased background systematics, and substantially improved nebular polarization maps. Applying "simultaneous fitting" to early IXPE Crab data extracts additional phase bins with significant polarization. These bins show interesting departures from the well-known optical polarization sweeps, although additional exposure will be needed for precise model confrontation.

The characteristic orbital period of the inner-most objects within the galactic census of planetary and satellite systems appears to be nearly universal, with $P$ on the order of a few days. This paper presents a theoretical framework that provides a simple explanation for this phenomenon. By considering the interplay between disk accretion, magnetic field generation by convective dynamos, and Kelvin-Helmholtz contraction, we derive an expression for the magnetospheric truncation radius in astrophysical disks, and find that the corresponding orbital frequency is independent of the mass of the host body. Our analysis demonstrates that this characteristic frequency corresponds to a period of $P\sim3$ days, although intrinsic variations in system parameters are expected to introduce a factor of $\sim2-3$ spread in this result. Standard theory of orbital migration further suggests that planets should stabilize at an orbital period that exceeds disk truncation by a small margin. Cumulatively, our findings predict that the periods of close-in bodies should span $P\sim2-12$ days - a range that is consistent with observations.

The connection between inner small planets and outer giant planets is crucial to our understanding of planet formation across a wide range of orbital separations. While Kepler provided a plethora of compact multi-planet systems at short separations ($\lesssim 1$ AU), relatively little is known about the occurrence of giant companions at larger separations and how they impact the architectures of the inner systems. Here, we use the catalog of systems from the Kepler Giant Planet Search (KGPS) to study how the architectures of the inner transiting planets correlate with the presence of outer giant planets. We find that for systems with at least three small transiting planets, the distribution of inner-system gap complexity ($\mathcal{C}$), a measure of the deviation from uniform spacings, appears to differ ($p \lesssim 0.02$) between those with an outer giant planet ($50 M_\oplus \leq M_p\sin{i} \leq 13 M_{\rm Jup}$) and those without any outer giants. All four inner systems (with 3+ transiting planets) with outer giant(s) have a higher gap complexity ($\mathcal{C} > 0.32$) than 79% (19/24) of the inner systems without any outer giants (median $\mathcal{C} \simeq 0.06$). This suggests that one can predict the occurrence of outer giant companions by selecting multi-transiting systems with highly irregular spacings. We do not find any correlation between outer giant occurrence and the size (similarity or ordering) patterns of the inner planets. The larger gap complexities of inner systems with an outer giant hints that massive external planets play an important role in the formation and/or disruption of the inner systems.

COCONUT is a global coronal magnetohydrodynamic model recently developed. In order to achieve robustness and fast convergence to steady-state, several assumptions have been made during its development, such as prescribing filtered photospheric magnetic maps for representing the magnetic field in the lower corona. This filtering leads to smoothing and lower magnetic field values at the inner boundary, resulting in an unrealistically high plasma beta.In this paper, we examine the effects of prescribing such filtered magnetograms and formulate a method for achieving more realistic plasma beta values without losing computational performance. We demonstrate the effects of the highly pre-processed magnetic maps and the resulting high plasma beta on the features in the domain. Then, in our new approach, we shift the inner boundary to 2 Rs and preserve the prescribed highly filtered magnetic map. This effectively reduces the prescribed plasma beta and leads to a more realistic setup. The method is applied on a magnetic dipole, a minimum (2008) and a maximum (2012) solar activity case, to demonstrate its effects. The results obtained with the proposed approach show significant improvements in the resolved density and radial velocity profiles, and far more realistic values of the plasma \{beta} at the boundary and inside the computational domain. This is also demonstrated via synthetic white light imaging and with the validation against tomography. The computational performance comparison shows similar convergence to a limit residual on the same grid when compared to the original setup. Considering that the grid can be further coarsened with this new setup, the operational performance can be additionally increased if needed. The newly developed method is thus deemed as a good potential replacement of the original setup for operational purposes.

We consider the polarized Sunyaev-Zel'dovich (pSZ) effect for a tomographic probe of cosmic birefringence, including all relevant terms of the pSZ effect in the CMB observables, some of which were ignored in the previous works. The pSZ effect produces late-time polarization signals from the scattering of the local temperature quadrupole seen by an electron. We forecast the expected constraints on cosmic birefringence at the late time of the universe with the pSZ effect. We find that the birefringence angles at $2\lesssim z\lesssim 5$ are constrained at a sub-degree level by the cross-correlations between CMB $E$- and $B$-modes or between CMB $B$-modes and remote quadrupole $E$-modes using data from LiteBIRD, CMB-S4, and LSST. In particular, the cross-correlation between large-scale CMB $B$-modes and remote-quadrupole $E$-modes has a much smaller bias from the Galactic foregrounds and is useful to cross-check the results from the $EB$ power spectrum.

Hydrogen-rich Type II supernovae (SNe II) are the most frequently observed class of core-collapse SNe (CCSNe). However, most studies that analyse large samples of SNe II lack events with absolute peak magnitudes brighter than -18.5 mag at rest-frame optical wavelengths. Thanks to modern surveys, the detected number of such luminous SNe II (LSNe II) is growing. There exist several mechanisms that could produce luminous SNe II. The most popular propose either the presence of a central engine (a magnetar gradually spinning down or a black hole accreting fallback material) or the interaction of supernova ejecta with circumstellar material (CSM) that turns kinetic energy into radiation energy. In this work, we study the light curves and spectral series of a small sample of six LSNe II that show peculiarities in their H$\alpha$ profile, to attempt to understand the underlying powering mechanism. We favour an interaction scenario with CSM that is not dense enough to be optically thick to electron scattering on large scales -- thus, no narrow emission lines are observed. This conclusion is based on the observed light curve (higher luminosity, fast decline, blue colours) and spectral features (lack of persistent narrow lines, broad H$\alpha$ emission, lack of H$\alpha$ absorption, weak or nonexistent metal lines) together with comparison to other luminous events available in the literature. We add to the growing evidence that transients powered by ejecta-CSM interaction do not necessarily display persistent narrow emission lines.

We provide a detailed analysis on tidal tails of the globular cluster M3 (NGC 5272). We first discover clear extra-tidal structures with slight S-shape near the cluster. This inspires us to examine the existence of its long tidal tails. We highlight potential stream stars using proper motions (PMs) of a model stream combined with the cluster's locus in a color-magnitude diagram (CMD). A 35 deg long leading tail and a 21 deg long trailing tail are successfully detected at the same time. Their corresponding overdensities can be recognized in CMD and PM space after subtracting background. We estimate stream width, star number density and surface brightness for both tails, as well as the distance variation along the entire stream. We then verify the connection of M3 and the Sv\"{o}l stream. Finally, we tabulate 11 member stars belonging to the M3 tidal stream with available spectroscopic observations.

The cosmic-ray database, CRDB, has been gathering cosmic-ray data for the community since 2013. We present a new release, CRDV v4.1, providing many new quantities and data sets, with several improvements made on the code and web interface, and with new visualisation tools. CRDB relies on the mySQL database management system, jquery and tsorter libraries for queries and sorting, and php web pages and ajax protocol for displays. A REST interface enables user queries from command line or scripts. A new (pip-installable) CRDB python library is developed and extensive jupyter notebook examples are provided. This release contains cosmic-ray dipole anisotropy data, high-energy $\bar{p}/p$ upper limits, some unpublished LEE and AESOP lepton time series, many more ultra-high energy data, and a few missing old data sets. It also includes high-precision data from the last three years, in particular the hundreds of thousands AMS-02 and PAMELA data time series (time-dependent plots are now enabled).All these data are shown in a gallery of plots, which can be easily reproduced from the public notebook examples. CRDB contains 314902 data points from 487 publications, in 4092 sub-experiments from 126 experiments.

Protoplanets accreting matter seem to feature a strong emission at $\text{H}\alpha$ wavelength and the contrast of the system is then weaker at these wavelengths and thus more accessible. Measuring this emission from a protoplanet allow to characterize the accretion processes hence to study the mechanisms of planetary formation. The Fibered Interferometer foR a Single Telescope (FIRST) instrument uses the concept of fibered pupil remapping in the visible, integrated on the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) telescope. Its principle is to sub-divide the entrance pupil into sub-pupils whose beams are injected into single mode optical fibers. The lattest have the double advantage of applying a spatial filtering on the wavefront thus suppressing the optical aberrations at the scale of the sub-pupils, at the same time as the remapping of the sub-pupils. FIRST has demonstrated that it can detect stellar binaries at a resolution below the diffraction limit of the telescope with this concept. The replica of a new version of this instrument (FIRSTv2) has been built in laboratory to allow the development of the integrated optics technology. My thesis work demonstrates its feasibility for HRA imaging, in order to improve the contrast performances of the instrument. The fringes of each pair-wise combination of sub-pupils are measured via a temporal sampling and their phase are estimated. Data analysis estimates the spectral differential phase, which is a self-calibrated observable of the atmospheric and instrumental perturbations. Thus, from a protoplanetary source simulator, I was able to demonstrate that FIRSTv2 could detect a protoplanetary companion at a separation equivalent to 0.7 l / B from the central source, with a contrast of about 0.5. Finally, FIRSTv2 were integrated on the SCExAO bench for its first light.

Characterizing the interior structure of exoplanets is essential for understanding their diversity, formation, and evolution. As the interior of exoplanets is inaccessible to observations, an inverse problem must be solved, where numerical structure models need to conform to observable parameters such as mass and radius. This is a highly degenerate problem whose solution often relies on computationally-expensive and time-consuming inference methods such as Markov Chain Monte Carlo. We present ExoMDN, a machine-learning model for the interior characterization of exoplanets based on Mixture Density Networks (MDN). The model is trained on a large dataset of more than 5.6 million synthetic planets below 25 Earth masses consisting of an iron core, a silicate mantle, a water and high-pressure ice layer, and a H/He atmosphere. We employ log-ratio transformations to convert the interior structure data into a form that the MDN can easily handle. Given mass, radius, and equilibrium temperature, we show that ExoMDN can deliver a full posterior distribution of mass fractions and thicknesses of each planetary layer in under a second on a standard Intel i5 CPU. Observational uncertainties can be easily accounted for through repeated predictions from within the uncertainties. We use ExoMDN to characterize the interior of 22 confirmed exoplanets with mass and radius uncertainties below 10% and 5% respectively, including the well studied GJ 1214 b, GJ 486 b, and the TRAPPIST-1 planets. We discuss the inclusion of the fluid Love number $k_2$ as an additional (potential) observable, showing how it can significantly reduce the degeneracy of interior structures. Utilizing the fast predictions of ExoMDN, we show that measuring $k_2$ with an accuracy of 10% can constrain the thickness of core and mantle of an Earth analog to $\approx13\%$ of the true values.

The powerful, high-energy magnetic activities of young stars play important roles in the magnetohydrodynamics in the innermost parts of the protoplanetary disks. In addition, the associated UV and X-ray emission dictates the photochemistry; moreover, the corona activities can affect the atmosphere of a newborn extra-solar planet. How the UV and X-ray photons are generated, and how they illuminate the disks, are not well understood. Here we report the analyses of the optical and infrared (OIR) photometric monitoring observations and the high angular-resolution centimeter band images of the low-mass (M1 type) pre-main sequence star, DM Tau. We found that the OIR photometric light curves present periodic variations, which is consistent with that the host young star is rotating in the same direction as the natal disk and is hosting at least one giant cold spot. In addition, we resolved that the ionized gas in the DM Tau disk is localized, and its spatial distribution is varying with time. All the present observations can be coherently interpreted, if the giant cold spot is the dominant anisotropic UV and/or X-ray source that illuminates the ambient cone-like region. These results indicate that a detailed theoretical model of the high-energy protostellar emission is essential in the understanding of the space weather around the extra-solar planets and the origin of life.

IceCube collaboration has previously reported an evidence for neutrino signal from a Seyfert galaxy NGC 1068. This may suggest that all Seyfert galaxies emit neutrinos. To test this hypothesis, we identify the best candidate neutrino sources among nearby Seyfert galaxies, based on their hard X-ray properties. Only two other sources, NGC 4151 and NGC 3079 are expected to be detectable in 10 years of IceCube data. We find an evidence for neutrino signal from both sources in publicly available ten-year IceCube dataset. The chance coincidence probability to find the observed neutrino count excesses in the directions of the two out of two expected sources, in addition to the previously reported brightest source, is p<2.6e-7. This corresponds to the detection of Seyfert galaxies as a neutrino source class.

We present a comprehensive analysis of the Galactic mid-infrared (MIR) bubble [HKS2019] E70 (E70) by adopting a multi-wavelength approach to understand the physical environment and star formation scenario around it. We identified a small (radius ~1.7 pc) stellar clustering inside the E70 bubble and its distance is estimated as 3.26 +/- 0.45 kpc. This cluster is embedded in the molecular cloud and hosts massive stars as well as young stellar objects (YSOs), suggesting active star formation in the region. The spectral type of the brightest star 'M1' of the E70 cluster is estimated as O9V and a circular ring/shell of gas and dust is found around it. The diffuse radio emission inside this ring/shell, the excess pressure exerted by the massive star 'M1' at the YSOs core, and the distribution of photo-dissociation regions (PDRs), a Class I YSO, and two ultra-compact (UC) H II regions on the rim of this ring/shell, clearly suggest positive feedback of the massive star 'M1' in the region. We also found a low-density shell-like structure in 12 CO(J=1-0) molecular emission along the perimeter of the E70 bubble. The velocity structure of the 12 CO emission suggests that the feedback from the massive star appears to have expelled the molecular material and subsequent swept-up material is what appears as the E70 bubble.

It is discussed in detail the complete mathematical model of gravitational lensing on a single cosmic string (CS) of general shape and position with respect to the line of sight. CS are one-dimensional extended objects assuredly predicted by modern cosmology. The presence of CS changes the global geometry of the Universe, could clarify the properties of the early Universe, including inflation models, and could serve as a unique proof of higher-dimensional theories. Despite the fact that CS have not yet been reliable detected, there are several strong independent indications of the existence of the CS, based of CMB analysis and search of gravitational lens chains with special properties. However, early considered models of straight CS presented only a small fraction of the general CS-configurations to be observed. Now we propose model which could significantly increase the possibilities of CS observational search. It is considered more realistic models have necessarily include the inclinations and bends of the CS. Besides, the recent analysis of observational data on the search for gravitational-lens candidates, shows a large number of pairs that could be explained by the complex geometry of the CS.

Many molecules observed in the interstellar medium are thought to result from thermal desorption of ices. Parameters such as desorption energy and pre-exponential frequency factor are essential to describe the desorption of molecules. Experimental determinations of these parameters are missing for many molecules, including those found in the interstellar medium. The objective of this work is to expand the number of molecules for which desorption parameters are available, by collecting and re-analysing experimental temperature programmed desorption data that are present in the literature. Transition State Theory (TST) is used in combination with the Redhead equation to determine desorption parameters. Experimental data and molecular constants (e.g., mass, moment of inertia) are collected and given as input. Using the Redhead-TST method, the desorption parameters for 133 molecules have been determined. The Redhead-TST method is found to provide reliable results that agree well with desorption parameters determined with more rigorous experimental methods. The importance of using accurately determined pre-exponential frequency factors to simulate desorption profiles is emphasised. The large amount of data allows to look for trends, the most important is the relationship log$_{10}$($\nu$) = 2.65ln($m$) + 8.07, where $\nu$ is the pre-exponential frequency factor and $m$ the mass of the molecule. The data collected in this work allow to model the thermal desorption of molecules and help understand changes in chemical and elemental composition of interstellar environments.

The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective stars exhibit very similar period-activity relations, hinting that dynamos generated by stars with and without radiative cores hold similar properties. Here, using kinematic ages, we discover an abrupt change in the stellar spin-down law across the fully convective boundary. We found that fully convective stars exhibit a higher angular momentum loss rate, corresponding to a torque that is $\sim$ 2.25 times higher for a given angular velocity than partially convective stars around the fully convective boundary. This requires a dipole field strength that is larger by a factor of $\sim$2.5, a mass loss rate that is $\sim$4.2 times larger, or some combination of both of those factors. Since stellar-wind torques depend primarily on large-scale magnetic fields and mass loss rates, both of which derive from magnetic activity, the observed abrupt change in spin-down law suggests that the dynamos of partially and fully convective stars may be fundamentally different

Dramatic recent progress in understanding galactic chemical evolution (GCE) has been driven partly by direct observations of the distant past with HST and JWST of stellar abundances from giant high-resolution spectroscopic surveys (APOGEE, GALAH) and the complementary power of Gaia astrometry and photometry. Focusing on archaeology, I give a rapid-fire, and I hope synthesizing, review of work my collaborators and I have done on theoretical modeling and observational interpretation. I discuss (1) the interleaved but distinguishable roles of stellar scale astrophysics and galactic scale astrophysics in governing GCE, (2) the use of abundance ratio trends to empirically infer nucleosynthetic yields, (3) the uncertainty in the overall scale of yields and its degeneracy with the importance of galactic outflows, (4) the emergence of equilibrium in GCE, (5) the dimensionality of the stellar distribution in chemical abundance space, and (6) insights from chemical abundances on the early history of the Milky Way, including measurements of the intrinsic scatter of abundance ratios in metal-poor stars (-2 < [Fe/H] < -1), suggesting that a typical halo star at this metallicity is enriched by the products of $N \sim 50$ supernovae mixed over $\sim 10^5 M_\odot$ of star-forming gas.

Recent observations with the \textit{James Webb} Space Telescope (JWST) have further refined the spectroscopic redshift of GN-z11, one of the most distant galaxies identified with the \textit{Hubble} Space Telescope (HST) at $z=10.603$. The presence of extremely dense gas ($>10^{10}$ cm$^{-3}$), the detection of high-ionisation lines and of CII*1335 emission, as well as the presence of an ionisation cone, indicate that GN-z11 also hosts an Active Galactic Nucleus (AGN). Further photometric and spectroscopic follow-up demonstrates that it lies in a large-scale, overdense structure with possible signatures of Population III (PopIII) stars in its halo. Surprisingly, Ly$\alpha$ has also been detected despite the expected largely neutral inter-galactic medium at such a redshift. We exploit recent JWST/NIRSpec IFU observations to demonstrate that the Ly$\alpha$ emission in GN-z11 is part of an extended halo with a minimum size of 0.8--3.2 kpc, depending on the definition used to derive the halo size. The surface brightness of the Ly$\alpha$ halo around GN-z11 appears consistent with Ly$\alpha$ halos observed around $z\sim6$ quasars. At the wavelength of Ly$\alpha$ at $z\sim$10.6, we identify three other emission line candidates within the IFU Field-of-View with no UV rest-frame counterpart visible in deep images from the JWST/NIRCam. If confirmed, this could be the first evidence that the local region of GN-z11 represents a candidate protocluster core, forming just 400 Myr after the Big Bang. We give a first estimate of the dark matter halo mass of this structure ($M_h$=2.96$^{+0.44}_{-0.39} \times$10$^{10}$ M$_{\odot}$), consistent with a Coma-like cluster progenitor.

Interplanetary solar radio type III bursts provide the means for remotely studying and tracking energetic electrons propagating in the interplanetary medium. Due to the lack of direct radio source imaging, several methods have been developed to determine the source positions from space-based observations. Moreover, none of the methods consider the propagation effects of anisotropic radio-wave scattering, which would strongly distort the trajectory of radio waves, delay their arrival times, and affect their apparent characteristics. We investigate the source positions and directivity of an interplanetary type III burst simultaneously observed by Parker Solar Probe, Solar Orbiter, STEREO, and Wind and compare the results of applying the intensity fit and timing methods with ray-tracing simulations of radio-wave propagation with anisotropic density fluctuations. The simulation calculates the trajectories of the rays, their time profiles at different viewing sites, and the apparent characteristics for various density fluctuation parameters. The results indicate that the observed source positions are displaced away from the locations where emission is produced, and their deduced radial distances are larger than expected from density models. This suggests that the apparent position is affected by anisotropic radio-wave scattering, which leads to an apparent position at a larger heliocentric distance from the Sun. The methods to determine the source positions may underestimate the apparent positions if they do not consider the path of radio-wave propagation and incomplete scattering at a viewing site close to the intrinsic source position.

NSV 14264 and NSV 14172 are suspected to be variable stars of RR Lyr type (Brun, 1964). They were observed during three nights in October 2018 with a 25cm diameter telescope. These observations completed by ASAS-SN survey data bring to the conclusion that these two stars are not RR Lyraes but constant stars in the limit of the precision of the present photometry. The analysis of GAIA data allows to say that NSV 14264 is a main sequence dwarf similar to the Sun but that NSV 14172 is a yellow giant star located in the HR diagram at the limit between RR Lyraes and CW cepheids; however, it does not pulsate with significant amplitude.

Radio galaxies are a key population to understand the importance of relativistic jets in AGN feedback. We present the results of a systematic, broadband X-ray spectral analysis of hard X-ray selected radio galaxies to investigate their nuclear structures. In this study, we focus on the seven most radio-loud, X-ray obscured narrow line radio galaxies in the \textit{Swift}/BAT 70 month AGN catalog. The spectra from 0.5 keV up to 66 keV obtained with \textit{Suzaku} and \textit{NuSTAR} of six objects are newly analyzed here by utilizing the X-ray clumpy torus model (XCLUMPY), whereas we refer to Ogawa et al. (2021) for the results of Centaurus A. We find that these radio galaxies have similar torus covering fractions compared with radio quiet AGNs at the same Eddington ratios ($-3 < \log \lambda_{\rm Edd} < -1$). This result implies that (1) the torus structure is not an important factor that determines the presence of jets and (2) AGN jets have physically little effect on the torus.

In radio astronomy, the challenge of reconstructing a sky map from time ordered data (TOD) is known as an inverse problem. Standard map-making techniques and gridding algorithms are commonly employed to address this problem, each offering its own benefits such as producing minimum-variance maps. However, these approaches also carry limitations such as computational inefficiency and numerical instability in map-making and the inability to remove beam effects in grid-based methods. To overcome these challenges, this study proposes a novel solution through the use of the conditional invertible neural network (cINN) for efficient sky map reconstruction. With the aid of forward modeling, where the simulated TODs are generated from a given sky model with a specific observation, the trained neural network can produce accurate reconstructed sky maps. Using the five-hundred-meter aperture spherical radio telescope (FAST) as an example, cINN demonstrates remarkable performance in map reconstruction from simulated TODs, achieving a mean squared error of $2.29\pm 2.14 \times 10^{-4}~\rm K^2$, a structural similarity index of $0.968\pm0.002$, and a peak signal-to-noise ratio of $26.13\pm5.22$ at the $1\sigma$ level. Furthermore, by sampling in the latent space of cINN, the reconstruction errors for each pixel can be accurately quantified.

YZ Ret is the first X-ray flash detected classical nova, and is also well observed in optical, X-ray, and gamma-ray. We propose a comprehensive model that explains the observational properties. The white dwarf mass is determined to be $\sim 1.33 ~M_\odot$ that reproduces multiwavelength light curves of YZ Ret, from optical, X-ray, and to gamma-ray. We show that a shock is naturally generated far outside the photosphere because winds collide with themselves. The derived lifetime of shock explains some of the temporal variations of emission lines. The shocked shell significantly contributes to the optical flux in the nebular phase. The decline trend of shell emission in the nebular phase is close to $\propto t^{-1.75}$ and the same as the universal decline law of classical novae, where $t$ is the time from the outburst.

Aims. We tested the new atomic model using atmospheric models of stars of different spectral types: the Sun (dG2), HD22049 (dK2, Epsilon Eridani), GJ 832 (dM2), and GJ 581 (dM3). Methods. We used new Breit-Pauli distorted-wave (DW) multiconfiguration calculations, which proved to be relevant for many transitions in the mid-infrared (MIR) range. The new atomic model of Mg I includes the following: i) recomputed ECS data through the DW method, including the superlevels. ii) For the nonlocal thermodynamic equilibrium (NLTE) population calculations, 5676 theoretical transitions were added (3001 term-to-term). iii) All of these improvements were studied in the Sun and the stars listed above. Results. The Mg distribution between ionization states for stars with different effective temperatures was compared. For the Sun and Epsilon Eridani, Mg II predominates with more than 95 %, while for GJ 832 and GJ 581, Mg I represents more than 72 % of the population. Moreover, in the latter stars, the amount of Mg forming molecules in their atmosphere is at least two orders of magnitude higher. Regarding the NLTE population, a noticeable lower variability in the departure coefficients was found, indicating a better population coupling for the new model. Comparing the synthetic spectrum calculated with the older and new Mg I atomic model, these results show minimal differences in the visible range but they are stronger in the IR for all of the stars. This aspect should be considered when using lines from this region as indicators. Nevertheless, some changes in the spectral type were found, also emphasizing the need to test the atomic models in different atmospheric conditions. The most noticeable changes occurred in the FUV and NUV, obtaining a higher flux for the new atomic model regardless of the spectral type.

In this work we terminate inflation during a phase of Constant Roll by means of a waterfall field coupled to the inflaton and a spectator field. The presence of a spectator field means that inflation does not end at a single point, $\phi_e$, but instead has some uncertainty resulting in a stochastic end of inflation. We find that even modestly coupled spectator fields can drastically increase the abundance of Primordial Black Holes (PBHs) formed by many orders of magnitude. The power spectrum created by the inflaton can be as little as $10^{-4}$ during a phase of Ultra Slow-Roll and still form a cosmologically relevant number of PBHs. We conclude that the presence of spectator fields, which very generically will alter the end of inflation, is an effect that cannot be ignored in realistic models of PBH formation.

Advanced LIGO and Virgo's third observing run brought another binary neutron star merger (BNS) and the first neutron-star black-hole (NSBH) mergers. While no confirmed kilonovae (KNe) was identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant ($H_0$), the Galactic enrichment from $r$-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network (IGWN) and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically-driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of $0.43^{+0.58}_{-0.26}$ ($1.97^{+2.68}_{-1.2}$) for the Zwicky Transient Facility (Rubin Observatory).

Measurements of 21cm intensity mapping (IM) during the dark ages can potentially provide us with an unprecedented window on high redshifts and small scales. One of the main advantages this can bring involves the possibility to probe the nature of dark matter. Tests of dark matter models with the large-scale structure of the Universe are limited by non-linearities and astrophysical effects, which are not present for IM measurements during the dark ages. In this paper we focus on constraining the model in which dark matter is comprised, totally on in part, by ultra-light axion-like particles around the $10^{-18}-10^{-22}$ eV mass scale. For this model, the angular power spectrum of 21cm brightness temperature fluctuations will exhibit a small-scale suppression. However, this effect is intertwined with the imprint of baryon-dark matter relative velocity at recombination, causing at the same time an enhancement at large-scales, which is affected by the mass and abundance of axion dark matter. In this work we forecast how future radio arrays will be able to constrain ultra-light axion mass through both these effects on the angular power spectrum.

We present our findings on the spectral analysis of seven magnetic white dwarfs that were presumed to be double degenerates. We obtained time-resolved spectroscopy at the Gemini Observatory to look for evidence of binarity or fast rotation. We find three of our targets have rotation periods of less than an hour based on the shifting positions of the Zeeman-split H$\alpha$ components: 13, 35, and 39 min, and we find one more target with a ~hour long period that is currently unconstrained. We use offset dipole models to determine the inclination, magnetic field strength, and dipole offset of each target. The average surface field strengths of our fast rotators vary by 1-2 MG between different spectra. In all cases, the observed absorption features are too shallow compared to our models. This could be due to extra flux from a companion for our three low-mass targets, but the majority of our sample likely requires an inhomogeneous surface composition. Including an additional magnetic white dwarf with similar properties presented in the literature, we find that 5 of the 8 targets in this sample show field variations on minute/hour timescales. A crystallization driven dynamo can potentially explain the magnetic fields in three of our targets with masses above $0.7~M_{\odot}$ but another mechanism is still needed to explain their rapid rotation. We suggest that rapid rotation or low-masses point to binary evolution as the likely source of magnetism in 7 of these 8 targets.

Compressive strength is a key to understanding the internal structure of dust aggregates in protoplanetary disks and their resultant bodies, such as comets and asteroids in the Solar System. Previous work has modeled the compressive strength of highly-porous dust aggregates with volume filling factors lower than 0.1. However, a comprehensive understanding of the compressive strength from low ($<0.1$) to high ($>0.1$) volume filling factors is lacking. In this paper, we investigate the compressive strength of dust aggregates by using aggregate compression simulations resolving constituent grains based on JKR theory to formulate the compressive strength comprehensively. We perform a series of numerical simulations with moving periodic boundaries mimicking the compression behavior. As a result, we find that the compressive strength becomes sharply harder when the volume filling factor exceeds 0.1. We succeed in formulating the compressive strength comprehensively by taking into account the rolling motion of aggregates for low volume filling factors and the closest packing of aggregates for high volume filling factors. We also find that the dominant compression mechanisms for high volume filling factors are sliding and twisting motions, while rolling motion dominates for low volume filling factors. We confirm that our results are in good agreement with previous numerical studies. We suggest that our analytical formula is consistent with the previous experimental results if we assume the surface energy of silicate is $\simeq210\pm90\mathrm{\ mJ\ m^{-2}}$. Now, we can apply our results to properties of small compact bodies, such as comets, asteroids, and pebbles.

White dwarf stars are the most common stellar fossils. When in binaries, they make up the dominant form of compact object binary within the Galaxy and can offer insight into different aspects of binary formation and evolution. One of the most remarkable white dwarf binary systems identified to date is AR Scorpii (henceforth AR Sco). AR Sco is composed of an M-dwarf star and a rapidly-spinning white dwarf in a 3.56-hour orbit. It shows pulsed emission with a period of 1.97 minutes over a broad range of wavelengths, which led to it being known as a white dwarf pulsar. Both the pulse mechanism and the evolutionary origin of AR Sco provide challenges to theoretical models. Here we report the discovery of the first sibling of AR Sco, J191213.72-441045.1 (henceforth J1912-4410), which harbours a white dwarf in a 4.03-hour orbit with an M-dwarf and exhibits pulsed emission with a period of 5.30 minutes. This discovery establishes binary white dwarf pulsars as a class and provides support for proposed formation models for white dwarf pulsars.

We constrain the growth index $\gamma$ by performing a full-shape analysis of the power spectrum multipoles measured from the BOSS DR12 data. We adopt a theoretical model based on the Effective Field theory of the Large Scale Structure (EFTofLSS) and focus on two different cosmologies: $\gamma$CDM and $\gamma \nu$CDM, where we also vary the total neutrino mass. We explore different choices for the priors on the primordial amplitude $A_s$ and spectral index $n_s$, finding that informative priors are necessary to alleviate degeneracies between the parameters and avoid strong projection effects in the posterior distributions. Our tightest constraints are obtained with 3$\sigma$ Planck priors on $A_s$ and $n_s$: we obtain $\gamma = 0.647 \pm 0.085$ for $\gamma$CDM and $\gamma = 0.612^{+0.075}_{-0.090}$, $M_\nu < 0.30$ for $\gamma \nu$CDM at 68\% c.l., in both cases $\sim 1\sigma$ consistent with the $\Lambda$CDM prediction $\gamma \simeq 0.55$. Additionally, we produce forecasts for a Stage-IV spectroscopic galaxy survey, focusing on a DESI-like sample. We fit synthetic data-vectors for three different galaxy samples generated at three different redshift bins, both individually and jointly. Focusing on the constraining power of the Large Scale Structure alone, we find that forthcoming data can give an improvement of up to $\sim 85\%$ in the measurement of $\gamma$ with respect to the BOSS dataset when no CMB priors are imposed. On the other hand, we find the neutrino mass constraints to be only marginally better than the current ones, with future data able to put an upper limit of $M_\nu < 0.27~{\rm eV}$. This result can be improved with the inclusion of Planck priors on the primordial parameters, which yield $M_\nu < 0.18~{\rm eV}$.

Radio and near-infrared observations have observed dozens of protoplanetary disks that host spiral arm features. Numerical simulations have shown that companions may excite spiral density waves in protoplanetary disks via companion-disk interaction. However, the lack of direct observational evidence for spiral-driving companions poses challenges to current theories of companion-disk interaction. Here we report multi-epoch observations of the binary system HD 100453 with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) facility at the Very Large Telescope. By recovering the spiral features via robustly removing starlight contamination, we measure spiral motion across 4 yr to perform dynamical motion analyses. The spiral pattern motion is consistent with the orbital motion of the eccentric companion. With this first observational evidence of a companion driving a spiral arm among protoplanetary disks, we directly and dynamically confirm the long-standing theory on the origin of spiral features in protoplanetary disks. With the pattern motion of companion-driven spirals being independent of companion mass, here we establish a feasible way of searching for hidden spiral-arm-driving planets that are beyond the detection of existing ground-based high-contrast imagers.

We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN2023ixf, obtained with the Submillimeter Array (SMA) at 2.6-18.6 days after explosion. The observations were obtained as part the SMA Large Program POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN2023ixf, with the deepest limits of $L_\nu(230\,{\rm GHz})\lesssim 8.6\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 2.7 and 7.7 days, and $L_\nu(230\,{\rm GHz})\lesssim 3.4\times 10^{25}$ erg s$^{-1}$ Hz$^{-1}$ at 18.6 days. These limits are about a factor of 2 times dimmer than the mm emission from SN2011dh (IIb), about an order of magnitude dimmer compared to SN1993J (IIb) and SN2018ivc (IIL), and about 30 times dimmer than the most luminous non-relativistic SNe in the mm-band (Type IIb/Ib/Ic). Using these limits in the context of analytical models that include synchrotron self-absorption and free-free absorption we place constraints on the proximate circumstellar medium around the progenitor star, to a scale of $\sim 2\times 10^{15}$ cm, excluding the range $\dot{M}\sim {\rm few}\times 10^{-6}-10^{-2}$ M$_\odot$ yr$^{-1}$ (for a wind velocity, $v_w=115$ km s$^{-1}$, and ejecta velocity, $v_{\rm eje}\sim (1-2)\times 10^4$ km s$^{-1}$). These results are consistent with an inference of the mass loss rate based on optical spectroscopy ($\sim 2\times 10^{-2}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$), but are in tension with the inference from hard X-rays ($\sim 7\times 10^{-4}$ M$_\odot$ yr$^{-1}$ for $v_w=115$ km s$^{-1}$). This tension may be alleviated by a non-homogeneous and confined CSM, consistent with results from high-resolution optical spectroscopy.

In the Milky Way, the distribution of stars in the $[\alpha/\mathrm{Fe}]$ vs. $[\mathrm{Fe/H}]$ and $[\mathrm{Fe/H}]$ vs. age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disk. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the APOGEE survey. We find that the low- and high-$\alpha$ disk stars have a clean separation in copula space and use this to provide an automated separation of the $\alpha$ sequences using a purely statistical approach. This separation reveals that the high-$\alpha$ disk ends at the same [$\alpha$/Fe] and age at high $[\mathrm{Fe/H}]$ as the low-$[\mathrm{Fe/H}]$ start of the low-$\alpha$ disk, thus supporting a sequential formation scenario for the high- and low-$\alpha$ disks. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age $\tau$ and metallicity $[\mathrm{Fe/H}]$ conditional on Galactocentric radius $R$ and height $z$ in the range $0 < R < 20$ kpc and $|z| < 2$ kpc. The resulting trends in the age-metallicity correlation with radius, height, and [$\alpha$/Fe] demonstrate a $\approx 0$ correlation wherever kinematically-cold orbits dominate, while the naively-expected negative correlation is present where kinematically-hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age-metallicity structure of the low-$\alpha$ disk.

In the presence of self-interactions, the post-inflationary evolution of the inflaton field is driven into the non-linear regime by the resonant growth of its fluctuations. The once spatially homogeneous coherent inflaton is converted into a collection of inflaton particles with non-vanishing momentum. Fragmentation significantly alters the energy transfer rate to the inflaton's offspring during the reheating epoch. In this work we introduce a formalism to quantify the effect of fragmentation on particle production rates, and determine the evolution of the inflaton and radiation energy densities, including the corresponding reheating temperatures. For an inflaton potential with a quartic minimum, we find that the efficiency of reheating is drastically diminished after backreaction, yet it can lead to temperatures above the big bang nucleosynthesis limit for sufficiently large couplings. In addition, we use a lattice simulation to estimate the spectrum of induced gravitational waves, sourced by the scalar inhomogeneities, and discuss detectability prospects. We find that a Boltzmann approach allows to accurately predict some of the main features of this spectrum.

Observational evidence suggests the existence of dark matter (DM), which comprises approximately $84.4\%$ of matter in the universe. Recent advances in tabletop quantum sensor technology have enabled searches for nongravitational interactions of DM. Our experiment named ChangE utilizes Coupled Hot Atom eNsembles to search for liGht dark mattEr and new physics. We identify a strongly-coupled hybrid spin-resonance (HSR) regime that enhances the bandwidth of $^{21}$Ne nuclear spin by three orders of magnitude while maintaining high sensitivity. In combination with a self-compensating mode (SC) for low frequencies, we present a comprehensive broadband search for axion-like dark matter with Compton frequencies in the range of $[0.01, 1000]$ Hz. We set new constraints on the DM interactions with neutrons and protons, accounting for the stochastic effect. For the axion-neutron coupling, our results reach a low value of $|g_{ann}|\le 3\times 10^{-10}$ in the frequency range $[0.02 , 4]$ Hz surpassing astrophysical limits and provide the strongest laboratory constraints in the $[10, 100]$ Hz range. For the axion-proton coupling, we offer the best terrestrial constraints for the frequency below 100 Hz.

Modern large-scale physics experiments create datasets with sizes and streaming rates that can exceed those from industry leaders such as Google Cloud and Netflix. Fully processing these datasets requires both sufficient compute power and efficient workflows. Recent advances in Machine Learning (ML) and Artificial Intelligence (AI) can either improve or replace existing domain-specific algorithms to increase workflow efficiency. Not only can these algorithms improve the physics performance of current algorithms, but they can often be executed more quickly, especially when run on coprocessors such as GPUs or FPGAs. In the winter of 2023, MIT hosted the Accelerating Physics with ML at MIT workshop, which brought together researchers from gravitational-wave physics, multi-messenger astrophysics, and particle physics to discuss and share current efforts to integrate ML tools into their workflows. The following white paper highlights examples of algorithms and computing frameworks discussed during this workshop and summarizes the expected computing needs for the immediate future of the involved fields.

We study the effect of the isovector-scalar meson $a_0$(980) on the properties of nuclear matter and the neutron star (NS) matter by constructing a parity doublet model with including the $a_0$ meson based on the chiral SU(2)$_L\times$SU(2)$_R$ symmetry. We also include the $\omega$-$\rho$ mixing contribution to adjust the slope parameter at the saturation. We find that, when the chiral invariant mass of nucleon $m_0$ is smaller than about 800 MeV, the existence of $a_0$(980) enlarges the symmetry energy by strengthening the repulsive $\rho$ meson coupling. On the other hand, for large $m_0$ where the Yukawa coupling of $a_0$(980) to nucleon is small, the symmetry energy is reduced by the effect of $\omega$-$\rho$ mixing. We then construct the equation of state (EoS) of a neutron star matter to obtain the mass-radius relation of NS. We find that, in most choices of $m_0$, the existence of $a_0$(980) stiffens the EoS and makes the radius of NS larger. We then constrain the chiral invariant mass of nucleon from the observational data of NS, and find that $580 \,\text{ MeV} \lesssim m_0 \lesssim 860 \,\text{ MeV} $ for $L_0=57.7$ MeV.

Given an Equation of State (EOS) for neutron star (NS) matter, there is a unique mass-radius sequence characterized by a maximum mass $M_{\rm{NS}}^{\max}$ at radius $R_{\max}$. We first show analytically that the $M_{\rm{NS}}^{\max}$ and $R_{\max}$ scale linearly with two different combinations of NS central pressure $P_{\rm{c}}$ and energy density $\varepsilon_{\rm{c}}$ by dissecting perturbatively the dimensionless Tolman-Oppenheimer-Volkoff (TOV) equations governing NS internal variables. The scaling relations are then verified via 87 widely used and rather diverse phenomenological as well as 17 microscopic NS EOSs with/without considering hadron-quark phase transitions and hyperons by solving numerically the original TOV equations. The EOS of densest NS matter allowed before it collapses into a black hole (BH) is then obtained. Using the universal $M_{\rm{NS}}^{\max}$ and $R_{\max}$ scalings and NICER (Neutron Star Interior Composition Explorer) and XMM-Newton mass-radius observational data for PSR J0740+6620, a very narrow constraining band on the NS central EOS is extracted directly from the data for the first time without using any specific input EOS model.

A characteristic observational signature of cosmic strings are short duration gravitational wave (GW) bursts. These have been searched for by the LIGO-Virgo-KAGRA (LVK) collaboration, and will be searched for with LISA. We point out that these burst signals are repeated, since cosmic string loops evolve quasi-periodically in time, and will always appear from essentially the same position in the sky. We estimate the number of GW repeaters for LVK and LISA, and show that the string tension that can be probed scales as detector sensitivity to the sixth power, which raises hope for detection in future GW detectors. The observation of repeated GW bursts from the same cosmic string loop helps distinguish between the GW waveform parameters and the sky-localization.

According to General Relativity, an isolated black hole in vacuum shall be described by the Kerr metric, whose geodesic equations are integrable. The violation of integrability leads to chaos for particles moving around the black hole. This chaotic dynamics could leave imprints on the associated gravitational waveform and could be tested with upcoming observations. In this paper, we investigate the chaotic orbital dynamics induced by the violation of a certain spacetime symmetry, the circularity. Specifically, we focus on the resonant orbits of a particular noncircular spacetime as an example and find that they form chains of Birkhoff islands on Poincar\'e surfaces of section. We compare the island structures with those generated in typical nonintegrable but circular spacetimes. The islands of stability induced by noncircularity appear asymmetric on the most common Poincar\'e surface of section at the equatorial plane. The asymmetric patterns of islands vary discontinuously when the spacetime parameters transit through integrable regions. The origin of such features is explained in the context of perturbation analysis by considering the orbits associated with stable fixed points on the section. Possible observational implications about testing circularity through gravitational wave detection are discussed.

We couple the issue of evolution in the laws of physics with that of violations of energy conservation. We define evolution in terms of time variables canonically dual to ``constants'' (such as $\Lambda$, the Planck mass or the gravitational coupling), mimicking a procedure associated with one formulation of unimodular gravity. We then introduce variability via a dependence of {\it other} fundamental ``constants'' on these clocks. Although this is not needed, sharper results are obtained if this procedure violates local Lorentz invariance, which we define in the spirit of Horava-Lifshitz theories (modifying a $3+1$ split action, so that a Lorentz invariant 4D reassembly is no longer possible). We find that variability in the ``laws of physics'' generically leads to violations of energy conservation if either a matter parameter varies as a function of a gravitational clock, or a gravity parameter depends on a matter clock, with the other combinations sterile. We illustrate this with a variety of clocks (associated with matter, the speed of light, the Ricci scalar, etc) and parameters (mainly the gravitational and matter speed of light, but also the cosmological constant). We can accommodate in this construction (and improve) several early Varying Speed of Light solutions to the flatness and cosmological constant problem, allowing for variability effects related to the spatial curvature and $\Lambda$ to cause creation of radiation and a Hot Big Bang. But we can also go well beyond, for example modelling signature change by a change in the sign of $c^2$, thereby obtaining a {\it classical} realization of the Hartle-Hawking no-boundary proposal, among other developments.

We investigate the production of non-thermal fermionic dark matter particles during the reheating era following slow roll inflation, driven by axion-like pseudo-Nambu-Goldstone boson, $\varphi$, that is non-minimally coupled to the curvature scalar, ${\cal R}$. We consider two types of non-minimal coupling, one of which is $\xi \varphi^2 {\cal R}$ (referred to as NM-N), and the other one is $\alpha\left( 1+ \cos\left(\frac{\varphi}{f_a}\right)\right)$ (referred to as NMP-N), where $\alpha$ and $\xi$ are dimensionless parameters and $f_a$ is an energy scale. We determine benchmark values for both inflationary scenarios satisfying current bounds from Cosmic Microwave Background (CMB) radiation measurement and find the mass of inflaton to be $\sim {\cal O}\left(10^{12}\right) \text{GeV}$ for both inflationary scenarios and tensor-to-scalar ratio, $r\sim 0.0177$ (for NM-N) and $\sim 0.0097$ (for NMP-N) which fall inside $1-\sigma$ contour on scalar spectral index versus $r$ plane of Planck 2018+Bicep 3+Keck Array 2018 joint analysis, and can be probed by future CMB observations e.g. Simons Observatory. We then show that dark matter particles produced from the decay of inflaton can fully match the present-day cold dark matter (CDM) yield, as well as other cosmological constraints, if the coupling value between inflaton and dark matter, $y_\chi$, and the dark matter mass, $m_\chi$, are within the range $10^{-6}\gtrsim y_\chi\gtrsim 10^{-17}$ and ${\cal O}(1\text{GeV})\lesssim m_\chi\lesssim m_\phi/2$, where $m_\phi$ denotes the mass of inflaton. The exact range of $y_\chi$ and $m_\chi$ varies with different benchmark values as well as parameters of inflation, like energy scale of inflation and $r$, some of which are within the reach of next generation CMB experiments.

We present measurements of electron and hole trapping in three COSI germanium cross-strip detectors. By characterizing the relative charge collection efficiency (CCE) as a function of interaction depth, we show that intrinsic trapping of both electrons and holes have significant effects on the spectroscopic performance of the detectors. We find that both the electron and hole trapping vary from detector to detector, demonstrating the need for empirical trapping measurements and corrections. Using our measurements of charge trapping, we develop a continuous depth-dependent second-order energy correction procedure. We show that applying this empirical trapping correction produces significant improvements to spectral resolution and to the accuracy of the energy reconstruction.

Atom interferometry detectors like AION, ZAIGA, and AEDGE will be able to detect gravitational waves (GWs) at dHz covering the band between large space-based laser interferometers LISA/TianQin/Taiji and ground-based facilities LIGO/Virgo/KAGRA. They will detect the late inspiral and merger of GW sources containing intermediate-mass black holes (IMBHs) in the mass range $10^2-10^5\,{\rm M_\odot}$. We study how accurately the parameters of an IMBH binary can be measured using the noise curve of AION. Furthermore, we propose a detection scheme where the early inspiral of the binary is detected using the regular broadband mode while the merger is detected using the resonant mode. We find that by using such a detection scheme the signal-to-noise ratio (SNR) of the detection as well as the detection accuracy of the parameters can be enhanced compared to the full detection of the signal using the broadband mode. We, further, assess the impact of the necessary detection gap while switching from broadband to resonant mode studying the case of a short ($30\,{\rm s}$) and a long ($600\,{\rm s}$) gap. We find that the improvement in SNR and detection accuracy is bigger for the shorter gap but that even in the case of the long gap such a scheme can be beneficial.

We review the motivation for axions, discuss benchmark axion models, and report on the ongoing and planned axion experiments in Hamburg and their discovery potential.

The 12C(p,{\gamma})13N reaction is the onset process of both the CNO and Hot CNO cycles that drive massive star, Red and Asymptotic Giant Branch star and novae nucleosynthesis. The 12C(p,{\gamma})13N rate affects the final abundances of the stable 12,13C nuclides, with ramifications for meteoritic carbon isotopic abundances and the s-process neutron source strength. Here, a new underground measurement of the 12C(p,{\gamma})13N cross-section is reported. The present data, obtained at the Felsenkeller shallow-underground laboratory in Dresden (Germany), encompass the 320-620 keV center of mass energy range to include the wide and poorly constrained E = 422 keV resonance that dominates the rate at high temperatures. This work S-factor results, lower than literature by 25%, are included in a new comprehensive R-matrix fit, and the energy of the 1+ first excited state of 13N is found to be 2369.6(4) keV, with radiative and proton width of 0.49(3) eV and 34.9(2) keV respectively. A new reaction rate, based on present R-matrix fit and extrapolation, is suggested.

This work provides, at lower order, general analytical solutions for the orbital separation, merging time, and orbital frequency of binary systems emitting gravitational waves while being submitted to mass variations. Specific features, depending on the exponent of the mass derivative, are investigated in details. Two phenomenologically interesting cases are explicitly considered : i) binaries formed by two light primordial black holes submitted to Hawking evaporation and ii) bodies driven by a Bondi accretion of phantom dark energy. It is shown that three different regimes arise, including an intricate non-monotonic behaviour of the system. We study subtle imprints that could be associated with those phenomena. A careful analysis of the conditions of validity of the different hypotheses performed is finally carried out.

In this paper, we examine chaotic inflation within the context of the energy-momentum squared gravity (EMSG) focusing on the energy-momentum powered gravity (EMPG) that incorporates the functional $f(\mathbb{T}^2)\propto (\mathbb{T}^2)^{\beta}$ in the Einstein-Hilbert action, in which $\beta$ is a constant and $\mathbb{T}^2\equiv T_{\mu \nu}T^{\mu \nu}$ where $T_{\mu \nu}$ is the energy-momentum tensor, which we consider to represent a single scalar field with a power-law potential. We demonstrate that the presence of EMSG terms allows the single-field monomial chaotic inflationary models to fall within current observational constraints, which are otherwise disfavored by Planck and BICEP/Keck findings. We show that the use of a non-canonical Lagrangian with chaotic potential in EMSG can lead to significantly larger values of the non-Gaussianity parameter, $f_{\rm Nl}^{\rm equi}$ whereas EMSG framework with canonical Lagrangian gives rise to results similar to those of the standard single-field model.

The non-radial oscillations of the neutron stars (NSs) have been suggested as a useful tool to probe the composition of neutron star matter (NSM). With this scope in mind, we consider a large number of equations of states (EOSs) that are consistent with nuclear matter properties and pure neutron matter EOS based on a chiral effective field theory (chEFT) calculation for the low densities and perturbative QCD EOS at very high densities. This ensemble of EOSs is also consistent with astronomical observations, gravitational waves in GW170817, mass and radius measurements from Neutron star Interior Composition ExploreR (NICER). We analyze the robustness of known universal relations (URs) among the quadrupolar $f$ mode frequencies, masses and radii with such a large number of EOSs and we find a new UR that results from a strong correlation between the $f$ mode frequencies and the radii of NSs. Such a correlation is very useful in accurately determining the radius from a measurement of $f$ mode frequencies in the near future. We also show that the quadrupolar $f$ mode frequencies of NS of masses 2.0 M$_\odot$ and above lie in the range $\sim$ 2-3 kHz in this ensemble of physically realistic EOSs. A NS of mass 2M$_{\odot}$ with a low $f$ mode frequency may indicate the existence of non-nucleonic degrees of freedom.