Papers for Wednesday, Sep 14 2022

We present an analysis of six independent on-sky datasets taken with the Keck-II/NIRC2 instrument. Using the off-axis point spread function (PSF) reconstruction software AIROPA, we extract stellar astrometry, photometry, and other fitting metrics in order to characterize the performance of this package. We test the effectiveness of AIROPA to reconstruct the PSF across the field of view in varying atmospheric conditions, number and location of PSF reference stars, stellar crowding and telescope position angle (PA). We compare the astrometric precision and fitting residuals between a static PSF model and a spatially varying PSF model that incorporates instrumental aberrations and atmospheric turbulence during exposures. Most of the fitting residuals we measure show little to no improvement in the variable-PSF mode over the single-PSF mode. For one of the data sets, we find photometric performance is significantly improved (by ${\sim}10\times$) by measuring the trend seen in photometry as a function of off-axis location. For nearly all other metrics we find comparable astrometric and photometric precision across both PSF modes, with a ${\sim}13$% smaller astrometric uncertainty in variable-PSF mode in the best case. We largely confirm that the spatially variable PSF does not significantly improve the astrometric and other PSF fitting residuals over the static PSF for on-sky observations. We attribute this to unaccounted instrumental aberrations that are not characterized through afternoon adaptive optics (AO) bench calibrations.

We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that had gone undetected until recently, making it difficult to interpret these events within the current multi-wavelength picture of the flaring process. The May 6 event is the smallest stellar millimeter flare detected to date. We compare the relationship between the soft X-ray and millimeter emission to that observed in solar flares. The X-ray and optical flare energies of 10$^{30.3\pm0.2}$ and 10$^{28.9\pm0.1}$ erg, respectively, the coronal temperature of T=11.0$\pm$2.1 MK, and the emission measure of 9.5$\pm$2.2 X 10$^{49}$ cm$^{-3}$ are consistent with M-X class solar flares. We find the soft X-ray and millimeter emission during quiescence are consistent with the Gudel-Benz Relation, but not during the flare. The millimeter luminosity is >100X higher than that of an equivalent X1 solar flare and lasts only seconds instead of minutes as seen for solar flares.

Ultra-Diffuse Galaxies are both extreme products of galaxy evolution and extreme environments in which to test our understanding of star formation. In this work, we contrast the spatially resolved star formation activity of a sample of 22 HI-selected UDGs and 35 low-mass galaxies from the NASA Sloan Atlas (NSA) within 120 Mpc. We employ a new joint SED fitting method to compute star formation rate and stellar mass surface density maps that leverage the high spatial resolution optical imaging data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the UV coverage of GALEX, along with HI radial profiles estimated from a subset of galaxies that have spatially resolved HI maps. We find that the UDGs have low star formation efficiencies as a function of their atomic gas down to scales of 500 pc. We additionally find that the stellar mass-weighted sizes of our UDG sample are unremarkable when considered as a function of their HI mass -- their stellar sizes are comparable to the NSA dwarfs at fixed HI mass. This is a natural result in the picture where UDGs are forming stars normally, but at low efficiencies. We compare our results to predictions from contemporary models of galaxy formation, and find in particular that our observations are difficult to reproduce in models where UDGs undergo stellar expansion due to vigorous star formation feedback should bursty star formation be required down to $z=0$.

The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, $M_\star \sim 10^9\,{\rm M}_\odot$) may be the primary building blocks of the stellar halo of our Galaxy. We use cosmological simulations of the $\Lambda$ Cold Dark Matter model to investigate LMC-mass galaxies at $z=1-2$ using a semi-analytic model of galaxy formation. We find that LMC analogues at $z=2$ evolve until the present day along three distinct pathways: (1) those that are destroyed in Milky Way-mass hosts; (2) those that are themselves the main progenitors of Milky Way-mass galaxies; and (3) those that survive until $z=0$, with stellar mass $\sim$1.0 dex lower than typical Milky Ways. Our model predicts that the properties of these galaxies at $z=2$ (stellar metallicities, sizes, gas content etc.) are more or less indistinguishable, irrespective of which of these pathways is eventually taken; a survey targeting such galaxies in this redshift range would struggle to tell apart a 'destroyed' stellar halo progenitor from a 'surviving' LMC analogue. The only factor that determines the eventual fate of these galaxies is their proximity to a neighbouring Milky Way main progenitor at $z=2$: while the typical separation to a 'surviving' galaxy is around 7 Mpc, it is only 670 kpc to a 'destroyed' galaxy. Our results suggest that the old stellar populations in the Milky Way may be essentially indistinguishable from the progenitors of its stellar halo.

With the aim of bringing substantial insight to the fundamental question of how galaxies acquire their material for star-formation, we present in this work the first comprehensive characterisation of the galaxy connectivity (i.e. the number of filamentary streams connected to a galaxy) in relation with the cosmic environment, and a statistical exploration of the impact of connectivity on the star-formation rate at z=2. We detect kpc-scale filaments directly connected to galaxies by applying the DisPerSE filament finder to the dark matter (DM) density around 2977 massive central galaxies ($M_* > 10^{8}$ $\mathrm{M}_\odot / h$) of the TNG50-1 simulation. Our results demonstrate that more than half of the galaxies are connected to two or three streams. About $15\%$ of them are connected to zero streams, i.e. these galaxies are completely disconnected from the local web. We examine a variety of factors that could influence the connectivity finding out that this parameter increases with mass, decreases with local density for low mass galaxies, and does not depend on local environment, estimated by the Delaunay tessellation, for high mass galaxies. In order to better understand the relation between the small-scale streams and the large-scale environment, galaxies are further classified according to their location in different cosmic web environments. On average, we find that galaxies are connected to two streams in filaments, three in cluster cores, but are not connected to streams in cluster outskirts. Finally, we show that, at fixed local density and cosmic web environment, the specific star-formation (sSFR) of low mass galaxies is up to $7 \, \sigma$ enhanced due to connectivity. A milder impact is found for high mass galaxies, hinting at different relative efficiencies of matter inflow via filamentary streams in galaxies of different masses.

The increasing number of distant galaxies observed with ALMA by the ALPINE and REBELS surveys, and the early release observations of the James Webb Space Telescope (JWST) promise to revolutionize our understanding of cosmic star formation and the assembly of normal, dusty galaxies. Here we introduce a new suite of cosmological simulations performed with \texttt{dustyGadget} to interpret high-redshift data. We investigate the comoving star formation history, the stellar mass density and a number of galaxy scaling relations such as the galaxy main sequence (MS), the stellar-to-halo mass and dust-to-stellar mass relations at $z > 4$. The predicted star formation rate and total stellar mass density rapidly increase in time with a remarkable agreement with available observations, including recent JWST ERO and DD-ERS data at $z \geq 8$. A well defined galaxy MS is found already at $z < 10$ following a non evolving power-law, in agreement with both JWST and REBELS data at the low/high-mass end respectively, and consistent with a star formation efficiently sustained by gas accretion and a specific star formation rate increasing with redshift, as established by recent observations. A population of low-mass galaxies ($\rm{Log(M_\star/M_\odot)} < 9$) at $z \leq 7$ exceeding present estimates in the stellar mass function is also responsible of a large scatter in the stellar-to-halo and dust-to-stellar mass relations. Future JWST observations will provide invaluable constraints on these low-mass galaxies, helping to shed light on their role in cosmic evolution.

The outskirts of galaxies have been studied from multiple perspectives for the past few decades. However, it is still unknown if all galaxies have clear-cut edges like everyday objects. We address this question by developing physically motivated criteria to define the edges of galaxies. Based on the gas density threshold required for star formation, we define the edge of a galaxy as the outermost radial location associated with a significant drop in either past or ongoing in-situ star formation. We explore $\sim$1000 low-inclination galaxies with a wide range in morphology (dwarfs to ellipticals) and stellar mass ($10^7 M_{\odot} < M_{\star} < 10^{12}M_{\odot}$). The location of the edges of these galaxies ($R_{\rm edge}$) are visually identified as the outermost cut-off or truncation in their radial profiles using deep multi-band optical imaging from the IAC Stripe82 Legacy Project. We find this characteristic feature at the following mean stellar mass density which varies with galaxy morphology: $2.9\pm0.10\,M_{\odot}$/pc$^2$ for ellipticals, $1.1\pm0.04\,M_{\odot}/$pc$^2$ for spirals and $0.6\pm0.03\,M_{\odot}/$pc$^2$ for present-day star forming dwarfs. Additionally, we find that $R_{\rm edge}$ depends on its age (colour) where bluer galaxies have larger $R_{\rm edge}$ at a fixed stellar mass. The resulting stellar mass--size plane using $R_{\rm edge}$ as a physically motivated galaxy size measure has a very narrow intrinsic scatter ($\lesssim 0.06$ dex). These results highlight the importance of new deep imaging surveys to explore the growth of galaxies and trace the limits of star formation in their outskirts.

Disks of gas accreting onto supermassive black holes, powering active galactic nuclei (AGN), can capture stars from nuclear star clusters or form stars in situ via gravitational instability. The density and thermal conditions of these disks can result in rapid accretion onto embedded stars, dramatically altering their evolution in comparison to stars in the interstellar medium. Theoretical models predict that, when subjected to sufficiently rapid accretion, fresh gas replenishes hydrogen in the cores of these stars as quickly as it is burned into helium, reaching a quasi-steady state. Such massive, long-lived ("immortal") stars may be capable of dramatically enriching AGN disks with helium, and would increase the helium abundance in AGN broad-line regions relative to that in the corresponding narrow-line regions and hosts. We investigate how the helium abundance of AGN disks alters the evolution of stars embedded therein. We find, in agreement with analytical arguments, that stars at a given mass are more luminous at higher helium mass fractions, and so undergo more radiation-driven mass-loss. We further find that embedded stars tend to be less massive in disks with higher helium mass fractions, and that immortal stars are less common in such disks. Thus, disk composition can alter the rates of electromagnetic and gravitational wave transients as well as further chemical enrichment by embedded stars.

In addition to occupying the extreme, diffuse tail of the dwarf galaxy population, Ultra-Diffuse Galaxies (UDGs) are themselves a key laboratory in which to study star formation in extreme low-density environments. In the second paper of this series, we compare the spatially resolved star formation activity of 22 HI-selected UDGs and 21 "normal" dwarf galaxies within 120 Mpc to predictions within the pressure-regulated, feedback-modulated (PRFM) theory of star formation. To do so, we employ a joint SED fitting method that allows us to estimate star formation rate and stellar mass surface density from UV-optical imaging. We find that the PRFM framework extends successfully to the UDG regime - although the UDGs in our sample show unusually low star formation rate surface densities given their HI content, this low star formation efficiency can be naturally explained by the diffuse structure of the UDGs. In fact, when cast in the PRFM framework, the relationship between midplane pressure and star formation in the UDG sample is in good agreement not only with the "normal" dwarf reference sample, but also with measurements from more massive galaxies. Our results suggest that despite their low star formation efficiencies, the HI-rich UDGs need not be forming stars in an exotic manner. We also find that the UDGs are likely H$_2$-poor compared even to the overall dwarf population.

We determine the mass functions (MFs) and the dynamical parameters of 15 nearby open clusters (OCs) using the unprecedented data set of the Gaia Early Data Release 3. We select the members of each cluster by combining the photometric (colour and magnitude) and astrometric (parallax and proper motions) parameters of stars, minimizing the contamination from Galactic field interlopers. By comparing the observed distribution of stars along the cluster main sequence with the best-fitting synthetic population, we find the present-day MF and the binary fraction of the OCs, along with their dynamical parameters like mass, half-mass radius, and half-mass relaxation time. We found that the global present-day MF of OCs are consistent with a single power-law function, $F(m)\propto m^\alpha$, with slopes $-3<\alpha<-0.6$ including both subsolar, $0.2<m/\text{M}_\odot<1$, and supersolar mass regimes. A significant correlation between the MF-slope and the ratio of age to half-mass relaxation time is evidenced, similarly to the same conclusion already observed among Galactic globular clusters. However, OCs evolve along different tracks in comparison with the globular clusters, possibly indicating primordial differences in their initial mass function (IMF). The comparison with Monte Carlo simulations suggests that all the analysed OCs could have been born with an IMF with slope $\alpha_{\text{IMF}}<-2.3$. We also show that the less evolved OCs have a MF consistent with that of the solar neighbourhood, indicating a possible connection between the dissolution of OCs and the formation of the Galactic disc.

The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Useful insights can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing with observed XLFs. We focus on high-mass XRBs (HMXBs) and study the effects on the XLF of various, poorly-constrained assumptions regarding physical processes such as the common-envelope phase, the core-collapse, and wind-fed accretion. We use the new binary population synthesis code POSYDON and generate 96 synthetic XRB populations corresponding to different combinations of model assumptions. The generated XLFs are feature-rich, deviating from the commonly assumed single power law. We find a break in our synthetic XLF at luminosity $\sim 10^{38}$erg/s, similarly to observed XLFs. However, we find also a general overabundance of XRBs (up to a factor of $\sim$10 for certain model parameter combinations) driven primarily by XRBs with black hole accretors. Assumptions about the transient behavior of Be-XRBs, asymmetric supernova kicks, and common-envelope physics can significantly affect the shape and normalization of our synthetic XLFs. We find that less well-studied assumptions regarding the orbit circularization at the onset of Roche-lobe overflow and criteria for the formation of a wind-fed X-ray emitting accretion disk around black holes can also impact our synthetic XLFs and reduce the discrepancy with observations. Due to model uncertainties, our synthetic XLFs do not always agree well with observations. However, different combinations of model parameters leave distinct imprints on the shape of the synthetic XLFs and can reduce this deviation, revealing the importance of large-scale parameter studies and highlighting the power of XRBs in constraining binary evolution theory.

We study the kinematics of the interstellar medium (ISM) viewed "down the barrel" in 20 gravitationally lensed galaxies during Cosmic Noon ($z=1.5 - 3.5$). We use moderate-resolution spectra ($R\sim4000$) from Keck/ESI and Magellan/MagE to spectrally resolve the ISM absorption in these galaxies into $\sim$10 independent elements and use double Gaussian fits to quantify the velocity structure of the gas. We find that the bulk motion of gas in this galaxy sample is outflowing, with average velocity centroid $\left<v_{cent}\right>=-148 $km$\,$s$^{-1}$ ($\pm109 $km$\,$s$^{-1}$ scatter) measured with respect to the systemic redshift. 16 out of the 20 galaxies exhibit a clear positive skewness, with a blueshifted tail extending to $\sim -500$ km$\,$s$^{-1}$. The velocity width is considerably larger in the lensed galaxy sample compared to strong absorption systems viewed in quasar spectra which probe larger impact parameters, suggesting that absorbing gas seen in our sample is in close proximity to the host galaxies ($\lesssim 10$s of kpc). We examine scaling relations in outflow velocities with galaxy stellar mass and star formation rate (SFR), finding correlations consistent with a momentum-driven wind scenario. Our measured outflow velocities are also comparable to those reported for FIRE-2 and TNG50 cosmological simulations at similar redshift and galaxy properties. We also consider implications for interpreting results from lower-resolution spectra. We demonstrate that while velocity centroids are accurately recovered at lower resolution, the skewness, velocity width and probes of high velocity gas (e.g., $v_{95}$) are biased at $R\lesssim2000$. This work represents the largest available sample of well-resolved outflow velocity structure at $z>2$, and highlights the need for good spectral resolution to recover accurate properties.

We present an extensive study on the X-ray intraday variability of ten TeV-emitting high synchrotron peaked blazars (HBLs): 1ES 0229+200, 1ES 0414+009, PKS 0548-322, 1ES 1101-232, 1H 1219+301, H 1426+428, Mrk 501, 1ES 1959+650, PKS 2005-489, and 1ES 2344+514 made with twenty-five XMM-Newton pointed observations during its operational period. Intraday variability has been estimated in three energy bands: soft (0.3--2 keV), hard (2--10 keV) and total (0.3--10 keV). Although seven out of these ten TeV HBLs exhibited some intraday variability at three-sigma levels no major variations exceeding six percent were detected. We explored the spectral properties of the sample by extracting the hardness ratio from the soft and hard bands; no significant variations in the hardness ratio were observed in any source. We performed power spectral density analyses on the variable light-curves by fitting power-laws, yielding slopes lying in the range from 1.11 to 2.93 for different HBLs. We briefly discuss possible emission mechanisms and carry out rough estimates for magnetic fields, electron Lorentz factors and emission region sizes for seven of these HBLs.

The third Gaia data release includes a catalog of exoplanets and exoplanet candidates identified via the star's astrometric motion. This paper reports on tests for consistency between the Gaia two-body orbital solutions and precise Doppler velocities, for the stars currently amenable to such a comparison. For BD 17-0063, HD 81040, and HD 132406, the Gaia orbital solution and the Doppler data were found to be consistent, and were fitted jointly to obtain the best possible constraints on the planets' orbits and masses. Inconsistencies were found for 4 stars: HD 111232, probably due to additional planets that were not included in the astrometric model; HD 175167 and HR 810, possibly due to inaccurate treatment of non-Gaussian uncertainties in the Gaia orbital solutions; and HIP 66074, for unknown reasons. Consistency tests were also performed for HD 114762, which was reported in 1989 to have a brown dwarf or exoplanet but has since been shown to be binary star. The joint Gaia-Doppler analysis shows the secondary mass to be 0.215 +/- 0.013 Msun and the orbital inclination to be 3.63 +/- 0.06 degrees.

We characterise the properties and evolution of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in overlapping regions of the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT), covering the redshift range $0.20<z<0.80$. Using this sample, we measure no change in the ICL's stellar content (between 50-300\,kpc) over this redshift range in clusters with log${_{10}}(M_{\rm 200m,SZ}$/M$_{\odot})>$14.4. We also measure the stellar mass - halo mass (SMHM) relation for the BCG+ICL system and find that the slope, $\beta$, which characterises the dependence of $M_{\rm 200m,SZ}$ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth ($z<2$) occurs beyond the BCG's core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG 300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation's slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core. Moreover, for the lower halo mass clusters (log$_{10}(M_{\rm 200m,SZ}$/M$_{\odot})<$14.59 ), we detect a modest change in the colour gradient's slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts.

We present TESS photometry of the asynchronous polar BY Cam, which undergoes a beat-cycle between the 199.384-min white dwarf (WD) spin period and the 201.244-min orbital period. This results in changes in the flow of matter onto the WD. The TESS light curve covers 92% of the beat cycle once and 71% of the beat cycle twice. The strongest photometric signal, at 197.560-min, is ascribed to a side-band period. During times of light-curve stability, the photometry modulates at the spin frequency, supporting our WD spin-period identification. Both one-pole and two-pole accretion configurations repeat from one beat cycle to the next with clear and repeatable beat-phase dependent intensity variations. To explain these, we propose the operation of a magnetic valve at L1. The magnetic valve modulates the mass-transfer rate, as evidenced by a factor of 5 variation in orbital-averaged intensity, over the course of the beat cycle in a repeatable manner. The accretion stream threading distance from the WD is also modulated at the beat-period, because of the variation of the WD magnetic field with respect to the stream and because of changes in the mass transfer rate due to the operation of the magnetic valve. Changes in the threading distance result in significant shifts in the position of accreting spots around the beat cycle. As a consequence, only the faintest photometric minima allow for an accurate ephemeris determination. Three regions on the white dwarf appear to receive most of the accretion flow, suggestive of a complex WD magnetic field.

We explore the kinematic scaling relations of 38 dwarf galaxies in the Fornax Cluster using observations from the SAMI integral field spectrograph. We focus on the Fundamental Plane (FP), defined by the physical properties of the objects (scale length, surface brightness and velocity dispersion) and the Stellar Mass (Fundamental) Plane, where surface brightness is replaced by stellar mass, and investigate their dynamical-to-stellar-mass ratio. We confirm earlier results that the Fornax dEs are significantly offset above the FP defined by massive, hot stellar systems. For the Stellar Mass (Fundamental) Plane, which shows much lower scatter, we find that young and old dwarf galaxies lie at about the same distance from the plane, all with comparable scatter. We introduce the perpendicular deviation of dwarf galaxies from the Stellar Mass Plane defined by giant early-types as a robust estimate of their DM fraction, and find that the faintest dwarfs are systematically offset above the plane, implying that they have a higher dark matter fraction. This result is confirmed when estimating the dynamical mass of our dEs using a virial mass estimator, tracing the onset of dark matter domination in low mass stellar systems. We find that the position of our galaxies on the Stellar Mass FP agrees with the galaxies in the Local Group. This seems to imply that the processes determining the position of dwarf galaxies on the FP depend on the environment in the same way, whether the galaxy is situated in the Local Group or in the Fornax Cluster.

We test the usefulness of the intermediate ionisation lines AlIII 1860 and CIII] 1909 as reliable virial mass estimators for quasars. We identify a sample of 309 quasars from the SDSS DR16 in the redshift range 1.2 < z < 1.4 to have [OII] 3728 recorded on the same spectrum of AlIII 1860, SiIII 1890, and CIII] 1909. We set the systemic quasar redshift using careful measurements of [OII]. We then classified the sources as Population A, extreme Population A (xA) and Population B, and analysed the 1900\AA\ blend using multi-component models to look for systematic line shifts of the AlIII and CIII] along the quasar main sequence. We do not find significant shifts of the AlIII line peak in Pop. B and the wide majority of Pop. A. For Pop. xA, a small median blueshift of -250 km/s was observed, motivating a decomposition of the AlIII line profile into a virialized component centred at rest-frame and a blueshifted component for an outflow emission. For Pop. B objects, we proved the empirical necessity to fit a redshifted very broad component (VBC), clearly seen in CIII], and analysed the physical implications on a Pop. B composite spectrum using CLOUDY simulations. We find consistent black hole mass estimations using AlIII and CIII] as virial estimators for the bulk of Population A. AlIII (and even CIII]) is a reliable virial black hole mass estimator for Pop. A and B objects. xA sources deserve special attention due to the significant blueshifted excess observed in the line profile of AlIII, although not as large as those observed in CIV 1549.

We present the analysis of three faint clusters of the Small Magellanic Cloud RZ82, HW42 and RZ158. We employed the SOAR telescope instrument SAM with adaptive optics, allowing us to reach to V~23-24 mag, unprecedentedly, a depth sufficient to measure ages of up to about 10-12Gyr. All three clusters are resolved to their centres, and the resulting colour-magnitude diagrams (CMDs) allow us to derive ages of 3.9, 2.6, and 4.8Gyr respectively. These results are significantly younger than previous determinations (7.1, 5.0, and 8.3Gyr, respectively), based on integrated photometry or shallower CMDs. We rule out older ages for these clusters based on deep photometry and statistical isochrone fitting. We also estimate metallicities for the three clusters of [Fe/H]=-0.68, -0.57 and -0.90, respectively. These updated ages and metallicities are in good agreement with the age-metallicity relation for the bulk of SMC clusters. Total cluster masses ranging from ~7-11x10^3Mo were estimated from integrated flux, consistent with masses estimated for other SMC clusters of similar ages. These results reduce the number of SMC clusters known to be older than about 5 Gyr and highlight the need of deep and spatially resolved photometry to determine accurate ages for older, low-luminosity SMC star clusters.

We present deep high-resolution ($\sim$50 mas, 8 au) ALMA 0.88 and 1.3 mm continuum observations of the LkCa 15 disk. The emission morphology shows an inner cavity and three dust rings at both wavelengths, but with slightly narrower rings at the longer wavelength. Along a faint ring at 42 au, we identify two excess emission features at $\sim$10$\sigma$ significance at both wavelengths: one as an unresolved clump and the other as an extended arc, separated by roughly 120 degrees in azimuth. The clump is unlikely to be a circumplanetary disk (CPD) as the emission peak shifts between the two wavelengths even after accounting for orbital motion. Instead, the morphology of the 42 au ring strongly resembles the characteristic horseshoe orbit produced in planet--disk interaction models, where the clump and the arc trace dust accumulation around Lagrangian points $L_{4}$ and $L_{5}$, respectively. The shape of the 42 au ring, dust trapping in the outer adjacent ring, and the coincidence of the horseshoe ring location with a gap in near-IR scattered light, are all consistent with the scenario of planet sculpting, with the planet likely having a mass between those of Neptune and Saturn. We do not detect point-like emission associated with a CPD around the putative planet location ($0.''27$ in projected separation from the central star at a position angle of $\sim$60\degr), with upper limits of 70 and 33 $\mu$Jy at 0.88 and 1.3 mm, respectively, corresponding to dust mass upper limits of 0.02--0.03 $M_{\oplus}$.

The time-scales associated with various stages of the star formation process represent major unknowns in our understanding of galactic evolution, as well as of star and planet formation. This is the second paper in a series aiming to establish a multi-tracer time-line of star formation in the Large Magellanic Cloud (LMC), focusing on the lifecycle of molecular clouds. We use a statistical method to determine a molecular cloud lifetime in the LMC of $t_{\text{CO}}=11.8^{+2.7}_{-2.2}$ Myr. This short time-scale is similar to the cloud dynamical time, and suggests that molecular clouds in the LMC are largely decoupled from the effects of galactic dynamics and have lifetimes set by internal processes. This provides a clear contrast to atomic clouds in the LMC, of which the lifetimes are correlated with galactic dynamical time-scales. We additionally derive the time-scale for which molecular clouds and HII regions co-exist as $t_{\text{fb}}=1.2^{+0.3}_{-0.2}$ Myr, implying an average feedback front expansion velocity of 12 km s$^{-1}$, consistent with expansion velocities of HII regions in the LMC observed directly using optical spectroscopy. Taken together, these results imply that the molecular cloud lifecycle in the LMC proceeds rapidly and is regulated by internal dynamics and stellar feedback. We conclude by discussing our measurements in the context of previous work in the literature, which reported considerably longer lifetimes for molecular clouds in the LMC, and find that these previous findings resulted from a subjective choice in timeline calibration that is avoided by our statistical methodology.

We present a systematic analysis of 191 stripped-envelope supernovae (SE SNe), aimed to compute their $^{56}$Ni masses from the luminosity in their radioactive tails ($M_\mathrm{Ni}^\mathrm{tail}$) and/or in their maximum light, and the mean $^{56}$Ni and iron yields of SE SNe and core-collapse SNe. Our sample consists of SNe of types IIb, Ib, and Ic from the literature and from the Zwicky Transient Facility Bright Transient Survey. We use color curves to infer host galaxy reddenings and the representative $R_V$ value for each SN type. To calculate luminosities from optical photometry, we compute bolometric corrections using 49 SE SNe with optical and near-IR photometry. We find that the equation of Khatami & Kasen relating peak time and luminosity is not a reliable estimator of the $^{56}$Ni masses of SE SNe. Instead, we find a correlation between $M_\mathrm{Ni}^\mathrm{tail}$, peak time, peak luminosity, and decline rate, which allows measuring individual $^{56}$Ni masses to a precision of 14%. Applying this method to the whole sample, we find, for SNe IIb, Ib, and Ic, mean $^{56}$Ni masses of $0.066\pm0.006$, $0.082\pm0.009$, and $0.132\pm0.011$ M$_{\odot}$, respectively. After accounting for their relative rates, for SE SNe as a whole we compute mean $^{56}$Ni and iron yields of $0.090\pm0.005$ and $0.097\pm0.006$ M$_{\odot}$, respectively. Combining these results with the recent Type II SN mean $^{56}$Ni mass derived by Rodr\'iguez et al., core-collapse SNe, as a whole, have mean $^{56}$Ni and iron yields of $0.055\pm0.006$ and $0.058\pm0.007$ M$_{\odot}$, respectively. We highlight that the Arnett model, Arnett's rule, and hydrodynamical models typically overestimate the $^{56}$Ni masses of SE SNe by 75, 90 and 65%, respectively.

The unknown cause of the correlation between Type Ia supernova (SN Ia) Hubble residuals and their host-galaxy masses (the "mass step") may bias cosmological parameter measurements. To better understand the mass step, we develop a SALT3 light-curve model for SN cosmology that uses the host-galaxy masses of 296 low-redshift SNe Ia to derive a spectral-energy distribution$-$host-galaxy mass relationship. The resulting model has larger average Ca H&K and Si II equivalent widths in low-mass host galaxies, at 2.3$\sigma$ and 2.2$\sigma$ significance, indicating higher explosion energies per unit mass. The model has phase-dependent changes in SN Ia colors as a function of host mass, indicating intrinsic differences in mean broad-band light curves. Although the model provides a better fit to the SN data overall, it does not substantially reduce data$-$model residuals for a typical light curve in our sample nor does it significantly reduce Hubble residual dispersion. This is because we find that previous SALT models parameterized most host-galaxy dependencies with its first principal component, although they failed to model some significant spectral variations. Our new model is luminosity- and cosmology-independent, and applying it to data reduces the mass step by $0.021\pm0.002$ mag (uncertainty accounts for correlated data sets); these results indicate that $\sim$35% of the mass step can be attributed to luminosity-independent effects. This SALT model version could be trained using alternative host-galaxy properties and at different redshifts, and therefore will be a tool for understanding redshift-dependent correlations between SNe Ia and their host properties as well as their impact on cosmological parameter measurements.

Abridged. We aim to better characterise the evolution and fates, and determine updated nucleosynthetic yields of intermediate-mass stars between primordial and EMP metallicity (Z=1e-10, 1e-8, 1e-7, 1e-6 and 1e-5). We also probed uncertainties in the nucleosynthesis of the oldest intermediate-mass stars during the asymptotic giant branch (AGB) phase. We analysed the evolution of models from their main sequence, through the thermally pulsing AGB (TP-AGB), to the latest stages of their evolution, using the Monash-Mount Stromlo stellar evolution code MONSTAR. The results were post-processed with the code MONSOON, which allowed for the determination of the nucleosynthetic yields of 77 species up to 62Ni. As reported in former works, we identified proton ingestion episodes (PIEs) in our lowest-mass lowest-Z models. Models of Z=1e-10 and Z=1e-8 in a narrow initial mass range around 5 Msun experience the cessation of thermal pulses, and their final fates as type-I1/2 supernovae cannot be discarded. All the models of initial mass of about 6-7 Msun experience a corrosive second dredge-up and undergo significant metal enrichment in their envelopes. This allows them to develop a solar-like TP-AGB or TP-super-AGB, ultimately becoming white dwarfs. Except for those undergoing the cessation of thermal pulses, all of our models show the nucleosynthetic signatures of both efficient third dredge-up and hot-bottom burning, with the activation of the NeNa cycle and the MgAlSi chains. This leads to the creation of vast amounts of CNO, with typical [N/Fe] > 4), and the characteristic abundance signature [N/Fe] > [C/Fe] > [O/Fe]. Due to differences in input physics (mostly related to convection and convective boundaries), our nucleosynthetic yields present dramatic differences with respect to recent results existing in the literature for intermediate-mass models of similar metallicities.

The European Space Agency Rosetta/Philae mission to Comet 67P/Churyumov-Gerasimenko in 2014-2016 is the most complete and diverse investigation of a comet carried out thus far. Yet, many physical and chemical properties of the comet remain uncertain or unknown, and cometary activity is still not a well-understood phenomenon. We here attempt to place constraints on the nucleus abundances and sublimation front depths of H2O and CO2 ice, and to reconstruct how the nucleus evolved throughout the perihelion passage. We employ the thermophysical modelling code 'Numerical Icy Minor Body evolUtion Simulator', or NIMBUS, to search for conditions under which the observed H2O and CO2 production rates are simultaneously reproduced before and after perihelion. We find that the refractories to water-ice mass ratio of relatively pristine nucleus material is mu~1, that airfall material has mu~2, and that the molar abundance of CO2 relative H2O near 30 per cent. The dust mantle thickness is typically < 2 cm. The average CO2 sublimation front depths near aphelion were ~3.8 m and ~1.9 m on the northern and southern hemispheres, respectively, but varied substantially with time. We propose that airfall material is subjected to substantial fragmentation and pulverisation due to thermal fatigue during the aphelion passage. Sub-surface compaction of material due to CO2 activity near perihelion seems to have reduced the diffusivity in a measurable way.

With the upgrade from GRAVITY to GRAVITY$^+$ the instrument will evolve into an all-sky interferometer that can observe faint targets, such as high redshift AGN. Observing the faintest targets requires reducing the noise sources in GRAVITY as much as possible. The dominant noise source, especially in the blue part of the spectrum, is the backscattering of the metrology laser light onto the detector. To reduce this noise we introduce two new metrology modes. With a combination of small hardware changes and software adaptations, we can dim the metrology laser during the observation without losing the phase referencing. For single beam targets, we can even turn off the metrology laser for the maximum SNR on the detector. These changes lead to an SNR improvement of over a factor of two averaged over the whole spectrum and up to a factor of eight in the part of the spectrum currently dominated by laser noise.

A large fraction of Type Ia supernova (SN Ia) observations over the next decade will be in the near-infrared (NIR), at wavelengths beyond the reach of the current standard light-curve model for SN Ia cosmology, SALT3 ($\sim 2800$--8700$A$ central filter wavelength). To harness this new SN Ia sample and reduce future light-curve standardization systematic uncertainties, we train SALT3 at NIR wavelengths (SALT3-NIR) up to 2 $\mu$m with the open-source model-training software SALTShaker, which can easily accommodate future observations. Using simulated data we show that the training process constrains the NIR model to $\sim 2$--3% across the phase range ($-20$ to $50$ days). We find that Hubble residual (HR) scatter is smaller using the NIR alone or optical+NIR compared to optical alone, by up to $\sim 30$% depending on filter choice (95% confidence). There is significant correlation between NIR light-curve stretch measurements and luminosity, with stretch and color corrections often improving HR scatter by up to $\sim20%$. For SN Ia observations expected from the \textit{Roman Space Telescope}, SALT3-NIR increases the amount of usable data in the SALT framework by $\sim 20$% at redshift $z\lesssim0.4$ and by $\sim 50$% at $z\lesssim0.15$. The SALT3-NIR model is part of the open-source {\tt SNCosmo} and {\tt SNANA} SN Ia cosmology packages.

Cluster substructure and ram pressure stripping in individual galaxies are among the primary evidence for the ongoing growth of galaxy clusters as they accrete galaxies and groups from their surroundings. We present a multi-wavelength study of the center of the Hydra I galaxy cluster, including exquisite new MeerKAT HI and DECam Halpha imaging which reveal conclusive evidence for ram pressure stripping in NGC 3312, NGC 3314a and NGC 3314b through compressed HI contours, well-defined HI tails, and ongoing star formation in the stripped gas. In particular, we quantify the stripped material in NGC 3312 and NGC 3314a, which makes up between 8% and 35% of the gas still in the disk, is forming stars at ~0.5 M_Sun yr^-1, and extends ~30-60 kpc from the main disk. The estimated stellar mass in the tails is an order of magnitude less than the HI mass. A fourth "ring" galaxy at the same velocity does not show signs of ram pressure in HI. In addition, we use the HI and stellar morphologies, combined with a Beta model of the hot intracluster medium, to constrain the real distances of the galaxies to the cluster center, and we use the chance alignment of NGC 3314b behind NGC 3314a to break the degeneracy between whether the galaxies are in front or in back of the cluster. The drag seen in the HI tails supports our preferred scenario that NGC 3312 and NGC 3314a are moving towards us as part of a foreground substructure which has already passed its pericenter and is on "out fall" from the cluster. The high surviving HI content of the galaxies may suggest that the substructure/intragroup medium can protect them from the harshest effects of ram pressure, or that in fact the galaxies are on more tangential orbits.

We present detailed C, O, Na, Mg, Si, Ca, Ti, V, Fe, Zr, Ba, and Eu abundance measurements for 20 red giant branch (RGB) stars in the LMC star cluster NGC 1846 ([Fe/H] = -0.59). This cluster is 1.95 Gyr old and lies just below the supposed lower age limit (2 Gyr) for the presence of multiple populations in massive star clusters. Our measurements are based on high and low-resolution VLT/FLAMES spectra combined with photometric data from HST. Corrections for non-local thermodynamic equilibrium effects are also included for O, Na, Mg, Si, Ca, Fe and Ba. Our results show that there is no evidence for multiple populations in this cluster based on the lack of any intrinsic star-to-star spread in the abundances of Na and O: we place 95 \% confidence limits on the intrinsic dispersion for these elements of $\leq 0.07$ and $\leq 0.09$ dex, respectively. However, we do detect a significant spread in the carbon abundances, indicating varying evolutionary mixing occurring on the RGB that increases with luminosity. Overall, the general abundance patterns for NGC 1846 are similar to those seen in previous studies of intermediate-age LMC star clusters and field stars.

Arrays of low-temperature microcalorimeters provide a promising technology for X-ray astrophysics: the imaging spectrometer. A camera with at least several thousand pixels, each of which has an energy-resolving power ($E/\Delta E\urss{FWHM}$) of a few thousand across a broad energy range (200~eV to 10~keV or higher), would be a revolutionary instrument for the study of energetic astrophysical objects and phenomena. Signal readout is a critical enabling technology. Multiplexed readout, in which signals from multiple pixels are combined into a single amplifier channel, allows a kilo pixel-scale microcalorimeter array to meet the stringent requirements for power consumption, mass, volume, and cooling capacity in orbit. This chapter describes three different multiplexed-readout technologies for transition-edge-sensor microcalorimeters: time-division multiplexing, frequency-domain multiplexing, and microwave-SQUID multiplexing. For each multiplexing technique, we present the basic method, discuss some design considerations and parameters, and show the state of the art. The chapter concludes with a brief discussion of future prospects.

We present a study of the compact blister HII region BFS 10 and its highly filamentary molecular cloud. We utilize 12CO observations from the Five College Radio Astronomy Observatory to determine the distance, size, mass, and velocity structure of the molecular cloud. Infrared observations obtained from the UKIRT Infrared Deep Sky Survey and the Spitzer Infrared Array Camera, as well as radio continuum observations from the Canadian Galactic Plane Survey, are used to extract information about the central HII region. This includes properties such as the ionizing photon rate and infrared luminosity, as well as identifying a rich embedded star cluster associated with the central O9 V star. Time-scales regarding the expansion rate of the HII region and lifetime of the ionizing star reveal a high likelihood that BFS 10 will develop into a bipolar HII region. Although the region is expected to become bipolar, we conclude from the clouds velocity structure that there is no evidence to support the idea that star formation at the location of BFS 10 was triggered by two colliding clouds. A search for embedded young stellar objects (YSOs) within the molecular cloud was performed. Two distinct regions of YSOs were identified; one region associated with the rich embedded cluster and another sparse group associated with an intermediate mass YSO.

We report the discovery of a candidate binary system consisting of a black hole (BH) and a red giant branch star from the Gaia DR3. This binary system is discovered from 64096 binary solutions for which both astrometric and spectroscopic data are available. For this system, the astrometric and spectroscopic solutions are consistent with each other, making this system a confident candidate of a BH binary. The primary (visible) star in this system, Gaia DR3 5870569352746779008, is a red giant branch whose mass is quite uncertain. Fortunately, albeit the uncertainty of the primary's mass, the secondary (dark) object in this system has a mass of $>5.25$ $M_\odot$ with a probability of $99$ %, based on the orbital parameters. The mass of the secondary object is much larger than the maximum neutron star mass ($\sim 2.0$ $M_\odot$), which indicates that the secondary object is likely a BH. We argue that, if this dark object is not a BH, this system must be a more exotic system, in which the primary red giant branch star orbits around a triple star system (or a higher-order multiple star system) whose total mass is more than $5.25$ $M_\odot$. Future deep photometric observations are awaited to rule out such an exotic possibility and to determine whether or not this system is a genuine BH binary. If this is a genuine BH binary, this has the longest period ($1352.25 \pm 45.50$ days) among discovered so far.

Interferometric techniques such as aperture masking have the potential to enhance spatial resolution capabilities when imaging moderate-contrast sources with small angular size, such as close-in exoplanets and circumstellar disks around distant young stars. The Slicer Combined with an Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument, currently under development, is a lenslet integral field spectrograph that will enable the W. M. Keck Observatory to carry out high-contrast direct imaging of exoplanets between 2 and 5 microns. We explore the potential benefit of aperture masking to SCALES by testing the contrast achievable by several mask designs. The scalessim software package was used to simulate observations at wavelength bins in the M, L, and K bands, with optical path difference (OPD) maps used to simulate realistic Keck adaptive optics performance. Noise from astrophysical and instrumental sources was also applied to simulated signals. Mask designs were assessed based on depth of the generated contrast curves.

Massive binary black holes (MBBHs) in nearby galactic centers, if any, may be nano-Hertz gravitational wave (GW) sources for pulsar timing arrays (PTAs) to detect. Normally the objective GWs for PTA experiments are approximated as plane waves because its sources are presumably located faraway. For nearby GW sources, however, this approximation may be inaccurate due to the curved GW wave front and the GW strength changes along the paths of PTA pulsar pulses. In this paper, we analyze the near-field effect in the PTA detection of nearby sources and find it is important if the source distance is less than a few tens Mpc, and ignoring this effect may lead to a significant signal-to-noise underestimation especially when the source distance is comparable to the pulsar distances. As examples, we assume a nano-Hertz MBBH source located at either the Galactic Center (GC) or the Large Magellanic Cloud (LMC) according to the observational constraints/hints on the MBBH parameter space, and estimate its detectability by current/future PTAs. We find that the GC MBBH may be detectable by the Square Kilometer Array (SKA) PTA. It is challenging for detecting the LMC MBBH; however, if a number ($N\gtrsim10$) of stable millisecond pulsars can be found in the LMC center, the MBBH may be detectable via a PTA formed by these pulsars. We further illustrate the near-field effects on the PTA detection of an isotropic GW background contributed mainly by nearby GW sources, and the resulting angular correlation is similar to the Hellings-Downs curve.

We consider the possibility of primordial black hole, PBH, formation sourced by a rise in the power spectrum. The power spectrum becomes large at late times due to decay of the inflaton into vectors through a $\phi F \tilde{F}$ coupling. Two background inflaton models which are well supported by current Planck data are considered, natural inflation and hilltop inflation. Many of the papers considering formation of PBHs have considered a peaked power spectrum where $P_{\zeta}$ gets small again at late times. This avoids overproducing miniature PBHs which would evaporate and could violate BBN and CMB bounds. This paper examines the other way of avoiding these bounds, producing PBHs from perturbations formed closer to the end of inflation such that the PBHs are too small to violate these bounds. This has the advantage of allowing for simpler models in that no additional feature is needed to be added to evade constraints. Although these black holes would have evaporated, they can be close to without exceeding current BBN bounds, making it possible the signature will be revealed in the future. We calculate how the various model parameters affect the mass and number of PBHs produced. Any evidence for PBHs sourced from an inflationary power spectrum would provide evidence for inflation on a drastically different energy scale from the CMB, and thus would be highly valuable in answering what occurred during inflation.

Both observational and theoretical studies of black-hole activity or active galactic nucleus (AGN) feedback have been ongoing since the first indication of supermassive black holes powering quasar activity in the 1960s. Although several crucial astrophysical questions have been answered in the following decades, a number of open problems remain, in particular how AGN feedback operates over nearly eight orders of magnitude - from scales of $\sim 10^{-3}\,{\rm pc}$ to the galaxy-cluster scales of a few hundred kiloparsecs. At the beginning of June 2022, about 50 junior as well as senior researchers met in Brno for the post-lockdown edition of the Cologne-Prague-Brno meeting to try to connect the dots.

We analyze the cosmic-ray variations during a significant Forbush decrease observed with world-wide networks of ground-based neutron monitors and muon detectors during November 3-5, 2021. Utilizing the difference between primary cosmic-ray rigidities monitored by neutron monitors and muon detectors, we deduce the rigidity spectra of the cosmic-ray density (or omnidirectional intensity) and the first- and second-order anisotropies separately, for each hour of data. A clear two-step decrease is seen in the cosmic-ray density with the first $\sim2\%$ decrease after the interplanetary shock arrival followed by the second $\sim5\%$ decrease inside the magnetic flux rope (MFR) at 15 GV. Most strikingly, a large bidirectional streaming along the magnetic field is observed in the MFR with a peak amplitude of $\sim5\%$ at 15 GV which is comparable to the total density decrease inside the MFR. The bidirectional streaming could be explained by adiabatic deceleration and/or focusing in the expanding MFR, which have stronger effects for pitch angles near 90$^\circ$, or by selective entry of GCRs along a leg of the MFR. The peak anisotropy and density depression in the flux rope both decrease with increasing rigidity. The spectra vary dynamically indicating that the temporal variations of density and anisotropy appear different in neutron monitor and muon detector data.

The activity of the Sun alternates between a solar minimum and a solar maximum, the former corresponding to a period of "quieter" status of the heliosphere. During solar minimum, it is in principle more straightforward to follow eruptive events and solar wind structures from their birth at the Sun throughout their interplanetary journey. In this paper, we report analysis of the origin, evolution, and heliospheric impact of a series of solar transient events that took place during the second half of August 2018, i.e. in the midst of the late declining phase of Solar Cycle 24. In particular, we focus on two successive coronal mass ejections (CMEs) and a following high-speed stream (HSS) on their way towards Earth and Mars. We find that the first CME impacted both planets, whilst the second caused a strong magnetic storm at Earth and went on to miss Mars, which nevertheless experienced space weather effects from the stream interacting region (SIR) preceding the HSS. Analysis of remote-sensing and in-situ data supported by heliospheric modelling suggests that CME--HSS interaction resulted in the second CME rotating and deflecting in interplanetary space, highlighting that accurately reproducing the ambient solar wind is crucial even during "simpler" solar minimum periods. Lastly, we discuss the upstream solar wind conditions and transient structures responsible for driving space weather effects at Earth and Mars.

We investigate the planetary transit detectability in the presence of stellar rotational activity from light curves for young M-dwarfs and estimate improvements of the detection at near-infrared (NIR) wavelengths. Making maps of the transit signal detection efficiency over the orbital period and planetary radius with light curves of members of four clusters, Hyades, Praesepe, Pleiades, and Upper Scorpius observed by the K2 mission, we evaluate the detectability for the rotation period and modulation semi-amplitude. We find that the detection efficiency remarkably decreases to about 20% for rapidly rotators with P_{rot} <= 1 d and the lack of planets in Pleiades is likely due to the high fraction of rapidly rotating M-dwarfs. We also evaluate the improvements of the planet detection with NIR photometry via tests using mock light curves assuming that the signal amplitude of stellar rotation decreases at NIR wavelengths. Our results suggest that NIR photometric monitoring would double relative detection efficiency for transiting planetary candidates with P_{rot} <= 1 d and find planets around M-dwarfs with approximately 100 Myr missing in the past transit surveys from the space.

Hierarchical mergers in bound gravitational environments can explain the presence of black holes with masses greater than $\sim 100 M_{\odot}$. Evidence for this process is found in the third LIGO-Virgo-KAGRA gravitational-wave transient catalog (GWTC-3). We study the efficiency with which hierarchical mergers can produce higher and higher masses using a simple model of forward evolution of binary black hole populations in gravitationally bound systems like stellar clusters. The model relies on pairing probability and initial mass functions for the black hole population, along with numerical relativity fitting formulas for the mass, spin and kick speed of the merger remnant. If unequal mass pairing is disfavored, we show that the retention probability decreases significantly with later generations of the binary black hole population. Our model also predicts the distribution of masses of each black hole merger generation. We find that two of the subdominant peaks in the GWTC-3 component mass spectrum are consistent with second and third generation mergers in dense environments. With more binary black hole detections, our model can be used to infer the black hole initial mass function and pairing probability exponent.

The overall understanding of cosmic dust particles is mainly inferred from the different Earth-based measurements of interplanetary dust particles and space missions such as Giotto, Stardust and Rosetta. The results from these measurements indicate the presence of a wide variety of morphologically significant dust particles. To interpret the light scattering and thermal emission observations arising due to dust in different regions of space, it is necessary to generate computer modelled realistic dust structures of various shape, size, porosity, bulk density, aspect ratio and material inhomogenity. The present work introduces a java package called Rough Ellipsoid Structure Tool (REST), which is a collection of multiple algorithms, that aims to craft realistic rough surface cosmic dust particles from spheres, super-ellipsoids and fractal aggregates depending on the measured bulk-density and porosity. Initially, spheres having $N_d$ dipoles or lattice points are crafted by selecting random material and space seed cells to generate strongly damaged structure, rough surface and poked structure. Similarly, REST generates rough surface super-ellipsoids and poked structure super-ellipsoids from initial super-ellipsoid structures. REST also generates rough fractal aggregates which are fractal aggregates having rough surface irregular grains. REST has been applied to create agglomerated debris, agglomerated debris super-ellipsoids and mixed morphology particles. Finally, the light scattering properties of the respective applied structures are studied to ensure their applicability. REST is a flexible structure tool, which shall be useful to generate various types of dust structures that can be applied to study the physical properties of dust in different regions of space.

Measurement of the universe expansion rate through the cosmic chronometers proves to be a novel approach to understanding cosmic history. Although it provides a direct determination of the Hubble parameters at different redshifts, it suffers from underlying systematic uncertainties. In this work, we analyze the recent cosmic chronometer data with and without systematic uncertainties and investigate how they affect the results. We perform our analysis in both model-dependent and independent methods to avoid any possible model bias. In the model-dependent approach, we consider the $\Lambda$CDM, wCDM and CPL models. On the Other hand, since the Gaussian process provides a unique tool to study data including a non-diagonal covariance matrix, our model-independent analysis is based on the Gaussian process.

Context: A challenge with radial-velocity (RV) data is disentangling the origin of signals either due to a planetary companion or to stellar activity. In fact, the existence of a planetary companion has been proposed, as well as contested, around the relatively bright, nearby M3.0V star AD Leo at the same period as the stellar rotation of 2.23d. Aims: We further investigate the nature of this signal. We introduce new CARMENES optical and near-IR RV data and an analysis in combination with archival data taken by HIRES and HARPS, along with more recent data from HARPS-N, GIANO-B, and HPF. Also, we address the confusion concerning the binarity of AD Leo. Methods: We consider possible correlations between the RVs and various stellar activity indicators accessible with CARMENES. We applied models within a Bayesian framework to determine whether a Keplerian model, a red-noise quasi-periodic model using a Gaussian process, or a mixed model would explain the observed data best. We also exclusively focus on spectral lines potentially associated with stellar activity. Results: The CARMENES RV data agree with the previously reported periodicity of 2.23d, correlate with some activity indicators, and exhibit chromaticity. However, when considering the entire RV data set, we find that a mixed model composed of a stable and a variable component performs best. Moreover, when recomputing the RVs using only spectral lines insensitive to activity, there appears to be some residual power at the period of interest. We therefore conclude that it is not possible to determinedly prove that there is no planet orbiting in synchronization with the stellar rotation given our data, current tools, machinery, and knowledge of how stellar activity affects RVs. We do rule out planets more massive than 27M_E (=0.084M_J). We also exclude any binary companion around AD Leo with Msini > 3-6M_J on orbital periods <14yr.

The elasticity of neutron star crust is important for adequate interpretation of observations. To describe elastic properties one should rely on theoretical models. The most widely used is Coulomb crystal model (system of point-like charges on neutralizing uniform background), in some works it is corrected for electron screening. These models neglect finite size of nuclei. This approximation is well justified except for the innermost crustal layers, where nuclei size becomes comparable with the inter-nuclear spacing. Still, even in those dense layers it seems reasonable to apply the Coulomb crystal result, if one assumes that nuclei are spherically symmetric: Coulomb interaction between them should be the same as interaction between point-like charges. This argument is indeed correct, however, as we point here, shear of crustal lattice generates (microscopic) quadrupole electrostatic potential in a vicinity of lattice cites, which induces deformation on the nuclei. We analyze this problem analytically within compressible liquid drop model, using ionic spheroid model (which is generalization of well known ion sphere model). In particular, for ground state crust composition the effective shear modulus is reduced for a factor of $1-u^{5/3}/(2+3\,u-4\,u^{1/3})$, where u is the filling factor (ratio of the nuclei volume to the volume of the cell). This result is universal and does not depend on the applied nucleon interaction model. For the innermost layers of inner crust u~0.2 leading to reduction of the shear modulus by ~25%, which can be important for correct interpretation of quasi-periodic oscillations in the tails of magnetar flares.

Context: The TESS team periodically issues a new list of transiting exoplanet candidates based on the analysis of the accumulating light curves obtained by the satellite. The list includes the estimated epochs, periods, and durations of the potential transits. As the point spread function (PSF) of TESS is relatively wide, follow-up photometric observations at higher spatial resolution are required in order to exclude apparent transits that are actually blended background eclipsing binaries (BEBs). Aims: The Gaia space mission, with its growing database of epoch photometry and high angular resolution, enables the production of distinct light curves for all sources included in the TESS PSF, up to the limiting magnitude of Gaia. This paper reports the results of an ongoing Gaia-TESS collaboration that uses the Gaia photometry to facilitate the identification of BEB candidates and even to confirm on-target candidates in some cases. Methods: We inspected the Gaia photometry of the individual sources included in the TESS PSF, searching for periodic dimming events compatible with their ephemerides and uncertainties, as published by TESS. The performance of the search depends mainly on the number of Gaia measurements during transit and their precision. Results: Since February 2021, the collaboration has been able to confirm 126 on-target candidates and exclude 124 as BEBs. Since June 2021, when our search methodology matured, we have been able to identify on the order of 5% as on-target candidates and another 5% as BEBs. Conclusions: This synergistic approach is combining the complementary capabilities of two of the astronomical space missions of NASA and ESA. It serves to optimize the process of detecting new planets by making better use of the resources of the astronomical community.

Reconstructing the large scale density and velocity fields from surveys of galaxy distances, is a major challenge for cosmography. The data is very noisy and sparse. Estimated distances, and thereby peculiar velocities, are strongly affected by the Malmquist-like lognormal bias. Two algorithms have been recently introduced to perform reconstructions from such data: the Bias Gaussian correction coupled with the Wiener filter (BGc/WF) and the HAMLET implementation of the Hamiltonian Monte Carlo forward modelling. The two methods are tested here against mock catalogs that mimic the Cosmicflows-3 data. Specifically the reconstructed cosmography and moments of the velocity field (monopole, dipole) are examined. A comparison is made to the ``exact'' wiener filter as well - namely the Wiener Filter in the unrealistic case of zero observational errors. This is to understand the limits of the WF method. The following is found. In the nearby regime ($d \lesssim 40 {\rm Mpc}/h$) the two methods perform roughly equally well. HAMLET does slightly better in the intermediate regime ($ 40 \lesssim d \lesssim 120 {\rm Mpc}/h$). The main differences between the two appear in the most distant regime ($d \gtrsim 120 {\rm Mpc}/h$), close to the edge of the data. The HAMLET outperforms the BGc/WF in terms of better and tighter correlations, yet in the distant regime the HAMLET yields a somewhat biased reconstruction. Such biases are missing from the BGc/WF reconstruction. In sum, both methods perform well and create reliable reconstructions with significant differences apparent when details are examined.

We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences ($<3\sigma$) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a $\Lambda$CDM cosmology. We obtain $\Omega_{\rm m} = 0.30\pm 0.04$, $\sigma_8 = 0.81\pm 0.07 $ and $S_8 = 0.81\pm 0.04$, which implies a $\sim 0.4\sigma$ shift in the $\Omega_{\rm}-S_8$ plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses.

The imaging and characterization of a larger range of exoplanets, down to young Jupiters and exo-Earths will require accessing very high contrasts at small angular separations with an increased robustness to aberrations, three constraints that drive current instrumentation development. This goal relies on efficient coronagraphs set up on extremely large diameter telescopes such as the Thirty Meter Telescope (TMT), the Giant Magellan Telescope (GMT), or the Extremely Large Telescope (ELT). However, they tend to be subject to specific aberrations that drastically deteriorate the coronagraph performance: their primary mirror segmentation implies phasing errors or even missing segments, and the size of the telescope imposes large spiders, generating low-wind effect as already observed on the Very Large Telescope (VLT)/SPHERE instrument or at the Subaru telescope, or adaptive-optics-due petaling, studied in simulations in the ELT case. The ongoing development of coronagraphs has then to take into account their sensitivity to such errors. We propose an innovative method to generate coronagraphs robust to primary mirror phasing errors and low-wind and adaptive-optics-due petaling effect. This method is based on the apodization of the segment or petal instead of the entire pupil, this apodization being then repeated to mimic the pupil redundancy. We validate this so-called Redundant Apodized Pupil (RAP) method on a James Webb Space Telescope-like pupil composed of 18 hexagonal segments segments to align, and on the VLT architecture in the case of residual low-wind effect.

The detected exoplanet population displays a dearth of planets with sizes of about two Earth radii, the so-called radius gap. This is interpreted as an evolutionary effect driven by a variety of possible atmospheric mass loss processes of exoplanets. For mass loss driven by an exoplanet's irradiation by stellar X-ray and extreme-UV photons, the time evolution of the stellar magnetic activity is important. It is known from observations of open stellar clusters that stars of the same age and mass do not all follow the same time evolution of activity-induced X-ray and extreme-UV luminosities. Here we explore how a realistic spread of different stellar activity tracks influences the mass loss and radius evolution of a simulated population of small exoplanets and the observable properties of the radius gap. Our results show qualitatively that different saturation time scales, i.e. the young age at which stellar high-energy emission starts to decline, and different activity decay tracks over moderate stellar ages can cause changes in the population density of planets in the gap, as well as in the observable width of the gap. We also find that while the first 100 million years of mass loss are highly important to shape the radius gap, significant evolution of the gap properties is expected to take place for at least the first 500-600 million years, i.e. the age of the Hyades cluster. Observations of exoplanet populations with defined ages will be able to shed more light on the radius gap evolution.

Oncoming exoplanet spectro-imagers like the Planetary Camera and Spectrograph (PCS) for the Extremely Large Telescope (ELT) will aim for a new class of exoplanets, including Earth-like planets evolving around M dwarfs i.e., closer than 0.1'' with contrasts around 10^-8. This can be achieved with coronagraphs to modulate the wavefront. Classical coronagraphs are not optimal: 1) they impose a planetary photon loss, which is particularly problematic when the instrument includes a high spectral-resolution spectrograph, 2) some aberrations such as the missing segments of the ELT are dynamic and not compatible with a static coronagraph design, 3) the coupling of the exoplanet image with a fiber for spectroscopy only requires the electric field to be controlled on a small region of the detector. Such instruments would benefit from an adaptive tool to modulate the wavefront in both amplitude and phase. We propose to combine in the pupil plane a deformable mirror (DM) to control the phase and a digital micro-mirror device (DMD) i.e., an array made of 1920*1080 micro-mirrors able to switch between two positions, to control its amplitude. If the DM is already well-known in the field in particular for adaptive optics applications, the DMD has so far not been fully considered. At IPAG, we are currently assembling a testbed called CIDRE (Coronagraphy for DiRect Imaging of Exoplanets) to develop, test, calibrate, and validate the combination of these two components with a Lyot coronagraph. Since March 2022, CIDRE is assembled albeit without the Lyot coronagraph yet. The first few months have been dedicated to the calibration of the DMD. Since May 2022, it is operational and used to test dynamic amplitude apodization coronagraphs (so-called Shaped Pupils). This proceeding presents the set up of the CIDRE testbench and the first experimental results on adaptive Shaped Pupils obtained with the DMD.

We present analysis of the far ultraviolet (FUV) emission of sources in the central region of the Coma cluster (z=0.023) using the data taken by the UVIT aboard the multi-wavelength satellite mission AstroSat. We find a good correlation between the UVIT FUV flux and the fluxes in both wavebands of the Galex mission, for the common sources. We detect stars and galaxies, amongst which the brightest (r <= 17 mag) galaxies in the field of view are mostly members of the Coma cluster. We also detect three quasars (z = 0.38, 0.51, 2.31), one of which is likely the farthest object observed by the UVIT so far. In almost all the optical and UV colour-colour and colour-magnitude planes explored in this work, the Coma galaxies, other galaxies and bright stars could be separately identified, but the fainter stars and quasars often coincide with the faint galaxies. We have also investigated galaxies with unusual FUV morphology which are likely to be galaxies experiencing ram-pressure stripping in the cluster. Amongst others, two confirmed cluster members which were not investigated in the literature earlier, have been found to show unusual FUV emission. All the distorted sources are likely to have fallen into the cluster recently, and hence have not virialised yet. A subset of our data have optical spectroscopic information available from the archives. For these sources (~ 10% of the sample), we find that 17 galaxies identify as star-forming, 18 as composite and 13 as host galaxies for active galactic nuclei, respectively on the emission-line diagnostic diagram.

The Cherenkov Telescope Array Observatory (CTAO) will be the largest and most advanced ground-based facility for gamma-ray astronomy. Several dozens of telescopes will be operated at both the Northern and Southern Hemisphere. With the advent of multi-messenger astronomy, many new large science infrastructures will start science operations and target-of-opportunity observations will play an important role in the operation of the CTAO. The Array Control and Data Acquisition (ACADA) system deployed on each CTAO site will feature a dedicated sub-system to manage external and internal scientific alerts: the Transients Handler. It will receive, validate, and process science alerts in order to determine if target-of-opportunity observations can be triggered or need to be updated. Various tasks defined by proposal-based configurations are processed by the Transients Handler. These tasks include, among others, the evaluation of observability of targets and their correlation with known sources or objects. This contribution will discuss the concepts and design of the Transients Handler and its integration in the ACADA system.

We report a galaxy group candidate HPC1001 at $z\approx3.7$ in the COSMOS field. This structure was selected as a high galaxy overdensity at $z>3$ in the COSMOS2020 catalog. It contains ten candidate members, of which eight are assembled in a $10''\times10''$ area with the highest sky density among known protoclusters and groups at $z>3$. Four out of ten sources were also detected at 1.2$~$mm with Atacama Large Millimeter Array continuum observations. Photometric redshifts, measured by four independent methods, fall within a narrow range of $3.5<z<3.9$ and with a weighted average of $z=3.65\pm0.07$. The integrated far-IR-to-radio spectral energy distribution yields a total UV and IR star formation rate ${\rm SFR}\approx 900~M_{\odot}~yr^{-1}$. We also estimated a halo mass of $\sim10^{13}~M_\odot$ for the structure, which at this redshift is consistent with potential cold gas inflow. Remarkably, the most massive member has a specific star formation rate and dust to stellar mass ratio of $M_{\rm dust}/M_{*}$ that are both significantly lower than that of star-forming galaxies at this redshift, suggesting that HPC1001 could be a $z\approx3.7$ galaxy group in maturing phase. If confirmed, this would be the earliest structure in maturing phase to date, and an ideal laboratory to study the formation of the earliest quiescent galaxies as well as cold gas accretion in dense environments.

High-contrast imaging instruments need extreme wavefront control to directly image exoplanets. This requires highly sensitive wavefront sensors which optimally make use of the available photons to sense the wavefront. Here, we propose to numerically optimize Fourier-filtering wavefront sensors using automatic differentiation. First, we optimize the sensitivity of the wavefront sensor for different apertures and wavefront distributions. We find sensors that are more sensitive than currently used sensors and close to the theoretical limit, under the assumption of monochromatic light. Subsequently, we directly minimize the residual wavefront error by jointly optimizing the sensing and reconstruction. This is done by connecting differentiable models of the wavefront sensor and reconstructor and alternatingly improving them using a gradient-based optimizer. We also allow for nonlinearities in the wavefront reconstruction using Convolutional Neural Networks, which extends the design space of the wavefront sensor. Our results show that optimization can lead to wavefront sensors that have improved performance over currently used wavefront sensors. The proposed approach is flexible, and can in principle be used for any wavefront sensor architecture with free design parameters.

Ice naturally forms in the disordered or ``amorphous'' state when accreted from vapor at temperatures and pressures found in the interstellar medium and in the frigid, low density outer regions of the Sun's protoplanetary disk. It is therefore the expected form of ice in comets and other primitive bodies that have escaped substantial heating since formation. Despite expectations, however, the observational evidence for amorphous ice in comets remains largely indirect. This is both because the spectral features of amorphous ice are subtle and because the solar system objects for which we possess high quality data are mostly too close to the Sun and too hot for amorphous ice to survive near the surface, where it can be detected. This chapter reviews the properties of amorphous ice, the evidence for its existence and its consequences for the behavior of comets.

The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chemical mechanisms of HD 80606 b. We first run three-dimensional general circulation models (GCMs) to establish the atmospheric thermal and dynamical structures for different atmospheric metallicities and internal heat. Based on the GCM output, we then adopted a 1D time-dependent photochemical model to investigate the compositional variation along the eccentric orbit. The transition of the circulation patterns of HD 80606 b matched the dynamics regimes in previous works. Our photochemical models show that efficient vertical mixing leads to deep quench levels of the major carbon and nitrogen species and the quenching behavior does not change throughout the eccentric orbit. Instead, photolysis is the main driver of the time-dependent chemistry. A transient state of [CO]/[CH$_4$] $>$ 1 after periastron is confirmed for all metallicity and internal heat cases. The upcoming JWST Cycle 1 GO program will be able to track this real-time CH$_4$--CO conversion and infer the chemical timescale. Furthermore, sulfur species initiated by sudden heating and photochemical forcing exhibit both short-term and long-term cycles, opening an interesting avenue for detecting sulfur on exoplanets.

We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above $10^{19}$ eV. Photons in the zenith angle range from 30$^\circ$ to 60$^\circ$ can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95\% CL are set to an $E^{-2}$ diffuse flux of ultra-high energy photons above $10^{19}$ eV, $2{\times}10^{19}$ eV and $4{\times}10^{19}$ eV amounting to $2.11{\times}10^{-3}$, $3.12{\times}10^{-4}$ and $1.72{\times}10^{-4}$ km$^{-2}$ sr$^{-1}$ yr$^{-1}$, respectively. While the sensitivity of the present search around $2 \times 10^{19}$ eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated.

Molecular line ratios, such as HCN(1-0)/HCO$^+$(1-0) and HCN(4-3)/CS(7-6) are routinely used to identify AGN activity in galaxies. Such ratios are however hard to interpret as they are highly dependent on the physics and energetics of the gas and hence can seldom be used a as a unique unambiguous diagnostic. We use the composite galaxy NGC 1068 as a "laboratory", to investigate whether molecular line ratios between HCN, HCO$^+$ and CS are useful tracers of AGN-dominated gas and determine the origin of the differences in such ratios across different types of gas. Such determination will allow a more rigorous use of such ratios. We first empirically examine the aforementioned ratios at different angular resolutions to quantify correlations. We then use LTE and non-LTE analyses coupled with Markov Chain Monte Carlo (MCMC) sampling in order to determine the origin of the underlying differences in ratios. We propose that at high spatial resolution (< 50 pc) the HCN(4-3)/CS(2-1) is a reliable tracer of AGN activity. Finally we find that the variations in ratios are not a consequence of different densities or temperature but of different fractional abundances yielding to the important result that it is essential to consider the what chemical processes are at play when drawing conclusions from radiative transfer calculations. Upon analysis at varying spatial scales previous proposed as well as a new molecular line ratio have been shown to have varying levels of consistency. We have also determined from investigation of radiative transfer modelling of our data that it is essential to consider the chemistry of the species when reaching conclusions from radiative transfer calculations.

Cygnus X-1, as the first discovered black hole binary, is a key source for understanding the mechanisms of state transitions, and the scenarios of accretion in extreme gravity fields. We present a spectral-timing analysis of observations taken with the Insight-HXMT mission, focusing on the spectral-state dependent timing properties in the broad energy range of 1--150 keV, thus extending previous RXTE-based studies to both lower and higher energies. Our main results are the following: a) We successfully use a simple empirical model to fit all spectra, confirming that the reflection component is stronger in the soft state than in the hard state; b) The evolution of the total fractional root mean square (rms) depends on the selected energy band and the spectral shape, which is a direct result of the evolution of the power spectral densities (PSDs); c) In the hard/intermediate state, we see clear short-term variability features and a positive correlation between central frequencies of the variability components and the soft photon index $\Gamma_1$, also at energies above 15 keV. The power spectrum is dominated by red noise in the soft state instead. These behaviors can be traced to at least 90 keV; d) The coherence and the phase-lag spectra show different behaviors dependent on different spectral shapes.

We have gathered near-infrared $zyJ$-band high resolution spectra of nearly 300 field red giant stars with known lithium abundances in order to survey their \species{He}{i} $\lambda$10830 absorption strengths. This transition is an indicator of chromospheric activity and/or mass loss in red giants. The majority of stars in our sample reside in the red clump or red horizontal branch based on their $V-J,M_V$ color-magnitude diagram and their Gaia \teff, \logg\ values. Most of our target stars are Li-poor in the sense of having normally low Li abundances, defined here as \eps{Li}~$<$~1.25. Over 90\% of these Li-poor stars have weak $\lambda$10830 features. But more than half of the 83 Li-rich stars (\eps{Li}~$>$~1.25) have strong $\lambda$10830 absorptions. These large $\lambda$10830 lines signal excess chromospheric activity in Li-rich stars; there is almost no indication of significant mass loss. The Li-rich giants also may have a higher binary fraction than do Li-poor stars, based on their astrometric data. It appears likely that both residence on the horizontal branch and present or past binary interaction play roles in the significant Li-He connection established in this survey.

Recent population studies have searched for a subpopulation of primordial black holes (PBHs) in the gravitational-wave (GW) events so far detected by LIGO/Virgo/KAGRA (LVK), in most cases adopting a phenomenological PBH mass distribution. When deriving such population from first principles in the standard scenario, however, the equation of state of the Universe at the time of PBH formation may strongly affect the PBH abundance and mass distribution, which ultimately depend on the power spectrum of cosmological perturbations. Here we improve on previous population studies on several aspects: (i) we adopt state-of-the-art PBH formation models describing the collapse of cosmological perturbations across the QCD epoch; (ii) we perform the first Bayesian multi-population inference on GW data including PBHs and directly using power spectrum parameters instead of phenomenological distributions; (iii) we critically confront the PBH scenario with LVK phenomenological models describing the GWTC-3 catalog both in the neutron-star and in the BH mass ranges, also considering PBHs as subpopulation of the total events. Our results confirm that LVK observations prevent the majority of the dark matter to be in the form of stellar mass PBHs. We find that the best fit PBH model can comprise a small fraction of the total events, in particular it can naturally explain events in the mass gaps. If the lower mass-gap event GW190814 is interpreted as a PBH binary, we predict that LVK should detect up to a few subsolar mergers and one to $\approx 30$ lower mass gap events during the upcoming O4 and O5 runs. Finally, mapping back the best-fit power spectrum into an ultra slow-roll inflationary scenario, we show that the latter predicts detectable PBH mergers in the LVK band, a stochastic GW background detectable by current and future instruments, and may include the entirety of dark matter in asteroid-mass PBHs.

Numerical simulations within a cold dark matter (DM) cosmology form halos whose density profiles have a steep inner slope (`cusp'), yet observations of galaxies often point towards a flat central `core'. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM halos. We train the network on simulated dwarf galaxies from the NIHAO and AURIGA projects, which include both DM cusps and cores: line-of-sight velocities and 2D spatial distributions of their stars are used as inputs to obtain a PDF representing the probability of predicting a specific inner slope. The model recovers accurately the expected DM profiles: $\sim$82$\%$ of the galaxies have a derived inner slope within $\pm$0.1 of their true value, while $\sim$98$\%$ within $\pm$0.3. We apply our model to four Local Group dwarf spheroidal galaxies and find results consistent with those obtained with the Jeans modelling based code GravSphere: the Fornax dSph has a strong indication of possessing a central DM core, Carina and Sextans have cusps (although the latter with large uncertainties), while Sculptor shows a double peaked PDF indicating that a cusp is preferred, but a core can not be ruled out. Our results show that simulation-based inference with neural networks provide a innovative and complementary method for the determination of the inner matter density profiles in galaxies, which in turn can help constrain the properties of the elusive DM.

Three- and five-minute oscillations are commonly found in any sunspot. As they are modulated by the internal thermal and magnetic structures of a sunspot, therehence, they could be used as an effective tool for sunspot seismology. In this paper, we investigate the properties of oscillations in sunspot groups with varying size and magnetic field, and aim to establish the relationships between sunspot oscillations and its internal structure comparatively. We selected three groups of unipolar sunspot with approximately axial-symmetric magnetic field and calculated their Fourier spectra based on the Ultraviolet(UV)/Extreme ultraviolet(EUV) emission intensity variations recorded by the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). We found that the distribution of three minute oscillation is defined by the joint effect of diverging magnetic field and the stratification of sunspot atmosphere. Its distribution could be modified by any invading magnetic structures in the umbra. Whereas the five minute oscillations are more prominent in small spots, it implies that five minute oscillation is very closely connected with umbral dynamics.

A new model atom of Sc II was constructed using the most up-to-date atomic data. For the testing purpose, the non-local thermodynamic equilibrium (non-LTE) calculations were carried out for three stars with reliably determined atmospheric parameters: the Sun, HD~61421 (Procyon), and HD~84937. Accounting for deviations from LTE leads to smaller abundance errors compared with the LTE case and consistent within the error bars abundances obtained from different Sc II lines. Solar non-LTE abundance log eps_Sun = 3.12+-0.05 exceeds the meteoritic abundance recommended by Lodders (2021), by 0.08~dex. But agreement within 0.02~dex with the meteoritic abundance is obtained for Procyon. Using high-resolution spectra, we determined the scandium LTE and non-LTE abundances for 56 stars in the metallicity range -2.62 <= [Fe/H] <= 0.24. The dependence of [Sc/Fe] on [Fe/H] demonstrates a similarity with the behavior of the alpha-process elements: scandium is enhanced relative to iron ([Sc/Fe] ~ 0.2) for [Fe/H] < -1, and [Sc/Fe] decreases with increasing [Fe/H] for the higher metallicity. There is a hint of a tight relation between abundances of scandium and titanium. The results obtained provide observational constraints to the scenarios of scandium origin.

The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep ground-based resolved-star data sets from Subaru's Hyper Suprime-Cam. We present one new ultrafaint dwarf galaxy in the M81 group and identify five additional extremely low surface brightness candidate ultrafaint dwarfs that reach deep into the ultrafaint regime to $M_V \sim -6$ (similar to current limits for Andromeda satellites). These candidates' luminosities and sizes are similar to known Local Group dwarf galaxies Tucana B, Canes Venatici I, Hercules, and Bo\"otes I. Most of these candidates are likely to be real, based on tests of our techniques on blank fields. Intriguingly, all of these candidates are spatially clustered around NGC 3077, which is itself an M81 group satellite in an advanced state of tidal disruption. This is somewhat surprising, as M81 itself and its largest satellite M82 are both substantially more massive than NGC 3077 and by virtue of their greater masses, would have been expected to host as many or more ultrafaint candidates. These results lend considerable support to the idea that satellites of satellites are an important contribution to the growth of satellite populations around Milky Way-mass galaxies.

Baryon acoustic oscillations (BAOs) are a powerful probe of the expansion history of our Universe and are typically measured in the two-point statistics of a galaxy survey, either in Fourier space or in configuration space. In this work, we report a first measurement of BAOs from a joint fit of power spectrum and correlation function multipoles. We tested our new framework with a set of 1000 mock catalogs and showed that our method yields smaller biases on BAO parameters than individually fitting power spectra or correlation functions, or when combining them with the Gaussian approximation method. Our estimated uncertainties are slightly larger than those from the Gaussian approximation, likely due to noise in our sample covariance matrix, the larger number of nuisance parameters, or the fact that our new framework does not rely on the assumption of Gaussian likelihoods for the BAO parameters. However, we argue that our uncertainties are more reliable since they rely on fewer assumptions, and because our method takes correlations between Fourier and configuration space at the level of the two-point statistics. We performed a joint analysis of the luminous red galaxy sample of the extended baryon oscillation spectroscopic survey (eBOSS) data release 16, obtaining $D_H/r_d = 19.27 \pm 0.48$ and $D_M/r_d = 17.77 \pm 0.37$, in excellent agreement with the official eBOSS consensus BAO-only results $D_H/r_d = 19.33 \pm 0.53$ and $D_M/r_d =17.86 \pm 0.33$.

In this work we present the spectrum modeling results for the newly discovered luminous blue variable (LBV) in the NGC 1156 galaxy. Extended atmosphere models were calculated using the non-LTE code CMFGEN. We have obtained the luminosity of the discovered LBV $L\simeq (1.6\pm0.2) \times 10^{5} L_{\odot}$, effective temperature $T_{\text{eff}}=7.9\pm0.4$ kK and mass-loss rate $\dot{M}f^{-0.5}= (8.2\pm1.0) \times 10^{-4} M_{\odot}\text{yr}^{-1}$. The hydrogen abundance in the wind is $\approx20\,$% for the metallicity $Z=0.5 Z_{\odot}$ of the host galaxy.

The DArk Matter Particle Explorer (DAMPE) is a space-based particle detector launched on December 17th, 2015 from the Jiuquan Satellite Launch Center (China). The main goals of the DAMPE mission are the study of galactic cosmic rays (CR), the electron-positron energy spectrum, gamma-ray astronomy, and indirect dark matter search. Among its sub-detectors, the deep calorimeter makes DAMPE able to measure electrons and gamma-ray spectra up to 10 TeV, and CR nuclei spectra up to hundreds of TeV, with unprecedented energy resolution. This high-energy region is important in order to search for electron-positron sources, for dark matter signatures in space, and to clarify CR acceleration and propagation mechanisms inside our galaxy. A general overview of the DAMPE experiment will be presented in this work, along with its main results and ongoing activities.

BlueMUSE is an integral field spectrograph in an early development stage for the ESO VLT. For our design of the data reduction software for this instrument, we are first reviewing capabilities and issues of the pipeline of the existing MUSE instrument. MUSE has been in operation at the VLT since 2014 and led to discoveries published in more than 600 refereed scientific papers. While BlueMUSE and MUSE have many common properties we briefly point out a few key differences between both instruments. We outline a first version of the flowchart for the science reduction, and discuss the necessary changes due to the blue wavelength range covered by BlueMUSE. We also detail specific new features, for example, how the pipeline and subsequent analysis will benefit from improved handling of the data covariance, and a more integrated approach to the line-spread function, as well as improvements regarding the wavelength calibration which is of extra importance in the blue optical range. We finally discuss how simulations of BlueMUSE datacubes are being implemented and how they will be used to prepare the science of the instrument.

We investigate the time evolution of bias of cosmic density fields. We perform numerical simulations of the evolution of the cosmic web for the conventional $\Lambda$ cold dark matter ($\Lambda$CDM) model. The simulations cover a wide range of box sizes $L=256 - 1024\Mpc$, and epochs from very early moments $z=30$ to the present moment $z=0$. We calculate spatial correlation functions of galaxies, $\xi(r)$, using dark matter particles of the biased $\Lambda$CDM simulation. We analyse how these functions describe biasing properties of the evolving cosmic web. We find that for all cosmic epochs the bias parameter, defined through the ratio of correlation functions of selected samples and matter, depends on two factors: the fraction of matter in voids and in the clustered population, and the luminosity (mass) of galaxy samples. Gravity cannot evacuate voids completely, thus there is always some unclustered matter in voids, thus the bias parameter of galaxies is always greater than unity, over the whole range of evolution epochs. We find that for all cosmic epochs bias parameter values form regular sequences, depending on galaxy luminosity (particle density limit), and decreasing with time.

The emission line spectra of WR stars are often formed completely in the optically thick stellar wind. Hence, any assumption on the wind velocity law in a spectral analysis has a profound impact on the determination of the stellar parameters. By comparing Potsdam Wolf-Rayet (PoWR) model spectra calculated with different $\beta$ laws, we show that the velocity field heavily influences the spectra: by using the appropriate $\beta$ laws, the entire range of late and early types can be covered with the same stellar model.

The kinematics of molecular gas are crucial for setting the stage for star formation. One key question related to the kinematic properties of gas is how they depend on the spatial scale. We aim to describe the CO spectra, velocity dispersions, and especially the linewidth-size relation, of molecular gas from cloud (parsec-) scales to kiloparsec scales in a complete region within the Milky Way disk. We utilise the census of molecular clouds within 2 kpc from our earlier work, together with CO emission data for them from the literature. We study the kinematics and the Larson's relations for the sample of individual clouds. We also mimic a face-on view of the Milky Way and analyse the kinematics of the clouds within apertures of 0.25-2 kpc in size. In this way, we describe the scale-dependency of the CO gas kinematics and Larson's relations. We describe the spectra of CO gas at cloud scales and in apertures between 0.25-2 kpc in our survey area. The spectra within the apertures are relatively symmetric but show non-Gaussian high-velocity wings. At cloud-scales, our sample shows a linewidth-size relation \sigma_v=1.5*R^{0.3\pm0.1} with a large scatter. The mass-size relation in the sample of clouds is M_{CO}= 794*R^{1.5\pm0.5}. The relations are also present for the apertures at kpc-scales. A suggestive dependency on galactic environment is seen, with apertures closer to the Galactic centre and the Sagittarius spiral arm having slightly higher velocity dispersions. We explore the possible effect of a diffuse component in the survey area, and find that such a component would widen the CO spectra and could flatten the linewidth-size relation. Understanding the nature of the possible diffuse CO component and its effects on observations is crucial for connecting Galactic and extragalactic data.

The collision of a primordial black hole with a neutron star results in the black hole eventually consuming the entire neutron star. However, if the black hole is magnetically charged, and therefore stable against decay by Hawking radiation, the consequences can be quite different. Upon colliding with a neutron star, a magnetic black hole very rapidly comes to a stop. For large enough magnetic charge, we show that this collision can be detected as a sudden change in the rotation period of the neutron star, a glitch or anti-glitch.We argue that the magnetic primordial black hole, which then settles to the core of the neutron star, does not necessarily devour the entire neutron star; the system can instead reach a long-lived, quasi-stable equilibrium. Because the black hole is microscopic compared to the neutron star, most stellar properties remain unchanged compared to before the collision. However, the neutron star will heat up and its surface magnetic field could potentially change, both effects potentially observable.

The nuclear region of Type 1 AGNs has only been partially resolved so far in the near-infrared (IR) where we expect to see the dust sublimation region and the nucleus directly without obscuration. Here we present the near-IR interferometric observation of the brightest Type 1 AGN NGC4151 at long baselines of ~250 m using the CHARA Array, reaching structures at hundred micro-arcsecond scales. The squared visibilities decrease down to as low as ~0.25, definitely showing that the structure is resolved. Furthermore, combining with the previous visibility measurements at shorter baselines but at different position angles, we show that the structure is elongated *perpendicular* to the polar axis of the nucleus, as defined by optical polarization and a linear radio jet. A thin-ring fit gives a minor/major axis ratio of ~0.7 at a radius ~0.5 mas (~0.03 pc). This is consistent with the case where the sublimating dust grains are distributed preferentially in an equatorial plane in a ring-like geometry, viewed at an inclination angle of ~40 deg. Recent mid-IR interferometric finding of polar-elongated geometry at a pc scale, together with a larger-scale polar outflow as spectrally resolved by the HST, would generally suggest a dusty, conical and hollow outflow being launched presumably in the dust sublimation region. This might potentially lead to a polar-elongated morphology in the near-IR, as opposed to the results here. We discuss a possible scenario where an episodic, one-off anisotropic acceleration formed a polar-fast and equatorially-slow velocity distribution, having lead to an effectively flaring geometry as we observe.

Low internal friction coatings are key components of advanced technologies such as optical atomic clocks and high-finesse optical cavity and often lie at the forefront of the most advanced experiments in Physics. Notably, increasing the sensitivity of gravitational-wave detectors depends in a very large part on developing new coatings, which entails developing more suitable methods and models to investigate their loss angle. In fact, the most sensitive region of the detection band in such detectors is limited by the coating thermal noise, which is related to the loss angle of the coating. Until now, models which describe only ideal physical properties have been adopted, wondering about the use of one or more loss angles to describe the mechanical properties of coatings. Here we show the presence of a systematic error ascribed to inhomogeneity of the sample at its edges in measuring the coating loss angle. We present a model for disk-shaped resonators, largely used in loss angle measurements, and we compare the theory with measurements showing how this systematic error impacts on the accuracy with which the loss model parameters are known.

Very-long-baseline interferometry (VLBI) is a challenging observational technique, which requires in-depth knowledge about radio telescope instrumentation, interferometry, and the handling of noisy data. The reduction of the raw data is mostly left to the scientists and demands the use of complex algorithms implemented in comprehensive software packages. The correct application of these algorithms necessitates a good understanding of the underlying techniques and physics that are at play. The verification of the processed data produced by the algorithms demands a thorough understanding of the underlying interferometric VLBI measurements. This review describes the latest techniques and algorithms that scientists should know about when analyzing VLBI data.

We investigate theoretical and observational aspects of a warm inflation scenario driven by the $\beta$-exponential potential, which generalizes the well-known power law inflation. In such a scenario, the decay of the inflaton field into radiation happens during the inflationary phase. In our study, we consider a dissipation coefficient ($\Gamma$) with cubic dependence on the temperature ($T$) and investigate the consequences in the inflationary dynamics, focusing on the impact on the spectral index $n_s$, its running $n_{run}$ and tensor-to-scalar ratio $r$. We find it possible to realize inflation in agreement with current cosmic microwave background data in weak and strong dissipation regimes. We also investigate theoretical aspects of the model in light of the swampland conjectures, as warm inflation in the strong dissipation regime has been known as a way to satisfy the three conditions currently discussed in the literature. We find that when $\Gamma\propto T^3$, the $\beta$-exponential model can be accommodated into the conjectures.

We present a simulation experiment of a pipeline based on machine learning algorithms for neutral hydrogen (HI) intensity mapping (IM) surveys with different telescopes. The simulation is conducted on HI signals, foreground emission, thermal noise from instruments, strong radio frequency interference (sRFI), and mild RFI (mRFI). We apply the Mini-Batch K-Means algorithm to identify sRFI, and Adam algorithm to remove foregrounds and mRFI. Results show that there exists a threshold of the sRFI amplitudes above which the performance of our pipeline enhances greatly. In removing foregrounds and mRFI, the performance of our pipeline is shown to have little dependence on the apertures of telescopes. In addition, the results show that there are thresholds of the signal amplitudes from which the performance of our pipeline begins to change rapidly. We consider all these thresholds as the edges of the signal amplitude ranges in which our pipeline can function well. Our work, for the first time, explores the feasibility of applying machine learning algorithms in the pipeline of IM surveys, especially for large surveys with the next-generation telescopes.

The probability of primordial black hole (PBH) formation is known to be boosted during the Quantum Chromodynamics (QCD) crossover due to a slight reduction of the equation of state. This induces a high peak and other features in the PBH mass distribution. But the impact of this variation during the PBH formation has been so far neglected. In this work we simulate for the first time the formation of PBHs by taking into account the varying equation of state at the QCD epoch, compute the over-density threshold using different curvature profiles and find that the resulting PBH mass distributions are significantly impacted. The expected merger rate distributions of early and late PBH binaries is comparable to the ones inferred from the GWTC-3 catalog for dark matter fractions in PBHs within $0.1 < f_{\rm PBH} <1 $. The distribution of gravitational-wave events estimated from the volume sensitivity could explain mergers around $30-50 M_\odot$, with asymmetric masses like GW190814, or in the pair-instability mass gap like GW190521. However, none of the considered cases leads to a multi-modal distribution with a secondary peak around $8-15 M_\odot$, as suggested by the GWTC-3 catalog, possibly pointing to a mixed population of astrophysical and primordial black holes.

Planetary atmospheres are commonly thought to result from the efficient outgassing of cooling magma oceans. During this stage, vigorous convective motions in the molten interior are believed to rapidly transport the dissolved volatiles to shallow depths where they exsolve and burst at the surface. This assumption of efficient degassing and atmosphere formation has important implications for planetary evolution, but has never been tested against fluid dynamics considerations. Yet, during a convective cycle, only a finite fraction of the magma ocean can reach the shallow depths where volatiles exsolution can occur, and a large-scale circulation may prevent a substantial magma ocean volume from rapidly reaching the planetary surface. Therefore, we conducted computational fluid dynamics experiments of vigorous 2D and 3D Rayleigh-B\'enard convection at Prandtl number of unity to characterize the ability of the convecting fluid to reach shallow depths at which volatiles are exsolved and extracted to the atmosphere. Outgassing efficiency is essentially a function of the magnitude of the convective velocities. This allows deriving simple expressions to predict the time evolution of the amount of outgassed volatiles as a function of the magma ocean governing parameters. For plausible cases, the time required to exsolve all oversaturated water can exceed the magma ocean lifetime in a given highly vigorous transient stage, leading to incomplete or even negligible outgassing. Furthermore, the planet size and the initial magma ocean water content, through the convective vigor and the exsolution depth, respectively, strongly affect magma oceans degassing efficiency, possibly leading to divergent planetary evolution paths and resulting surface conditions. Overall, despite vigorous convection, for a significant range of parameters, convective degassing appears not as efficient as previously thought.

We calculate the primordial power spectrum of tensor perturbations, within the emergent universe scenario, incorporating a version of the Continuous Spontaneous Localization (CSL) model as a mechanism capable of: breaking the initial symmetries of the system, generating the perturbations, and also achieving the quantum-to-classical transition of such perturbations. We analyze how the CSL model modifies the characteristics of the B-mode CMB polarization power spectrum, and we explore their differences with current predictions from the standard concordance cosmological model. We have found that, regardless of the CSL mechanism, a confirmed detection of primordial B-modes that fits to a high degree of precision the shape of the spectrum predicted from the concordance $\Lambda$CDM model, would rule out one of the distinguishing features of the emergent universe. Namely, achieving a best fit to the data consistent with the suppression observed in the low multipoles of the angular power spectrum of the temperature anisotropy of the CMB. On the contrary, a confirmed detection that accurately exhibits a suppression of the low multipoles in the B-modes, would be a new feature that could be considered as a favorable evidence for the emergent scenario.

We construct a set of hyperonic equations of state (EoS) by assuming SU(3) symmetry within the baryon octet and by using a covariant density functional (CDF) theory approach. The low-density regions of our EoS are constrained by terrestrial experiments, while the high-density regime is modeled by systematically varying the nuclear matter skewness coefficient $Q_{\rm sat}$ and the symmetry energy slope $L_{\rm sym}$. The sensitivity of the EoS predictions is explored in terms of $z$ parameter of the SU(3) symmetric model that modifies the meson-hyperon coupling constants away from their SU(6) symmetric values. Our results show that model EoS based on our approach can support static Tolman-Oppenheimer-Volkof (TOV) masses in the range $2.3$-$2.5\,M_{\odot}$ in the large-$Q_{\rm sat}$ and small-$z$ regime, however, such stars contain only a trace amount of hyperons compared to SU(6) models. We also construct uniformly rotating Keplerian configurations for our model EoS for which the masses of stellar sequences may reach up to $3.0\,M_{\odot}$. These results are used to explore the systematic dependence of the ratio of maximum masses of rotating and static stars, the lower bound on the rotational frequency of the models that will allow secondary masses in the gravitational waves events to be compact stars with $M_2 \lesssim 3.0\,M_{\odot}$ and the strangeness fraction on the model parameters. We conclude that very massive stellar models can be, in principle, constructed within the SU(3) symmetric model, however, they are nucleonic-like as their strangeness fraction drops below 3\%.

We study the response function of the Einstein Telescope to kinematic Doppler anisotropies, which represent one of the guaranteed properties of the stochastic gravitational wave background. If the frequency dependence of the stochastic background changes slope within the detector frequency band, the Doppler anisotropic contribution to the signal can not be factorized in a part depending on frequency, and a part depending on direction. For the first time, we study the detector response function to Doppler anisotropies without making any factorizable Ansatz. Moreover, we do not assume that kinematic effects are small, and we derive general formulas valid for any relative velocity among frames. We apply our findings to three well-motivated examples of background profiles: power-law, broken power-law, and models with a resonance motivated by primordial black hole scenarios. We derive the signal-to-noise ratio associated with an optimal estimator for the detection of non-factorizable kinematic anisotropies, and we study it for representative examples.

The Lunar Lander Neutron and Dosimetry (LND) Experiment aboard the Chang$'$E-4 Lander on the lunar-far side measures energetic charged and neutral particles and monitors the corresponding radiation levels. During solar quiet times, galactic cosmic rays (GCRs) are the dominating component of charged particles on the lunar surface. Moreover, the interaction of GCRs with the lunar regolith also results in upward directed albedo protons which are measured by the LND. In this work, we used calibrated LND data to study the GCR primary and albedo protons. We calculate the averaged GCR proton spectrum in the range of 9 368 MeV and the averaged albedo proton flux between 64.7 and 76.7 MeV from June 2019 (the 7th lunar day after Chang$'$E-4$'$s landing) to July 2020 (the 20th lunar day). We compare the primary proton measurements of LND with the Electron Proton Helium INstrument (EPHIN) on SOHO. The comparison shows a reasonable agreement of the GCR proton spectra among different instruments and illustrates the capability of LND. Likewise, the albedo proton measurements of LND are also comparable with measurements by the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) during solar minimum. Our measurements confirm predictions from the Radiation Environment and Dose at the Moon (REDMoon) model. Finally, we provide the ratio of albedo protons to primary protons for measurements in the energy range of 64.7-76.7 MeV which confirms simulations over a broader energy range.

Neutron stars are the densest objects in the Universe. In this paper we consider so-called inner crust - the layer, where neutron-excess nuclei are immersed into degenerate gas of electrons and sea of quasi-free neutrons. It was generally believed that spherical nuclei become unstable with respect to quadrupole deformations at high densities and here we consider this instability. Within perturbative approach we show that spherical nuclei with equilibrium number density are, in fact, stable with respect to infinitesimal quadrupole deformation. This is due to background of degenerate electrons and associated electrostatic potential which maintain stability of spherical nuclei. However, if the number of atomic nuclei per unit volume is much less than the equilibrium value, instability can arise. To avoid confusion we stress that our results are limited to infinitesimal deformations and do not guaranty strict thermodynamic stability of spherical nuclei. In particular, they does not exclude that substantially non-spherical nuclei (so-called pasta phase) represent thermodynamic equilibrium state of the densest layers of neutron star crust. Rather our results points that spherical nuclei can be metastable even if they are not energetically favourable and the timescale of transformation of spherical nuclei to the pasta phases should be estimated subsequently.

Owed to their compactness, neutron stars involve strong gravity and extreme density physics. Nevertheless, at present, there are a variety of problems where progress (at least conceptually) can be made in the context of weak gravity. Motivated by this we examine how accurately one can model neutron stars using the post-Newtonian approximation to general relativity. In general, we find there is a significant degree of freedom in how the post-Newtonian equations of stellar structure can be formulated. We discuss this flexibility in the formulation and provide examples to demonstrate the impact on stellar models. We also consider the (closely related) problem of building neutron stars using isotropic coordinates. In this context, we provide a new strategy for solving the equations (based on a scaling argument) which significantly simplifies the problem.

The effects of strong magnetic fields on the deconfinement phase transition expected to take place in the interior of massive neutron stars is studied in detail for the first time. For hadronic matter, the very general density-dependent relativistic mean-field (DD-RMF) model is employed, while the simple, but effective Vector-Enhanced Bag model (vBag) model is used to study quark matter. Magnetic-field effects are incorporated into the matter equation of state and in the general-relativity solutions, which also satisfy Maxwell's equations. We find that, for large values of magnetic dipole moment, the maximum mass, canonical-mass radius, and dimensionless tidal deformability obtained for stars using spherically-symmetric TOV equations and axisymmetric solutions attained through the LORENE library differ considerably. The deviations depend on the stiffness of the equation of state and on the star mass being analyzed. This points to the fact that, unlike what was assumed previously in the literature, magnetic field thresholds for the correct assumption of isotropic stars and the proper use of TOV equations depend on the matter composition and interactions.

Equations of state for a cold neutron star's interior are presented in three-column tables that relate the baryonic density, the energy density, and the pressure. A few analytical expressions for those tables have been established these past two decades, as a convenient way to present a large number of nuclear models for neutron star matter. Some of those analytical representations are based on nonunified equations of state, in the sense that the high and the low density part of the star are not computed with the same nuclear model. Fits of equations of state based on a piecewise polytropic representation are revised by using unified tables of equations of state, that is to say models which have been calculated consistently for the core and the crust. A set of 52 unified equations of state is chosen. Each one is divided in seven polytropes via an adaptive segmentation, and two parameters per polytrope are fitted to the tabulated equation of state. The total mass, radius, tidal deformability and moment of inertia of neutron stars are modelled from the fits and compared with the quantities calculated from the original tables to ensure the accuracy of the fits on macroscopic parameters. We provide the polytropes parameters for 15 nucleonic relativistic mean field models, seven hyperonic relativistic mean field models, five hybrid relativistic mean-field models, 24 nucleonic Skyrme models, and one ab initio model. The fit error on the macroscopic parameters of neutron stars is small and well within the estimated measurement accuracy from current and next generation telescopes.

High resolution polarization maps of the Cosmic Microwave Background (CMB) are on high demand, since the discovery of primordial B-Modes in the polarization patterns would confirm the inflationary phase of the Universe that would have taken place before the emission of the CMB. Transition Edge Sensors (TES) and Microwave Kinetic Inductance Detectors (MKID) are the predominant detector technologies of cryogenic detector array based CMB instruments that search for primordial B-Modes. In this paper we propose another type of cryogenic detector to be used for CMB survey: A magnetic microbolometer (MMB) that is based on a paramagnetic temperature sensor. It is an adaption of state-of-the-art metallic magnetic calorimeters (MMCs) that are meanwhile a key technology for high resolution $\alpha$, $\beta$, $\gamma$ and X-ray spectroscopy as well as the study of neutrino mass. A complete simulation framework was developed that accounts for the electrical and thermal properties of the bolometer and that can be used to obtain its responsivity and bandwidth, as well as estimating noise. A brief proof of concept case study is analyzed, taking into account typical constraints in CMB measurements and reliable microfabrication processes, to assess the suitability of metallic magnetic sensors in CMB experiments. The results show that MMBs provide a promising technology for CMB polarization measurements as their sensitivity can be tuned for background limited detection of the sky while simultaneously maintaining a low time response to avoid degradation of the point-source response of the telescope. As the sensor technology and its fabrication techniques are compatible with TES based bolometric detector arrays, a change of detector technology would even come with very low cost.

Expressions for the BER of M-ary PPM & biphase DPSK modulations in the presence of noise are derived using analytical, statistical methods. The PPM expression is verified via Poisson statistics based simulation. BER expressions are then applied to a representative set of receiving telescope & sky spectral radiance parameters in order to assess performance of PPM & DPSK relative to one another. Finally, efficiency & additional considerations are discussed.