Papers for Friday, Nov 18 2022

With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of Simulation-Based inference (SBI) and amortized Neural Posterior Estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; and b) amortized calculations of posterior distributions of said galaxy parameters at the modest cost of a few galaxy fits with MCMC methods. We utilise the SED generator and inference framework Prospector to generate simulated spectra, and train a dataset of 2$\times$10$^6$ spectra (corresponding to a 5-parameter SED model) with NPE. We show that SBI -- with its combination of fast and amortized posterior estimations -- is capable of inferring accurate galaxy stellar masses and metallicities. Our uncertainty constraints are comparable to or moderately weaker than traditional inverse-modeling with Bayesian MCMC methods (e.g., 0.17 and 0.26 dex in stellar mass and metallicity for a given galaxy, respectively). We also find that our inference framework conducts rapid SED inference (0.9-1.2$\times$10$^5$ galaxy spectra via SBI/SNPE at the cost of 1 MCMC-based fit). With this work, we set the stage for further work that focuses of SED fitting of galaxy spectra with SBI, in the era of JWST galaxy survey programs and the wide-field Roman Space Telescope spectroscopic surveys.

We present re-processed flux calibrated spectra of 406 stars from the UVES-POP stellar library in the wavelength range 320-1025 nm, which can be used for stellar population synthesis. The spectra are provided in the two versions having spectral resolving power R=20,000 and R=80,000. Raw spectra from the ESO data archive were re-reduced using the latest version of the UVES data reduction pipeline with some additional algorithms that we developed. The most significant improvements in comparison with the original UVES-POP release are: (i) an updated Echelle order merging, which eliminates "ripples" present in the published spectra, (ii) a full telluric correction, (iii) merging of non-overlapping UVES spectral setups taking into account the global continuum shape, (iv) a spectrophotometric correction and absolute flux calibration, and (v) estimates of the interstellar extinction. For 364 stars from our sample, we computed atmospheric parameters $T_\mathrm{eff}$, surface gravity log $g$, metallicity [Fe/H], and $\alpha$-element enhancement [$\alpha$/Fe] by using a full spectrum fitting technique based on a grid of synthetic stellar atmospheres and a novel minimization algorithm. We also provide projected rotational velocity $v\sin i$ and radial velocity $v_{rad}$ estimates. The overall absolute flux uncertainty in the re-processed dataset is better than 2% with sub-% accuracy for about half of the stars. A comparison of the recalibrated UVES-POP spectra with other spectral libraries shows a very good agreement in flux; at the same time, $Gaia$ DR3 BP/RP spectra are often discrepant with our data, which we attribute to spectrophotometric calibration issues in $Gaia$ DR3.

The gravitational waves from the binary neutron star merger GW170817 were accompanied by a multi-wavelength electromagnetic counterpart, which confirms the association of the merger with a short gamma-ray burst (sGRB). The afterglow observations implied that the event was accompanied by a narrow, $\sim 5~$deg, and powerful, $\sim 10^{50}$ erg, jet. We study the propagation of a Poynting flux-dominated jet within the merger ejecta (kinematic, neutrino-driven and MRI turbulence-driven) of a neutrino-radiation-GR-MHD simulation of two coalescing neutron stars. We find that the presence of a post-merger low-density/low-pressure polar cavity, that arose due to angular momentum conservation, is crucial to let the jet break out. At the same time the ejecta collimates the jet to a narrow opening angle. The collimated jet has a narrow opening angle of $\sim 4$-$7$ deg and an energy of $10^{49}$-$10^{50}~$erg, in line with the observations of GW170817 and other sGRBs.

Subhalo abundance matching (SHAM) has played an important role in improving our understanding of how galaxies populate their host dark matter halos. In essence, the SHAM framework is to find a dark matter halo property that best correlates with an attribute of galaxies, such as stellar mass. The peak value of the maximum circular velocity ($V_{\rm max}$) a halo/subhalo has ever attained throughout its lifetime, $V_{\rm peak}$, has been a popular choice for SHAM. A recent study by Tonnesen & Ostriker (2021) suggested that quantity $\phi$, which combines the present-day $V_{\rm max}$ and the peak value of halo dark matter mass, performs better in predicting stellar mass than $V_{\rm peak}$. Inspired by their approach, in this work, we find that further improvement can be achieved by a quantity $\psi_5$ that combines the 90th percentile of $V_{\rm max}$ a halo/subhalo has ever achieved with the 60th percentile of the dark matter halo time variation rate. Tests based on the simulation IllustrisTNG300 show that our new SHAM scheme, with just three free parameters, can improve the stellar mass prediction and mass-dependent clustering by 15% and 16% from $\phi$, respectively, over the redshift range $z=0-2$.

We present results of a Transiting Exoplanet Survey Satellite (TESS) search for short-period pulsations in compact stellar objects observed in years 1 and 3 of the TESS mission, during which the southern ecliptic hemisphere was targeted. We describe the TESS data used and the details of the search method. For many of the targets, we use unpublished spectroscopic observations to classify the objects. From the TESS photometry, we clearly identify 43 short-period hot-subdwarf pulsators, including 32 sdB stars, eight sdOB stars, two sdO stars, and, significantly, one He-sdOB star, which is the first of this kind to show short-period pulsations. Eight stars show signals at both low and high frequencies, and are therefore ``hybrid'' pulsators. We report the list of prewhitened frequencies and we show the amplitude spectra calculated from the TESS data. We make an attempt to identify possible multiplets caused by stellar rotation, and we select four candidates with rotation periods between 1 and 12.9d. The most interesting targets discovered in this survey should be observed throughout the remainder of the TESS mission and from the ground. Asteroseismic investigations of these data sets will be invaluable in revealing the interior structure of these stars and will boost our understanding of their evolutionary history. We find three additional new variable stars but their spectral and variability types remain to be constrained.

Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing with thermo-turbulent star formation and mechanical supernova feedback, we gradually increase the complexity of physics incorporated leading to full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate combinations of these processes, comparing them with observations to constrain what are the main mechanisms determining dwarf galaxy properties. We find hydrodynamical `SN feedback-only' simulations struggle to produce realistic dwarf galaxies, leading either to overquenched or too centrally concentrated, dispersion dominated systems when compared to observed field dwarfs. Accounting for radiation with cosmic rays results in extended and rotationally-supported systems. Spatially `distributed' feedback leads to realistic stellar and HI masses as well as kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, HI and HII after star formation events. Episodic star formation combined with its associated feedback induces more core-like dark matter central profiles, which our `SN feedback-only' models struggle to achieve. Our results demonstrate the complexity of physical processes required to capture realistic dwarf galaxy properties, making tangible predictions for integral field unit surveys, radio synchrotron emission, and for galaxy and multi-phase interstellar medium properties that JWST will probe.

We aim to investigate the AGN/ICM interplay in ZwCl 235, a galaxy cluster with high X-ray flux, an extended central radio galaxy, and evidence of multi-phase gas at its center. Using archival data from the Chandra telescope, the VLASS survey, the LOTSS survey and the VLBA telescope, we perform a complete study of ZwCl 235, dissecting the dynamics of the ICM, the thermodynamic state of the central gas, and the properties of the BCG. By means of radial profiles and 2D spectral maps, we measure the temperature, entropy and cooling time of the ICM, and we compare the morphology of the central radio galaxy with the surrounding medium. We find evidence that ZwCl 235 is a sloshing cool core cluster in which the activity of the central galaxy has excavated a pair of cavities and possibly uplifted enriched gas to an altitude of ~30 kpc. In the cluster core, the lowest entropy ICM is preferentially found in a 20 kpc-long filament tangential to the southern radio lobe of the AGN. We argue that the observed cool (~1.3 keV) filament is likely produced by a combination of sloshing and stimulated ICM cooling, that may be fueling the central supermassive black hole. Additionally, we determine that the X-ray emission of the BCG originates from a ~1.4 keV plasma kernel which extends for 5 kpc in radius and has a short cooling time (~240 Myr), and could represent the thermal corona of the BCG. Overall, we propose that several sources (the large scale ICM, the low entropy filament and the ~1.4 keV kernel) of cold material are currently feeding the central AGN, and that the ICM cooling cycle expectations are met from the macro scales (between 5 - 100 kpc) to the meso scales (<5 kpc) of AGN feedback.

We forward-model mass-weighted ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of $\sim 1$ Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the UniverseMachine model. We find total quenching times of $t_\mathrm{Q}=3.7\pm 0.4$ Gyr and $t_\mathrm{Q}=4.0\pm 0.2$ Gyr after first pericentre, for $9<\log(M_{\star}/\mathrm{M}_{\odot})<10$ and $10<\log(M_{\star}/\mathrm{M}_{\odot})<10.5$ galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and timescale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for our lower and higher stellar mass bins, respectively, and exponential SFR suppression timescales are $\tau_{\mathrm{env}}\leq 1.0$ Gyr and $\tau_{\mathrm{env}}\sim 2.3$ Gyr for our lower and higher stellar mass bins, respectively. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the timescale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage.

Massive ($\geq$8 $M_\odot$) stars perish via one of two fates: core-collapse supernovae (CCSNe), which release synthesized heavy elements, or failed supernovae, thereby forming black holes. In the conventional Galactic chemical evolution (GCE) scheme, a substantial portion of massive stars, e.g., all stars in the mass range of 8-100 $M_{\odot},$ are assumed to enrich the Galaxy with their nucleosynthetic products. However, this hypothesis conflicts with the observations, namely, few CCSNe whose progenitor stars are more massive than $\sim$18 $M_{\odot}.$ Here, we show that the chemical characteristics shaped by local thin disk stars are compatible with the predictions by enrichment via CCSNe limited to less massive progenitors in the new paradigm of Galactic dynamics that allows stars to migrate from the inner disk. This renewed GCE model predicts that the bursting star formation events$-$which are considered to take place in the Galactic bulge as well as in the thick disk$-$generate more numerous low-mass CCSNe than those expected from the locally determined canonical initial mass function. This finding suggests a high rate of CCSNe in early-type galaxies, which reflects a unique cosmic history of the CCSN rate. With considerable contributions from these galaxies to the cosmic star formation rates in the early Universe, we predict a more steeply increasing slope of the CCSN rate with increasing redshift than that in proportion to cosmic star formation. This predicted redshift evolution agrees well with the measured rates for 0 $\lesssim$ z $\lesssim$ 0.8; however, its predicted CCSN rate for higher-$z$ calls for more precise data from future surveys.

In galaxy clusters, the hot intracluster medium (ICM) can develop a striking multi-phase structure around the brightest cluster galaxy. Much work has been done on understanding the origin of this central nebula, but less work has studied its eventual fate after the originally filamentary structure is broken into individual cold clumps. In this paper we perform a suite of 30 (magneto-)hydrodynamical simulations of kpc-scale cold clouds with typical parameters as found by galaxy cluster simulations, to understand whether clouds are mixed back into the hot ICM or can persist. We investigate the effects of radiative cooling, small-scale heating, magnetic fields, and (anisotropic) thermal conduction on the long-term evolution of clouds. We find that filament fragments cool on timescales shorter than the crushing timescale, fall out of pressure equilibrium with the hot medium, and shatter, forming smaller clumplets. These act as nucleation sites for further condensation, and mixing via Kelvin-Helmholtz instability, causing cold gas mass to double within 75 Myr. Cloud growth depends on density, as well as on local heating processes, which determine whether clouds undergo ablation- or shattering-driven evolution. Magnetic fields slow down but don't prevent cloud growth, with the evolution of both cold and warm phase sensitive to the field topology. Counter-intuitively, anisotropic thermal conduction increases the cold gas growth rate compared to non-conductive clouds, leading to larger amounts of warm phase as well. We conclude that dense clumps on scales of $500$ pc or more cannot be ignored when studying the long-term cooling flow evolution of galaxy clusters.

The Sloan Digital Sky Survey IV (SDSS-IV) APOGEE-2 primary science goal was to observe red giant stars throughout the Galaxy to study its dynamics, morphology, and chemical evolution. The APOGEE instrument, a high-resolution 300 fiber H-band (1.55-1.71 micron) spectrograph, is also ideal to study other stellar populations in the Galaxy, among which are a number of star forming regions and young open clusters. We present the results of the determination of six stellar properties ($T_{eff}$, $\log{g}$, [Fe/H], $L/L_\odot$, $M/M_\odot$, and ages) for a sample that is composed of 3360 young stars, of sub-solar to super-solar types, in sixteen Galactic star formation and young open cluster regions. Those sources were selected by using a clustering method that removes most of the field contamination. Samples were also refined by removing targets affected by various systematic effects of the parameter determination. The final samples are presented in a comprehensive catalog that includes all six estimated parameters. This overview study also includes parameter spatial distribution maps for all regions and Hertzprung-Russell ($L/L_\odot$ vs. $T_{eff}$) diagrams. This study serves as a guide for detailed studies on individual regions, and paves the way for the future studies on the global properties of stars in the pre-main sequence phase of stellar evolution using more robust samples.

According to the current experimental data, the Higgs vacuum appears to be metastable due to the development of a second lower ground state in its potential. Consequently, vacuum decay would induce the nucleation of true vacuum bubbles with catastrophic consequences for our Universe and therefore we are motivated to study possible stabilising mechanisms in the early universe. In our latest investigation (2207.00696), we studied the electroweak metastability in the context of the observationally favoured model of Starobinsky inflation. Following the motivation and techniques from our first study (2011.037633), we obtained constraints on the Higgs curvature coupling $\xi$, while embedding the SM on the modified gravity scenario $R+R^2$, which introduces Starobinsky inflation naturally. This had significant repercussions for the effective Higgs potential in the form of additional negative terms that destabilize the false vacuum. Another important aspect lay in the definition for the end of inflation, as bubble nucleation is most prominent during its very last moments. Our results dictated that these stronger lower $\xi$-bounds are very sensitive to the final moments of inflation, where spacetime deviates increasingly from de Sitter.

We analyze jet features found in the coma of Comet 21P/Giacobini-Zinner (21P/GZ) during its 2018 perihelion passage using narrowband CN photometric imaging in order to determine the comet's rotational period, and constrain the CN gas outflow velocity and rotational state through the analysis of azimuthally enhanced morphological features. We find that 21P/GZ has a periodicity of either 7.4 $\pm$ 0.1 or 10.7 $\pm$ 0.2 hours. We measure a lower limit to the outflow velocity for the Northern Jet of 730 $\pm$ 30 m s$^{-1}$ and for the Southern Jet of 740 $\pm$ 30 m s$^{-1}$. We analyze the morphologies of the jet features and determine that the Northern pointing jet possesses a corkscrew pattern, and utilize that knowledge to determine a rotational pole position at a Right Ascension of ${169^{+28}_{-23}}^{\circ}$ and a Declination of ${73^{+5}_{-11}}^{\circ}$, with undetermined sense of rotation.

We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prior simulation evidence, the proposed mechanism reproduces the known scaling of the viscous $\alpha$-parameter, $\alpha \propto \beta_z^{-1/2}$. Here, $\beta_z = 8\pi p_g/B_{z0}^2$ is the initial plasma $\beta$-parameter on the disc midplane, $B_{z0}$ is the net field, and $p_g $ is the midplane gas pressure. We generalize the argument to discs with strong suprathermal toroidal magnetic fields, where the MRI growth rate is modified from the weak-field limit. Additional sources of turbulence are required if such discs are to become magnetically elevated, with the increased scale heights near the midplane that are seen in simulations. We speculate that tearing modes, associated with current sheets broadened by the effective resistivity, are a possible source of enhanced turbulence in elevated discs.

We report on Chandra X-ray observations of ASASSN-18tb/SN 2018fhw, a low luminosity Type Ia supernova that showed a H line in its optical spectrum. No X-ray emission was detected at the location of the SN. Upper limits to the luminosity of up to 3 $\times 10^{39}$ erg s$^{-1}$ are calculated, depending on the assumed spectral model, temperature and column density. These are compared to two Type Ia-CSM SNe, SN 2005gj and SN 2002ic, that have been observed with Chandra in the past. The upper limits are lower than the X-ray luminosity found for the Type Ia-CSM SN 2012ca, the only Type Ia SN to have been detected in X-rays. Consideration of various scenarios for the H$\alpha$ line suggests that the density of the surrounding medium at the time of H$\alpha$ line detection could have been as high as 10$^8$ cm$^{-3}$, but must have decreased below 5 $\times\, 10^6$ cm$^{-3}$ at the time of X-ray observation. Continual X-ray observations of SNe which show a H line in their spectrum are necessary in order to establish Type Ia SNe as an X-ray emitting class.

The Mg II h&k lines are amongst the best diagnostic tools of the upper solar chromosphere. This region of the atmosphere is of particular interest as it is the lowest region of the Sun's atmosphere where the magnetic field is dominant in the energetics and dynamics, defining its structure. While highly successful in the photosphere and lower to mid chromosphere, numerical models have produced syntheticMg II lines that do not match the observations well. We present a number of large scale models with magnetic field topologies representative of the quiet Sun, ephemeral flux regions and plage, and also models where the numerical resolution is high and where we go beyond the MHD paradigm. The results of this study show models with a much improved correspondence with \iris\ observations both in terms of intensities and widths, especially underscoring the importance of chromospheric mass loading and of capturing the magnetic field topology and evolution in simulations. This comes in addition to the importance of capturing the generation of small scale velocity fields and including non-equilibrium ionization and ion-neutral interaction effects. Understanding and modeling all these effects and their relative importance is necessary in order to reproduce observed spectral features.

The POLISH2 optical polarimeter has been in operation at the Lick Observatory 3-m Shane telescope since 2011, and it was commissioned at the Gemini North 8-m in 2016. This instrument primarily targets exoplanets, asteroids, and the Crab pulsar, but it has also been used for a wide variety of planetary, galactic, and supernova science. POLISH2's photoelastic modulators, employed instead of rotating waveplates or ferro-electric liquid crystal modulators, offer the unprecedented ability to achieve sensitivity and accuracy of order 1 ppm (0.0001%), which are difficult to obtain with conventional polarimeters. Additionally, POLISH2 simultaneously measures intensity (Stokes I), linear polarization (Stokes Q and U), and circular polarization (Stokes V), which fully describe the polarization state of incident light. We document our laboratory and on-sky calibration methodology, our archival on-sky database, and we demonstrate conclusive detection of circular polarization of certain objects.

Tadpole galaxies are metal-poor dwarfs with typically one dominant star-forming region, giving them a head-tail structure when inclined. A metallicity drop in the head suggests that gas accretion with even lower metallicity stimulated the star formation. Here we present multiband HST WFC3 and ACS images of four nearby (<25 Mpc) tadpoles, SBS0, SBS1, Kiso 3867, and UM461, selected for their clear metallicity drops shown in previous spectroscopic studies. Properties of the star complexes and compact clusters are measured. Each galaxy contains from 3 to 10 young stellar complexes with 10^3-10^5 Msun of stars ~3-10 Myr old. Between the complexes, the disk has a typical age of ~3 Gyr. Numerous star clusters cover the galaxies, both inside and outside the complexes. The combined cluster mass function, made by normalizing the masses and counts before stacking, is a power law with a slope of -1.12+-0.14 on a log-log plot and the combined distribution function of cluster lifetime decays with age as t^{-0.65+-0.24}. A comparison between the summed theoretical Lyman continuum (LyC) emission from all the clusters, given their masses and ages, is comparable to or exceeds the LyC needed to excite the observed Halpha in some galaxies, suggesting LyC absorption by dust or undetected gas in the halo, or perhaps galaxy escape.

Coadded astronomical images are created by stacking multiple single-exposure images. Because coadded images are smaller in terms of data size than the single-exposure images they summarize, loading and processing them is less computationally expensive. However, image coaddition introduces additional dependence among pixels, which complicates principled statistical analysis of them. We present a principled Bayesian approach for performing light source parameter inference with coadded astronomical images. Our method implicitly marginalizes over the single-exposure pixel intensities that contribute to the coadded images, giving it the computational efficiency necessary to scale to next-generation astronomical surveys. As a proof of concept, we show that our method for estimating the locations and fluxes of stars using simulated coadds outperforms a method trained on single-exposure images.

The Terahertz Intensity Mapper (TIM) is designed to probe the star formation history in dust-obscured star-forming galaxies around the peak of cosmic star formation. This will be done via measurements of the redshifted 157.7 um line of singly ionized carbon ([CII]). TIM employs two R $\sim 250$ long-slit grating spectrometers covering 240-420 um. Each is equipped with a focal plane unit containing 4 wafer-sized subarrays of horn-coupled aluminum kinetic inductance detectors (KIDs). We present the design and performance of a prototype focal plane assembly for one of TIM's KID-based subarrays. Our design strictly maintain high optical efficiency and a suitable electromagnetic environment for the KIDs. The prototype detector housing in combination with the first flight-like quadrant are tested at 250 mK. Initial frequency scan shows that many resonances are affected by collisions and/or very shallow transmission dips as a result of a degraded internal quality factor (Q factor). This is attributed to the presence of an external magnetic field during cooldown. We report on a study of magnetic field dependence of the Q factor of our quadrant array. We implement a Helmholtz coil to vary the magnetic field at the detectors by (partially) nulling earth's. Our investigation shows that the earth magnetic field can significantly affect our KIDs' performance by degrading the Q factor by a factor of 2-5, well below those expected from the operational temperature or optical loading. We find that we can sufficiently recover our detectors' quality factor by tuning the current in the coils to generate a field that matches earth's magnetic field in magnitude to within a few uT. Therefore, it is necessary to employ a properly designed magnetic shield enclosing the TIM focal plane unit. Based on the results presented in this paper, we set a shielding requirement of |B| < 3 uT.

A comprehensive summary of the measurements made to characterize test mass charging due to the space environment during the LISA Pathfinder mission is presented. Measurements of the residual charge of the test mass after release by the grabbing and positioning mechanism, show that the initial charge of the test masses was negative after all releases, leaving the test mass with a potential in the range $-12$ mV to $-512$ mV. Variations in the neutral test mass charging rate between $21.7$ e s$^{-1}$ and $30.7$ e s$^{-1}$ were observed over the course of the 17-month science operations produced by cosmic ray flux changes including a Forbush decrease associated with a small solar energetic particle event. A dependence of the cosmic ray charging rate on the test mass potential between $-30.2$ e s$^{-1}$ V$^{-1}$ and $-40.3$ e s$^{-1}$ V$^{-1}$ was observed and this is attributed to a contribution to charging from low-energy electrons emitted from the gold surfaces of the gravitational reference sensor. Data from the on-board particle detector show a reliable correlation with the charging rate and with other environmental monitors of the cosmic ray flux. This correlation is exploited to extrapolate test mass charging rates to a 20-year period giving useful insight into the expected range of charging rate that may be observed in the LISA mission.

Neutrino-neutrino forward scatterings potentially induce collective neutrino oscillation in dense neutrino gases in astrophysical sites such as core-collapse supernovae (CCSN) and binary neutron star mergers (BNSM). In this paper, we present a detailed study of fast neutrino-flavor conversion (FFC), paying special attention to asymptotic states, by linear stability analysis and local simulations with a periodic boundary condition. We find that asymptotic states can be characterized by two key properties of FFC: (1) the conservation of lepton number for each flavor of neutrinos and (2) the disappearance of ELN(electron neutrino-lepton number)-XLN(heavy-leptonic one) angular crossings in the spatial- or time-averaged distributions. The system which initially has the positive (negative) ELN-XLN density reaches a flavor equipartition in the negative (positive) ELN-XLN angular directions, and the other part compensates it to preserve the conservation laws. These properties of FFCs offer an approximate scheme determining the survival probability of neutrinos in asymptotic states without solving quantum kinetic equations. We also demonstrate that the total amount of flavor conversion can vary with species-dependent neutrino distributions for identical ELN-XLN ones. Our results suggest that even shallow or narrow ELN angular crossings have the ability to drive large flavor conversion, exhibiting the need for including the effects of FFCs in the modeling of CCSN and BNSM.

We present the estimation of the solar observation with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). For both the quite Sun and the Sun with radio bursts, when pointing directly to the Sun, the total power received by FAST would be out of the safe operational range of the signal chain, even resulting in the damage to the receiver. As a conclusion, the Sun should be kept at least $\sim 2^{\circ}$ away from the main beam during the observing at $\sim 1.25 {\ \rm GHz}$. The separation for lower frequency should be larger. For simplicity, the angular separation between the FAST beam and the Sun is suggested to be $\sim 5^{\circ}$ for observations on 200 MHz or higher bands.

We report on the kinetic inductance detector (KID) array focal plane assembly design for the Terahertz Intensity Mapper (TIM). Each of the 2 arrays consists of 4 wafer-sized dies (quadrants), and the overall assembly must satisfy thermal and mechanical requirements, while maintaining high optical efficiency and a suitable electromagnetic environment for the KIDs. In particular, our design manages to strictly maintain a 50 $\mathrm{\mu m}$ air gap between the array and the horn block. We have prototyped and are now testing a sub-scale assembly which houses a single quadrant for characterization before integration into the full array. The initial test result shows a $>$95\% yield, indicating a good performance of our TIM detector packaging design.

The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high signal-to-noise, with a narrow central component, and a red wing, confirming previous results. The OH 1720 MHz line is clearly detected in (maser) emission, peaking at a velocity shifted by -10 to -15 km/s with respect to the HI peak. The 1612 MHz line is lost due to radio interferences. The OH 1667 MHz main line is tentatively detected in absorption, but not the 1665 MHz one. Over 30 years, a high variability is observed in optical depths, due to the rapid changes of the line of sight, caused by the superluminal motions of the radio knots. The HI line has varied by 20 per cent in depth, while the OH-1720 MHz depth has varied by a factor 4. The position of the central velocity and the widths also varied. The absorbing galaxy is an early-type spiral (maybe S0) seen edge-on, with a prominent dust lane, covering the whole disk. Given the measured mass concentration, and the radio continuum size at centimeter wavelengths (100 mas corresponding to 400 pc at z = 0.25), the width of absorption lines from the nuclear regions are expected up to 250 km/S. The narrowness of the observed lines (< 15 km/s) suggest that the absorption comes from an outer gas ring, as frequently observed in S0 galaxies. The millimetric lines are even narrower (< 1 km/s), which corresponds to the continuum size restricted to the core. The core source is covered by individual 1 pc molecular clouds, of column density a few 10^22 cm-2, which is compatible with the gas screen detected in X-rays.

The ponderomotive force is suggested to be the main mechanism to produce the so-called first ionization potential (FIP) effect - the enrichment of low FIP elements observed in the outer solar atmosphere. It is well known that the ionization of these elements occurs within the chromosphere. Therefore, this phenomenon is intimately tied to the plasma state in the chromosphere and the corona. In addition, the chromosphere is a highly complex region with a large variation in the ion-neutral collision frequencies, and hydrogen and helium ionization is largely out of equilibrium. For this study, we combine IRIS observations, a single fluid 2.5D radiative magnetohydrodynamics (MHD) model of the solar atmosphere, including ion-neutral interaction effects and non-equilibrium (NEQ) ionization effects, and a novel multi-fluid multi-species (MFMS) numerical code (Ebysus). Nonthermal velocities of Si IV measured from IRIS spectra can provide an upper limit for the strength of any high-frequency Alfven waves. With the single-fluid model, we investigate the possible impact of NEQ ionization within the region where the FIP may occur and the plasma properties in those regions. These models suggest that regions with strong enhanced network and type II spicules are associated with the presence of large ponderomotive forces. From the plasma properties from the single-fluid MHD model and IRIS observations, we initialize our multi-fluid models to investigate the multi-fluid effects on the ponderomotive force associated with Alfven waves. Our multi-fluid analysis reveals that collisions and NEQ ionization effects dramatically impact the behavior of the ponderomotive force in the chromosphere, and existing theories may need to be revisited.

We report on the first observation of the diffuse cosmic neutrino flux with the Baikal-GVD neutrino telescope. Using cascade-like events collected by Baikal-GVD in 2018--2021, a significant excess of events over the expected atmospheric background is observed. This excess is consistent with the high-energy diffuse cosmic neutrino flux observed by IceCube. The null cosmic flux assumption is rejected with a significance of 3.05$\sigma$. Assuming a single power law model of the astrophysical neutrino flux with identical contribution from each neutrino flavor, the following best-fit parameter values are found: the spectral index $\gamma_{astro}$ = $2.58^{+0.27}_{-0.33}$ and the flux normalization $\phi_{astro}$ = 3.04$^{+1.52}_{-1.21}$ per one flavor at 100 TeV.

As technology has improved, binary neutron star systems have been observed more frequently, in fact, the first gravitational wave to have an electromagnetic counterpart originated from the merger of two neutron stars (GW170817). Detecting these systems prior to merger may help recover essential data for developing an Equation of State for neutron stars. This paper examines the observability of detached eclipsing binary neutron stars prior to merger by simulating the potential observability of neutron star systems in the optical. It is found that it is not likely considering current instruments due to low visibility and inadequate time resolution, however, improvements in the future or a wide field X-ray instrument may offer a viable option for detecting these systems.

Assuming the common-spectrum process in the NANOGrav 12.5-year dataset has an origin of scalar induced gravitational waves, we study the enhancement of primordial curvature perturbations and the mass function of primordial black holes, by performing the Bayesian parameter inference for the first time. We obtain lower limits on the spectral amplitude, i.e. $\mathcal{A}\gtrsim10^{-2}$ at 95\% confidence level, when assuming the power spectrum of primordial curvature perturbations to follow a log-normal distribution function with width $\sigma$. In the limit of $\sigma\rightarrow0$, we find that the primordial black holes with $10^{-3}-10^{-2}$ solar mass are allowed to compose at least a fraction of dark matter. Such a mass range would be shifted to more massive regimes for larger values of $\sigma$, e.g. to a regime of $10^{-2}-10^{-1}$ solar mass for $\sigma=1$. We expect that planned gravitational-wave experiments are capable to reach at least $\mathcal{A}\sim10^{-4}$ and search for the primordial black holes over the whole parameter space. In addition, importance of multi-band detector networks is emphasized to accomplish our theoretical expectation.

A balanced ratio of ocean to land is believed to be essential for an Earth-like biosphere and one may conjecture that plate-tectonics planets should be similar in geological properties. After all, the volume of continental crust evolves towards an equilibrium between production and erosion. If the interior thermal states of Earth-sized exoplanets are similar to the Earth's, one might expect a similar equilibrium between continental production and erosion to establish and, hence, a similar land fraction. We will show that this conjecture is not likely to be true. Positive feedback associated with the coupled mantle water - continental crust cycle may rather lead to a manifold of three possible planets, depending on their early history: a land planet, an ocean planet and a balanced Earth-like planet. In addition, thermal blanketing of the interior by the continents enhances the sensitivity of continental growth to its history and, eventually, to initial conditions. Much of the blanketing effect is however compensated by mantle depletion in radioactive elements. A model of the long-term carbonate-silicate cycle shows the land and the ocean planet to differ by about 5 K in average surface temperature. A larger continental surface fraction results both in higher weathering rates and enhanced outgassing, partly compensating each other. Still, the land planet is expected to have a substantially dryer, colder and harsher climate possibly with extended cold deserts in comparison with the ocean planet and with the present-day Earth. Using a model of balancing water availability and nutrients from continental crust weathering, we find the bioproductivity and the biomass of both the land and ocean planet to be reduced by a third to half of Earth's. The biosphere on these planets might not be substantial enough to produce a supply of free oxygen.

We study the relation between stellar dynamo-wave propagation and the structure of the stellar magnetic field. Modeling dynamo waves by the well-known Parker migratory dynamo, we vary the intensity of dynamo drivers in order to obtain activity-wave propagation toward the Equator (as in the solar-activity cycle) or towards the Poles. We match the magnetic field in the dynamo active shell with that in the surrounding stellar material, using a simple dissipativ magnetohydrodynamic system for the transition region. Introducing a weak asymmetry between the stellar hemispheres, we study phase shifts of the dipole, quadrupole, and octupole magnetic components at various distances from the star to demonstrate that several-percent asymmetry in dynamo drivers are sufficient to obtain a realistic relation between solar dipole and quadrupole moments. We study the behavior of the stellar current sheets and show that for the poleward propagating activity it is substantially different from solar ones. In particular, we demonstrate conditions in which the conical current sheets propagate opposite to the solar directions.

We report results obtained from the optical and X-ray studies of the Be/X-ray binary 1A 0535+262/HD 245770 during the 2020 October giant X-ray outburst, using the 1.2 m telescope at Mount Abu Infrared observatory and AstroSat, respectively. The peak flux of the outburst was recorded to be around 11 Crab in the 15-50 keV range, the highest ever observed from the pulsar. We performed optical observations in the 6000-7200 angstroms band before, during, and after the outburst to investigate the evolution of the circumstellar disc of the Be star between 2020 February and 2022 February. Our optical spectra exhibit prominent emission lines at 6563 angstroms (H I), 6678 angstroms (He I), and 7065 angstroms (He I). We found a significantly variable Halpha line in the spectra. The single-peaked line profile appeared asymmetric with broad red- & blue-wings in the data before and during the outburst. The post-outburst observations, however, resulted in a double-peaked profile with asymmetry in the blue-wing. Our observations before the outburst confirmed a larger Be disc that decreased in size as the outburst progressed. Furthermore, the observed variabilities in the Halpha line profile and parameters suggest the presence of a highly misaligned, precessing, and warped Be disc. AstroSat observation of the pulsar detected pulsations at around 103.55 s in the light curve up to 110 keV. We found strongly energy-dependent pulse profiles with increasing contribution of the pulsing component in hard X-rays. The broadband spectral fitting in the 0.7-90.0 keV range confirmed the presence of the known cyclotron resonance scattering feature at around 46.3 keV.

The accurate estimation of photometric redshifts plays a crucial role in accomplishing science objectives of the large survey projects. The template-fitting and machine learning are the two main types of methods applied currently. Based on the training set obtained by cross-correlating the DESI Legacy Imaging Surveys DR9 galaxy catalogue and SDSS DR16 galaxy catalogue, the two kinds of methods are used and optimized, such as EAZY for template-fitting approach and CATBOOST for machine learning. Then the created models are tested by the cross-matched samples of the DESI Legacy Imaging SurveysDR9 galaxy catalogue with LAMOST DR7, GAMA DR3 and WiggleZ galaxy catalogues. Moreover three machine learning methods (CATBOOST, Multi-Layer Perceptron and Random Forest) are compared, CATBOOST shows its superiority for our case. By feature selection and optimization of model parameters, CATBOOST can obtain higher accuracy with optical and infrared photometric information, the best performance ($MSE=0.0032$, $\sigma_{NMAD}=0.0156$ and $O=0.88$ per cent) with $g \le 24.0$, $r \le 23.4$ and $z \le 22.5$ is achieved. But EAZY can provide more accurate photometric redshift estimation for high redshift galaxies, especially beyond the redhisft range of training sample. Finally, we finish the redshift estimation of all DESI DR9 galaxies with CATBOOST and EAZY, which will contribute to the further study of galaxies and their properties.

Exoplanets in protoplanetary disks cause localized deviations from Keplerian velocity in channel maps of molecular line emission. Current methods of characterizing these deviations are time consuming, and there is no unified standard approach. We demonstrate that machine learning can quickly and accurately detect the presence of planets. We train our model on synthetic images generated from simulations and apply it to real observations to identify forming planets in real systems. Machine learning methods, based on computer vision, are not only capable of correctly identifying the presence of one or more planets, but they can also correctly constrain the location of those planets.

Context: NOTT (formerly Hi-5) is a new high-contrast L' band (3.5-4.0 \textmu m) beam combiner for the VLTI with the ambitious goal to be sensitive to young giant exoplanets down to 5 mas separation around nearby stars. The performance of nulling interferometers in these wavelengths is affected both by fundamental noise from the background and by the contributions of instrumental noises. This motivates the development of end-to-end simulations to optimize these instruments. Aims: To enable the performance evaluation and inform the design of such instruments on the current and future infrastructures, taking into account the different sources of noise, and their correlation. Methods: SCIFYsim is an end-to-end simulator for single mode filtered beam combiners, with an emphasis on nulling interferometers. It is used to compute a covariance matrix of the errors. Statistical detection tests based on likelihood ratios are then used to compute compound detection limits for the instrument. Results: With the current assumptions on the performance of the wavefront correction systems, the errors are dominated by correlated instrumental errors down to stars of magnitude 6-7 in the L band, beyond which thermal background from the telescopes and relay system becomes dominant. Conclusions: SCIFYsim is suited to anticipate some of the challenges of design, tuning, operation and signal processing for integrated optics beam combiners. The detection limits found for this early version of NOTT simulation with the unit telescopes are compatible with detections at contrasts up to $10^5$ in the L band at separations of 5 to 80 mas around bright stars.

Solar modelling has long been split into ''internal'' and ''surface'' modelling, because of the lack of tools to connect the very different scales in space and time, as well as the widely different environments and dominating physical effects involved. Significant efforts have recently been put into resolving this disconnect. We address the outstanding bottlenecks in connecting internal convection zone and dynamo simulations to the surface of the Sun, and conduct a proof-of-concept high resolution global simulation of the convection zone of the Sun, using the task-based DISPATCH code framework. We present a new `volleyball' mesh decomposition, which has Cartesian patches tessellated on a sphere with no singularities. We use our new entropy based HLLS approximate Riemann solver to model magneto-hydrodynamics in a global simulation, ranging between 0.655 -- 0.995 R$_\odot$, with an initial ambient magnetic field set to 0.1 Gauss. The simulations develop convective motions with complex, turbulent structures. Small-scale dynamo action twists the ambient magnetic field and locally amplifies magnetic field magnitudes by more than two orders of magnitude within the initial run-time.

We investigate the formation of infant globular cluster (GC) candidates in high-resolution cosmological simulations from the First Billion Years (FiBY) project. By analysing the evolution of the systems in the energy and angular momentum plane, we identify the redshift at which the infant GCs first became gravitationally bound, and we find evidence of radial infall of their gaseous and stellar components. The collapse appears to be driven by internal self-gravity, however, the initial trigger is sourced from the external environment. The phase space behaviour of the infant GCs also allows us to identify some characteristic groupings of objects. Such a classification based on internal properties appears to be reflected in the formation environment: GC candidates that belong to the same class are found in host galaxies of similar morphology, with the majority of the infant GCs located in clumpy, irregular proto-galaxies. Finally, through the inspection of two GC candidates that contain only stars by z = 6, we find that supernova feedback is the main physical mechanism behind their dearth of gas and that the systems subsequently respond with an approximately adiabatic expansion. Such infant GC candidates already resemble the GCs we currently observe in the local Universe.

With the onset of large-scale astronomical surveys capturing millions of images, there is an increasing need to develop fast and accurate deconvolution algorithms that generalize well to different images. A powerful and accessible deconvolution method would allow for the reconstruction of a cleaner estimation of the sky. The deconvolved images would be helpful to perform photometric measurements to help make progress in the fields of galaxy formation and evolution. We propose a new deconvolution method based on the Learnlet transform. Eventually, we investigate and compare the performance of different Unet architectures and Learnlet for image deconvolution in the astrophysical domain by following a two-step approach: a Tikhonov deconvolution with a closed-form solution, followed by post-processing with a neural network. To generate our training dataset, we extract HST cutouts from the CANDELS survey in the F606W filter (V-band) and corrupt these images to simulate their blurred-noisy versions. Our numerical results based on these simulations show a detailed comparison between the considered methods for different noise levels.

While the probability density function (PDF) of the cosmic microwave background (CMB) convergence field approximately follows a Gaussian distribution, small contributions from structures at low redshifts make the overall distribution slightly non-Gaussian. Some of this late-time component can be modelled using the distribution of galaxies and subtracted off from the original CMB lensing map to produce a map of matter distribution at high redshifts. Using this high-redshift mass map, we are able to directly study the early phases of structure formation and look for deviations from our standard model. In this work, we forecast the detectability of signatures of non-Gaussianity due to nonlinear structure formation at $z>1.2$. Although we find that detecting such signatures using ongoing surveys will be challenging, we forecast that future experiments such as the CMB-S4 will be able to make detections of $\sim$ 7$\sigma$.

Recent progress in high-resolution transmission spectroscopy has offered new avenues in which to characterise the atmospheres of transiting exoplanets. High-resolution cross-correlation spectroscopy allows for the unambiguous detection of molecules/atoms. It has also been used to map both atmospheric dynamics and longitudinal variations in the abundance of species across the morning and evening limbs. We present multiple VLT/ESPRESSO observations of the ultra-hot Jupiter WASP-121b, from which we constrain relative abundances of various neutral metals consistently across all observations, whilst accounting for the distortion of the exoplanet's signal caused by traditional data processing techniques. We also constrain planetary orbital velocities and $T$-$P$ profiles. We compare our abundance constraints with previous constraints using VLT/UVES transmission spectroscopy of WASP-121b, and find our results to be consistent between observations, and also in agreement with stellar values for species previously detected in the atmosphere of WASP-121b. Our retrieval framework can also be used to identify potential exospheric species, resulting in extended absorption features beyond the transit equivalent Roche limit of WASP-121b ($R_{\rm eqRL}$ $\sim$ 1.3 $R_{\rm p}$). H$\alpha$, Fe II, and Ca II were found to extend to high altitudes ($1.54\pm0.04$ $R_{\rm p}$, $1.17\pm0.01$ $R_{\rm p}$, and $2.52\pm0.34$ $R_{\rm p}$, respectively), which are broadly consistent with literature values. The consistency of our constraints across multiple high-resolution observations is a strong validation of our model filtering and retrieval framework, as well as the stability of the atmosphere over the timescales of months/years, and could allow for planet formation processes to be inferred from future ground-based observations of exoplanetary atmospheres.

Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $\tau_N$, of the $9.8 \mu$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $\tau_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $\tau_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution.

Evidence for bright radio blazars being high-energy neutrino sources was found in recent years. Specifics of how and where these particles get produced still remain not fully determined. In this paper, we add 14 new IceCube detections from 2020-2022 to update our analysis of the neutrino-blazars connection. We test and refine earlier findings by utilizing the total of 71 track-like high-energy IceCube events from 2009-2022. We correlate them with the complete sample of 3412 extragalactic radio sources selected by their compact radio emission. We demonstrate that neutrinos are statistically associated with radio-bright blazars with the post-trial p-value of 3*10^-4. In addition to this statistical study, we confirm previous individual neutrino-blazar associations, find and discuss several new ones. Notably, PKS 1741-038 was selected earlier and had the second neutrino detected from its direction in 2022; PKS 0735+168 has experienced a major flare across the whole electromagnetic spectrum coincidently with a neutrino arrival from that direction in 2021.

Gamma-ray bursts (GRBs) are bright extragalactic flashes of gamma-ray radiation and briefly the most energetic explosions in the Universe. Their catastrophic origin (the merger of compact objects or the collapse of massive stars) drives the formation of a newborn compact remnant (black hole or magnetar) that powers two highly relativistic jets. To distinguish between magnetized and baryonic jet models and ultimately determine the power source for these energetic explosions, our team studies the polarization of the light during the first minutes after the explosion (using novel instruments on fully autonomous telescopes around the globe) to directly probe the magnetic field properties in these extragalactic jets. This technology allowed the detection of highly polarized optical light in GRB 120308A and confirmed the presence of mildly magnetized jets with large-scale primordial magnetic fields in a reduced sample of GRBs (e.g. GRB 090102, GRB 110205A, GRB 101112A, GRB 160625B). Here we discuss the observations of the most energetic and first GRB detected at very high TeV energies, GRB 190114C, which opens a new frontier in GRB magnetic field studies suggesting that some jets can be launched highly magnetized and that the collapse and destruction of these magnetic fields at very early times may have powered the explosion itself. Additionally, our most recent polarimetric observations of the jet of GRB 141220A indicate that, when the jetted ejected material is decelerated by the surrounding environment, the magnetic field amplification mechanisms at the front shock (needed to generate the observed synchrotron emission) produce small magnetic domains. These measurements validate theoretical expectations and contrast with previous observations that suggest large magnetic domains in collisionless shocks (i.e. GRB 091208B).

We use Functional Renormalisation Group (FRG) techniques to analyse the behaviour of a spectator field, $\sigma$, during inflation that obeys an overdamped Langevin equation. We briefly review how a derivative expansion of the FRG can be used to obtain Effective Equations of Motion (EEOM) for the one- and two-point function and derive the EEOM for the three-point function. We show how to compute quantities like the amplitude of the power spectrum and the spectral tilt from the FRG. We do this explicitly for a potential with multiple barriers and show that in general many different potentials will give identical predictions for the spectral tilt suggesting that observations are agnostic to localised features in the potential. Finally we use the EEOM to compute first-passage time (FPT) quantities for the spectator field. The EEOM for the one- and two-point function are enough to accurately predict the average time taken $\left\langle \mathcal{N}\right\rangle$ to travel between two field values with a barrier in between and the variation in that time $\delta \mathcal{N}^2$. It can also accurately resolve the full PDF for time taken $\rho (\mathcal{N})$, predicting the correct exponential tail. This suggests that an extension of this analysis to the inflaton can correctly capture the exponential tail that is expected in models producing Primordial Black Holes.

Solar wind heating rates have often been calculated by fitting plasma and magnetic field data with a set of model functions. In this letter, we show that the rates obtained by such an approach strongly depend on the rather arbitrary choice one makes for these model functions. An alternative approach, consisting in monitoring the radial evolution of the adiabatic invariants, based on locally and consistently measured plasma and magnetic field parameters, is free from such a flaw. We apply this technique to a recently released Helios proton dataset, and confirm the existence of a clear perpendicular heating of solar wind's protons. On the other hand, no significant change in the parallel adiabatic invariant is visible in the data. We conclude that to date, and in the distance range of 0.3 to 1 AU, no clear observation of a deviation of solar wind's protons from parallel adiabaticity has ever been made.

Euclid is a European Space Agency (ESA) mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey in visible and near-infrared bands, covering approximately 15 000 deg2 of the extragalactic sky in six years. Euclid will be equipped with a 1.2 m diameter Silicon Carbide (SiC) mirror telescope feeding two instruments built by the Euclid Consortium: a high-quality panoramic visible imager and a near-infrared photometer and spectrograph. These proceedings briefly describe the satellite and its instruments, which are optimised for pristine point spread function and reduced stray light, producing very crisp images. Furthermore, we summarise the survey strategy, the global scheduling, and the preparations for the satellite commissioning and the Science Data Centers to produce scientific data.

We present a novel, few-body computational framework designed to shed light on the likelihood of forming intermediate-mass (IM) and supermassive (SM) black holes (BHs) in nuclear star clusters (NSCs) through successive BH mergers, initiated with a single BH seed. Using observationally motivated NSC profiles, we find that the probability of a ${\sim}100 \, M_\odot$ BH to grow beyond ${\sim}1000 \, M_\odot$ through successive mergers ranges from ${\sim}0.1\%$ in low-density, low-mass clusters to nearly $90\%$ in high-mass, high-density clusters. However, in the most massive NSCs, the growth timescale can be very long ($\gtrsim 1\,$Gyr); vice versa, while growth is least likely in less massive NSCs, it is faster there, requiring as little as ${\sim}0.1\,$Gyr. The increased gravitational focusing in systems with lower velocity dispersions is the primary contributor to this behavior. We find that there is a simple "7-strikes-and-you're-in" rule governing the growth of BHs: our results suggest that if the seed survives 7 to 10 successive mergers without being ejected (primarily through gravitational wave recoil kicks), the growing BH will most likely remain in the cluster and will then undergo runaway, continuous growth all the way to the formation of an SMBH (under the simplifying assumption adopted here of a fixed background NSC). Furthermore, we find that rapid mergers enforce a dynamically-mediated "mass gap" between about ${50-300 \, M_\odot}$ in an NSC.

We present and analyze a series of synthetic galaxy survey fields based on the IllustrisTNG Simulation suite. With the Illustris public data release and JupyterLab service, we generated a set of twelve lightcone catalogs covering areas from 5 to 365 square arcminutes, similar to several JWST Cycle 1 programs, including JADES, CEERS, PRIMER, and NGDEEP. From these catalogs, we queried the public API to generate simple mock images in a series of broadband filters used by JWST-NIRCam and the Hubble Space Telescope cameras. This procedure generates wide-area simulated mosaic images that can support investigating the predicted evolution of galaxies alongside real data. Using these mocks, we demonstrate a few simple science cases, including morphological evolution and close pair selection. We publicly release the catalogs and mock images through MAST, along with the code used to generate these projects, so that the astrophysics community can make use of these products in their scientific analyses of JWST deep field observations.

Spatially resolved IFS maps in a sample of $532$ S0 galaxies from the MaNGA survey have unveiled the existence of inner rings ($\langle R\rangle\sim 1\,R_\mathrm{e}$) betraying ongoing star formation in a number of these objects. Activity gradients averaged over bins of galactocentric radius up to $\sim 1.5\,R_\mathrm{e}$ have been measured in the subspace defined by the first two principal components of the optical spectra of these galaxies. We find that the gradients sign is closely related to the presence of such rings in the spectral maps, which are specially conspicuous in the equivalent width of the H$\alpha$ emission line, EW(H$\alpha$), with a fractional abundance $\unicode{x2014}21\unicode{x2013}34\%\unicode{x2014}$ notably larger than that inferred from optical images. While the numbers of S0s with positive, negative, and flat activity gradients are comparable, star-forming rings are largely found in objects for which quenching proceeds from the inside-out, in good agreement with predictions from cosmological simulations studying S0 buildup. Assessment of these ringed structures indicates that their frequency increases with the mass of their hosts, that they have shorter lifetimes in galaxies with ongoing star formation, that they may feed on gas from the disks, and that the local environment does not play a relevant role in their formation. We conclude that the presence of inner rings in the EW(H$\alpha$) is a common phenomenon in fully formed S0s, possibly associated with annular disk resonances driven by weakly disruptive mergers preferentially involving a relatively massive primary galaxy and a tiny satellite strongly bound to the former.

The Earth's ionosphere introduces systematic effects that limit the performance of a radio interferometer at low frequencies ($\lesssim 1$\,GHz). These effects become more pronounced for severe geomagnetic activities or observations involving longer baselines of the interferometer. The uGMRT, a pathfinder for the Square Kilometre Array (SKA), is located in between the northern crest of the Equatorial Ionisation Anomaly (EIA) and the magnetic equator. Hence, this telescope is more prone to severe ionospheric conditions and is a unique radio interferometer for studying the ionosphere. Here, we present 235\,MHz observations with the GMRT, showing significant ionospheric activities over a solar minimum. In this work, we have characterised the ionospheric disturbances observed with the GMRT and compared them with ionospheric studies and observations with other telescopes like the VLA, MWA and LOFAR situated at different magnetic latitudes. We have estimated the ionospheric total electron content (TEC) gradient over the full GMRT array which shows an order of magnitude higher sensitivity compared to the Global Navigation Satellite System (GNSS). Furthermore, this article uses the ionospheric characteristics estimated from the observations with uGMRT, VLA, LOFAR and MWA to forecast the effects on the low-frequency observations with the SKA1-MID and SKA1-LOW in future.

Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfv\'enic, high plasma-beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, versus and direct numerical simulation, DNS). Using an energy transfer analysis framework we calculate the effective, numerical viscosities and resistivities and demonstrate and that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that -- contrary to hydrodynamic turbulence -- the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at $2{,}048^3$ cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it is remarkably successful in reproducing the baryon cycles of the simulated galaxies and their CGM. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$, in agreement with the simulations. We discuss the ways in which our model is complementary to existing approaches for modeling the CGM--galaxy connection and possible future directions.

I have used high-precision photometry and astrometry from the third data release of Gaia (DR3) to perform a survey for members of the Taurus star-forming region and young associations in its vicinity. This work has produced a new catalog of 532 adopted members of Taurus, which has only minor changes relative to the previous catalog from Esplin & Luhman 2019. I have used the Gaia astrometry to divide the Taurus members into 13 groups that have distinct kinematics. Meanwhile, I have identified 1378 candidate members of seven associations near Taurus. All of these associations have histograms of spectral types that peak near M5 (~0.15 Msun), resembling other young populations in the solar neighborhood. For the Taurus groups and neighboring associations, I have estimated ages from their sequences of low-mass stars in Hertzsprung-Russell diagrams. Most of the Taurus groups have median ages of ~1-3 Myr while the associations have ages ranging from 13 to 56 Myr. I have used mid-infrared photometry from the Wide-field Infrared Survey Explorer to search for excess emission from circumstellar disks among the candidate members of the associations. Disks are detected for 51 stars, 20 of which are reported for the first time in this work. Some recent studies have proposed that samples of older stars (>=10 Myr) found in the vicinity of Taurus represent a distributed population that is associated with the Taurus cloud complex. However, I find that most of those stars have kinematics that are inconsistent with any relationship with Taurus.

We consider sampling and detection strategies for solar illuminated space debris. We argue that the lowest detectable debris cross section may be reduced by 10-100x by analysis of phase-space-pixels rather than single frame data. The phase-space-pixel is a weighted stacking of pixels corresponding to a test debris trajectory within the very wide camera field-of-view (FOV). To isolate debris signals from background, exposure time is set to match the time it takes a debris to transit through the instantaneous field of view. Debris signatures are detected though a generalized Hough transform of the data cube. Radiometric analysis of line integrals shows that that sub-cm objects in Low Earth Orbit can be detected and assigned full orbital parameters by this approach

The origin of the so-called Galactic Center Excess in GeV gamma rays has been debated for more than 10 years. What makes this excess so interesting is the possibility of interpreting it as additional radiation consistent with that expected from dark matter annihilation. Alternatively, the excess could come from undetected point sources. In this work, we examine the following questions: Since the majority of the previously reported interpretations of this excess are highly dependent on the simulation, how does the model used for the simulation affect the interpretations? Are such uncertainties taken into account? When different models lead to different conclusions, there may be a general gap between these simulations and reality that influences our conclusions. To investigate these questions, we build an ultra-fast and powerful inference pipeline based on convolutional deep ensemble networks and test the interpretations with a wide range of different models to simulate the excess. We find that our conclusions (dark matter or not) strongly depend on the type of simulation and that this is not revealed by systematic uncertainties. Furthermore, we measure whether reality lies in the simulation parameter space and conclude that there is a gap to reality in all simulated models. Our approach offers a means to assess the severity of the reality gap in future works. Our work questions the validity of conclusions (dark matter) drawn about the GCE in other works: Has the reality gap been closed and at the same time is the model correct?

Monitoring changes in tree cover for rapid assessment of deforestation is considered the critical component of any climate mitigation policy for reducing carbon. Here, we map tropical tree cover and deforestation between 2015 and 2022 using 5 m spatial resolution Planet NICFI satellite images over the state of Mato Grosso (MT) in Brazil and a U-net deep learning model. The tree cover for the state was 556510.8 km$^2$ in 2015 (58.1 % of the MT State) and was reduced to 141598.5 km$^2$ (14.8 % of total area) at the end of 2021. After reaching a minimum deforested area in December 2016 with 6632.05 km$^2$, the bi-annual deforestation area only showed a slight increase between December 2016 and December 2019. A year after, the areas of deforestation almost doubled from 9944.5 km$^2$ in December 2019 to 19817.8 km$^2$ in December 2021. The high-resolution data product showed relatively consistent agreement with the official deforestation map from Brazil (67.2%) but deviated significantly from year of forest cover loss estimates from the Global Forest change (GFC) product, mainly due to large area of fire degradation observed in the GFC data. High-resolution imagery from Planet NICFI associated with deep learning technics can significantly improve mapping deforestation extent in tropics.

One of the unique features of quantum gravity is the lack of local observables and the completeness of boundary observables. We show that the existence of boundary observables in scalar field cosmologies where $a(t)\sim t^{p}$ is equivalent to TCC, which implies $p\leq 1$. Moreover, the mass of weakly-coupled particles must decay like $m\lesssim t^{1-2p}$ to ensure that they yield non-trivial boundary observables. This condition can be expressed in terms of the scalar field that drives the cosmology as $m\lesssim\exp(-c\phi)$ where $c$ depends on the scalar potential. The strongest bound we find is achieved for $V\sim \exp(-2\phi/\sqrt{d-2})$ where $c=1/\sqrt{d-2}$. These results connect some of the most phenomenologically interesting Swampland conjectures to the most basic version of holography.

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above ${\cal O}(10^{-3}\,{\rm eV}$)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which, is to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below ${\cal O}(10^{-5}\,{\rm eV}$)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, leading to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

It is well-known that first order phase transitions in the early universe can be a powerful source of observable stochastic gravitational wave backgrounds. Any such gravitational wave background must exhibit large-scale anisotropies at least as large as those seen in the CMB $\sim 10^{-5}$, providing a valuable new window onto the (inflationary) origins of primordial fluctuations. While significantly larger fractional anisotropies are possible (for example, in multi-field inflation) and would be easier to interpret, it has been argued that these can only be consistent with CMB bounds if the gravitational wave signal is correspondingly smaller. In this paper, we show that this argument, which relies on assuming radiation dominance of the very early universe, can be evaded if there is an era of early matter dominance of a certain robust type. This allows large gravitational wave anisotropies to be consistent with observable signals at proposed future gravitational wave detectors. Constraints from the CMB on large scales, as well as primordial black hole and mini-cluster formation on small scales, and secondary scalar-induced gravitational waves are all taken into account.

The $(n,\alpha)$ reactions play an important role for the energy generation and the synthesis of chemical elements in the stars, as well as for nuclear engineering and medical applications. The aim of this study is to explore the evolution of $(n,\alpha)$ reactions in Fe and Sn isotope chains in order to assess their properties with the increase of neutrons in target nucleus, and compare with other relevant neutron induced reactions. Model calculations of the cross sections are based on the statistical Hauser-Feshbach model in TALYS implementation, using global optical model potential that is additionally adjusted by the $(n,\alpha)$ cross section data for $^{54}$Fe and $^{118}$Sn. The calculations of $(n,\alpha)$ reactions in Fe and Sn isotopes provide the insight into their isospin dependence and properties over the complete relevant range of neutron energies. The results show the evolution of the cross sections with pronounced maxima at low-mass isotopes, and rather strong decrease for neutron-rich nuclei consistent with the reduction of the reaction $Q$-value and increased contributions from other exit channels from compound nucleus. The analysis of the Maxwellian averaged cross sections at temperatures in stellar environment shows that while the $(n,\alpha)$ reactions contribute for the low-mass isotopes, in neutron induced reactions with nuclei with neutron excess, $\gamma$ and neutron emission dominate.

Bayesian inference with nested sampling requires a likelihood-restricted prior sampling method, which draws samples from the prior distribution that exceed a likelihood threshold. For high-dimensional problems, Markov Chain Monte Carlo derivatives have been proposed. We numerically study ten algorithms based on slice sampling, hit-and-run and differential evolution algorithms in ellipsoidal, non-ellipsoidal and non-convex problems from 2 to 100 dimensions. Mixing capabilities are evaluated with the nested sampling shrinkage test. This makes our results valid independent of how heavy-tailed the posteriors are. Given the same number of steps, slice sampling is outperformed by hit-and-run and whitened slice sampling, while whitened hit-and-run does not provide as good results. Proposing along differential vectors of live point pairs also leads to the highest efficiencies, and appears promising for multi-modal problems. The tested proposals are implemented in the UltraNest nested sampling package, enabling efficient low and high-dimensional inference of a large class of practical inference problems relevant to astronomy, cosmology, particle physics and astronomy.

We study the homogeneous and anisotropic dynamics of pseudoscalar inflation coupled to an SU($N$) gauge field. To see how the initially anisotropic universe is isotropized in such an inflation model, we derive the equations to obtain axisymmetric SU($N$) gauge field configurations in Bianchi type-I geometry and discuss a method to identify their isotropic subsets which are the candidates of their late-time attractor. Each isotropic solution is characterized by the corresponding SU(2) subalgebra of the SU($N$) algebra. It is shown numerically that the isotropic universe is a universal late-time attractor in the case of the SU(3) gauge field. Interestingly, we find that a transition between the two distinct gauge-field configurations characterized by different SU(2) subalgebras can occur during inflation. We clarify the conditions for this to occur. This transition could leave an observable imprint on the CMB and the primordial gravitational wave background.

Operating satellites in Very Low Earth Orbit (VLEO) benefits the already expanding New Space industry in applications including Earth Observation and beyond. However, long-term operations at such low altitudes require propulsion systems which compensate the large atmospheric drag forces, so that, when using conventional propulsion systems, the amount of storable propellant limits the maximum mission lifetime. The latter can be avoided by employing Atmosphere-Breathing Electric Propulsion (ABEP) system, which collects the residual atmospheric particles and use them as propellant for an electric thruster. Thus, the requirement of on-board propellant storage can ideally be nullified. At the Institute of Space Systems (IRS) of the University of Stuttgart, an intake and a RF Helicon-based Plasma Thruster (IPT) for ABEP system are developed within the Horizons 2020 funded DISCOVERER project. In order to assess possible future use cases, this paper proposes and analyzes several novel ABEP based mission scenarios. Beginning with technology demonstration mission in VLEO, more complex mission scenarios are derived and discussed in detail. These include, amongst others, orbit maintenance around Mars as well as refuelling and space tug missions. The results show that the potential of ABEP system can exceed drag compensation missions and show that a multitude of different future mission applications exist.

The structure and evolution of Primordial Antimatter domains and Dark matter objects are analysed. Relativistic low-density antimatter domains are described. The Relativistic FRW perfect-fluid solution is found for the characterization of i) ultra-high density antimatter domains, ii) high-density antimatter domains, and iii) dense anti-matter domains. The possible sub-domains structures is analyzed. The structures evolved to the time of galaxy formation are outlined. Comparison is given with other primordial celestial objects. The features of antistars are outlined. In the case of WIMP dark matter clumps, the mechanisms of their survival to the present time are discussed. The cosmological features of neutrino clumping due to fifth force are examined.

In a previous paper we have shown that superluminal particles are allowed by the general relativistic theory of gravity provided that the metric is locally Euclidean. Here we calculate the probability density function of a canonical ensemble of superluminal particles as function of temperature. This is done for both space-times invariant under Lorentz symmetry group, and for space times invariant under an Euclidean symmetry group. Although only the Lorentzian metric is stable for normal matter density, an Euclidian metric can be created under special gravitational circumstances and persist in a limited region of space-time consisting of the very early universe which is characterized by extremely high densities and temperatures. Superluminal particles also allow attaining thermodynamic equilibrium at a shorter duration and also suggest a rapid expansion of the matter density, thus making mechanism such as inflation (which demands invoking and ad-hoc scalar field) redundant. This is in accordance with Occam's razor.

We argue that there is a lower bound of order $10^{-18}$ eV on dark matter mass if it is produced after inflation via a process with finite correlation length.