Papers for Wednesday, Sep 06 2023

Numerical relativity simulations are crucial for studying black holes and have been instrumental in the detection of gravitational waves by the LVK. However, these simulations produce vast amounts of data that must be processed in order to perform studies, create models, and use them with gravitational wave detection pipelines. This paper introduces mayawaves, an open-source python library for processing, studying, and exporting numerical relativity simulations performed using the Einstein Toolkit and MAYA. Mayawaves streamlines the process of analyzing simulations with an intuitive interface, greatly reducing the learning curve for numerical relativity.

This paper describes the overview of the FLASHING (Finest Legacy Acquisitions of SiO-/ H$_2$O-maser Ignitions by the Nobeyama Generation) project promoted using the 45 m telescope of Nobeyama Radio Observatory, which aims to intensively monitor H$_2$O (22 GHz) and SiO (43 GHz) masers associated with so-called "water fountain" sources. Here we show scientific results on the basis of the data taken in for the first five seasons of FLASHING, from 2018 December to 2023 April). We have found the evolution of the H$_2$O maser spectra, such as new spectral components breaking the record of the jet's top speed and/or systematic velocity drifts in the spectrum indicating acceleration or deceleration of the maser gas clumps. For the 43 GHz SiO maser emission, we have found its new detection in a source while its permanent disappearance in other source. Our finding may imply that the jets from these water fountains can be accelerated or decelerated, and show how cicumstellar envelopes have been destroyed.

The equations of motion for a Lagrangian mainly refer to the acceleration equations, which can be obtained by the Euler--Lagrange equations. In the post-Newtonian Lagrangian form of general relativity, the Lagrangian systems can only maintain a certain post-Newtonian order and are incoherent Lagrangians since the higher-order terms are omitted. This truncation can cause some changes in the constant of motion. However, in celestial mechanics, Hamiltonians are more commonly used than Lagrangians. The conversion from Lagrangian to Hamiltonian can be achieved through the Legendre transformation. The coordinate momentum separable Hamiltonian can be computed by the symplectic algorithm, whereas the inseparable Hamiltonian can be used to compute the evolution of motion by the phase-space expansion method. Our recent work involves the design of a multi-factor correction map for the phase-space expansion method, known as the correction map method. In this paper, we compare the performance of the implicit algorithm in post-Newtonian Lagrangians and the correction map method in post-Newtonian Hamiltonians. Specifically, we investigate the extent to which both methods can uphold invariance of the motion's constants, such as energy conservation and angular momentum preservation. Ultimately, the results of numerical simulations demonstrate the superior performance of the correction map method, particularly with respect to angular momentum conservation.

The recent release of the final, complete survey of Lyman-alpha baryonic acoustic oscillation measurements provides the most significant and accurate data base for studying cosmic geometry at an effective redshift z_eff=2.334, which is inaccessible to other sources. In this Letter, we use these data to select among four distinct cosmologies: Planck LCDM, the R_h=ct universe, the Milne universe and Einstein-de Sitter. Given the breadth and depth of the Lyman-alpha study, this BAO measurement alone provides a strong model comparison, complementary to previous studies that combined Lyman-$\alpha$ data with measurements at lower redshifts. Though both approaches are useful, the latter tends to dilute the disparity between model predictions and the observations. We therefore examine how the models compare to each other strictly based on the BAO scale measured in the Lyman-alpha forest and background quasars. We find that Milne and Einstein-de Sitter are strongly ruled out by these data. There is also strong evidence disfavoring the standard model. The Lyman-alpha measurements are completely consistent with the cosmic geometry predicted by R_h=ct. As such, evidence continues to grow that the zero active mass condition from general relativity ought to be an essential ingredient in LCDM.

We propose a self-consistent explanation of Rieger-type periodicities, the Schwabe cycle, and the Suess-de Vries cycle in terms of resonances of various wave phenomena with gravitational forces exerted by the orbiting planets. Starting on the high-frequency side, we show that the two-planet spring tides of Venus, Earth and Jupiter are able to excite magneto-Rossby waves which can be linked with typical Rieger-type periods. We argue then that the 11.07-year beat period of those magneto-Rossby waves synchronizes an underlying conventional $\alpha-\Omega$-dynamo, by periodically changing either the field storage capacity in the tachocline or some portion of the $\alpha$-effect therein. We also strengthen the argument that the Suess-de Vries cycle appears as an 193-year beat period between the 22.14-year Hale cycle and a spin-orbit coupling effect related with the 19.86-year rosette-like motion of the Sun around the barycenter.

One of the most important questions in astrophysics is what causes galaxies to stop forming stars. Previous studies have shown a tight link between quiescence and black hole mass. Other studies have revealed that quiescence is also associated with 'starvation', the halting of gas inflows, which results in the remaining gas being used up rapidly by star formation and in rapid chemical enrichment. In this work we find the final missing link between these two findings. Using a large sample of galaxies, we uncover the intrinsic dependencies of the stellar metallicity on galaxy properties. In the case of the star-forming galaxies, the stellar metallicity is driven by stellar mass. However, for passive galaxies the stellar metallicity is primarily driven by the black hole mass, as traced by velocity dispersion. This result finally reveals the connection between previous studies, where the integrated effect of black hole feedback prevents gas inflows, starving the galaxy, which is seen by the rapid increase in the stellar metallicity, leading to the galaxy becoming passive.

Using the full six years of imaging data from the Dark Energy Survey, we study the surface brightness profiles of galaxy cluster central galaxies and intra-cluster light. We apply a ``stacking'' method to over four thousand galaxy clusters identified by the redMaPPer cluster finding algorithm in the redshift range of 0.2 to 0.5. This yields high signal-to-noise radial profile measurements of the central galaxy and intra-cluster light out to 1 Mpc from the cluster center. Using redMaPPer richness as a cluster mass indicator, we find that the intra-cluster light brightness has a strong mass dependence throughout the 0.2 to 0.5 redshift range, and the dependence grows stronger at a larger radius. In terms of redshift evolution, we find some evidence that the central galaxy, as well as the diffuse light within the transition region between the cluster central galaxy and intra-cluster light within 80 kpc from the center, may be growing over time. At larger radii, more than 80 kpc away from the cluster center, we do not find evidence of additional redshift evolution beyond the cluster mass dependence, which is consistent with the findings from the IllustrisTNG hydrodynamic simulation. We speculate that the major driver of intra-cluster light growth, especially at large radii, is associated with cluster mass growth. Finally, we find that the color of the cluster central galaxy and intra-cluster light displays a radial gradient that becomes bluer at a larger radius, which is consistent with a stellar stripping and disruption origin of intra-cluster light as suggested by simulation studies.

We explore the implications of the light curve of the early TeV gamma-ray afterglow of GRB221009A reported by the LHAASO collaboration. We show that the reported offset of the reference time, $T_*$, allows the determination of the relativistic jet activation time, which occurs approximately $200\,\mathrm{s}$ after the GBM trigger time and closely precedes the moment at which GBM was saturated. We find that while the LHAASO data do not exclude the homogeneous circumburst medium scenario, the progenitor wind scenario looks preferable, finding excellent agreement with the expected size of the stellar bubble. We conclude that the initial growth of the light curve is dominated by processes internal to the jet or by gamma-gamma attenuation on the photons emitted during the prompt phase. Namely, either the activation of the acceleration process or the decrease of internal gamma-gamma absorption can naturally explain the initial rapid flux increase. The subsequent slow flux growth phase observed up to $T_*+18\,\mathrm{s}$ is explained by the build-up of the synchrotron radiation -- the target for inverse Compton scattering, which is also supported by a softer TeV spectrum measured during this period. The duration of this phase allows an almost parameter-independent determination of the jet's initial Lorentz factor, $\Gamma_0\approx600$, and magnetic field strength, $B'\sim0.3\,\mathrm{G}$. These values appear to match well those previously revealed through spectral modeling of the GRB emission.

Measurement of the largest angular scale ($\ell < 30$) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization, $\tau$, and search for the signature of inflation through the detection of primordial $B$-modes. We present an analysis of maps covering nearly 75% of the sky made from the ground-based $40\,\mathrm{GHz}$ channel of the Cosmology Large Angular Scale Surveyor (CLASS) from August 2016 to May 2022. Using fast front-end polarization modulation from the Atacama Desert in Chile, we show this channel achieves higher sensitivity than the analogous frequencies from satellite measurements in the range $10 < \ell < 100$. After a final calibration adjustment, noise simulations show the CLASS linear (circular) polarization maps have a white noise level of $125 \,(130)\,\mathrm{\mu K\, arcmin}$. We measure the Galaxy-masked $EE$ and $BB$ spectra of diffuse synchrotron radiation and compare to space-based measurements at similar frequencies. In combination with external data, we expand measurements of the spatial variations of the synchrotron spectral energy density (SED) to include new regions of the sky and measure the faint diffuse SED in the harmonic domain. We place a new upper limit on a background of circular polarization in the range $5 < \ell < 125$ with the first bin showing $D_\ell < 0.023$ $\mathrm{\mu K^2_{CMB}}$ at 95% confidence. These results establish a new standard for recovery of the largest-scale CMB polarization from the ground and signal exciting possibilities when the higher sensitivity and higher frequency CLASS channels are included in the analysis.

Theory of the physics of the early hot universe leads to a prediction of baryon acoustic oscillations that has received confirmation from the pair-wise separations of galaxies in samples of hundreds of thousands of objects. Evidence is presented here for the discovery of a remarkably strong individual contribution to the baryon acoustic oscillation (BAO) signal at z=0.068, an entity that is given the name Ho'oleilana. The radius of the 3D structure is 155/h_{75} Mpc. At its core is the Bootes supercluster. The Sloan Great Wall, CfA Great Wall, and Hercules complex all lie within the BAO shell. The interpretation of Ho'oleilana as a BAO structure with our preferred analysis implies a value of the Hubble constant of 76.9+8.2-4.8 km/s/Mpc.

We have selected 337 intermediate and high-mass YSOs ($1.5$ to $20$ M$_{\odot}$) well-characterised with spectroscopy. By means of the clustering algorithm HDBSCAN, we study their clustering and association properties in the Gaia DR3 catalogue as a function of stellar mass. We find that the lower mass YSOs ($1.5-4$ M$_{\odot}$) have clustering rates of $55-60\%$ in Gaia astrometric space, a percentage similar to the one found in the T Tauri regime. However, intermediate-mass YSOs in the range $4-10$ M$_{\odot}$ show a decreasing clustering rate with stellar mass, down to $27\%$. We find tentative evidence suggesting that massive YSOs ($>10$ M$_{\odot}$) often appear $-$yet not always$-$ clustered. We put forward the idea that most massive YSOs form via a mechanism that demands many low-mass stars around them. However, intermediate-mass YSOs form in a classical core-collapse T Tauri way, yet they do not appear often in the clusters around massive YSOs. We also find that intermediate and high-mass YSOs become less clustered with decreasing disk emission and accretion rate. This points towards an evolution with time. For those sources that appear clustered, no major correlation is found between their stellar properties and the cluster sizes, number of cluster members, cluster densities, or distance to cluster centres. In doing this analysis, we report the identification of 55 new clusters. We present tabulated all the derived cluster parameters for the considered intermediate and high-mass YSOs.

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA.

The estimates of the delivery of icy planetesimals from the feeding zone of Proxima Centauri c (with mass equal to 7mE, mE is the mass of the Earth) to inner planets b and d were made. They included the studies of the total mass of planetesimals in the feeding zone of planet c and the probabilities of collisions of such planetesimals with inner planets. This total mass could be about 10-15mE. It was estimated based on studies of the ratio of the mass of planetesimals ejected into hyperbolic orbits to the mass of planetesimals collided with forming planet c. At integration of the motion of planetesimals, the gravitational influence of planets c and b and the star was taken into account. In most series of calculations, planetesimals collided with planets were excluded from integrations. Based on estimates of the mass of planetesimals ejected into hyperbolic orbits, it was concluded that during the growth of the mass of planet c the semi-major axis of its orbit could decrease by at least a factor of 1.5. Depending on possible gravitational scattering due to mutual encounters of planetesimals, the total mass of material delivered by planetesimals from the feeding zone of planet c to planet b was estimated to be between 0.002mE and 0.015mE. Probably, the amount of water delivered to Proxima Centauri b exceeded the mass of water in Earth's oceans. The amount of material delivered to planet d could be a little less than that delivered to planet b.

We present a comprehensive set of forecasts for the cross-correlation signal between 21cm intensity mapping and galaxy redshift surveys. We focus on the data sets that will be provided by the SKAO for the 21cm signal, DESI and Euclid for galaxy clustering. We build a likelihood which takes into account the effect of the beam for the radio observations, the Alcock-Paczynski effect, a simple parameterization of astrophysical nuisances, and fully exploit the tomographic power of such observations in the range $z=0.7-1.8$ at linear and mildly non-linear scales ($k<0.25 h/$Mpc). The forecasted constraints, obtained with Monte Carlo Markov Chains techniques in a Bayesian framework, in terms of the six base parameters of the standard $\Lambda$CDM model, are promising. The predicted signal-to-noise ratio for the cross-correlation can reach $\sim 50$ for $z\sim 1$ and $k\sim 0.1 h/$ Mpc. When the cross-correlation signal is combined with current Cosmic Microwave Background (CMB) data from Planck, the error bar on $\Omega_{\rm c}\,h^2$ and $H_0$ is reduced by a factor 3 and 6, respectively, compared to CMB only data, due to the measurement of matter clustering provided by the two observables. The cross-correlation signal has a constraining power that is comparable to the auto-correlation one and combining all the clustering measurements a sub-percent error bar of 0.33% on $H_0$ can be achieved, which is about a factor 2 better than CMB only measurement. Finally, as a proof-of-concept, we test the full pipeline on the real data measured by the MeerKat collaboration (Cunnington et al. 2022) presenting some (weak) constraints on cosmological parameters.

The Wide Field Survey Telescope (WFST) is a dedicated time-domain multi-band ($u$, $g$, $r$, $i$, and $z$) photometric survey facility under construction. In this paper, we present a preliminary study that assesses the quality of photometric redshifts based on WFST by utilizing mock observations derived with the galaxy catalog in the COSMOS/UltraVISTA field. We apply the template fitting technique to estimate photometric redshifts by using the ZEBRA photometric-redshift code and adopting a modified set of adaptive templates. We evaluate the bias (median relative offset between the output photometric redshifts and input redshifts), normalized median absolute deviation ($\sigma_{\rm NMAD}$) and outlier fraction ($f_{\rm outlier}$) of photometric redshifts in two typical WFST observational cases, the single 30-second exposure observations (hereafter shallow mode) and co-added 50-minute exposure observations (hereafter deep mode). We find bias$\la0.006$, $\sigma_{\rm NMAD}\la0.03$, and $f_{\rm outlier}\la5\%$ in the shallow mode and bias$\approx 0.005$, $\sigma_{\rm NMAD}\approx 0.06$, and $f_{\rm outlier}\approx 17\%$--$27\%$ in the deep mode, respectively, under various lunar phases. Combining the WFST mock observational data with that from the upcoming CSST and Euclid surveys, we demonstrate that the $z_{\rm phot}$ results can be significantly improved, with $f_{\rm outlier}\approx 1\%$ and $\sigma_{\rm NMAD}\approx 0.02$.

We discuss problems of planetesimal migration in the emerging Solar System and exoplanetary systems. Protoplanetary disk evolution models and the formation of planets are considered. The formation of the Moon and of the asteroid and trans-Neptunian belts is studied. We show that Earth and Venus could acquire more than half of their mass in 5 million years, and their outer layers could accumulate the same material from different parts of the feeding zone of these planets. The migration of small bodies toward the terrestrial planets from various regions of the Solar System is simulated numerically. Based on these computations, we conclude that the mass of water delivered to the Earth by planetesimals, comets, and carbonaceous chondrite asteroids from beyond the ice line could be comparable to the mass of Earth's oceans. The processes of dust migration in the Solar System and sources of the zodiacal cloud are considered.

AXIS is a Probe-class mission concept that will provide high-throughput, high-spatial-resolution X-ray spectral imaging, enabling transformative studies of high-energy astrophysical phenomena. To take advantage of the advanced optics and avoid photon pile-up, the AXIS focal plane requires detectors with readout rates at least 20 times faster than previous soft X-ray imaging spectrometers flying aboard missions such as Chandra and Suzaku, while retaining the low noise, excellent spectral performance, and low power requirements of those instruments. We present the design of the AXIS high-speed X-ray camera, which baselines large-format MIT Lincoln Laboratory CCDs employing low-noise pJFET output amplifiers and a single-layer polysilicon gate structure that allows fast, low-power clocking. These detectors are combined with an integrated high-speed, low-noise ASIC readout chip from Stanford University that provides better performance than conventional discrete solutions at a fraction of their power consumption and footprint. Our complementary front-end electronics concept employs state of the art digital video waveform capture and advanced signal processing to deliver low noise at high speed. We review the current performance of this technology, highlighting recent improvements on prototype devices that achieve excellent noise characteristics at the required readout rate. We present measurements of the CCD spectral response across the AXIS energy band, augmenting lab measurements with detector simulations that help us understand sources of charge loss and evaluate the quality of the CCD backside passivation technique. We show that our technology is on a path that will meet our requirements and enable AXIS to achieve world-class science.

The presence of a population of a large number ($\sim$400) of almost coeval (100--300 Myr) super star clusters (SSCs) in the disk of M82 offers an opportunity to construct the Cluster Initial Mass Function (CIMF) from the observed present-day Cluster Mass Function (CMF). We carry out the dynamical and photometric evolution of the CMF assuming the clusters move in circular orbits under the gravitational potential of the host galaxy using the semi-analytical simulation code EMACSS. We explore power-law and log-normal functions for the CIMFs, and populate the clusters in the disk assuming uniform, power-law, and exponential radial distribution functions. We find that the observed CMF is best produced by a CIMF that is power-law in form with an index of 1.8, for a power-law radial distribution function. More importantly, we establish that the observed turn-over in the present-day CMF is the result of observational incompleteness rather than due to dynamically induced effects, or an intrinsically log-normal CIMF, as was proposed for the fossil starburst region B of this galaxy. Our simulations naturally reproduce the mass-radius relation observed for a sub-sample of M82 SSCs.

Looking to the future of exo-Earth imaging from the ground, core technology developments are required in visible extreme adaptive optics (ExAO) to enable the observation of atmospheric features such as oxygen on rocky planets in visible light. UNDERGROUND (Ultra-fast AO techNology Determination for Exoplanet imageRs from the GROUND), a collaboration built in Feb. 2023 at the Optimal Exoplanet Imagers Lorentz Workshop, aims to (1) motivate oxygen detection in Proxima Centauri b and analogs as an informative science case for high-contrast imaging and direct spectroscopy, (2) overview the state of the field with respect to visible exoplanet imagers, and (3) set the instrumental requirements to achieve this goal and identify what key technologies require further development.

As we look to the next generation of adaptive optics systems, now is the time to develop and explore the technologies that will allow us to image rocky Earth-like planets; wavefront control algorithms are not only a crucial component of these systems, but can benefit our adaptive optics systems without requiring increased detector speed and sensitivity or more effective and efficient deformable mirrors. To date, most observatories run the workhorse of their wavefront control as a classic integral controller, which estimates a correction from wavefront sensor residuals, and attempts to apply that correction as fast as possible in closed-loop. An integrator of this nature fails to address temporal lag errors that evolve over scales faster than the correction time, as well as vibrations or dynamic errors within the system that are not encapsulated in the wavefront sensor residuals; these errors impact high contrast imaging systems with complex coronagraphs. With the rise in popularity of machine learning, many are investigating applying modern machine learning methods to wavefront control. Furthermore, many linear implementations of machine learning methods (under varying aliases) have been in development for wavefront control for the last 30-odd years. With this work we define machine learning in its simplest terms, explore the most common machine learning methods applied in the context of this problem, and present a review of the literature concerning novel machine learning approaches to wavefront control.

It has long been claimed that novae reaching the highest luminosity at the peak of their eruptions appear to fade the fastest from maximum light. The relationship between peak brightness and fade rate is known as the Maximum-Magnitude, Rate-of-Decline (MMRD) relation. Lightcurve parameters for the most recent sample of M31 recurrent novae are presented and used to buttress the case that the observed MMRD relation can be explained as a consequence of observational selection effects coupled with expectations from standard nova models.

The surface component of the IceCube Neutrino Observatory, IceTop, consists of an array of ice-Cherenkov tanks measuring the electromagnetic signal as well as low-energy ($\sim\rm{GeV}$) muons from cosmic-ray air showers. In addition, accompanying high-energy (above a few $100\,\rm{GeV}$) muons can be observed in coincidence in the deep in-ice detector. A combined measurement of the low- and high-energy muon content is of particular interest for tests of hadronic interaction models as well as for cosmic-ray mass discrimination. However, since IceTop does not feature dedicated muon detectors, an estimation of the low-energy muon component of individual air showers is challenging. In this work, a two-component lateral distribution function (LDF), using separate descriptions for the electromagnetic and muon lateral distributions of the detector signals, is introduced as a new approach for the estimation of low-energy muons in air showers on an event-by-event basis. The principle of the air-shower reconstruction using the two-component LDF, as well as its reconstruction performance with respect to primary energy and number of low-energy muons will be discussed.

A search for pulse signals in a area with declinations of $+52\degr <\delta <+55\degr$ was carried out on the LPA LPI radio telescope. When processing ten months of observations recorded in six frequency channels with a channel width of 415 kHz and a total bandwidth of 2.5 MHz, 22 thousand events were found with a pronounced dispersion delay of signals over frequency channels, i.e. having signs of pulsar pulses. It turned out that the found pulses belong to four known pulsars and two new rotating radio transients (RRATs). An additional pulse search conducted in 32-channel data with a channel width of 78 kHz revealed 8 pulses for the transient J0249+52 and 7 pulses for the transient J0744+55. Periodic radiation of transients was not detected. The analysis of observations shows that the found RRATs are most likely pulsars with nullings, where the proportion of nulling is greater than 99.9\%.

In the manuscript, we check properties of electron densities $n_e$ traced by flux ratio $R_{sii}$ of [S~{\sc ii}]$\lambda6716$\AA~ to [S~{\sc ii}]$\lambda6731$\AA~ in narrow emission line regions (NLRs) between Type-1 AGN and Type-2 AGN in SDSS DR12. Under the framework of Unified Model considering kpc-scale structures, similar $n_e$ in NLRs should be expected between Type-1 AGN and Type-2 AGN. Based on reliable measurements of [S~{\sc ii}] doublet with measured parameters at least five times larger than corresponding uncertainties, there are 6039 Type-1 AGN and 8725 Type-2 AGN (excluding the Type-2 LINERs and the composite galaxies) collected from SDSS DR12. Then, lower $R_{sii}$ (higher $n_e$) in NLRs can be well confirmed in Type-1 AGN than in Type-2 AGN, with confidence level higher than 5$\sigma$, even after considering necessary effects including effects of electron temperatures traced by [O~{\sc iii}]$\lambda4364,4959,5007$\AA~ on estimating $n_e$ in NLRs. Two probable methods are proposed to explain the higher $n_e$ in NLRs in Type-1 AGN. First, the higher $n_e$ in NLRs of Type-1 AGN could indicate longer time durations of AGN activities in Type-1 AGN than in Type-2 AGN, if AGN activities triggering galactic-scale outflows leading to more electrons injecting into NLRs were accepted to explain the higher $n_e$ in NLRs of Type-2 AGN than HII galaxies. Second, the lower $n_e$ in NLRs of Type-2 AGN could be explained by stronger star-forming contributions in Type-2 AGN, considering lower $n_e$ in HII regions. The results provide interesting challenges to the commonly and widely accepted Unified Model of AGN.

The electron cyclotron maser instability (ECMI) of extraordinary mode waves was investigated with the parameters observed in Saturn's kilometric radiation (SKR) sources. Previous studies employed simplified dispersion relations, and did not consider the excitation of the relativistic (R) mode. This mode is introduced by considering the relativistic effect in plasmas consisting of both cold and hot electrons. Using particle-in-cell simulations, we investigated the excitation of R and X modes based on the measured data. Using the reported value of the density ratio of energetic to total electrons $n_e/n_0=24\%$, the most unstable mode is the R mode. The escaping X-mode emissions are amplified only if the energetic electrons are dominant with $n_e/n_0 \ge 90\%$. For these cases, only the X mode is excited and the R mode disappears due to its strong coupling. The results are well in line with the linear kinetic theory of ECMI. The properties of both the R and X modes are consistent with the observed SKR emissions. This raises questions about the nature of the measured electric field fluctuations within ``presumed'' SKR sources. The study provides new insights into the ECMI process relevant to SKR emission mechanisms.

Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450{\mu}$m-selected galaxies from deep JCMT surveys will also be strongly lensed. Follow-up ALMA campaigns have so far found one potential lens candidate, but without clear compelling evidence e.g. from lensing arcs. Here we report the discovery of a compelling gravitational lens system confirming the lensing population predictions, with a $z_{s} = 3.4 {\pm} 0.4$ submm source lensed by a $z_{spec} = 0.360$ foreground galaxy within the COSMOS field, identified through public JWST imaging of a $450{\mu}$m source in the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) catalogue. These systems will typically be well within the detectable range of future wide-field surveys such as Euclid and Roman, and since sub-millimetre galaxies are predominantly very red at optical/near-infrared wavelengths, they will tend to appear in near-infrared channels only. Extrapolating to the Euclid-Wide survey, we predict tens of thousands of strongly lensed near-infrared galaxies. This will be transformative for the study of dusty star-forming galaxies at cosmic noon, but will be a contaminant population in searches for strongly lensed ultra-high-redshift galaxies in Euclid and Roman.

The merging of a binary system involving two neutron stars (NSs), or a black hole (BH) and a NS, often results in the emission of an electromagnetic (EM) transient. One component of this EM transient is the epic explosion known as a kilonova (KN). The characteristics of the KN emission can be used to probe the equation of state (EoS) of NS matter responsible for its formation. We predict KN light curves from computationally simulated BH-NS mergers, by using the 3D radiative transfer code \texttt{POSSIS}. We investigate two EoSs spanning most of the allowed range of the mass-radius diagram. We also consider a soft EoS compatible with the observational data within the so-called 2-families scenario in which hadronic stars coexist with strange stars. Computed results show that the 2-families scenario, characterized by a soft EoS, should not produce a KN unless the mass of the binary components are small ($M_{\rm BH} \leq 6M_{\odot}$, $M_{\rm NS} \leq 1.4M_{\odot}$) and the BH is rapidly spinning ($\chi_{\rm BH} \geq 0.3$). In contrast, a strong KN signal potentially observable from future surveys (e.g. VRO/LSST) is produced in the 1-family scenario for a wider region of the parameter space, and even for non-rotating BHs ($\chi_{\rm BH} = 0$) when $M_{\rm BH} = 4M_{\odot}$ and $M_{\rm NS} = 1.2M_{\odot}$. We also provide a fit that allows for the calculation of the unbound mass from the observed KN magnitude, without running timely and costly radiative transfer simulations. Findings presented in this paper will be used to interpret light curves anticipated during the fourth observing run (O4), of the advanced LIGO, advanced Virgo and KAGRA interferometers and thus to constrain the EoS of NS matter.

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at non-zero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve proto-neutron stars born in supernova explosions, binary neutron star mergers and accreting neutron stars. The temperature-density range studied covers the transition from a non-degenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (a) the dynamical screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (b) the correlations of the ionic component in a one-component plasma, and (c) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropies introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviation of relaxation time and components of the conductivity tensor from the average are below $\le 10\%$ except close to crust-core transition, where non-spherical nuclear structures are expected. Our results can be used in dissipative magneto-hydrodynamics (MHD) simulations of warm compact stars.

We conducted research on the classification and physical properties of 10 objects from the HASH (Hong Kong/Australian Astronomical Observatory/Strasbourg Observatory H-alpha Planetary Nebula (PN)) database with small angular sizes (< 8\arcsec) in the northern hemisphere. The sample consisted of 6 Likely PNe, 2 new candidates, one emission-line star, and one object of unknown nature. Among them, we observed 4 objects for the first time using the medium-resolution TFOSC spectrograph located on the RTT150 cm of the T\"UB\.ITAK National Observatory (TUG). To investigate the classification of the observed objects, we utilized the emission line ratios of [O III]/H$_{\gamma}$, [O III]/H$_{\beta}$, [N II]/H$_{\alpha}$ and [S II]/H$_{\alpha}$ and diagnostic diagrams such as the Sabbadin-Minello-Bianchini (SMB) and Baldwin-Phillips-Terlevich (BPT). When considering a broader range of diagnostic criteria compared to those provided in the literature, our analyses resulted in the reclassification of 4 objects from Likely PNe to True PNe and the retention of the previous classification for the remaining 6 objects. In addition, we obtained various physical conditions such as electron temperatures, electron densities, logarithmic extinction coefficients, and excitation classes for the 10 objects under study. Our analysis revealed that the ionic abundances of the majority of these objects were in agreement with Galactic PNe. Our spectral observations have led to the updating of 10 PNe in the HASH database.

Spatially compact objects with extremely red color in the rest-frame optical to near-infrared (0.4--3.0 ${\rm \mu m}$) and blue color in the rest-frame ultraviolet (UV; 0.2--0.4 ${\rm \mu m}$) have been discovered at $5 < z < 9$ using the James Webb Space Telescope (JWST). These extremely red objects (JWST-EROs) exhibit spectral energy distributions (SEDs) that are difficult to explain using a single component of either star-forming galaxies or quasars, leading to two-component models in which the blue UV and extremely red optical are explained using less-dusty and dusty spectra of galaxies or quasars, respectively. Here, we report the remarkable similarity in SEDs between JWST-EROs and blue-excess dust-obscured galaxies (BluDOGs) identified at $2 < z < 3$. BluDOGs are a population of active galactic nuclei (AGNs) with blackhole masses of $\sim10^{8-9}$ M$_\odot$, which are one order of magnitude larger than those in some JWST-EROs. The Eddington accretion rates of BluDOGs are one or higher, whereas those of JWST-EROs are in the range of 0.1--1. Therefore, JWST-EROs are less massive, less active, and more common counterparts in higher-$z$ of BluDOGs in cosmic noon. Conversely, JWST-EROs have a significantly higher fraction of those with blue-excess than DOGs. We present the average UV spectra of BluDOGs as a comparison to JWST-EROs and discuss a coherent evolutionary scenario for dusty AGN populations.

Multi-Messenger Astrophysics (MMA) has produced a wealth of data with much more to come in the future. This enormous data set will reveal new insights into the physics of Core Collapse SuperNovae (CCSN), Binary Neutron Star Mergers (BNSM), and many other objects where it is actually possible, if not probable, that new physics is in operation. To tease out different possibilities, we will need to analyze signals from photons, neutrinos, gravitational waves, and chemical elements. This task is made all the more difficult when it is necessary to evolve the neutrino component of the radiation field and associated quantum-mechanical property of flavor in order to model the astrophysical system of interest -- a numerical challenge that has not been addressed to this day. In this work, we take a step in this direction by adopting the technique of angular-integrated moments with a truncated tower of dynamical equations and a closure, convolving a flavor-transformation with spatial transport to evolve the neutrino radiation quantum field. We show that moments capture the dynamical features of Fast Flavor Instabilities (FFI) and provide comparable results to a more precise particle-in-cell method. We propose areas for improvement in the future.

Deploying \textit{multiple sharp transitions} (MSTs) under a unified framework, we investigate the formation of Primordial Black Holes (PBHs) and the production of Scalar Induced Gravitational Waves (SIGWs) by incorporating one-loop corrected renormalized-resummed scalar power spectrum. With effective sound speed parameter, $1 \leq c_s \leq 1.17$, the direct consequence is the generation of PBH masses spanning $M_{\rm PBH}\sim{\cal O}(10^{-31}M_{\odot}- 10^{4}M_{\odot})$, thus evading well known \textit{No-go theorem} on PBH mass. Our results align coherently with the extensive NANOGrav 15-year data and the sensitivities outlined by other terrestrial and space-based experiments (e.g.: LISA, HLVK, BBO, HLV(O3), etc.).

We study the evolution of photon injections with a power-law type spectrum inserted at various epochs of the universe, and obtain constraints on their parameter space from multiple different cosmological probes. Our work is motivated by the realistic possibility of having extended photon spectra from astrophysical and exotic sources. Going beyond a $\delta$-function like approximation, the physics becomes richer and the constraining power of cosmological probes starts to depend on the photon injection history in a complex way. As a toy model, we first consider a decaying particle scenario, and then generalize to a more model independent power law type injection in redshift. Different combinations of our parameters can be mapped to a wide variety of realistic astrophysical and exotic sources, providing useful benchmarks for study in future work.

This paper introduces a novel variability report generator developed for the Large Array Survey Telescope (LAST), a cost-effective multi-purpose telescope array conducting a wide survey of the variable sky in the visible-light spectrum. Designed to automate variability detection, the report generator identifies candidate variable stars by employing adjustable thresholds to detect periodic and non-periodic variables. The program outputs a visual and tabular photometric report for each candidate variable source from a given LAST sub-image. Functioning as a whitepaper, this document also provides a concise overview of LAST, discussing its design, data workflow, and variability search performance.

We present a set of isochrone-tailored spectral libraries for analyzing composite spectra of low-metallicity galaxies. Specifically, we have computed synthetic spectra for stars of all initial masses for isochrones at metallicities Z=0.002 and Z=0.0004, with and without considering rotation, constructed by the Geneva group (Ekstr\"{o}m et al., 2011; Georgy et al.. 2013; Groh et al., 2019). We also present a Python program for integrating the individual spectra with a given initial mass function.

We propose a novel method to quantify the assembly histories of dark matter halos with the redshift evolution of the mass-weighted spatial variance of their progenitor halos, i.e. the protohalo size history. We find that the protohalo size history for each individual halo at z~0 can be described by a double power-law function. The amplitude of the fitting function strongly correlates to the central-to-total stellar mass ratios of descendant halos. The variation of the amplitude of the protohalo size history can induce a strong halo assembly bias effect for massive halos. This effect is detectable in observation using the central-to-total stellar mass ratio as a proxy of the protohalo size. The correlation to the descendant central-to-total stellar mass ratio and the halo assembly bias effect seen in the protohalo size are much stronger than that seen in the commonly adopted half-mass formation time derived from the mass accretion history. This indicates that the information loss caused by the compression of halo merger trees to mass accretion histories can be captured by the protohalo size history. Protohalo size thus provides a useful quantity to connect protoclusters across cosmic time and to link protoclusters with their descendant clusters in observations.

We present results from wide-field imaging of the resolved stellar populations of the dwarf spheroidal galaxies Cassiopeia III (And XXXII) and Perseus I (And XXXIII), two satellites in the outer stellar halo of the Andromeda galaxy (M31). Our WIYN pODI photometry traces the red giant star population in each galaxy to ~2.5-3 half-light radii from the galaxy center. We use the Tip of the Red Giant Branch (TRGB) method to derive distances of (m-M)_0 = 24.62+/-0.12 mag (839 (+48,-450) kpc, or 156 (+16,-13) kpc from M31) for Cas III and 24.47+/-0.13 mag (738 (+48,-45) kpc, or 351 (+17,-16) kpc from M31) for Per I. These values are consistent within the errors with TRGB distances derived from a deeper Hubble Space Telescope study of the galaxies' inner regions. For each galaxy, we derive structural parameters, total magnitude, and central surface brightness. We also place upper limits on the ratio of neutral hydrogen gas mass to optical luminosity, confirming the gas-poor nature of both galaxies. We combine our data set with corresponding data for the M31 satellite galaxy Lacerta I (And XXXI) from earlier work, and search for substructure within the RGB star populations of Cas III, Per I, and Lac I. We find an overdense region on the west side of Lac I at a significance level of 2.5-3-sigma and a low-significance filament extending in the direction of M31. In Cas III, we identify two modestly significant overdensities near the center of the galaxy and another at two half-light radii. Per I shows no evidence for substructure in its RGB star population, which may reflect this galaxy's isolated nature.

We present a comprehensive analysis of a Transiting Exoplanet Survey Satellite (TESS) high-quality light curve for a young brown dwarf, MHO~4 having spectral type M7.0, in the Taurus star-forming region. We investigate the rotation periods and characterize the BD's dynamic atmosphere and surface features. We present light curve analysis of MHO~4, and estimate the rotation period to be around 2.224~d. Remarkably, MHO~4 exhibits two significant flaring events. Furthermore, we also estimated bolometric flare energies to be within the energy range of $10^{34}$ to $10^{35}$ erg, which sits in the superflare category.

We investigate the influence of Toomre's $Q$ parameter on the bar-forming dynamics of Maclaurin disk using $N$-body simulations. According to the Toomre's criterion, local velocity dispersion parametrized by $Q\geq 1$ is required to suppress the local axisymmetric instability but, in turn, it deviates particle orbits from nearly circular limit in which particle natural frequencies are calculated. We resolve this by including the effect of velocity dispersion, as the pressure potential, into the effective potential with the gravitational potential. With this formulation, circular orbit approximation is retrieved. The effective potential hypothesis can describe the $Q$-dependences of angular and epicyclic motions of the bar-forming processes and the established bars reasonably well provided that $Q\geq 1$. This indicates the influence of initial $Q$ that is imprinted in the entire disk dynamics, not only that $Q$ serves as the stability indicator. In addition, we perform the stability test for the disk-in-halo systems. With the presence of halo, disks are more susceptible to the bar formation as seen by the elevated critical $Q$ than that for the isolated disk. This is attributed to the differential rotation that builds the unstable non-axisymmetric spiral modes more efficiently which are the ingredients of bar instability.

Using a set of high-resolution N-body simulations, we extend the unified distribution model of cold dark matter (CDM) subhaloes to the warm dark matter(WDM) case. The same model framework combining the unevolved mass function, unevolved radial distribution, and tidal stripping can predict the mass function and spatial distribution of subhaloes in both CDM and WDM simulations. The dependence of the model on the DM particle property is universally parameterized through the half-mode mass of the initial power spectrum. Compared with the CDM model, the WDM model differs most notably in two aspects. 1) In contrast to the power-law form in CDM, the unevolved subhalo mass function for WDM is scale-dependent at the low mass end due to the cut-off in the initial power spectrum. 2) WDM subhaloes are more vulnerable to tidal stripping and disruption due to their lower concentrations at accretion time. Their survival rate is also found to depend on the infall mass. Accounting for these differences, the model predicts a final WDM subhalo mass function that is also proportional to the unevolved subhalo mass function. The radial distribution of WDM subhaloes is predicted to be mass-dependent. For low mass subhaloes, the radial distribution is flatter in the inner halo and steeper in the outer halo compared to the CDM counterpart, due to the scale-dependent unevolved mass function and the enhanced tidal stripping. The code for sampling subhaloes according to our generalized model is available at https://github.com/fhtouma/subgen2 .

The electromagnetic transient following a binary neutron star merger is known as a kilonova (KN). Owing to rapid expansion velocities and small ejecta masses, KNe rapidly transition into the Non-Local Thermodynamic Equilibrium (NLTE) regime. In this study, we present synthetic NLTE spectra of KNe from 5 to 20 days after merger using the \texttt{SUMO} spectral synthesis code. We study three homogeneous composition, 1D multi-zone models with characteristic electron fractions of $Y_e \sim 0.35, 0.25$ and $0.15$. We find that emission features in the spectra tend to emerge in windows of reduced line blocking, as the ejecta are still only partially transparent even at 20 days. For the $Y_e \sim 0.35$ (lanthanide-free) ejecta, we find that the neutral and singly ionised species of Rb, Sr, Y and Zr dominate the spectra, all with good potential for identification. We directly test and confirm an impact of Sr on the 10000 angstrom spectral region in lanthanide-free ejecta, but also see that its signatures may be complex. We suggest the Rb I $\rm{5p^{1}}$- $\rm{5s^{1}}$ 7900 angstrom transition as a candidate for the $\lambda_0 \sim$ 7500--7900 angstrom P-Cygni feature in AT2017gfo. For the $Y_e \sim 0.25$ and $0.15$ compositions, lanthanides are dominant in the spectral formation, in particular Nd, Sm, and Dy. We identify key processes in KN spectral formation, notably that scattering and fluorescence play important roles even up to 20 days after merger, implying that the KN ejecta are not yet optically thin at this time.

We analyze the longitude-velocity diagram (LVD) of the CO-line emission from archival data and use the most accurate rotation curve (RC) of the Milky Way to transform radial velocity to face-on position in the Galactic plane. We point out that the face-on transformation is highly sensitive to the adopted RC, especially in the inner Milky Way, in the sense that deviations of the RC from the true rotation velocity yield an artifact hole or overcrowded concentration along the tangent circle for over- or under-estimated RC. Even if the RC is sufficiently accurate, non-circular motion such as with the 3 kpc expanding ring causes significant artifacts in the resulting face-on-map as long as a circular rotation is assumed. On the other hand, if we properly take into account the non-circular motion, it can be used to solve the near-far degeneracy problem of determination of kinematic distance. We thus propose a new method to solve the degeneracy by incorporating the expanding motion of a ring or arms. We apply the method to the LVD of the 3-kpc expanding ring and present its face-on map projected onto the galactic plane for the first time.

Many modern millimeter and submillimeter (``mm-wave'') telescopes for astronomy are deploying more detectors by increasing detector pixel density, and with the rise of lithographed detector architectures and high-throughput readout techniques, it is becoming increasingly practical to overfill the focal plane. However, when the pixel pitch $p_{\rm pix}$ is small compared to the product of the wavelength $\lambda$ and the focal ratio $F$, or $p_{\mathrm{pix}} \lesssim 1.2 F \lambda$, the Bose term of the photon noise correlates between neighboring detector pixels due to the Hanbury Brown & Twiss (HBT) effect. When this HBT effect is non-negligible, the array-averaged sensitivity scales with detector count $N_{\mathrm{det}}$ less favorably than the uncorrelated limit of $N_{\mathrm{det}}^{-1/2}$. In this paper, we present a general prescription to calculate this HBT correlation based on a quantum optics formalism and extend it to polarization-sensitive detectors. We then estimate the impact of HBT correlations on the sensitivity of a model mm-wave telescope and discuss the implications for focal-plane design.

Kilonova spectra provide us with information of r-process nucleosynthesis in neutron star mergers. However, it is still challenging to identify individual elements in the spectra mainly due to lack of experimentally accurate atomic data for heavy elements in the near-infrared wavelengths. Recently, Domoto et al. (2022) proposed the absorption features around 14500 A in the observed spectra of GW170817/AT2017gfo as Ce III lines. But they used theoretical transition probabilities (gf-values) whose accuracy is uncertain. In this paper, we derive the astrophysical gf-values of the three Ce III lines, aiming at verification of this identification. We model high resolution H-band spectra of four F-type supergiants showing the clear Ce III absorption features by assuming stellar parameters derived from optical spectra in literatures. We also test the validity of the derived astrophysical gf-values by estimating Ce III abundances in Ap stars. We find that the derived astrophysical gf-values of the Ce III lines are systematically lower by about 0.25 dex than those used in previous work of kilonovae, while they are still compatible within the uncertainty range. By performing radiative transfer simulations of kilonovae with the derived gf-values, we find that the identification of Ce III as a source of the absorption features in the observed kilonova spectra still stands, even considering the uncertainties in the astrophysical gf-values. This supports identification of Ce in the spectra of GW170817/AT2017gfo.

It is recognized that some core-collapse supernovae (SNe) show a double-peaked radio light curve within a few years since the explosion. A shell of circumstellar medium (CSM) detached from the SN progenitor has been considered to play a viable role in characterizing such a re-brightening of radio emission. Here, we propose another mechanism that can give rise to the double-peaked radio light curve in core-collapse SNe. The key ingredient in the present work is to expand the model for the evolution of the synchrotron spectral energy distribution (SED) to a generic form, including fast and slow cooling regimes, as guided by the widely-accepted modeling scheme of gamma-ray burst afterglows. We show that even without introducing an additional CSM shell, the radio light curve would show a double-peaked morphology when the system becomes optically thin to synchrotron self-absorption at the observational frequency during the fast cooling regime. We can observe this double-peaked feature if the transition from fast cooling to slow cooling regime occurs during the typical observational timescale of SNe. This situation is realized when the minimum Lorentz factor of injected electrons is initially large enough for the non-thermal electrons' SED to be discrete from the thermal distribution. We propose SN 2007bg as a special case of double-peaked radio SNe that can be explained by the presented scenario. Our model can serve as a potential diagnostic for electron acceleration properties in SNe.

Recently, Primordial Black Holes (PBHs) as a remarkable candidate for Dark Matter (DM) content of the universe have attracted a lot of interests from the scientific community. So, in this study generation of PBHs from Higgs potential in the presence of a tiny bump in non-canonical inflationary model has been inquired. It is demonstrated that, a viable inflationary era can be driven through the Higgs potential with self coupling constant $\lambda=0.13$, in non-canonical framework with a power-law Lagrangian density. Furthermore, setting a suitable function of inflaton field as a correction term (like a bump) to the Higgs potential, causes the inflaton to slow down for a while. In such a short time span, the amplitude of the scalar perturbations power spectrum on small scales grows up sufficiently versus CMB scales. In addition to the bump feature, the enhancing effect of the $\alpha$ parameter of the Lagrangian on the amplitude of the scalar power spectrum has been shown. Fine tuning of three parameter Cases of our model, gives rise to generate three Cases of PBHs with masses ${\cal O}(10)M_{\odot}$, ${\cal O}(10^{-6})M_{\odot}$ and ${\cal O}(10^{-13})M_{\odot}$, which can be appropriate to expound LIGO-VIRGO events, microlensing events in OGLE data and nearly the entirety of DM content, respectively. Generation of these PBHs accompanies by propagating of secondary Gravitational Waves (GWs). It is illustrated that, the contemporary density parameter spectra $(\Omega_{\rm GW_0})$ of these GWs can be tracked down by GWs detectors. Lastly, a power-law behavior for the spectra of $\Omega_{\rm GW_0}$ with regard to frequency has been pointed out. The power index in the infrared domain has the log-reliant form as $n=3-2/\ln(f_c/f)$.

Ultra-high energy cosmic rays (UHECRs) have been studied with the data of the Pierre Auger Observatory for more than fifteen years. An essential feature of the Observatory is its hybrid design: UHECRs are detected through the observation of the associated extensive air showers (EASs) with different and complementary techniques. The analyses of the multi-detector data have enabled high-statistics and high-precision studies of the UHECR energy spectrum, mass composition and distribution of arrival directions. The resulting science of UHECRs is summarized in this contribution. While no discrete source of UHECRs has been identified so far, the extragalactic origin of the particles has been recently confirmed from the arrival directions above 8~EeV, and the ring is closing around nearby astrophysical sites at higher energies. Also, the established upper limits on fluxes of UHE neutrinos and photons have implications on multi-messenger studies and beyond-the-Standard-Model (BSM) physics.

Understanding the interplay of stellar feedback and turbulence in the interstellar medium (ISM) is essential to modeling the evolution of galaxies. To determine the timescales over which stellar feedback drives turbulence in the ISM, we performed a spatially resolved, multi-wavelength study of the nearby star-forming dwarf galaxy UGC 4305 (aka Holmberg II). As indicators of turbulence on local scales (400 pc), we utilized ionized gas velocity dispersion derived from IFU H$\alpha$ observations and atomic gas velocity dispersion and energy surface densities derived from HI synthesis observations with the Very Large Array. These indicators of turbulence were tested against star formation histories over the past 560 Myr derived from Color-Magnitude Diagrams (CMD) using Spearman's rank correlation coefficient. The strongest correlation identified at the 400 pc scale is between measures of HI turbulence and star formation 70-140 Myr ago. We repeated our analysis of UGC 4305's current turbulence and past star formation activity on multiple physical scales ($\sim$560, and 800 pc) to determine if there are indications of changes in the correlation timescale with changes to the physical scale. No notable correlations were found at larger physical scales emphasizing the importance of analyzing star formation driven turbulence as a local phenomenon.

Environmental effects on the evolution of galaxies have been one of the leading questions in galaxy studies for decades. In this work, we investigate the relationship between the star formation activity of galaxies and their environmental matter density using the cosmological hydrodynamic simulation Simba. The star formation activity indicators we explore include the star formation efficiency (SFE), specific star formation rate (sSFR) and molecular hydrogen mass fraction ($f^*_{H_2}$) and the environment is considered as the large-scale environmental matter density, calculated based on the stellar mass of nearby galaxies on a 1 Mpc/h grid using the cloud in cell (CIC) method. Our sample includes galaxies with $9<\log(M_*/M_{\odot})$ at $0<z<4$, divided into three mass bins to disentangle the effects of mass and environment on the galactic star formation activity. For low- to intermediate-mass galaxies at low-redshifts ($z<1.5$), we find that the star formation efficiency of those in high-density regions are $\sim 0.3$ dex lower than those in low-density regions. However, there is no significant environmental dependence of the star formation efficiency for massive galaxies over all our redshift range, and low- to intermediate-mass galaxies at high redshifts ($z > 1.5$). We present a scaling relation for the depletion time of molecular hydrogen (${t_{depl}}=1/SFE$) as a function of galaxy parameters including environmental density. Our findings provide a framework for quantifying the environmental effects on the star formation activities of galaxies as a function of stellar mass and redshift. The most significant environmental dependence is seen at later cosmic times ($z<1.5$) and towards lower stellar masses ($9<\log(M_*/M_{\odot})<10$). Future large galaxy surveys can use this framework to look for the environmental dependence of the star formation activity and examine our predictions.

We investigate the previous microlensing data collected by the KMTNet survey in search of anomalous events for which no precise interpretations of the anomalies have been suggested. From this investigation, we find that the anomaly in the lensing light curve of the event KMT-2021-BLG-1547 is approximately described by a binary-lens (2L1S) model with a lens possessing a giant planet, but the model leaves unexplained residuals. We investigate the origin of the residuals by testing more sophisticated models that include either an extra lens component (3L1S model) or an extra source star (2L2S model) to the 2L1S configuration of the lens system. From these analyses, we find that the residuals from the 2L1S model originate from the existence of a faint companion to the source. The 2L2S solution substantially reduces the residuals and improves the model fit by $\Delta\chi^2=67.1$ with respect to the 2L1S solution. The 3L1S solution also improves the fit, but its fit is worse than that of the 2L2S solution by $\Delta\chi^2=24.7$. According to the 2L2S solution, the lens of the event is a planetary system with planet and host masses $(M_{\rm p}/M_{\rm J}, M_{\rm h}/M_\odot)=\left( 1.47^{+0.64}_{-0.77}, 0.72^{+0.32}_{-0.38}\right)$ lying at a distance $\D_{\rm L} =5.07^{+0.98}_{-1.50}$~kpc, and the source is a binary composed of a subgiant primary of a late G or an early K spectral type and a main-sequence companion of a K spectral type. The event demonstrates the need of sophisticated modeling for unexplained anomalies for the construction of a complete microlensing planet sample.

GRB 220101A is the most distant gamma-ray burst detected by the Fermi-LAT to date, at a redshift z = 4.618. It is also a very energetic event, with an equivalent isotropic energy of $3.6\times10^{54}$ erg. We jointly analyzed the Fermi/GBM and LAT observations of GRB 220101A with two independent approaches, and found a significant spectral break at sub-100 MeV energies during the prompt emission. The fast variability of the emission suggests that this spectral attenuation is caused by internal opacity to pair creation. Regardless of the nature of the emission processes assumed in the spectral analysis, we infer a moderate value for the jet Lorentz factor, $\Gamma\sim110$, and find that all of the high-energy emission was produced above and near the photosphere, at a distance of $\sim10^{14}$ cm from the central engine. We compare these results with the four other LAT-detected gamma-ray bursts with similar properties.

Solar eruptive activities could occur in weak magnetic field environments and over large spatial scales, especially relevant to eruptions involving intermediate or quiescent solar filaments. To handle the large scales, we implement and apply a flux rope embedding method using regularized Biot-Savart laws in the spherical coordinate system. Combined with a potential field source surface model and a magneto-frictional method, a nonlinear force-free field comprising a flux rope embedded in a potential field is constructed. Using the combined nonlinear force-free field as the initial condition, we then perform a zero-$\beta$ data-constrained magnetohydrodynamic (MHD) simulation for an M8.7 flare at 03:38 UT on 2012 January 23. The MHD model reproduces the eruption process, flare ribbon evolution (represented by the quasi-separatrix layer evolution) and kinematics of the flux rope. This approach could potentially model global-scale eruptions from weak field regions.

Combining cosmological probes has consolidated the standard cosmological model with percent precision, but some tensions have recently emerged when certain parameters are estimated from the local or primordial Universe. Whether this is due to hidden systematics or it points towards new physics is still under debate, however, it is crucial to study as many probes as possible to cross-check the results with independent methods and provide additional pieces of information to the cosmological puzzle. In this work, by combining several late-Universe probes sampling the redshift range $0<z<10$, namely, Type Ia SuperNovae, Baryon Acoustic Oscillations, Cosmic Chronometers and Gamma-Ray Bursts, we aim to derive cosmological constraints independently of local or early-Universe anchors. To test the standard cosmological model and its various extensions, considering an evolving Dark Energy Equation of State and the curvature as a free parameter, we analyse each probe individually and all their possible permutations. Assuming a flat $\Lambda$CDM model, the full combination of probes provides $H_0=67.2^{+3.4}_{-3.2}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_m=0.325\pm0.015$ ($68\%$ C.L.). Considering a flat $w$CDM model, we measure $w_0=-0.91^{+0.07}_{-0.08}$ ($68\%$ C.L.), while by relaxing the flatness assumption ($\Lambda$CDM model, $95\%$ C.L.) we obtain $\Omega_k=0.125^{+0.167}_{-0.165}$. Finally, we analytically characterize the degeneracy directions and the relative orientation of the probes' contours. By calculating the Figure--of--Merit, we quantify the synergies among independent methods, estimate the constraining power of each probe and identify which provides the best contribution to the inference process. Pending the new cosmological surveys, this study confirms the exigency for new emerging cosmological probes in the landscape of modern cosmology.

The non-local thermodynamical equilibrium (NLTE) line formation of Y I and Y II is considered in 1D LTE model atmospheres of F-G-K-type stars. The model atom was constructed with the most up-to-date atomic data, including quantum cross sections and rate coefficients for transitions in inelastic collisions of Y I and Y II with hydrogen atoms. For seven reference stars, we obtained an agreement between NLTE abundances inferred from the two ionization stages, while the difference in LTE abundance (Y I - Y II) can reach up to -0.31 dex. In the atmospheres of F-G-K-type stars, for both Y I and Y II lines, the NLTE abundance corrections are positive. In solar metallicity stars, the NLTE abundance corrections for Y II lines do not exceed 0.12 dex, while in atmospheres of metal-poor stars they do not exceed 0.21 dex. For Y I lines, the NLTE abundance corrections can reach up to 0.5 dex. We determined the yttrium NLTE abundances for a sample of 65 F and G dwarfs and subgiants in the -2.62~$\leq$~[Fe/H]~$\leq$~+0.24 metallicity range, using high-resolution spectra. For stars with [Fe/H]~$\leq$~-1.5, [Y/Fe] versus [Fe/H] diagram reveals positive trend with an average value of [Y/Fe]~$\simeq$~0. For metal-poor stars, among Sr, Y, and Zr, the arrangement [Sr/Fe] < [Y/Fe] < [Zr/Fe] remains consistent. The current study is useful for the Galactic chemical evolution research. The model atom will be applied for NLTE yttrium abundance determination in very metal-poor stars studied with LAMOST and Subaru.

We review what could be astronomy from the Moon in the next decades in the visible domain. After a short review observational approaches, from photometry to high contrast and high angular resolution imaging, We essentially focus on some promising scientific objectives, from Solar System to the extragalactic domain. At the end, I add a proposal to use the Earth-Moon system to test fundamental physics. Since this meeting is dedicated to the next decades of Astronomy from the Moon, we consider projects and science objectives for several decades from now.

For the past several decades, numerous attempts have been made to model the climate of Mars with extensive studies focusing on the planet's dynamics and the understanding of its climate. While physical modeling and data assimilation approaches have made significant progress, uncertainties persist in comprehensively capturing and modeling the complexities of Martian climate. In this work, we propose a novel approach to Martian climate modeling by leveraging machine learning techniques that have shown remarkable success in Earth climate modeling. Our study presents a deep neural network designed to accurately model relative humidity in Gale Crater, as measured by NASA's Mars Science Laboratory ``Curiosity'' rover. By utilizing simulated meteorological variables produced by the Mars Planetary Climate Model, a robust Global Circulation Model, our model accurately predicts relative humidity with a mean error of 3\% and an $R^2$ score of 0.92. Furthermore, we present an approach to predict quantile ranges of relative humidity, catering to applications that require a range of values. To address the challenge of interpretability associated with machine learning models, we utilize an interpretable model architecture and conduct an in-depth analysis of its internal mechanisms and decision making processes. We find that our neural network can effectively model relative humidity at Gale crater using a few meteorological variables, with the monthly mean surface H$_2$O layer, planetary boundary layer height, convective wind speed, and solar zenith angle being the primary contributors to the model predictions. In addition to providing a fast and efficient method to modeling climate variables on Mars, this modeling approach can also be used to expand on current datasets by filling spatial and temporal gaps in observations.

The challenges inherent to time-domain multi-messenger astronomy require strategic actions so that adapted, optimized follow-up observations are performed efficiently. In particular, poorly localized events require dedicated tiling and/or targeted, follow-up campaigns so that the region in which the source really is can be efficiently covered, increasing the chances to detect the multi-wavelength counterpart. We have developed the python package "tilepy" to rapidly derive the observation scheduling of large uncertainty localization events by small/mid-FoV instruments. We will describe several mature follow-up scheduling strategies. These range from an option to use of low-resolution grids, to the full integration of sky regions and targeted observations using galaxy catalogs. The algorithms consider the visibility constraints of customisable observatories and allow to schedule observations in both astronomical darkness and in moonlight conditions. Developed initially to provide a rapid response to gravitational wave (GW) alerts by Imaging Atmospheric Cherenkov Telescopes (IACTs), they have been proven successful, as shown by the GW follow-up during O2 and O3 with the H.E.S.S. telescopes, and particularly in the follow-up of GW170817, the first binary neutron star (BNS) merger ever detected. Here we will present a generalisation of these rapid strategies to other alerts showing large uncertainties in the localization, like Gamma-Ray Burst (GRB) alerts from Fermi-GBM. We will also demonstrate the flexibility of {\it tilepy} in scheduling observations for a large variety of observatories. We will conclude by describing the latest developments of these algorithms that are able to derive optimised follow-up schedules across multiple observatories and networks of telescopes.

We used high resolution spectra acquired at the Magellan Telescope to measure radial and rotational velocities of approximately 200 stars in the Galactic globular cluster NGC 3201. The surveyed sample includes Blue Stragglers Stars (BSSs) and reference stars in different evolutionary stages (main sequence turn-off, sub-giant, red giant and asymptotic giant branches). The average radial velocity value ($\langle V_r\rangle = 494.5 \pm 0.5$ km s$^{-1}$) confirms a large systemic velocity for this cluster and was used to distinguish 33 residual field interlopers. The final sample of member stars counts 67 BSSs and 114 reference stars. Similarly to what is found in other clusters, the totality of the reference stars has negligible rotation ($<20$ km s$^{-1}$), while the BSS rotational velocity distribution shows a long tail extending up to $\sim 200$ km s$^{-1}$, with 19 BSSs (out of 67) spinning faster than 40 km s$^{-1}$. This sets the percentage of fast rotating BSSs to $\sim 28\%$. Such a percentage is roughly comparable to that measured in other loose systems ($\omega$ Centauri, M4 and M55) and significantly larger than that measured in high-density clusters (as 47 Tucanae, NGC 6397, NGC 6752 and M30). This evidence supports a scenario where recent BSS formation (mainly from the evolution of binary systems) is occurring in low-density environments. We also find that the BSS rotational velocity tends to decrease for decreasing luminosity and surface temperature, similarly to what is observed in main sequence stars. Hence, further investigations are needed to understand the impact of BSS internal structure on the observed rotational velocities.

Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study a complete sample of massive clusters beyond 7 (5) $\times$ 10$^{14}$ M$_{\odot}$ at redshift 0 (1) on a $(1.48 \ \mathrm{Gpc})^3$ volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 $\times$ 10$^9$ M$_{\odot}$) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark matter only version of these simulations, we demonstrate that baryonic physics preserves significantly low mass subhalos (galaxies) as have also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark matter only simulations of the same objects and taking the same limit in subhalo mass we find significantly more cuspy galaxy density profiles towards the center of the clusters, where the low mass substructures would tend to concentrate. We conclude that using dark matter only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high resolution hydro simulations we conclude that a mass resolution of 1.8 $\times$ 10$^8$ h$^{-1}$ M$_{\odot}$ will be needed for The Three Hundred simulations to approach the expected magnitude limits for the Euclid survey.

This paper examines the energetics of a convective flow subject to an oscillation with a period $t_{\rm osc}$ much smaller than the convective timescale $t_{\rm conv}$, allowing for compressibility and uniform rotation. We show that the energy of the oscillation is exchanged with the kinetic energy of the convective flow at a rate $D_R$ that couples the Reynolds stress of the oscillation with the convective velocity gradient. For the equilibrium tide and inertial waves, this is the only energy exchange term, whereas for p modes there are also exchanges with the potential and internal energy of the convective flow. Locally, $\left| D_R \right| \sim u^{\prime 2} /t_{\rm conv}$, where $u^{\prime}$ is the oscillating velocity. If $t_{\rm conv} \ll t_{\rm osc}$ and assuming mixing length theory, $\left| D_R \right|$ is $\left( \lambda_{\rm conv}/ \lambda_{\rm osc} \right)^2$ smaller, where $\lambda_{\rm conv}$ and $\lambda_{\rm osc}$ are the characteristic scales of convection and the oscillation. Assuming local dissipation, we show that the equilibrium tide lags behind the tidal potential by a phase $\delta \left( r \right) \sim r \omega_{\rm osc}/ \left( g \left( r \right) t_{\rm conv} \left( r \right)\right)$, where $g$ is the gravitational acceleration. The equilibrium tide can be described locally as a harmonic oscillator with natural frequency $\left( g / r \right)^{1/2}$ and subject to a damping force $-u^{\prime}/t_{\rm conv}$. Although $\delta \left( r \right) $ varies by orders of magnitude through the flow, it is possible to define an average phase shift $\overline{\delta}$ which is in good agreement with observations for Jupiter and some of the moons of Saturn. Finally, $1 / \overline{\delta}$ is shown to be equal to the standard tidal dissipation factor.

The brightness of faculae and network depends on the angle at which they are observed and the magnetic flux density. Close to the limb, assessment of this relationship has until now been hindered by the increasingly lower signal in magnetograms. This preliminary study aims at highlighting the potential of using simultaneous observations from different vantage points to better determine the properties of faculae close to the limb. We use data from the Solar Orbiter/Polarimetric and Helioseismic Imager (SO/PHI), and the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), recorded at $\sim60^\circ$ angular separation of their lines of sight at the Sun. We use continuum intensity observed close to the limb by SO/PHI and complement it with the co-observed $B_{\rm LOS}$ from SDO/HMI, originating closer to disc centre (as seen by SDO/HMI), thus avoiding the degradation of the magnetic field signal near the limb. We derived the dependence of facular brightness in the continuum on disc position and magnetic flux density from the combined observations of SO/PHI and SDO/HMI. Compared with a single point of view, we were able to obtain contrast values reaching closer to the limb and to lower field strengths. We find the general dependence of the limb distance at which the contrast is maximum on the flux density to be at large in line with single viewpoint observations, in that the higher the flux density is, the closer the turning point lies to the limb. There is a tendency, however, for the maximum to be reached closer to the limb when determined from two vantage points. We note that due to the preliminary nature of this study, these results must be taken with caution. Our analysis shows that studies involving two viewpoints can significantly improve the detection of faculae near the solar limb and the determination of their brightness contrast relative to the quiet Sun.

To overcome the high optical extinction, near-infrared observations are needed for probing the microlensing events toward the Galactic center. The 2015-2019 UKIRT microlensing survey toward the Galactic center is the first dedicated precursor near-infrared (NIR) survey for the Nancy Grace Roman Space Telescope. We here analyze the online data from the UKIRT microlensing survey, reaching $l=b=0^\circ$. Using the event-finder algorithm of KMTNet with the $\Delta \chi^2$ threshold of 250, we find 522 clear events, 436 possible events, and 27 possible anomalous events. We fit a point-source point-lens (PSPL) model to all the clear events and derive the PSPL parameters with uncertainties using a Markov chain Monte Carlo method. Assuming perfect detection efficiency, we compute the uncorrected event rates, which should serve as the lower limits on the true event rate. We find that the uncorrected NIR event rates are likely rising toward the Galactic center and higher than the optical event rates.

Since their discoveries in 1967, Gamma-Ray Bursts (GRBs) continue to be one of the most researched objects in astrophysics. Multi-messenger observations are key to gaining a deeper understanding of these events. In order to facilitate such measurements, fast and accurate localization of the gamma-ray prompt emission is required. As traditional localization techniques are often time consuming or prone to significant systematic errors, here we present a novel method which can be applied on the POLAR-2 observatory. POLAR-2 is a dedicated GRB polarimeter, which will be launched towards the China Space Station (CSS) in 2025. The CSS provides POLAR-2 access to a GPU, which makes it possible and advantageous to run a Deep Learning model on it. In this work, we explore the possibility to identify GRBs in real time and to infer their location and spectra with deep learning models. Using POLAR simulations and data, a feasibility experiment was performed to implement this method on POLAR-2. Our results indicate that using this method, in combination with real time data downlinking capabilities, POLAR-2 will be able to provide accurate localization alerts within 2 minutes of the GRB onset.

We report on the chemical composition of the very metal-poor ([Fe/H]=-2.9) star LAMOST J1645+4357 that is identified to be a red giant having peculiar abundance ratios by Li et al. (2022). The standard abundance analysis is carried out for this object and the well studied metal-poor star HD~122563 that has similar atmospheric parameters. LAMOST J1645+4357 has a remarkable abundance set, highlighted by these features: (1) Nitrogen is significantly enhanced ([N/Fe]=+1.4) and the total abundance of C and N is also very high ([(C+N)/Fe]=+0.9); (2) alpha-elements are over-abundant with respect to iron as generally found in very metal-poor stars; (3) Ti, Sc, Co and Zn are significantly deficient; (4) Cr and Mn are enhanced compared to most of very metal-poor stars; (5) Sr and Ba are deficient and the Sr/Ba ratio ([Sr/Ba]=-1.0) is significantly lower than the value expected for the r-process. The overall abundance pattern of this object from C to Zn is well reproduced by a faint supernova model assuming spherical explosion, except for the excess of Cr and Mn which requires enhancement of incomplete Si burning or small contributions of a type Ia supernova or a pair-instability supernova. There remains, however, a question why the abundance pattern of this star is so unique among very metal-poor stars.

Ly$\alpha$ nebulae ubiquitously found around z>2 quasars can supply unique constraints on the properties of the Circumgalactic Medium, such as its density distribution, provided the quasar halo mass is known. We present a new method to constrain quasar halo masses based on the line-of-sight velocity dispersion maps of Ly$\alpha$ nebulae. By using MUSE-like mock observations obtained from cosmological hydrodynamic simulations under the assumption of maximal quasar fluorescence, we show that the velocity dispersion radial profiles of Ly$\alpha$-emitting gas are strongly determined by gravity and that they are thus self-similar with respect to halo mass when rescaled by the virial radius. Through simple analytical arguments and by exploiting the kinematics of HeII1640\.A emission for a set of observed nebulae, we show that Ly$\alpha$ radiative transfer effects plausibly do not change the shape of the velocity dispersion profiles but only their normalisation without breaking their self-similarity. Taking advantage of these results, we define the variable $\eta^{140-200}_{40-100}$ as the ratio of the median velocity dispersion in two specifically selected annuli and derive an analytical relation between $\eta^{140-200}_{40-100}$ and the halo mass which can be directly applied to observations. We apply our method to 37 observed quasar Ly$\alpha$ nebulae at 3<z<4.7 and find that their associated quasars are typically hosted by ~$10^{12.16 \pm 0.14}$ M$_{\odot}$ haloes independent of redshift within the explored range. This measurement, which is completely independent of clustering methods, is consistent with the lowest mass estimates based on quasar auto-correlation clustering at z~3 and with quasar-galaxies cross-correlation results.

Proto-neutron stars (PNS) are born hot, with temperatures exceeding a few times $10^{10}$ K. In these conditions, the PNS crust is expected to be made of a Coulomb liquid composed of an ensemble of different nuclear species. We perform a study of the beta-equilibrated PNS crust in the liquid phase in a self-consistent multi-component plasma (MCP) approach, thus allowing us to consistently calculate the impurity parameter, often taken as a free parameter in cooling simulations. We developed a self-consistent MCP approach at finite temperature using a compressible liquid-drop description of the ions, with surface parameters adjusted to reproduce experimental masses. The treatment of the ion centre-of-mass motion was included through a translational free-energy term accounting for in-medium effects. The results of self-consistent MCP calculations are systematically compared with those performed in a perturbative and in the one-component plasma treatment. We show that the inclusion of non-linear mixing terms arising from the ion centre-of-mass motion leads to a breakdown of the ensemble equivalence between the one-component and MCP approach. Our findings illustrate that the abundance of light nuclei becomes important, eventually dominating the distribution at higher density and temperature. This is reflected in the impurity parameter, which, in turn, may have a potential impact on NS cooling. For practical applications, we also provide a fitting formula for the impurity parameter in the PNS inner crust. Our results obtained within a self-consistent MCP approach show important differences in the prediction of the PNS composition with respect to those obtained with a one-component or a perturbative MCP approximation, particularly in the deeper region of the crust. This highlights the importance of a full, self-consistent MCP calculation for reliable predictions of the PNS crust composition.

Using the absolute detection calibration and abundant detections of the OSSOS (Outer Solar System Origin Survey) project, we provide population measurements for the main Kuiper Belt. For absolute magnitude $H_r<8.3$, there are 30,000 non-resonant main-belt objects, with twice as many hot-component objects than cold, and with total mass of 0.014 $M_\Earth$, only 1/7 of which is in the cold belt (assuming a cold-object albedo about half that of hot component objects). We show that transneptunian objects with $5.5 < H_r < 8.3$ (rough diameters 400--100~km) have indistinguishable absolute magnitude (size) distributions, regardless of being in the cold classical Kuiper belt (thought to be primordial) or the `hot' population (believed to be implanted after having been formed elsewhere). We discuss how this result was not apparent in previous examinations of the size distribution due to the complications of fitting assumed power-law functional forms to the detections at differing depths. This shared size distribution is surprising in light of the common paradigm that the hot population planetesimals formed in a higher density environment much closer to the Sun, in an environment that also (probably later) formed larger (dwarf planet and bigger) objects. If this paradigm is correct, % << BG added clause our result implies that planetesimal formation was relatively insensitive to the local disk conditions and that the subsequent planet-building process in the hot population did not modify the shape of the planetesimal size distribution in this 50--300~km range.

The well-known anti-correlation between the optical/ultraviolet (UV) emission line equivalent widths of active galactic nuclei and the continuum luminosity (the so-called Baldwin effect) is a long-standing puzzle. One common hypothesis is that more luminous sources have softer spectral energy distribution (SED) in the extreme UV (EUV), as revealed by some observational studies. In this work we revisit this issue through cross-matching SDSS quasars with GALEX far-UV/near-UV catalogs and correcting the effect of a severe observational bias of significant UV detection incompleteness, i.e., the more luminous in observed-frame optical, the more likely detected in observed-frame UV. We find that, for GALEX detected quasars at 1.8 < z < 2.2, the rest-frame mean UV SED (~ 500 -- 3000 Angstrom) bewilderingly shows no luminosity dependence at log(\nu L_\nu(2200 Angstrom)) > 45 (up to 47.3), contrary to the standard thin disc model predictions and the observed Baldwin effect in this luminosity range. Probably, the universal mean UV SED is the result of a local atomic-originated process, and in fainter quasars stronger disk turbulence launching more clouds is the main origin of the Baldwin effect. After correcting for the absorption of the intergalactic medium, a rest-frame intrinsic mean EUV SED is derived from a sub-sample of bright quasars and is found to be much redder in the EUV than all previous quasar composite spectra, highlighting the significance of properly accounting for the sample incompleteness. Interestingly, the global consistence between our extremely red mean EUV SED and the line-driven wind model again supports an origin of a local physical process.

Ammonia is one of the most widely observed molecules in space, and many observations are able to resolve the hyperfine structure due to the electric quadrupole moment of the nitrogen nucleus. The observed spectra often display anomalies in the satellite components of the lines, which indicate substantial deviations from the local thermodynamic equilibrium. The interpretation of the spectra thus requires the knowledge of the rate coefficients for the hyperfine excitation of NH$_3$ induced by collisions with H$_2$ molecules, the dominant collider in the cold interstellar medium. In this paper we present the first such calculations using a recoupling approach. The rate coefficients are obtained for all hyperfine levels within rotation-inversion levels up to $j=4$ and temperatures up to 100 K by means of quantum scattering close-coupling calculations on an accurate, five-dimensional, potential energy surface. We show that the rate coefficients depart significantly from those obtained with the statistical approach and that they do not conform to any simple propensity rules. Finally, we perform radiative transfer calculations to illustrate the impact of our new rate coefficients by modelling the hyperfine line intensities of the inversion transition in ground state para-NH$_3$ ($j_k=1_1$) and of the rotational transition $1_0\rightarrow 0_0$ in ortho-NH$_3$.

The virial equation is used to clarify the nature of the dynamic evolution of a stellar system. Compared to the kinetic equation, it gives a deeper but incomplete description of the process of relaxation to a quasi-stationary state, which here means the fulfillment of the virial theorem. Analysis shows that the time to reach the virial equlibrium state $T_v$ is about two to three dozen dynamic time periods $T_d$. Namely, during $T_v$ the virial ratio, the mean harmonic radius, and the root-mean-square radius of the system fluctuate, and then the first two characteristics stabilize near their equilibrium values, while the root-mean-square radius continues to grow (possibly ad infinitum). This indicates a fundamentally different behavior of the moment of inertia of the system relative to the center of gravity and its potential energy, leading to the formation of a relatively small equilibrium core and an extended halo.

We conduct a line-by-line differential analysis of a sample of 125 co-moving pairs of stars (dwarfs and subgiants near solar metallicity). We obtain high precision stellar parameters with average uncertainties in effective temperature, surface gravity and metallicity of 16.5 K, 0.033 dex and 0.014 dex, respectively. We classify the co-moving pairs of stars into two groups, chemically homogeneous (conatal; |Delta[Fe/H]| $\le$ 0.04 dex) and inhomogeneous (non-conatal), and examine the fraction of chemically homogeneous pairs as a function of separation and effective temperature. The four main conclusions from this study are: (1) A spatial separation of \ds = 10$^6$ AU is an approximate boundary between homogeneous and inhomogeneous pairs of stars, and we restrict our conclusions to only consider the 91 pairs with \ds $\le$ 10$^6$ AU; (2) There is no trend between velocity separation and the fraction of chemically homogeneous pairs in the range \dv $\le$ 4 \kms; (3) We confirm that the fraction of chemically inhomogeneous pairs increases with increasing \teff\ and the trend matches a toy model of that expected from planet ingestion; (4) Atomic diffusion is not the main cause of the chemical inhomogeneity. A major outcome from this study is a sample of 56 bright co-moving pairs of stars with chemical abundance differences $\leq$ 0.02 dex (5\%) which is a level of chemical homogeneity comparable to that of the Hyades open cluster. These important objects can be used, in conjunction with star clusters and the \gaia\ ``benchmark'' stars, to calibrate stellar abundances from large-scale spectroscopic surveys.

Our previous investigation of neutron-star crusts, based on the functional BSk24, led to a substantial reduction of the pasta mantle when Strutinsky integral and pairing corrections were added on top of the fourth-order extended Thomas-Fermi method (ETF). Here, our earlier calculations are widened to a larger set of functionals within the same family, and we find that the microscopic corrections weaken significantly the influence of the symmetry energy. In particular, the correlation observed at the pure ETF level between the density for the onset of pasta formation and the symmetry energy vanishes, not only for the $L$ coefficient but also for the symmetry-energy values at the relevant densities. Moreover, the inclusion of microscopic corrections results in a much lower abundance of pasta for all functionals.

We present a comprehensive study of the excitation of CI fine-structure levels along 57 sight lines in the Large and Small Magellanic Clouds. The sight lines were selected by the detection of H$_2$ in FUSE spectra. Using archival HST/COS and HST/STIS spectra we detected absorption of CI fine-structure levels and measured their populations for 29 and 22 sight lines in the LMC and SMC, respectively. The CI column density ranges from $10^{13}$ to $10^{14}\,{\rm cm}^{-2}$ for the LMC and $10^{13}$ to $10^{15.4}\,{\rm cm}^{-2}$ for the SMC. We found excitation of CI fine-structure levels in the LMC and SMC to be 2-3 times higher than typical values in local diffuse ISM. Comparing excitation of both CI fine-structure levels and H$_2$ rotational levels with a grid of PDR Meudon models we find that neutral cold gas in the LMC and SMC is illuminated by stronger UV field than in local ISM ($\chi=5^{+7}_{-3}$ units of Mathis field for the LMC and $2^{+4}_{-1}$ for the SMC) and has on average higher thermal pressure ($\log p/k =4.2\pm0.4$ and $4.3\pm0.5$, respectively). Magellanic Clouds sight lines likely probe region near star-formation sites, which also affects the thermal state and CI/H$_2$ relative abundances. At the same time such high measurements of UV field are consistent with some values obtained at high redshifts. Together with low metallicities this make Magellanic Clouds to be an interesting test case to study of the central parts of high redshift galaxies.

HD molecule is one the most abundant molecule in the Universe and due to its sensibility to the conditions in the medium, it can be used to constrain physical parameters in the medium where HD resides. Lately we have shown that HD abundance can be enhanced in the low metallicity medium. Large and Small Magellanic Clouds give us an opportunity to study low metallicity galaxies in details towards different sightlines due to their proximity to our Galaxy. We revisited FUSE space telescope archival spectra towards bright stars in Magellanic Clouds to search for HD molecules, associated with the medium of these galaxies. We reanalysed H$_2$ absorption lines and constrained HD column density at the positions of H$_2$ components. We detected HD towards 24 sightlines (including 19 new detections). We try to measure cosmic ray ionization rate for several systems using measured $N({\rm HD})/N({\rm H_2})$, and in most cases get loose constraints due to insufficient quality of the FUSE spectra.

Compact obscured nuclei (CONs) are relatively common in the centers of local (U)LIRGs, yet their nature remains unknown. Both AGN activity and extreme nuclear starbursts have been suggested as plausible nuclear power sources. The prevalence of outflows in these systems suggest that CONs represent a key phase in the nuclear feedback cycle, in which material is ejected from the central regions of the galaxy. Here, we present results from MUSE for the confirmed local CON galaxy NGC4418. For the first time we spatially map the spectral features and kinematics of the galaxy in the optical, revealing several previously unknown structures. In particular, we discover a bilateral outflow along the minor axis, an outflowing bubble, several knot structures and a receding outflow partially obscured by the galactic disk. Based on the properties of these features, we conclude that the CON in NGC4418 is most likely powered by an AGN.

450 years after the explosion of the Type Ia SN1572, the dynamics of the Tycho supernova remnant can give us keys to understand the explosion mechanism and the interaction of the remnant with the interstellar medium. To probe the asymmetries and the evolution of the SNR, we track the ejecta dynamics using new methods applied to the deep X-ray observations available in the Chandra space telescope archive. For the line of sight velocity measurement Vz, we use the Doppler effect focused on the bright Si line in the 1.6-2.1 keV band. Using the component separation tool General Morphological Component Analysis (GMCA), we successfully disentangle the red and blueshifted Si ejecta emission. This allows us to reconstruct a map of the peak energy of the Si line with a total coverage of the SNR at a 2'' resolution and a proxy of the velocity in the line of sight. For the proper motions in the plane of the sky Vxy, we develop a new method, named Poisson Optical Flow, to measure the displacement of 2D features between the observations of 2003 and 2009. The result is a field of 1700 velocity vectors covering the entire SNR. These exhaustive 3D velocity measurements reveal the complex and patchy dynamics of the SNR. At the large-scale, an asymmetry with the North being dominantly blueshifted and the South redshifted is observed. The proper motion vector field Vxy highlights different dynamics between the East and the West parts of the SNR. The eastern velocity field is more disturbed by external inhomogeneities and the South-East ejecta knot. In particular, a slow-down is observed in the North-East which could be due to the interaction with higher densities as seen in other wavelengths. The vector field is also used to backtrace the center of the explosion which is then compared with potential stellar progenitors distances from the latest Gaia DR3, leaving only stars B and E as possible candidates.

In this study, we investigate a cosmological model involving a negative cosmological constant (AdS vacua in the dark energy sector). We consider a quintessence field on top of a negative cosmological constant and study its impact on cosmological evolution and structure formation. We use the power spectrum of the redshifted HI 21 cm brightness temperature maps from the post-reionization epoch as a cosmological probe. The signature of baryon acoustic oscillations (BAO) on the multipoles of the power spectrum is used to extract measurements of the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$. The projected errors on these are then subsequently employed to forecast the constraints on the model parameters ($\Omega_\Lambda, w_0, w_a$) using Markov Chain Monte Carlo techniques. We find that a negative cosmological constant with a phantom dark energy equation of state (EoS) and a higher value of $H_0$ is viable from BAO distance measurements data derived from galaxy samples. We also find that BAO imprints on the 21cm power spectrum obtained from a futuristic SKA-mid like experiment yield a $1-\sigma$ error on a negative cosmological constant and the quintessence dark energy EoS parameters to be $\Omega_\Lambda=-0.883^{0.978}_{-2.987} $ and $w_0=-1.030^{0.023}_{-0.082}$, $w_a=-0.088^{0.162}_{-0.343}$ respectively, which is competitive with other probes reported in the literature.

The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. We have previously reported the discovery of a small molecular outflow from the edge-on T Tauri star in the Bok globule CB26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the disk. We report new, high-resolution mm-interferometric observations of CB26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface. The IRAM PdBI was used to observe CO(2-1) at 1.3mm with a resolution of 0.5". Using a physical model of the disk, which was derived from the dust emission, we employed chemo-dynamical modeling combined with line radiative transfer to constrain kinematic parameters and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB26. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ~7deg, which can be traced out to vertical heights of ~500au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, which we trace down to vertical heights of 40au, where it is launched from the surface of the flared disk at radii of 20-45au. The observed CO outflow has a total momentum flux of 1e-5 Msun km/s/yr, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly an MHD disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk.

This paper presents a comprehensive analysis of four triply eclipsing triple star systems, namely TIC 88206187, TIC 14839347, TIC 298714297, and TIC 66893949. The four systems with third-body eclipses were found in the TESS lightcurves from among a sample of ~400 matches between known eclipsing binaries and the Gaia DR3 Non-Single Star (NSS; Gaia Collaboration 2022; Pourbaix et al. 2022) solution database. We combined photometric lightcurves, eclipse timing variations, archival spectral energy distributions, and theoretical evolution tracks in a robust photodynamical analysis to determine the orbital and system parameters. The triples have outer periods of 52.9, 85.5, 117, and 471 days, respectively. All dozen stars have masses $\lesssim$ 2.6 M$_\odot$. The systems are quite flat with mutual inclination angles between the inner and outer orbital planes that are all $\lesssim 4^\circ$. The outer mass ratios range from 0.39--0.76, consistent with our earlier collection of compact triply eclipsing triples. TIC 88206187 exhibits a fractional radius of the outer tertiary component $(r_B \equiv R_B/a_{\rm out})$ exceeding 0.1 (only the third such system known), and we consider its future evolution. Finally, we compare our photodynamical analysis results and the orbital parameters given in the Gaia DR3 NSS solutions, indicating decent agreement, but with the photodynamical results being more accurate.

Significant evidence for a stochastic gravitational-wave background has recently been reported by several Pulsar Timing Array observations. These studies have shown that, in addition to astrophysical explanations based on supermassive black hole binaries (SMBHBs), cosmological origins are considered equally important sources for these signals. To further explore these cosmological sources, in this study, we discuss the anisotropies in the cosmological gravitational wave background (CGWB) in a model-independent way. Taking the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 15-year dataset as a benchmark, we estimate the angular power spectra of the CGWB and their cross-correlations with cosmic microwave background (CMB) fluctuations and weak gravitational lensing. We find that the NANOGrav 15-year data implies suppressed Sachs-Wolf (SW) effects in the CGBW spectrum, leading to a marginally negative cross-correlation with the CMB at large scales. This procedure is applicable to signals introduced by different early universe processes and is potentially useful for identifying unique features about anisotropies of CGWB from future space-based interferometers and astrometric measurements.

The space sector is experiencing a flourishing growth and evidence is mounting that the near future will be characterized by a large amount of deep-space missions. In the last decade, CubeSats have granted affordable access to space due to their reduced manufacturing costs compared to traditional missions. At the present-day, most miniaturized spacecraft have thus far been deployed into near-Earth orbits, but soon a multitude of interplanetary CubeSats will be employed for deep-space missions as well. Nevertheless, the current paradigm for deep-space missions strongly relies on ground-based operations. Although reliable, this approach will rapidly cause saturation of ground slots, thereby hampering the current momentum in space exploration. At the actual pace, human-in-the-loop, flight-related operations for deep-space missions will soon become unsustainable. Self-driving spacecraft are challenging the current paradigm under which spacecraft are piloted in interplanetary space. They are intended as machines capable of traveling in deep space and autonomously reaching their destination. In EXTREMA, these systems are used to engineer ballistic capture (BC), thereby proving the effectiveness of autonomy in a complex scenario. The key is to accomplish low-thrust orbits culminating in BC. For this, a bundle of BC orbits named ballistic capture corridor (BCC) can be targeted far away from a planet. To achieve BC at Mars without any a priori instruction, an inexpensive and accurate method to construct BCC directly on board is required. Therefore, granting spacecraft the capability to manipulate stable sets in order to self-compute a BCC is crucial. The goal of the paper is to numerically synthesize a corridor exploiting the polynomial chaos expansion (PCE) technique, thereby applying a suited uncertainty propagation technique to BC orbit propagation.

Water is a very abundant molecule in star-forming regions. Its deuterium fractionation is an important tool for understanding its formation and evolution during the star and planet formation processes. While the HDO/H$_2$O ratio has been determined toward several Class 0 protostars and comets, the number of studies toward Class I protostars is limited. We aim to study the water deuteration toward the Class I binary protostar L1551 IRS5 and to investigate the effect of evolutionary stage and environment on variations in the water D/H ratio. Observations were made using the NOEMA interferometer. The HDO 3$_{1,2}$-2$_{2,1}$ transition at 225.9 GHz and the H$_2^{18}$O 3$_{1,3}$-2$_{2,0}$ transition at 203.4 GHz were covered with a spatial resolution of 0.5'' $\times$ 0.8'', while the HDO 4$_{2,2}$-4$_{2,3}$ transition at 143.7 GHz was observed with a resolution of 2.0'' $\times$ 2.5''. We used both LTE and non-LTE models. The three transitions are detected. The line profiles display two peaks, one at $\sim$6 km s$^{-1}$ and one at $\sim$9 km s$^{-1}$. We derive an HDO/H$_2$O ratio of (2.1 $\pm$ 0.8) $\times$ 10$^{-3}$ for the redshifted component and a lower limit of $>$ 0.3 $\times$ 10$^{-3}$ for the blueshifted component due to the blending with the redshifted CH$_3$OCH$_3$ emission. The HDO/H$_2$O in L1551 IRS5 is similar to the ratios in isolated Class 0 sources and to the Class I V883 Ori, while it is significantly higher than in the clustered Class 0 sources and the comets. This suggests that the chemistry of protostars in low source densities clouds share more similarities with the isolated sources than the protostars of very dense clusters. If Class 0 protostars with few sources around and isolated Class 0 objects are comparable in the HDO/H$_2$O ratio, it would mean that there is little water reprocessing from the Class 0 to Class I protostellar stage.

An outstanding question is whether the $\alpha$/Fe bimodality exists in disk galaxies other than in the Milky Way. Here we present a bimodality using our state-of-the-art galactic chemical evolution models that can explain various observations in the Andromeda Galaxy (M31) disks, namely, elemental abundances both of planetary nebulae, and of red-giant branch stars recently observed with the James Webb Space Telescope. We find that in M31 a high-$\alpha$ thicker-disk population out to 30 kpc formed by more intense initial star burst than in the Milky Way. We also find a young low-$\alpha$ thin disk within 14 kpc, which is formed by a secondary star formation M31 underwent about 2-4.5 Gyr ago, probably triggered by a wet merger. In the outer disk, however, the planetary nebula observations indicate a slightly higher-$\alpha$ young ($\sim$2.5 Gyr) population at a given metallicity, possibly formed by secondary star formation from almost pristine gas. Therefore, an $\alpha$/Fe bimodality is seen in the inner disk ($<$14 kpc), while only a slight $\alpha$/Fe offset of the young population is seen in the outer disk ($>$18 kpc). The appearance of the $\alpha$/Fe bimodality depends on the merging history at various galactocentric radii, and wide-field multi-object spectroscopy is required for unveiling the history of M31.

We draw attention to the bright galaxies that do not show a bar in their structure but have a flocculent spiral structure. Using the THINGS' and HERACLES' kinematic data for four barless galaxies (NGC~2841, NGC~3512, NGC~5055, NGC~7331) we built their mass models including dark halos. We concluded that the fraction of the dark matter does not exceed 50\% within the optical radii of the galaxies. This is too little to explain the lack of a bar in these galaxies. In an attempt to understand the featureless structure of these galaxies we constructed several $N$-body models with an initially reduced content of dark matter. We concluded that, in addition to the low mass of the dark halo, the decisive factor that leads to a barless disc is the start from an initially unstable state. An isolated dynamically cold disc (with the Toomre parameter $Q < 0.5$) settled into rotational equilibrium passes trough the short stage of violent instability with fragmentation and formation of stellar clumps. After that, it evolves passively and ends up with a featureless structure. We assume that the barless flocculent galaxies studied in the present work may be descendants of galaxies at high redshifts with rotation curves which are consistent with the high mass fraction of baryons relative to the total dark matter halo.

Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range 7, 3, and 2 mm. We have modelled observations assuming LTE and LVG approximation. We have also used the code Nautilus to study the time evolution of sulphur species. We have detected 20 sulphur species with a total gas-phase S abundance similar to that found in the envelopes of other Class 0 objects, but with significant differences in the abundances between sulphur carbon chains and sulphur molecules containing oxygen and nitrogen. Our results highlight the nature of B335 as a source especially rich in sulphur carbon chains unlike other Class 0 sources. The low presence or absence of some molecules, such as SO and SO+, suggests a chemistry not particularly influenced by shocks. We, however, detect a large presence of HCS+ that, together with the low rotational temperatures obtained for all the S species (<15 K), reveals the moderate or low density of the envelope of B335. We also find that observations are better reproduced by models with a sulphur depletion factor of 10 with respect to the sulphur cosmic elemental abundance. The comparison between our model and observational results for B335 reveals an age of 10$^4$$<$t$<$10$^5$ yr, which highlights the particularly early evolutionary stage of this source. B335 presents a different chemistry compared to other young protostars that have formed in dense molecular clouds, which could be the result of accretion of surrounding material from the diffuse cloud onto the protostellar envelope of B335. In addition, the analysis of the SO2/C2S, SO/CS, and HCS+/CS ratios within a sample of prestellar cores and Class 0 objects show that they could be used as good chemical evolutionary indicators of the prestellar to protostellar transition.

We present deep \textit{XMM-Newton}, Karl Jansky Very Large Array, and upgraded Giant Metrewave Radio Telescope observations of Abell 746, a cluster that hosts a plethora of diffuse emission sources that provide evidence for the acceleration of relativistic particles. Our new \textit{XMM-Newton} images reveal a complex morphology of the thermal gas with several substructures. We observe an asymmetric temperature distribution across the cluster: the southern regions exhibit higher temperatures, reaching $\sim$9\,keV, while the northern regions have lower temperatures ($\rm \leq4\,keV$), likely due to a complex merger. We find evidence of four surface brightness edges, of which three are merger-driven shock fronts. Combining our new data with the published LOw-Frequency ARray observations has unveiled the nature of diffuse sources in this system. The bright northwest relic shows thin filaments and high degree of polarization with aligned magnetic field vectors. We detect a density jump, aligned with the fainter relic to the north. To the south, we detect high-temperature regions, consistent with shock-heated regions and density jump coincident with the northern tip of the southern radio structure. Its integrated spectrum shows a high-frequency steepening. Lastly, we find that the cluster hosts large-scale radio halo emission. The comparison of the thermal and nonthermal emission reveals an anticorrelation between the bright radio and X-ray features at the center. Our findings suggest that Abell 746 is a complex system that involves multiple mergers.

Airborne dust plays an active role in determining the thermal structure and chemical composition of the present-day atmosphere of Mars and possibly the planet's climate evolution over time through radiative--convective and cloud microphysics processes. Thus, accurate measurements of the distribution and variability of dust are required. Observations from the Mars Global Surveyor/Thermal Emission Spectrometer Mars Mars Reconnaissance Orbiter/Mars Climate Sounder and Mars Express/Fourier Transform Spectrometer and the Curiosity Rover have limited capability to measure dust. We show that spacecraft occultation of the Martian atmosphere at far-infrared frequencies between 1 and 10 THz can provide the needed global and temporal data on atmospheric dust by providing co-located measurements of temperature and dust opacity from the top of the atmosphere all the way down to the surface. In addition, spacecraft occultation by a small-satellite constellation could provide global measurements of the development of dust storms.

The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models. These emulators are defined with artificial neural networks (ANNs) designed and trained to address the specificities inherent in ISM models. Indeed, such models often predict many observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities. We propose applying five strategies to address these characteristics: 1) an outlier removal procedure; 2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; 3) a dimension reduction technique that enables to adequately size the network architecture; 4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and 5) a dense architecture to ease the learning of simple relations. We compare the proposed ANNs with standard classes of interpolation methods to emulate the Meudon PDR code, a representative ISM numerical model. Combinations of the proposed strategies outperform all interpolation methods by a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. This work will enable efficient inferences on wide-field multiline observations of the ISM.

In the outskirts of galaxy clusters, entropy profiles measured from X-ray observations of the hot intracluster medium (ICM) drops off unexpectedly. One possible explanation for this effect is gas clumping, where pockets of cooler and denser structures within the ICM are present. Current observatories are unable to directly detect these hypothetical gas clumps. One of the science drivers of the proposed STAR-X observatory is to resolve these or similar structures. Its high spatial resolution, large effective area, and low instrumental background make STAR-X ideal for directly detecting and characterizing clumps and diffuse emission in cluster outskirts. The aim of this work is to simulate observations of clumping in clusters to determine how well STAR-X will be able to detect clumps, as well as what clumping properties reproduce observed entropy profiles. This is achieved by using yt, pyXSIM, SOXS, and other tools to inject ideally modeled clumps into three-dimensional models derived from actual clusters using their observed profiles from other X-ray missions. Radial temperature and surface brightness profiles are then extracted from mock observations using concentric annuli. We find that in simulated observations for STAR-X, a parameter space of clump properties exists where gas clumps can be successfully identified using wavdetect and masked, and are able to recover the true cluster profiles. This demonstrates that STAR-X could be capable of detecting substructure in the outskirts of nearby clusters and that the properties of both the outskirts and the clumps will be revealed.

We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities ($\sim10^{17-18}$ erg s$^{-1}$ Hz$^{-1}$) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant ($>$ 50 pc) radio emitting M dwarfs.

We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-matter particle mass $164\,M_{\odot}$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5-9.5}\,M_{\odot}$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modeled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/M_{\odot}\simeq7.6$) causes the inner cusp in the haloes density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major $\langle q \rangle$ and intermediate-to-major $\langle s \rangle$ axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at $z=6.14$, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass.

We present a study of the physical and chemical properties of the Extended Green Object (EGO) G19.01$-$0.03 using sub-arcsecond angular resolution Atacama Large Millimeter/submillimeter Array (ALMA) 1.05mm and Karl G. Jansky Very Large Array (VLA) 1.21cm data. G19.01$-$0.03 MM1, the millimetre source associated with the central massive young stellar object (MYSO), appeared isolated and potentially chemically young in previous Submillimeter Array observations. In our $\sim0.4''$-resolution ALMA data, MM1 has four low-mass millimetre companions within 0.12pc, all lacking maser or outflow emission, indicating they may be prestellar cores. With a rich ALMA spectrum full of complex organic molecules, MM1 does not appear chemically young, but has molecular abundances typical of high-mass hot cores in the literature. At the 1.05mm continuum peak of MM1, $\mathrm{N}(\mathrm{CH}_{3}\mathrm{OH})=(2.22\pm0.01)\times10^{18}$cm$^{-2}$ and $T_{\mathrm{ex}} = 162.7\substack{+0.3 \\ -0.5}$K based on pixel-by-pixel Bayesian analysis of LTE synthetic methanol spectra across MM1. Intriguingly, the peak CH$_{3}$OH $T_{\mathrm{ex}}=165.5\pm0.6$ K is offset from MM1's millimetre continuum peak by $0.22''\sim880$au, and a region of elevated CH$_{3}$OH $T_{\mathrm{ex}}$ coincides with free-free VLA 5.01cm continuum, adding to the tentative evidence for a possible unresolved high-mass binary in MM1. In our VLA 1.21cm data, we report the first NH$_{3}$(3,3) maser detections towards G19.01$-$0.03, along with candidate 25GHz CH$_{3}$OH $5(2,3)-5(1,4)$ maser emission; both are spatially and kinematically coincident with 44GHz Class I CH$_{3}$OH masers in the MM1 outflow. We also report the ALMA detection of candidate 278.3GHz Class I CH$_{3}$OH maser emission towards this outflow, strengthening the connection of these three maser types to MYSO outflows.

Fifth forces are ubiquitous in modified theories of gravity. In this paper, we analyze their effect on the Cepheid-calibrated cosmic distance ladder, specifically with respect to the inferred value of the Hubble constant ($H_0$). We consider a variety of effective models where the strength, or amount of screening, of the fifth force is estimated using proxy fields related to the large-scale structure of the Universe. For all models considered, the local distance ladder and the Planck value for $H_0$ agrees with a probability $\gtrsim 20 \, \%$, relieving the tension compared to the concordance model with data being excluded at $99 \, \%$ confidence. The alleviated discrepancy comes partially at the cost of an increased tension between distance estimates from Cepheids and the tip of the red-giant branch (TRGB). Demanding also that the consistency between Cepheid and TRGB distance estimates is not impaired, some fifth force models can still accommodate the data with a probability $\gtrsim 20 \, \%$. This provides incentive for more detailed investigations of fundamental theories on which the effective models are based, and their effect on the Hubble tension.

We present the first optical study of the supernova remnant (SNR) G7.7-3.7, with the aim of determining its evolutionary phase since it has been suggested to be the remnant of SN 386 AD. We obtained narrow-band images in the filters H$\alpha$ + [NII], H$\beta$, [OIII], [SII] that revealed faint optical emission in the southern region of the SNR consisting of two filaments elongated in the east-west direction aligned with the X-ray emitting region of the remnant. The filaments were seen in H$\alpha$ + [NII], [OIII] images and marginally in the [SII] images, with a non-detection in H$\beta$. Long-slit spectroscopy of three regions along one filament revealed large ratios of [SII] / H$\alpha$ = (1.6-2.5), consistent with that expected for a shock-heated SNR. The [SII] doublet ratio observed in two of the regions implies an upper limit for the electron density of the gas, with estimates falling below 400 cm$^{-3}$ and 600 cm$^{-3}$ in the respective areas. We discuss potential physical mechanisms that formed the observed optical filaments and we suggest that most likely they resulted by a collision of the SNR with a dense circumstellar shell lying at the southern region of the remnant.

We study the coupling of hydrodynamics and reactions in simulations of the double detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator / Strang splitting vs.\ the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and the nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation timestep to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

Most of the Active Galactic Nuclei (AGN) are radio-quiet (RQ) and, differently from radio-loud (RL) AGN, do not show signature of large-scale and powerful jets. The physical origin of their radio emission remains then broadly unclear. The observation of flat/inverted radio spectra at GHz frequencies seems to support however the presence of an unresolved synchrotron self-absorbed region in the close environment of the supermassive black hole. Its size could be as small as that of the X-ray corona. Since synchrotron self absorption decreases strongly with frequency, these sources need to be observed in the millimetric (mm) domain. We report here a 12h simultaneous mm-X-ray observation of the RQ AGN MCG+08-11-11 by NOEMA and NuSTAR. The mm flux shows a weak but clear increase along the pointing with a fractional variability of $2.0\pm0.1$\%. The 3-10 keV flux of NuSTAR also increases and shows a fractional variability of $7.0\pm1.5$\%. A structure function analysis shows a local maximum in the mm light curve corresponding to 2-3\% of variability on timescale of $\sim2\times10^4$ seconds (100-300 $R_g$ light crossing time). Assuming an optically thick mm emitting medium, this translates into an upper limit of its size of $\sim$1300 $R_g$. The observation of fast variability in radio-mm and X-ray wavelengths, as well as a similar variability trend, well support the mm emission to be emitted by a region close, and potentially related to, the X-ray corona like an outflow/weak jet.

We aim to understand what drives the IRX-\beta dust attenuation relation at intermediate redshift (0.5 < z < 0.8) in star-forming galaxies. We investigate the role of various galaxy properties in shaping this observed relation. We use robust [O ii] {\lambda}3727, [O iii] {\lambda}{\lambda}4959, 5007, and H\beta line detections of our statistical sample of 1049 galaxies to estimate the gas-phase metallicities. We derive key physical properties that are necessary to study galaxy evolution, such as the stellar masses and the star formation rates, using the spectral energy distribution fitting tool CIGALE. Equivalently, we study the effect of galaxy morphology (mainly the S\'ersic index n and galaxy inclination) on the observed IRX-\beta scatter. We also investigate the role of the environment in shaping dust attenuation in our sample. We find a strong correlation of the IRX-\beta relation on gas-phase metallicity in our sample, and also strong correlation with galaxy compactness characterized by the S\'ersic indexes. Correlations are also seen with stellar masses, specific star formation rates and the stellar ages of our sources. Metallicity strongly correlates with the IRX-\beta scatter, this also results from the older stars and higher masses at higher beta values. Galaxies with higher metallicities show higher IRX and higher beta values. The correlation with specific dust mass strongly shifts the galaxies away from the IRX-\beta relation towards lower \b{eta} values. We find that more compact galaxies witness a larger amount of attenuation than less compact galaxies. There is a subtle variation in the dust attenuation scatter between edge-on and face-on galaxies, but the difference is not statistically significant. Galaxy environments do not significantly affect dust attenuation in our sample of star-forming galaxies at intermediate redshift.

In this paper, we investigate the nonprincipal axis (NPA) rotational state of 1I/`Oumuamua -- the first interstellar object discovered traversing the inner Solar System -- from its photometric light curve. Building upon Mashchenko (2019), we develop a model which incorporates NPA rotation and {Sun-induced, time-varying} outgassing torques to generate synthetic light curves of the object. The model neglects tidal forces, which are negligible compared to outgassing torques over the distances that `Oumuamua was observed. We implement an optimization scheme that incorporates the NPA rotation model to calculate the initial rotation state of the object. We find that an NPA rotation state with an average period of $\langle P \rangle\simeq7.34$ hr best reproduces the photometric data. The discrepancy between this period and previous estimates is due to continuous period modulation induced by outgassing torques in the rotational model, {as well as different periods being used}. The best fit to the October 2017 data does not reproduce the November 2017 data (although the later measurements are too sparse to fit). The light curve is consistent with no secular evolution of the angular momentum, somewhat in tension with the empirical correlations between nuclear spin-up and cometary outgassing. The complex rotation of `Oumuamua may be {the result of primordial rotation about the smallest principal axis} if (i) the object experienced hypervolatile outgassing and (ii) our idealized outgassing model is accurate.

Theoretical models indicate that photoevaporative and magnetothermal winds play a crucial role in the evolution and dispersal of protoplanetary disks and affect the formation of planetary systems. However, it is still unclear what wind-driving mechanism is dominant or if both are at work, perhaps at different stages of disk evolution. Recent spatially resolved observations by Fang et al. (2023) of the [OI] 6300 Angstrom spectral line, a common disk wind tracer, in TW Hya revealed that about 80% of the emission is confined to the inner few au of the disk. In this work, we show that state-of-the-art X-ray driven photoevaporation models can reproduce the compact emission and the line profile of the [OI] 6300 Angstrom line. Furthermore, we show that the models also simultaneously reproduce the observed line luminosities and detailed spectral profiles of both the [OI] 6300 Angstrom and the [NeII] 12.8 micron lines. While MHD wind models can also reproduce the compact radial emission of the [OI] 6300 Angstrom line, they fail to match the observed spectral profile of the [OI] 6300 Angstrom line and underestimate the luminosity of the [NeII] 12.8 micron line by a factor of three. We conclude that, while we cannot exclude the presence of an MHD wind component, the bulk of the wind structure of TW Hya is predominantly shaped by a photoevaporative flow.

Fast radio bursts (FRBs) are millisecond-duration extragalactic radio transients. They apparently fall into repeaters and non-repeaters. However, such a classification has lacked a motivation on the physical picture. Here we propose a unified geometric model to distinguish between the repeaters and non-repeaters, in which the quasi-tangential (QT) propagation effect within the magnetospheric polar cap of a neutron star is considered. In this model, the non-repeaters arise from the sources whose emitting region has a smaller impact angle with respect to the magnetic axis, while the repeaters come from the sources whose emitting region has a larger impact angle. The observational discriminant polarization properties between the repeaters and non-repeaters are an important clue to verifying this unified geometric model since the polarization is sensitive to the QT propagation effect. Moreover, our model effectively explains all of the other discriminant properties, including bandwidth, duration, peak luminosity, energy, brightness temperature, time-frequency downward drifting, and repetition rate, providing compelling evidence for the magnetospheric origin of FRBs.

Young supernova remnants (SNRs) strongly modify surrounding magnetic fields, which in turn play an essential role in accelerating cosmic rays (CRs). X-ray polarization measurements probe magnetic field morphology and turbulence at the immediate acceleration site. We report the X-ray polarization distribution in the northeastern shell of SN1006 from a 1 Ms observation with the Imaging X-ray Polarimetry Explorer (IXPE). We found an average polarization degree of $22.4\pm 3.5\%$ and an average polarization angle of $-45.4\pm 4.5^\circ$ (measured on the plane of the sky from north to east). The X-ray polarization angle distribution reveals that the magnetic fields immediately behind the shock in the northeastern shell of SN 1006 are nearly parallel to the shock normal or radially distributed, similar to that in the radio observations, and consistent with the quasi-parallel CR acceleration scenario. The X-ray emission is marginally more polarized than that in the radio band. The X-ray polarization degree of SN 1006 is much larger than that in Cas A and Tycho, together with the relatively tenuous and smooth ambient medium of the remnant, favoring that CR-induced instabilities set the turbulence in SN 1006 and CR acceleration is environment-dependent.

There is a desire to observe the sun's poles to further deepen our understanding of solar activity. However, because of the large speeds needed to perform out-of-ecliptic plane maneuvers, chemical and electric rocket propulsion mechanisms have been proven to be costly and impractical, leaving alternative space technology systems like solar sails to be considered for these applications. In this paper, we study the possibility of using a solar sail as a probe observing the sun. We design and optimize the trajectories of the solar sail probe through the surface constraint approach, with the assumption that the sail moves on a displaced spherical surface. We first review the surface constraint approach, focusing on its important assumptions and limitations. Then, we solve and obtain a family of radial and azimuthal trajectory equations by choosing the correct constraint equation. We characterize the trajectories based on the functional dependence of the sail's orientation with the polar angle. Finally, we determine the trajectories of the probe that will give us the minimum flight time. Results show that increasing the number of mission stages decreases the total flight time, with minimal changes in the sail's radial and polar velocities. Furthermore, changing the functional dependence of the clock angle resets the azimuthal velocity, making the sail not reverse its direction and directly approach the sun along the spherical surface.

In this letter, we report the discovery of a fast neutral hydrogen outflow in SDSS J145239.38+062738.0, a merging radio galaxy containing an optical type I active galactic nuclei (AGN). This discovery was made through observations conducted by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) using redshifted 21-cm absorption. The outflow exhibits a blueshifted velocity likely up to $\sim-1000\,\rm km\,s^{-1}$ with respect to the systemic velocity of the host galaxy with an absorption strength of $\sim -0.6\,\rm mJy\,beam^{-1}$ corresponding to an optical depth of 0.002 at $v=-500\,\rm km\,s^{-1}$. The mass outflow rate ranges between $2.8\times10^{-2}$ and $3.6\, \rm M_\odot \, yr^{-1}$, implying an energy outflow rate ranging between $4.2\times10^{39}$ and $9.7\times10^{40}\rm\,erg\,s^{-1}$, assuming 100 K $<T_{\rm s}<$ 1000 K. Plausible drivers of the outflow include the star bursts, the AGN radiation, and the radio jet, the last of which is considered the most likely culprit according to the kinematics. By analysing the properties of the outflow, the AGN, and the jet, we find that if the HI outflow is driven by the AGN radiation, the AGN radiation seems not powerful enough to provide negative feedback whereas the radio jet shows the potential to provide negative feedback. Our observations contribute another example of a fast outflow detected in neutral hydrogen, as well as demonstrate the capability of FAST in detecting such outflows.

Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.07" (10 au), and in the 12CO, C18O, 13CO (J=2-1), and SO (J_N = 6_5-5_4) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks (eDisk). The continuum emission traces a non-axisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position-velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ~0.12-0.17 Msun and ~13-19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r~28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very low-mass protostars.

Comet 157P is a faint object with a history of being prone to unfortunate situations, circumstances, and/or coincidences. Several weeks after its 1978 discovery the comet disappeared and remained lost nonstop for twenty five years. Rediscovered in 2003 as a new comet, it was about 500 times brighter than in 1978, caught apparently in one of its outbursts. The comet was not detected 200 days after its 2016 perihelion, being fainter than mag 20, but 80 days later it was mag 16 and gradually fading back to mag 20 over a period of four months. The comet did not miss the opportunity to have a close encounter with Jupiter, having approached it to less than 0.3 AU on 2020 February 10. The 2017 outburst or surge of activity appears to have accompanied an event of nuclear fragmentation. The birth of a second companion is dated to the months following the Jupiter encounter. The series of weird episodes culminated near the 2022 perihelion, when one companion brightened to become observable for two weeks and after another two weeks the other flared up to be seen for the next two weeks. Unnoticed, this incredible coincidence fooled some experts into believing that a single object, designated 157P-B, was involved, even though its orbit left large residuals. I now offer representative fragmentation solutions for the two companions, the mean residuals amounting to +/-0".4 and +/-1".0, respectively.

The ubiquity of information transmission via molecular communication between cells is comprehensively documented on Earth; this phenomenon might even have played a vital role in the origin(s) and early evolution of life. Motivated by these considerations, a simple model for molecular communication entailing the diffusion of signaling molecules from transmitter to receiver is elucidated. The channel capacity $C$ (maximal rate of information transmission) and an optimistic heuristic estimate of the actual information transmission rate $\mathcal{I}$ are derived for this communication system; the two quantities, especially the latter, are demonstrated to be broadly consistent with laboratory experiments and more sophisticated theoretical models. The channel capacity exhibits a potentially weak dependence on environmental parameters, whereas the actual information transmission rate may scale with the intercellular distance $d$ as $\mathcal{I} \propto d^{-4}$ and could vary substantially across settings. These two variables are roughly calculated for diverse astrobiological environments, ranging from Earth's upper oceans ($C \sim 3.1 \times 10^3$ bits/s; $\mathcal{I} \sim 4.7 \times 10^{-2}$ bits/s) and deep sea hydrothermal vents ($C \sim 4.2 \times 10^3$ bits/s; $\mathcal{I} \sim 1.2 \times 10^{-1}$ bits/s) to the hydrocarbon lakes and seas of Titan ($C \sim 3.8 \times 10^3$ bits/s; $\mathcal{I} \sim 2.6 \times 10^{-1}$ bits/s).

One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more than 27,000 galaxy groups and clusters at $0.1 < z_{\rm cl} < 1.4$ with more than 1 million member galaxies selected from the Subaru Hyper Suprime-Cam. By combining various AGN selection methods based on infrared (IR), radio, and X-ray data, we identify 2,688 AGN. We find that (i) $f_{\rm AGN}$ increases with $z_{\rm cl}$ and (ii) $f_{\rm AGN}$ decreases with $R/R_{\rm 200}$. The main contributors to the rapid increase of $f_{\rm AGN}$ towards high-$z$ and cluster center are IR- and radio-selected AGN, respectively. Those results indicate that the emergence of the AGN population depends on the environment and redshift, and galaxy groups and clusters at high-$z$ play an important role in AGN evolution. We also find that cluster-cluster mergers may not drive AGN activity in at least the cluster center, while we have tentative evidence that cluster-cluster mergers would enhance AGN activity in the outskirts of (particularly massive) galaxy clusters.

Thermal protection systems are a critical component of planetary exploration, enabling probes to enter the atmosphere and perform in-situ measurements. The aero-thermal conditions encountered during entry are destination and vehicle dependent, ranging from relatively benign conditions at Mars and Titan, to extreme conditions at Venus and Jupiter. The thermal protection system is a single-point-of-failure for both entry probe and aerocapture missions, and hence must be qualified using ground tests to ensure mission success. The high density Carbon-Phenolic which was used in the Galileo and the Pioneer Venus missions is no longer available due to the lack of the manufacturing base for its raw materials. To address the need for Venus and outer planet missions, NASA has developed the Heatshield for Extreme Environment Entry Technology (HEEET). The present study uses the Aerocapture Mission Analysis Tool (AMAT) to perform a comparative study of the thermal protection system requirements for various planetary destinations and the applicability of HEEET for future entry and aerocapture missions. The heat rate and stagnation pressure for aerocapture is significantly less compared to probe entry. The large heat loads during aerocapture present a challenge, but HEEET is capable of sustaining large heat loads within a reasonable TPS mass fraction.

We test the consistency of dark matter velocity distributions obtained from dark matter-only numerical simulations with analytic predictions, using the publicly available Via Lactea 2 dataset as an example. We find that, well inside the scale radius, the velocity distribution obtained from numerical simulation is consistent with a function of a single integral of motion -- the energy -- and moreover is consistent with the result obtained from Eddington inversion. This indicates that the assumptions underlying the analytic result, namely, spherical symmetry, isotropy, and a static potential, are sufficiently accurate to govern the coarse properties of the velocity distribution in the inner regions of the halo. We discuss implications for the behavior of the high-velocity tail of the distribution, which can dominate dark matter annihilation from a $p$- or $d$-wave state.

Type II supernovae represent the most common stellar explosions in the Universe, for which the final stage evolution of their hydrogen-rich massive progenitors towards core-collapse explosion are elusive. The recent explosion of SN 2023ixf in a very nearby galaxy, Messier 101, provides a rare opportunity to explore this longstanding issue. With the timely high-cadence flash spectra taken within 1-5 days after the explosion, we can put stringent constraints on the properties of the surrounding circumstellar material around this supernova. Based on the rapid fading of the narrow emission lines and luminosity/profile of $\rm H\alpha$ emission at very early times, we estimate that the progenitor of SN 2023ixf lost material at a mass-loss rate $\dot{\rm M} \approx 6 \times 10^{-4}\, \rm M_{\odot}\,yr^{-1}$ over the last 2-3 years before explosion. This close-by material, moving at a velocity $v_{\rm w} \approx 55\rm \, km\,s^{-1}$, accumulates a compact CSM shell at the radius smaller than $7 \times 10^{14}$ cm from the progenitor. Given the high mass-loss rate and relatively large wind velocity presented here, together with the pre-explosion observations made about two decades ago, the progenitor of SN 2023ixf could be a short-lived yellow hypergiant that evolved from a red supergiant shortly before the explosion.

We report on the serendipitous discovery of an unprecedented interaction between the radio lobe of a radio galaxy and a spiral galaxy. The discovery was made thanks to LOFAR observations at 144 MHz of the galaxy cluster Abell 160 ($z=0.04317$) provided by the LOFAR Two-metre Sky Survey. The new low-frequency observations revealed that one of the radio plumes of the central galaxy GIN 049 overlaps with the spiral galaxy JO36. Previous studies carried out with MUSE revealed that the warm ionized gas in the disk of JO36, traced by the H$\alpha$ emission, is severely truncated with respect to the stellar disk. We further explore this unique system by including new uGMRT observations at 675 MHz to map the spectral index. The emerging scenario is that JO36 has interacted with the radio plume in the past 200-500 Myr. The encounter resulted in a positive feedback event for JO36 in the form of a star formation rate burst of $\sim14$ $M_\odot$ yr$^{-1}$. In turn, the galaxy passage left a trace in the radio-old plasma by re-shaping the old relativistic plasma via magnetic draping.

The existing matched filtering method for gravitational wave (GW) search relies on a template bank. The computational efficiency of this method scales with the size of the templates within the bank. Higher-order modes and eccentricity will play an important role when third-generation detectors operate in the future. In this case, traditional GW search methods will hit computational limits. To speed up the computational efficiency of GW search, we propose the utilization of a deep learning (DL) model bank as a substitute for the template bank. This model bank predicts the latent templates embedded in the strain data. Combining an envelope extraction network and an astrophysical origin discrimination network, we realize a novel GW search framework. The framework can predict the GW signal's matched filtering signal-to-noise ratio (SNR). Unlike the end-to-end DL-based GW search method, our statistical SNR holds greater physical interpretability than the $p_{score}$ metric. Moreover, the intermediate results generated by our approach, including the predicted template, offer valuable assistance in subsequent GW data processing tasks such as parameter estimation and source localization. Compared to the traditional matched filtering method, the proposed method can realize real-time analysis. The minor improvements in the future, the proposed method may expand to other scopes of GW search, such as GW emitted by the supernova explosion.

Magnetic fields are fundamental to the evolution of galaxies, playing a key role in the astrophysics of the interstellar medium and star formation. Large-scale ordered magnetic fields have been mapped in the Milky Way and nearby galaxies, but it is not known how early in the Universe such structures form. Here we report the detection of linearly polarized thermal emission from dust grains in a strongly lensed, intrinsically luminous galaxy that is forming stars at a rate more than a thousand times that of the Milky Way at redshift 2.6, within 2.5 Gyr of the Big Bang. The polarized emission arises from the alignment of dust grains with the local magnetic field. The median polarization fraction is of order one per cent, similar to nearby spiral galaxies. Our observations support the presence of a 5 kiloparsec-scale ordered magnetic field with a strength of around 500uG or lower, orientated parallel to the molecular gas disk. This confirms that such structures can be rapidly formed in galaxies, early in cosmic history.

For cosmic microwave background (CMB) polarization observations, calibration of detector polarization angles is essential. We have developed a fully remote controlled calibration system with a sparse wire grid that reflects linearly polarized light along the wire direction. The new feature is a remote-controlled system for regular calibration, which has not been possible in sparse wire grid calibrators in past experiments. The remote control can be achieved by two electric linear actuators that load or unload the sparse wire grid into a position centered on the optical axis of a telescope between the calibration time and CMB observation. Furthermore, the sparse wire grid can be rotated by a motor. A rotary encoder and a gravity sensor are installed on the sparse wire grid to monitor the wire direction. They allow us to achieve detector angle calibration with expected systematic error of $0.08^{\circ}$. The calibration system will be installed in small-aperture telescopes at Simons Observatory.

The Coronagraph Instrument of the Nancy Grace Roman Space Telescope (Roman Coronagraph) will be capable of both total intensity and polarization measurements of circumstellar disks. The polarimetric performance is impacted by polarization effects introduced by all mirrors before the Wollaston prisms. In this paper, we aim to characterize these effects for the Roman Coronagraph in bands 1 and 4 using the FALCO and PROPER packages. We simulate the effect of polarization aberrations that impact the polarimetric contrast and the instrumental polarization effects to study the polarimetric accuracy. We include spacecraft rolls, but leave out systematic camera noise. We find that polarimetric differential imaging (PDI) improves the contrast by a factor of six. The PDI contrast of $\sim 8 \times 10^{-11}$ is limited by polarized speckles from instrumental polarization effects and polarization aberrations. By injecting polarized companions with at various contrast levels and demodulating their polarimetric signal, we recover their source Stokes vector within 2%.

The future Habitable Worlds Observatory aims to characterize the atmospheres of rocky exoplanets around solar-type stars. The vector vortex coronagraph (VVC) is a main candidate to reach the required contrast of $10^{-10}$. However, the VVC requires polarization filtering and every observing band requires a different VVC. The triple-grating vector vortex coronagraph (tgVVC) aims to mitigate these limitations by combining multiple gratings that minimize the polarization leakage over a large spectral bandwidth. In this paper, we present laboratory results of a tgVVC prototype using the In-Air Coronagraphic Testbed (IACT) facility at NASA's Jet Propulsion Laboratory and the Space Coronagraph Optical Bench (SCoOB) at the University of Arizona Space Astrophysics Lab (UASAL). We study the coronagraphic performance with polarization filtering at 633 nm and reach a similar average contrast of $2 \times 10^{-8}$ between 3-18 $\lambda/D$ at the IACT, and $6 \times 10^{-8}$ between 3-14 $\lambda/D$ at SCoOB. We explore the limitations of the tgVVC by comparing the testbed results. We report on other manufacturing errors and ways to mitigate their impact.

We present a new gas-grain chemical model for the combined isotopic fractionation of carbon and nitrogen in molecular clouds, in which the isotope chemistry of carbon and nitrogen is coupled with a time-dependent description of spin-state chemistry. We updated the rate coefficients of some isotopic exchange reactions considered in the literature, and present here a set of new exchange reactions involving molecules substituted in $\rm ^{13}C$ and $\rm ^{15}N$ simultaneously. We apply the model to a series of zero-dimensional simulations representing a set of physical conditions across a prototypical prestellar core, exploring the deviations of the isotopic abundance ratios in the various molecules from the elemental isotopic ratios as a function of physical conditions and time. We find that the $\rm ^{12}C/^{13}C$ ratio can deviate from the elemental ratio by up to a factor of several depending on the molecule, and that there are highly time-dependent variations in the ratios. The $\rm HCN/H^{13}CN$ ratio, for example, can obtain values of less than 10 depending on the simulation time. The $\rm ^{14}N/^{15}N$ ratios tend to remain close to the assumed elemental ratio within $\sim$ ten per cent, with no clear trends as a function of the physical conditions. Abundance ratios between $\rm ^{13}C$-containing molecules and $\rm ^{13}C$+$\rm ^{15}N$-containing molecules show somewhat increased levels of fractionation due to the newly included exchange reactions, though still remaining within a few tens of per cent of the elemental $\rm ^{14}N/^{15}N$ ratio. Our results imply the existence of gradients in isotopic abundance ratios across prestellar cores, suggesting that detailed simulations are required to interpret observations of isotopically substituted molecules correctly, especially given that the various isotopic forms of a given molecule do not necessarily trace the same gas layers.

Temperature profiles of the hot galaxy cluster intracluster medium (ICM) have a complex non-linear structure that traditional parametric modelling may fail to fully approximate. For this study, we made use of neural networks, for the first time, to construct a data-driven non-parametric model of ICM temperature profiles. A new deconvolution algorithm was then introduced to uncover the true (3D) temperature profiles from the observed projected (2D) temperature profiles. An auto-encoder-inspired neural network was first trained by learning a non-linear interpolatory scheme to build the underlying model of 3D temperature profiles in the radial range of [0.02-2] R$_{500}$, using a sparse set of hydrodynamical simulations from the THREE HUNDRED PROJECT. A deconvolution algorithm using a learning-based regularisation scheme was then developed. The model was tested using high and low resolution input temperature profiles, such as those expected from simulations and observations, respectively. We find that the proposed deconvolution and deprojection algorithm is robust with respect to the quality of the data, the morphology of the cluster, and the deprojection scheme used. The algorithm can recover unbiased 3D radial temperature profiles with a precision of around 5\% over most of the fitting range. We apply the method to the first sample of temperature profiles obtained with XMM{\it -Newton} for the CHEX-MATE project and compared it to parametric deprojection and deconvolution techniques. Our work sets the stage for future studies that focus on the deconvolution of the thermal profiles (temperature, density, pressure) of the ICM and the dark matter profiles in galaxy clusters, using deep learning techniques in conjunction with X-ray, Sunyaev Zel'Dovich (SZ) and optical datasets.

We present MeerKAT HI line observations of the nearby interacting galaxy pair NGC 1512/1510. The MeerKAT data yield high-fidelity image sets characterised by an excellent combination of high angular resolution (~20") and and sensitivity (~0.08 Msun/pc^2), thereby offering the most detailed view of this well-studied system's neutral atomic hydrogen content, especially the HI co-located with the optical components of the galaxies. The stellar bulge and bar of NGC 1512 are located within a central HI depression where surface densities fall below 1 Msun/pc^2, while the galaxy's starburst ring coincides with a well-defined HI annulus delimited by a surface density of 3 Msun/pc^2. In stark contrast, the star-bursting companion, NGC 1510, has its young stellar population precisely matched to the highest HI over-densities we measure (~12.5 Msun/pc^2). The improved quality of the MeerKAT data warrants the first detailed measurements of the lengths and masses of the system's tidally-induced HI arms. We measure the longest of the two prominent HI arms to extend over ~27 kpc and to contain more than 30% of the system's total HI mass. We quantitatively explore the spatial correlation between HI and far-ultraviolet flux over a large range of HI mass surface densities spanning the outer disk. The results indicate the system's HI content to play an important role in setting the pre-conditions required for wide-spread, high-mass star formation. This work serves as a demonstration of the remarkable efficiency and accuracy with which MeerKAT can image nearby systems in HI line emission.

We perform a new general-relativistic viscous-radiation hydrodynamics simulation for supernova-like explosion associated with stellar core collapse of rotating massive stars to a system of a black hole and a massive torus paying particular attention to large-mass progenitor stars with the zero-age main-sequence mass of $M_\mathrm{ZAMS}=$20, 35, and 45$M_\odot$ of Ref.~\cite{Aguilera-Dena2020oct}. Assuming that a black hole is formed in a short timescale after the onset of the stellar collapse, the new simulations are started from initial data of a spinning black hole and infalling matter that self-consistently satisfy the constraint equations of general relativity. It is found that with a reasonable size of the viscous parameter, the supernova-like explosion is driven by the viscous heating effect in the torus around the black hole irrespective of the progenitor mass. The typical explosion energy and ejecta mass for the large-mass cases ($M_\mathrm{ZAMS}=35$ and $45M_\odot$) are $\sim 10^{52}$ erg and $\sim 5M_\odot$, respectively, with $^{56}$Ni mass larger than $0.15M_\odot$. These are consistent with the observational data of stripped-envelope and high-energy supernovae such as broad-lined type Ic supernovae. This indicates that rotating stellar collapses of massive stars to a black hole surrounded by a massive torus can be a central engine for high-energy supernovae. By artificially varying the angular velocity of the initial data, we explore the dependence of the explosion energy and ejecta mass on the initial angular momentum and find that the large explosion energy $\sim 10^{52}$ erg and large $^{56}$Ni mass $\geq 0.15M_\odot$ are possible only when a large-mass compact torus with mass $\gtrsim 1M_\odot$ is formed.

Axion-like particles (ALPs) could mix with photons in the presence of astrophysical magnetic fields, and result in oscillations in the high energy $\gamma$-ray spectra observed by experiments. In this work, we investigate the ALP-photon oscillation effect through the blazar Mrk 421 spectra of 15 periods observed by Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) and Fermi Large Area Telescope (Fermi-LAT). Compared with previous studies, we generate the mock data under the ALP hypothesis and apply the ${\rm CL_s}$ method to set constraints on the ALP parameters. This method is widely employed in high energy experiments and could avoid the possibility of excluding some parameter regions due to the fluctuation. We find that the ALP-photon coupling $g_{a\gamma}$ is constrained to be smaller than $\sim 2\times10^{-11}$ GeV$^{-1}$ for ALP mass ranging from $10^{-9}$ eV to $10^{-7}$ eV at a 95\% confidence level. The constraints obtained with the method based on the TS distribution under the null hypothesis, which is adopted in many previous astrophysical ALP studies, are also shown. Our results demonstrate that the joint constraints of all the periods from both methods are consistent. However, the latter method fails to provide constraints for some observation periods, whereas the ${\rm CL_s}$ method remains effective in such cases.

Over its 13 years of operation (1990 -- 2002), the Faint Object Camera (FOC) on board the Hubble Space Telescope (HST) observed 26 individual active galactic nuclei (AGNs) in ultraviolet (UV) imaging polarimetry. However, not all of the observations have been reduced and analyzed or set within a standardized framework. We plan to reduce and analyze the AGN observations that have been neglected in the FOC archives using a consistent, novel, and open-access reduction pipeline of our own. We then extend the method to the full AGN sample, thus leading to potential discoveries in the near future. We developed a new pipeline in Python that will be able to reduce all the FOC observations in imaging polarimetry in a homogeneous way. Most of the previously published reduced observations are dispersed throughout the literature, with the range of different analyses and approaches making it difficult to fully interpret the FOC AGN sample. By standardizing the method, we have enabled a coherent comparison among the different observational sets. In this first paper of a series exploring the full HST/FOC AGN sample, we present an exhaustively detailed account of how to properly reduce the observational data. Current progress in data-analysis is implemented in and has provided state-of-the-art UV polarimetric maps. We compare our new maps to the benchmark AGN case of NGC~1068 and successfully reproduce the main results previously published, while pushing the polarimetric exploration of this AGN futher, thanks to a finer resolution and a higher signal-to-noise ratio (S/N) than previously reported. We also present, for the first time, an optical polarimetric map of the radio-loud AGN IC~5063 and we examine the complex interactions between the AGN outflows and the surrounding interstellar medium (ISM).

In this paper, we investigate the process in which axion dark matter undergoes thermal friction, resulting in energy injection into dark radiation, with the aim of mitigating the Hubble tension and large-scale structure tension. In the early universe, this scenario led to a rapid increase in the energy density of dark radiation; in the late universe, the evolution of axion dark matter is similar to that of cold dark matter, with this scenario resembling decaying dark matter and serving to ease the large-scale structure tension. We employ cosmological observational data, including cosmic microwave background (CMB), baryon acoustic oscillation (BAO), supernova data (SNIa), $H_0$ measurement from SH0ES, and $S_8$ from the Dark Energy Survey Year-3 (DES), to study and analyze this model. Our results indicate that the thermal friction model offers partial alleviation of the large-scale structure tension, while its contribution on alleviating Hubble tension can be ignored. The new model constrained by the complete dataset yields the value of $S_8$ is $0.795\pm 0.011$ at 68\% confidence level, while the $\Lambda$CDM model yields a result of $0.8023\pm 0.0085$. In addition, when constrained solely by CMB, BAO, and SNIa data, the $\Lambda$CDM model exhibits a smaller $\chi^2_\mathrm{tot}$ statistical value. However, upon incorporating SH0ES and DES data, the new model exhibits a lower $\chi^2_\mathrm{tot}$ value, with a difference of -2.60 compared to the $\Lambda$CDM model.

We present tests of a new method to simultaneously estimate stellar population and emission line (EL) properties of galaxies out of S-PLUS photometry. The technique uses the AlStar code, updated with an empirical prior which greatly improves its ability to estimate ELs using only the survey's 12 bands. The tests compare the output of (noise-perturbed) synthetic photometry of SDSS galaxies to properties derived from previous full spectral fitting and detailed EL analysis. For realistic signal-to-noise ratios, stellar population properties are recovered to better than 0.2 dex in masses, mean ages, metallicities and $\pm 0.2$ mag for the extinction. More importantly, ELs are recovered remarkably well for a photometric survey. We obtain input $-$ output dispersions of 0.05--0.2 dex for the equivalent widths of $[\mathrm{O}\,\rm{II}]$, $[\mathrm{O}\,\rm{III}]$, H$\beta$, H$\alpha$, $[\mathrm{N}\,\rm{II}]$, and $[\mathrm{S}\,\rm{II}]$, and even better for lines stronger than $\sim 5$ $\mathring{A}$. These excellent results are achieved by combining two empirical facts into a prior which restricts the EL space available for the fits: (1) Because, for the redshifts explored here, H$\alpha$ and $[\mathrm{N}\,\rm{II}]$ fall in a single narrow band (J0660), their combined equivalent width is always well recovered, even when $[\mathrm{N}\,\rm{II}]$/H$\alpha$ is not. (2) We know from SDSS that $W_{H\alpha+[\mathrm{N}\,\rm{II}]}$ correlates with $[\mathrm{N}\,\rm{II}]$/H$\alpha$, which can be used to tell if a galaxy belongs to the left or right wings in the classical BPT diagnostic diagram. Example applications to integrated light and spatially resolved data are also presented, including a comparison with independent results obtained with MUSE-based integral field spectroscopy.

Edukoi is a software that aims to make interactive sonification suitable to convey and extract information. The program design is a modification of the software Herakoi, which sonifies images in real time mapping pitch to colour using a motion-aware approach for allowing users to interact with images through sound. The pitch-colour association of Hearkoi, albeit pleasing from the entertainment side, is not efficient for communicating specific information regarding colours and hues to listeners. Hence we modified it to create an instrument to be used by visually impaired and sighted children to explore images through sound and extract accurate information. We aim at building a flexible software that can be used in middle-schools for both art and science teaching. We tested its effectiveness using astronomical images, given the great fascination that astronomy always has on kids of all ages and backgrounds. Astronomy is also considered a very visual science, a characteristic that prevents students from learning this subject and having a related career. With this project we aim to challenge this belief and give to students the possibility to explore astronomical data through sound. Here we discuss our experiment, the choices we made regarding sound mappings, and what psychophysiological aspects we aim to evaluate to validate and improve Edukoi.

The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction ($f_{esc}$). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their $f_{esc}$ can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at $6 \leq z \leq 9$ from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average $6<z<9$ star-forming galaxies are compact in the rest-frame UV ($r_e \sim $ 0.4 kpc), are blue sources (UV-$\beta$ slope $\sim $ -2.17), and have a predicted $f_{esc}$ of about 0.13. A comparison of our results to models and predictions as well as an estimation of the ionizing budget suggests that low-mass galaxies with UV magnitudes fainter than $M_{1500} = -18$ that we currently do not characterize with JWST observations probably played a key role in the process of reionization.

The rotation and orientation of Mars is commonly described using two different sets of angles, namely the Euler angles wrt the Mars orbit plane and the right ascension, declination, and prime meridian location angles wrt the Earth equator at J2000 (as adopted by the IAU). We propose a formulation for both these sets of angles, which consists of the sum of a second degree polynomial and of periodic and Poisson series. Such a formulation is shown here to enable accurate (and physically sound) transformation from one set of angles to the other. The transformation formulas are provided and discussed in this paper. In particular, we point that the quadratic and Poisson terms are key ingredients to reach a transformation precision of 0.1 mas, even 30 years away from the reference epoch of the rotation model (e.g. J2000). Such a precision is required to accurately determine the smaller and smaller geophysical signals observed in the high-accuracy data acquired from the surface of Mars. In addition, we present good practices to build an accurate Martian rotation model over a long time span (30 years around J2000) or over a shorter one (e.g. lifetime of a space mission). We recommend to consider the J2000 mean orbit of Mars as the reference plane for Euler angles. An accurate rotation model should make use of up-to-date models for the rigid and liquid nutations, relativistic corrections in rotation, and polar motion induced by the external torque. Our transformation model and recommendations can be used to define the future IAU solution for the rotation and orientation of Mars using right ascension, declination, and prime meridian location. In particular, thanks to its quadratic terms, our transformation model does not introduce arbitrary and non-physical terms of very long period and large amplitudes, thus providing unbiased values of the rates and epoch values of the angles.

Ciii]{\lambda}{\lambda}1907, 1909 doublet emission line is a valuable tool for exploring early star-forming galaxies. It has been proposed as a potential alternative to the diminishing Ly{\alpha}{\lambda}1215.7 line for tracing galaxies during the Epoch of Reionization. In this study, we investigate the utility of the Ciii] line as a proxy for Ly{\alpha} in the reionization era by comparing the equivalent widths (EW) and velocity offsets of both emission lines. Our analysis focuses on star-forming galaxies at z \sim 3 - 4 from the VANDELS survey. We examined the spectra of 773 objects, identifying the rest-frame UV line Ciii]. Subsequently, we measured the EW of Ciii], Ly{\alpha}, and Heii. For objects displaying both Ciii] and Ly{\alpha} emission, the Ly{\alpha} velocity offsets was calculated. After removing 10 potential AGNs from our analysis, we detected Ciii] emission in 280/773 galaxies, with 139 receiving the highest confidence rating. The EW(Ciii]) had an average \sim 6 {\AA}, while EW(Ly{\alpha}) had an average \sim 18 {\AA}. Among the subset that showed both Ciii] and Ly{\alpha} (52/139), the average EW(Ciii]) was \sim 5 {\AA}, higher than those without Ly{\alpha} (EW(Ciii]) \sim 3 {\AA}). Additionally, all 52 galaxies in the Ciii] and Ly{\alpha} subset displayed a velocity shift ({\Delta}v_{Ly{\alpha}}), with average offset 533 km/s. This expanded dataset provides valuable insights, including a positive correlation between EW(Ciii]) and EW(Ly{\alpha}), confirming earlier findings. Furthermore, we report a promising anti-correlation between EW(Ciii]) and {\Delta}v_{Ly{\alpha}}, which may serve as a tool for inferring Ly{\alpha} properties and potentially detecting ionized bubbles at z \gt 6.

Jets of active galactic nuclei are potential accelerators of ultra high-energy cosmic rays. Supernovae can occur inside these jets and contribute to cosmic ray acceleration, particularly of heavy nuclei, but that contribution has been hardly investigated so far. We carried out a first dedicated exploration of the role of supernovae inside extragalactic jets in the production of ultra high-energy cosmic rays. We characterized the energy budget of supernova-jet interactions, and the maximum possible energies of the particles accelerated in those events, likely dominated by heavy nuclei. This allowed us to assess whether these interactions can be potential acceleration sites of ultra high-energy cosmic rays, or at least of their seeds. For that, we estimated the cosmic ray luminosity for different galaxy types, and compared the injection rate of cosmic ray seeds into the jet with that due to galactic cosmic ray entrainment. Since the supernova is fueled for a long time by the luminosity of the jet, the energy of a supernova-jet interaction can be several orders of magnitude greater than that of an isolated supernova. Thus, despite the low rate of supernovae expected to occur in the jet, they could still provide more seeds for accelerating ultra high-energy particles than cosmic ray entrainment from the host galaxy. Moreover, these interactions can create sufficiently efficient accelerators to be a source of cosmic rays with energies $\gtrsim 10$~EeV. Supernova-jet interactions can contribute significantly to the production of ultra high-energy cosmic rays, either directly by accelerating these particles themselves or indirectly by providing pre-accelerated seeds.

Based on the expected population of core collapse supernova remnants and the huge number of detected pulsars in the Galaxy, still representing only a fraction of the real population, pulsar wind nebulae are likely to constitute one of the largest classes of {extended} Galactic sources in many energy bands. For simple evolutionary reasons, the majority of the population is made of evolved systems, whose detection and identification are complicated by their reduced luminosity, the possible lack of X-ray emission (that fades progressively away with the age of the pulsar), and by their modified morphology with respect to young systems. Nevertheless they have gained renewed attention in recent years, following the detection of misaligned X-ray tails protruding from an increasing number of nebulae created by fast moving pulsars, and of extended TeV halos surrounding aged systems. Both these features are clear signs of an efficient escape of particles, with energy close to the maximum acceleration limit of the pulsar. Here we discuss the properties of those evolved systems and what we have understood about the process of particle escape, and the formation of observed features.

Recent radio observations and coincident neutrino detections suggest that some tidal disruption events (TDEs) exhibit late-time activities, relative to the optical emission peak, and these may be due to delayed outflows launched from the central supermassive black hole. We investigate the possibility that jets launched with a time delay of days to months, interact with a debris that may expand outwards. We discuss the effects of the time delay and expansion velocity on the outcomes of jet breakout and collimation. We find that a jet with an isotropic-equivalent luminosity of $\lesssim 5 \times 10^{45}\,{\rm erg/s}$ is likely to be choked for a delay time of $\sim 3$ months. We also study the observational signatures of such delayed choked jets. The jet-debris interaction preceding the breakout would lead to particle acceleration and the resulting synchrotron emission can be detected by current and near-future radio, optical and X-ray telescopes, and the expanding jet-driven debris could explain late-time radio emission. We discuss high-energy neutrino production in delayed choked jets, and the time delay can significantly alleviate the difficulty of the hidden jet scenario in explaining neutrino coincidences.

We use kinematic data of proper motions from Gaia of forty-two globular and open clusters from Large Magellanic Cloud (LMC) to explore the possibility of them having extragalactic origins. We find the difference between the proper motions of cluster stars and a surrounding patch of young LMC stars in each case. We find five globular clusters towards the north-east showing a high difference (> 0.11 mas/yr, or > 25 km/s). We also examine the statistical significance of this difference taking into account both measurement errors of cluster and surrounding stars as well as inherent dispersion of stellar motions in the local galactic environment. The five globular clusters (NGC 2005, NGC 2210, NGC 1978, Hodge 3 and Hodge 11) have mean proper motions that lie outside the 85% confidence interval of the mean of surrounding young stars, with a clear outlier (NGC 1978 outside 99.96% confidence) whose difference cannot be accounted for by statistical noise. A young cluster (NGC 2100) also fitting the criteria is ruled out owing to contrary evidence from literature. This indicates a possible interaction with a dwarf galaxy resulting in the accretion/disruption in path of the five globular clusters, or possibly one or more past merger(s) of smaller galaxy/galaxies with LMC from its north-eastern region. This direction also coincides with the location of Tarantula Nebula, suggesting the possibility of the interaction event or merger having triggered its star formation activity.

We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a quasar and a radio mode, with the radio mode tied to the jet energy production; (iii) a model tracking the development of black hole (BH) spins; (iv) more sophisticated modelling of environmental effects on satellite galaxies; and (v) automatic parameter exploration using Particle Swarm Optimisation. We focus on two timely research topics: the structural properties of galaxies and the quenching of massive galaxies. For the former, SHARK v2.0 is capable of producing a more realistic stellar size-mass relation with a plateau marking the transition from disk- to bulge-dominated galaxies, and scaling relations between specific AM and mass that agree well with observations. For the quenching of massive galaxies, SHARK v2.0 produces massive galaxies that are more quenched than the previous version, reproducing well the observed relations between star formation rate (SFR) and stellar mass, and specific SFR and BH mass at $z=0$. SHARK v2.0 produces a number density of massive-quiescent galaxies >1dex higher than the previous version, in good agreement with JWST observations at $z\le 5$; predicts a stellar mass function of passive galaxies in reasonably good agreement with observations at $0.5<z<5$; and environmental quenching to already be effective at $z=5$.

Optical spectra contain a wealth of information about the physical properties and formation histories of galaxies. Often though, spectra are too noisy for this information to be accurately retrieved. In this study, we explore how machine learning methods can be used to de-noise spectra and increase the amount of information we can gain without having to turn to sample averaging methods such as spectral stacking. Using machine learning methods trained on noise-added spectra - SDSS spectra with Gaussian noise added - we investigate methods of maximising the information we can gain from these spectra, in particular from emission lines, such that more detailed analysis can be performed. We produce a variational autoencoder (VAE) model, and apply it on a sample of noise-added spectra. Compared to the flux measured in the original SDSS spectra, the model values are accurate within 0.3-0.5 dex, depending on the specific spectral line and S/N. Overall, the VAE performs better than a principle component analysis (PCA) method, in terms of reconstruction loss and accuracy of the recovered line fluxes. To demonstrate the applicability and usefulness of the method in the context of large optical spectroscopy surveys, we simulate a population of spectra with noise similar to that in galaxies at $z = 0.1$ observed by the Dark Energy Spectroscopic Instrument (DESI). We show that we can recover the shape and scatter of the MZR in this "DESI-like" sample, in a way that is not possible without the VAE-assisted de-noising.

Pulse profiles of accreting millisecond pulsars can be used to determine neutron star (NS) parameters, such as their masses and radii, and therefore provide constraints on the equation of state of cold dense matter. Information obtained by the Imaging X-ray Polarimetry Explorer (IXPE) can be used to decipher pulsar inclination and magnetic obliquity, providing ever tighter constraints on other parameters. In this paper, we develop a new emission model for accretion-powered millisecond pulsars based on thermal Comptonization in an accretion shock above the NS surface. The shock structure was approximated by an isothermal plane-parallel slab and the Stokes parameters of the emergent radiation were computed as a function of the zenith angle and energy for different values of the electron temperature, the Thomson optical depth of the slab, and the temperature of the seed blackbody photons. We show that our Compton scattering model leads to a significantly lower polarization degree of the emitted radiation compared to the previously used Thomson scattering model. We computed a large grid of shock models, which can be combined with pulse profile modeling techniques both with and without polarization included. In this work, we used the relativistic rotating vector model for the oblate NS in order to produce the observed Stokes parameters as a function of the pulsar phase. Furthermore, we simulated the data to be produced by IXPE and obtained constraints on model parameters using nested sampling. The developed methods can also be used in the analysis of the data from future satellites, such as the enhanced X-ray Timing and Polarimetry mission.

In this work, we investigate the discovery potential of low-mass Galactic dark matter (DM) subhaloes for indirect searches of DM. We use data from the Via Lactea II (VL-II) N-body cosmological simulation, which resolves subhaloes down to $\mathcal{O}(10^4)$ solar masses and it is thus ideal for this purpose. First, we characterize the abundance, distribution and structural properties of the VL-II subhalo population in terms of both subhalo masses and maximum circular velocities. Then, we repopulate the original simulation with millions of subhaloes of masses down to about five orders of magnitude below the minimum VL-II subhalo mass (more than one order of magnitude in velocities). We compute subhalo DM annihilation astrophysical "J-factors" and angular sizes for the entire subhalo population, by placing the Earth at a random position but at the right galactocentric distance in the simulation. Thousands of these realizations are generated in order to obtain statistically meaningful results. We find that some nearby low-mass Galactic subhaloes, not massive enough to retain stars or gas, may indeed yield DM annihilation fluxes comparable to those expected from other, more massive and acknowledgeable DM targets like dwarf satellite galaxies. Typical angular sizes are of the order of the degree, thus subhaloes potentially appearing as extended sources in gamma-ray telescopes, depending on instrument angular resolution and sensitivity. Our work shows that low-mass Galactic subhaloes with no visible counterparts are expected to play a relevant role in current and future indirect DM search searches and should indeed be considered as excellent DM targets.

Cool astrophysical objects, such as (exo)planets, brown dwarfs, or asymptotic giant branch stars, can be strongly affected by condensation. Condensation does not only directly affect the chemical composition of the gas phase by removing elements but the condensed material also influences other chemical and physical processes in these object. This includes, for example, the formation of clouds in planetary atmospheres and brown dwarfs or the dust-driven winds of evolved stars. In this study we introduce FastChem Cond, a new version of the FastChem equilibrium chemistry code that adds a treatment of equilibrium condensation. Determining the equilibrium composition under the impact of condensation is complicated by the fact that the number of condensates that can exist in equilibrium with the gas phase is limited by a phase rule. However, this phase rule does not directly provide information on which condensates are stable. As a major advantage of FastChem Cond is able to automatically select the set stable condensates satisfying the phase rule. Besides the normal equilibrium condensation, FastChem Cond can also be used with the rainout approximation that is commonly employed in atmospheres of brown dwarfs or (exo)planets. FastChem Cond is available as open-source code, released under the GPLv3 licence. In addition to the C++ code, FastChem Cond also offers a Python interface. Together with the code update we also add about 290 liquid and solid condensate species to FastChem.

The growth of mega-constellations is rapidly increasing the number of rocket launches required to place new satellites in space. While Low Earth Orbit (LEO) broadband satellites help to connect unconnected communities and achieve the Sustainable Development Goals, there are also a range of negative environmental externalities, from the burning of rocket fuels and resulting environmental emissions. We present sustainability analytics for phase 1 of the three main LEO constellations including Amazon Kuiper (3,236 satellites), OneWeb (648 satellites), and SpaceX Starlink (4,425 satellites). In baseline scenarios over five years, we find a per subscriber carbon dioxide equivalent (CO$_2$eq) of 0.70$\pm$0.34 tonnes for Kuiper, 1.41$\pm$0.71 tonnes for OneWeb and 0.47$\pm$0.15 tonnes CO$_2$eq/subscriber for Starlink. However, in the worst-case emissions scenario these values increase to 3.02$\pm$1.48 tonnes for Kuiper, 1.7$\pm$0.71 tonnes for OneWeb and 1.04$\pm$0.33 tonnes CO$_2$eq/subscriber for Starlink, more than 31-91 times higher than equivalent terrestrial mobile broadband. Importantly, phase 2 constellations propose to increase the number of satellites by an order-of-magnitude higher, highlighting the pressing need to mitigate negative environmental impacts. Strategic choices in rocket design and fuel options can help to substantially mitigate negative sustainability impacts.

Measurements of neutron star masses, radii, and tidal deformability have direct connections to nuclear physics via the equation of state (EoS), which for the cold, catalyzed matter in neutron star cores is commonly represented as the pressure as a function of energy density. Microscopic models with exotic degrees of freedom display nontrivial structure in the speed of sound ($c_s$) in the form of first-order phase transitions and bumps, oscillations, and plateaus in the case of crossovers and higher-order phase transitions. We present a procedure based on Gaussian processes to generate an ensemble of EoSs that include nontrivial features. Using a Bayesian analysis incorporating measurements from X-ray sources, gravitational wave observations, and perturbative QCD results, we show that these features are compatible with current constraints. We investigate the possibility of a global maximum in $c_s$ that occurs within the densities realized in neutron stars -- implying a softening of the EoS and possibly an exotic phase in the core of massive stars -- and find that such a global maximum is consistent with, but not required by, current constraints.

The paper presents results of studying the kinematics of the ionized gas in the galaxy of the Large Magellanic Cloud type NGC 7292 obtained with the 2.5-m telescope of the Caucasian Mountain Observatory (CMO SAI MSU) and the 6-m BTA telescope of the Special Astrophysical Observatory (SAO RAS). Analysis of the velocity fields of the ionized and neutral hydrogen showed that the kinematic center of NGC 7292 located at the center of the bar, northwest of the photometric center of the galaxy (the southeastern end of the bar) previously taken as the center of NGC 7292. In addition to the circular rotation of the gas, the radial motions associated with the bar play a significant role in the kinematics of the disk. The observed perturbations of the gaseous-disk kinematics induced by the ongoing star formation do not exceed those caused by the bar. It is possible that part of the non-circular motions (at the southeastern end of the bar which is the brightest HII region) may be related to the effects of the capture of a dwarf companion or a gaseous cloud.

The recent detection of a stochastic gravitational wave background (SGWB) at nanohertz frequencies by pulsar timing arrays (PTAs) has sparked a flurry of interest. Beyond the standard interpretation that the progenitor is a network of supermassive black hole binaries, many exotic models have also been proposed, some of which can potentially offer a better fit to the data. We explore how the various connections between gravitational waves and CMB spectral distortions can be leveraged to help determine whether a SGWB was generated primordially or astrophysically. To this end, we present updated $k$-space window functions which can be used for distortion parameter estimation on enhancements to the primordial scalar power spectrum. These same enhancements can also source gravitational waves (GWs) directly at second order in perturbation theory, so-called scalar-induced GWs (SIGWs), and indirectly through the formation of primordial black holes (PBHs). We perform a mapping of scalar power spectrum constraints into limits on the GW parameter space of SIGWs for $\delta$-function features. We highlight that broader features in the scalar spectrum can explain the PTA results while simultaneously producing a spectral distortion (SD) within reach of future experiments. We additionally update PBH constraints from $\mu$- and $y$-type spectral distortions. Refined treatments of the distortion window functions widen existing SD constraints, and we find that a future CMB spectrometer could play a pivotal role in unraveling the origin of GWs imprinted at or below CMB anisotropy scales.

We identify a new mechanism that can lead to the destruction of small, close-in planetary satellites. If a small moon close to the planet has a sizable eccentricity and inclination, its ejecta that escape to planetocentric orbit would often re-impact with much higher velocity due to the satellite's and the fragment's orbits precessing out of alignment. If the impacts of returning ejecta result in net erosion, a runaway process can occur which may end in disruption of the satellite, and we term this process ``sesquinary catastrophe''. We expect the moon to re-accrete, but on an orbit with significantly lower eccentricity and inclination. We find that the large majority of small close-in moons in the Solar System, have orbits that are immune to sesquinary catastrophe. The exceptions include a number of resonant moonlets of Saturn for which resonances may affect the velocities of re-impact of their own debris. Additionally, we find that Neptune's moon Naiad (and to a lesser degree, Jupiter's Thebe) must have substantial internal strength, in line with prior estimates based on Roche limit stability. We also find that sesquinary instability puts important constraints on the plausible past orbits of Phobos and Deimos or their progenitors.

Bars are important drivers of galaxy evolution, influencing many physical processes and properties. Characterising bars is a difficult task, especially in large-scale surveys. In this work, we propose a novel morphological segmentation technique for determining bar lengths based on deep learning. We develop U-Nets capable of decomposing galaxy images into pixel masks highlighting the regions corresponding to bars and spiral arms. We demonstrate the versatility of this technique through applying our models to galaxy images from two different observational datasets with different source imagery, and to RGB colour and monochromatic galaxy imaging. We apply our models to analyse SDSS and Subaru HSC imaging of barred galaxies from the NA10 and SAMI catalogues in order to determine the dependence of bar length on stellar mass, morphology, redshift and the spin parameter proxy $\lambda_{R_e}$. Based on the predicted bar masks, we show that the relative bar scale length varies with morphology, with early type galaxies hosting longer bars. While bars are longer in more massive galaxies in absolute terms, relative to the galaxy disc they are actually shorter. We also find that the normalised bar length decreases with increasing redshift, with bars in early-type galaxies exhibiting the strongest rate of decline. We show that it is possible to distinguish spiral arms and bars in monochrome imaging, although for a given galaxy the estimated length in monochrome tends to be longer than in colour imaging. Our morphological segmentation technique can be efficiently applied to study bars in large-scale surveys and even in cosmological simulations.

Coronal fan loops rooted in sunspot umbra show outward propagating waves with subsonic phase speed and period around 3-min. However, their source region in the lower atmosphere is still ambiguous. We performed multi-wavelength observations of a clean fan loop system rooted in sunspot observed by Interface Region Imaging Spectrograph (IRIS) and Solar Dynamics Observatory (SDO). We utilised less explored property of frequency modulation of these 3-min waves from the photosphere to corona, and found them to be periodic with the ranges in 14-20 min, and 24-35 min. Based on our findings, we interpret that 3-min slow waves observed in the coronal fan loops are driven by 3-min oscillations observed at the photospheric footpoints of these fan loops in the umbral region. We also explored any connection between 3-min and 5-min oscillations observed at the photosphere, and found them to be poorly understood. Results provide clear evidence of magnetic coupling of the solar umbral atmosphere through propagation of 3-min waves along the fan loops at different atmospheric heights.

A significant part of barred disc galaxies exhibits boxy/peanut-shaped structures (B/PS bulges) at high inclinations. Another structure also associated with the bar is a barlens, often observed in galaxies in a position close to face-on. At this viewing angle, special kinematic tests are required to detect a 3D extension of the bars in the vertical direction (B/PS bulges). We use four pure $N$-body models of galaxies with B/PS bulges, which have different bar morphology from bars with barlenses to the so-called face-on peanut bars. We analyse the kinematics of our models to establish how the structural features of B/PS bulges manifest themselves in the kinematics for galaxies at intermediate inclinations and whether these features are related to the barlenses. We apply the dissection of the bar into different orbital groups to determine which of them are responsible for the features of the LOSVD (line-of-sight velocity distribution), i.e., for the deep minima of the $h_4$ parameter along the major axis of the bar. As a result, we claim that for our models at the face-on position, the kinematic signatures of a `peanut' indeed track the vertical density distribution features. We conclude that orbits responsible for such kinematic signatures differ from model to model. We pay special attention to the barlens model. We show that orbits assembled into barlens are not responsible for the kinematic signatures of B/PS bulges. The results presented in this work are applicable to the interpretation of IFU observations of real galaxies.

The generation and evolution of dust in galaxies are important tracers for star formation, and can characterize the rest-frame ultraviolet to infrared emission from the galaxies. In particular understanding dust in high-redshift galaxies are important for observational cosmology, as they would be necessary to extract information on star formation in the early universe. We update the public semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) to model the evolution of dust, focusing on high-redshift star-forming galaxies with stellar masses of $\sim 10^8$--$10^{10}M_\odot$ observed by ALMA ($z\approx 7$) and JWST ($z\approx 11$). We find that these galaxies should qualitatively differ in their star formation properties; while the samples in ALMA are explained by dust growth in normal star-forming galaxies, the lack of dust in the samples by JWST requires dust ejection by radiation pressure due to recent highly efficient star-formation within a few 10 Myr, with order 100 times higher efficiency than normal galaxies calibrated by A-SLOTH. Depending on where the JWST galaxies locate on the luminosity function, their bursty star formation histories inferred from our model can have impacts for rates of star formation, supernova explosion, stellar feedback, and detectability of dusty, mature galaxies in the very early universe.

Signatures of dark matter in celestial objects have become of increasing interest due to their powerful detection prospects. To test any of these signatures, the fundamental quantity needed is the rate in which dark matter is captured by celestial objects. Depending on whether dark matter is light, heavy, or comparable in mass to the celestial-body scattering targets, there are different considerations when calculating the capture rate. Furthermore, if dark matter has strong or weak interactions, the physical behaviour important for capture varies. Using both analytic approximations and simulations, we demonstrate how to treat dark matter capture in a range of celestial objects for arbitrary dark matter mass and interaction strength. We release our calculation framework as a public package available in both Python and Mathematica versions, called Asteria.

A large primordial lepton asymmetry is capable of explaining the baryon asymmetry of the Universe (BAU) through suppression of the electroweak sphaleron rates (``sphaleron freeze-in") which can lead to a first-order cosmic QCD transition with an observable gravitational wave (GW) signal. With next-to-leading order dimensional reduction and the exact 1-loop fluctuation determinant, we accurately compute the lepton asymmetry needed to realize this paradigm, finding it to be an order of magnitude smaller than previous estimates. Further, we apply an improved QCD equation of state capable of describing the phase transition line together with the critical endpoint leading to better agreement with lattice and functional QCD results. Based on this, we identify the range of lepton flavor asymmetries inducing a first-order cosmic QCD transition. We then extract the parameters relevant to the prediction of GW signal from a first-order cosmic QCD transition. Our result showcases the possibility of probing the sphaleron freeze-in paradigm as an explanation of BAU by future gravitational wave experiments like $\mu$Ares.

Since the derivation of a well-defined $D\rightarrow 4$ limit for 4D Einstein Gauss-Bonnet (4DEGB) gravity coupled to a scalar field, there has been interest in testing it as an alternative to Einstein's general theory of relativity. Using the Tolman-Oppenheimer-Volkoff (TOV) equations modified for 4DEGB gravity, we model the stellar structure of quark stars using a novel interacting quark matter equation of state. We find that increasing the Gauss-Bonnet coupling constant $\alpha$ or the interaction parameter $\lambda$ both tend to increase the mass-radius profiles of quark stars described by this theory, allowing a given central pressure to support larger quark stars in general. These results logically extend to cases where $\lambda < 0$, in which increasing the magnitude of the interaction effects instead diminishes masses and radii. We also analytically identify a critical central pressure in both regimes, below which no quark star solutions exist due to the pressure function having no roots. Most interestingly, we find that quark stars can exist below the general relativistic Buchdahl bound and Schwarzschild radius $R=2M$, due to the lack of a mass gap between black holes and compact stars in 4DEGB. Even for small $\alpha$ well within current observational constraints, we find that quark star solutions in this theory can describe Extreme Compact Objects (ECOs), objects whose radii are smaller than what is allowed by general relativity.

We derive quantum kinetic equations for mixing neutrinos including consistent forward scattering terms and collision integrals for coherent neutrino states. In practice, we reduce the general Kadanoff--Baym equations in a few clearly justified steps to a generalized density matrix equation that describes both the flavour- and particle-antiparticle coherences and is valid for arbitrary neutrino masses and kinematics. We then reduce this equation to a simpler particle-antiparticle diagonal limit and eventually to the ultra-relativistic limit. Our derivation includes simple Feynman rules for computing collision integrals with the coherence information. We also expose a novel spectral shell structure underlying the mixing phenomenon and quantify how the prior information on the system impacts on the QKE's, leading to a direct effect on its evolution. Our results can be used for example to accurately model neutrino distributions in hot and dense environments and to study the production and decay of mixing heavy neutrinos in colliders.

The history of the seemingly simple problem of straight line fitting in the presence of both $x$ and $y$ errors has been fraught with misadventure, with statistically ad hoc and poorly tested methods abounding in the literature. The problem stems from the emergence of latent variables describing the "true" values of the independent variables, the priors on which have a significant impact on the regression result. By analytic calculation of maximum a posteriori values and biases, and comprehensive numerical mock tests, we assess the quality of possible priors. In the presence of intrinsic scatter, the only prior that we find to give reliably unbiased results in general is a mixture of one or more Gaussians with means and variances determined as part of the inference. We find that a single Gaussian is typically sufficient and dub this model Marginalised Normal Regression (MNR). We illustrate the necessity for MNR by comparing it to alternative methods on an important linear relation in cosmology, and extend it to nonlinear regression and an arbitrary covariance matrix linking $x$ and $y$. We publicly release a Python/Jax implementation of MNR and its Gaussian mixture model extension that is coupled to Hamiltonian Monte Carlo for efficient sampling, which we call ROXY (Regression and Optimisation with X and Y errors).

Pre-merger localization of binary neutron stars (BNSs) is one of the most important scientific goals for the third generation (3G) gravitational wave (GW) detectors. It will enable the electromagnetic observation of the whole process of BNS coalescence, especially for the pre-merger and merger phases which have not been observed yet, opening a window for deeper understandings of compact objects. To reach this goal, we describe a novel combination of multi-band matched filtering and semi-analytical localization algorithms to achieve early-warning localization of long BNS signals in 3G detectors. Using our method we are able to efficiently simulate one month of observations with a three-detector 3G network, and show that it is possible to provide accurate sky localizations more than 30 minutes before the merger. Our simulation shows that there could be ~ 10 (~ 100) BNS events localized within 100 deg2, 20 (6) minutes before merger, per month of observation.

We consider a formulation of the Brans-Dicke theory in Jordan's frame for Bianchi-I spacetime within the framework of loop quantum gravity. The robustness of singularity resolutions due to the quantum effects is explicitly verified. We also present an exploration of the effects of quantum geometry on the background dynamics. In particular, starting with a detailed constraint analysis, we present the off-shell Lie algebra of the full theory in the geometrodynamic phase space. This is followed by casting the theory in terms of the Yang-Mills phase space variables suitable for invoking loop quantization and show that the system is highly nonlinear. Finally, we present the quantum-corrected effective dynamics that resolve the initial singularity

Being born hot from core-collapse supernova, the crust of the proto-neutron star is expected to be made of a Coulomb liquid and composed of an ensemble of different nuclear species. In this work, we study the beta-equilibrated proto-neutron-star crust in the liquid phase in a self-consistent multi-component approach, employing a compressible liquid-drop description of the ions including the ion centre-of-mass motion. Particular care is also devoted to the calculation of the rearrangement term, thus ensuring thermodynamic consistency. We compare the results of the multi-component plasma calculations with those obtained within a one-component (single-nucleus) approach, showing that important differences arise between the predictions of the two treatments. In particular, the abundances of helium clusters become important using a complete multi-component plasma approach, and eventually dominate the whole distribution at higher temperature in the crust.

Complexity will be more and more essential in high-energy physics. It is naturally extended into the very early universe. Considering the universe as a quantum chaotic system, the curvature perturbation of the scalar field is identified with the two-mode squeezed state. By solving the Schr$\ddot{o}$dinger equation, one can obtain the numerical solutions of the angle parameter and squeezing parameter. The solution of the squeezing parameter mainly determines the evolution of complexity. Our numeric indicates that the complexity of the modified dispersion relation will have a non-linear pattern after the horizon exits. Meanwhile, its corresponding Lyapunov index is also larger compared with the standard case. During the inflationary period, the complexity will irregularly oscillate and its scrambling time is also shorter compared with the standard case. Since the modified dispersion relation can be dubbed as the consequences of various frameworks of quantum gravity, it could be applicable to these frameworks. Finally, one can expect the framework of quantum gravity will lead to the fruitful evolution of complexity, which guides us in distinguishing various inflationary models.

We use 120 three dimensional direct numerical simulations to study fingering convection in non-rotating spherical shells. We investigate the scaling behaviour of the flow lengthscale, mean velocity and transport of chemical composition over the fingering convection instability domain defined by $1 \leq R_\rho \leq Le$, $R_\rho$ being the ratio of density perturbations of thermal and compositional origins. We show that the horizontal size of the fingers is accurately described by a scaling law of the form $\mathcal{L}_h/d \sim |Ra_T|^{-1/4} (1-\gamma)^{-1/4}/\gamma^{-1/4}$, where $d$ is the shell depth, $Ra_T$ the thermal Rayleigh number and $\gamma$ the flux ratio. Scaling laws for mean velocity and chemical transport are derived in two asymptotic regimes close to the two edges of the instability domain, namely $R_\rho \lesssim Le$ and $R_\rho \gtrsim 1$. For the former, we show that the transport follows power laws of a small parameter $\epsilon^\star$ measuring the distance to onset. For the latter, we find that the Sherwood number $Sh$, which quantities the chemical transport, gradually approaches a scaling $Sh\sim Ra_\xi^{1/3}$ when $Ra_\xi \gg 1$; and that the P\'eclet number accordingly follows $Pe \sim Ra_\xi^{2/3} |Ra_T|^{-1/4}$, $Ra_\xi$ being the chemical Rayleigh number. When the Reynolds number exceeds a few tens, a secondary instability may occur taking the form of large-scale toroidal jets. Jets distort the fingers resulting in Reynolds stress correlations, which in turn feed the jet growth until saturation. This nonlinear phenomenon can yield relaxation oscillation cycles.

Guided by previous non-perturbative lattice simulations of a two-step electroweak phase transition, we reformulate the perturbative analysis of equilibrium thermodynamics for generic cosmological phase transitions in terms of effective field theory (EFT) expansions. Based on thermal scale hierarchies, we argue that the scale of many interesting phase transitions is in-between the soft and ultrasoft energy scales, which have been the focus of studies utilising high-temperature dimensional reduction. The corresponding EFT expansions provide a handle to control the perturbative expansion, and allow us to avoid spurious infrared divergences, imaginary parts, gauge dependence and renormalisation scale dependence that have plagued previous studies. As a direct application, we present a novel approach to two-step electroweak phase transitions, by constructing separate effective descriptions for two consecutive transitions. Our approach provides simple expressions for effective potentials separately in different phases, a numerically inexpensive method to determine thermodynamics, and significantly improves agreement with the non-perturbative lattice simulations.

Observation of gravitational waves from inspiralling binary black holes has offered a unique opportunity to study the physical parameters of the component black holes. To infer these parameters, Bayesian methods are employed in conjunction with general relativistic waveform models that describe the source's inspiral, merger, and ringdown. The results depend not only on the accuracy of the waveform models but also on the underlying fiducial prior distribution used for the analysis. In particular, when the pre-merger phase of the signal is barely observable within the detectors' bandwidth, as is currently the case with intermediate-mass black hole binary signals in ground-based gravitational wave detectors, different prior assumptions can lead to different interpretations. In this study, we utilise the gravitational-wave inference library, $\texttt{Parallel Bilby}$, to evaluate the impact of mass prior choices on the parameter estimation of intermediate-mass black hole binary signals. While previous studies focused primarily on analysing event data, we offer a broader, more controlled study by using simulations. Our findings suggest that the posteriors in total mass, mass ratio and luminosity distance are contingent on the assumed mass prior distribution used during the inference process. This is especially true when the signal lacks sufficient pre-merger information and/or has inadequate power in the higher-order radiation multipoles. In conclusion, our study underscores the importance of thoroughly investigating similarly heavy events in current detector sensitivity using a diverse choice of priors. Absent such an approach, adopting a flat prior on the binary's redshifted total mass and mass ratio emerges as a reasonable choice, preventing biases in the detector-frame mass posteriors.

The conventional theory of magnetic field generation in a turbulent flow considers time-reversible random flows. However, real turbulent flows are known to be time irreversible: the presence of energy cascade is an intrinsic property of turbulence. We generalize the 'standard' model to account for the irreversibility. We show that even small time asymmetry leads to significant suppression of the dynamo effect at low magnetic Prandtl numbers, increases the generation threshold and may even make generation impossible for any magnetic Reynolds number. We calculate the magnetic energy increment as a function of the parameters of the flow.

We present a novel mechanism for the primordial black hole (PBH) production within the QCD axion framework. We take the case where the Peccei-Quinn symmetry breaks during inflation, resulting in a $N_{\rm DW}=1$ string-wall network that re-enters horizon sufficiently late. Therefore, closed axion domain walls naturally arising in the network are sufficiently large to collapse into PBHs. Our numerical simulation shows that $0.8\%$ of the total wall area is in the form of closed walls. Notably, this fraction is independent of any axion parameters, as its determination is firmly grounded in the principles of percolation theory. In addition, the relic abundance of dark matter is accounted for by free axions from the collapse of open walls bounded by strings. This framework yields a calculated PBH fraction of dark matter as 0.0256. The PBHs uniformly share the same mass, which spans from about $10^{-10}$ to $10^{-1}$ solar masses, corresponding to the classical QCD axion mass window $10^{-5}-10^{-2}$ eV and the re-entering horizon temperature $300-1$ MeV. Intriguingly, PBHs in this mechanism can naturally account for the gravitational-lensing events observed by the OGLE collaboration.

Invited contribution to the Encyclopedia of Mathematical Physics (2nd edition), providing an overview over some main ideas and results in quantum cosmology. Key points: Canonical quantisation of homogeneous, isotropic cosmology; discussion of ambiguities in this quantisation; Construction of explicit solutions, attempts at physical interpretation; conceptual issues (time evolution, unitarity, probability interpretation); Introduction of scalar fields or perfect fluids; Discussion of whether classical singularities are resolved; Addition of inhomogeneities, possible predictions for primordial cosmology; Connections to string theory and loop quantum gravity.

In Starobinsky inflation with a Weyl squared Lagrangian $-\alpha C^2$, where $\alpha$ is a coupling constant, we study the linear stability of cosmological perturbations on a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker background. In this theory, there are two dynamical vector modes propagating as ghosts for $\alpha>0$, whose condition is required to avoid tachyonic instabilities of vector perturbations during inflation. The tensor sector has four propagating degrees of freedom, among which two of them correspond to ghost modes. However, tensor perturbations approach constants after the Hubble radius crossing during inflation, and hence the classical instabilities are absent. In the scalar sector, the Weyl curvature gives rise to a ghost mode coupled to the scalaron arising from the squared Ricci scalar. We show that two gauge-invariant gravitational potentials, which are both dynamical in our theory, are subject to exponential growth after the Hubble radius crossing. There are particular gauge-invariant combinations like the curvature perturbations whose growth is suppressed, but it is not possible to remove the instability of other propagating degrees of freedom present in the perturbed metric. This violent and purely classical instability present in the scalar sector makes the background unviable. Furthermore, the presence of such classical instability makes the quantization of the modes irrelevant, and the homogeneous inflationary background is spoiled by the Weyl curvature term.

We investigate the weak-interaction-driven bulk-viscous transport properties of $npe$ matter in the neutrino transparent regime. Previous works assumed that the induced bulk viscosity correction to pressure, near beta equilibrium, is linear in deviations from the equilibrium charge fraction. We show that this is not always true for (some) realistic equations of state at densities between one and three times saturation density. This nonlinear nature of the perturbation around equilibrium motivates a far-from-beta-equilibrium description of bulk-viscous transport in neutron star mergers, which can be precisely achieved using a new Israel-Stewart formulation with resummed bulk and relaxation time transport coefficients. The computation of these transport coefficients depends on out-of-beta-equilibrium pressure corrections, which can be computed for a given equation of state. We calculate these coefficients for equations of state that satisfy the latest constraints from multi-messenger observations from LIGO/VIRGO and NICER. We show that varying the nuclear symmetry energy $J$ and its slope $L$ can significantly affect the transport coefficients and the nonlinear behavior of the out-of-equilibrium pressure corrections. Therefore, having better constraints on $J$ and $L$ will directly impact our understanding of bulk-viscous processes in neutron star mergers.

The lightest $Z_2$ odd particle in the scotogenic model, referred to as scotogenic dark matter (DM), is a widely studied candidate for DM. This scotogenic DM is generated through well-known thermal processes as well as via the evaporation of primordial black holes (PBHs). Recent reports suggested that the curvature fluctuations of PBHs during an epoch dominated by these entities in the early universe can serve as the source of so-called induced gravitational waves (GWs). In this study, we demonstrate that stringent constraints on the mass of scotogenic DM and PBHs can be obtained through the detection of induced GWs using future detectors.

Secondary messengers such as neutrinos and photons are expected to be produced in interactions of ultra-high-energy cosmic rays (UHECRs) with extragalactic background photons. Their propagation could be altered by the effects of Lorentz invariance violation. In this work, we have developed an extension of the SimProp code that includes some Lorentz-violating scenarios affecting the propagation of neutrinos. We present the corresponding expected cosmogenic neutrino fluxes for three different astrophysical scenarios for the production of UHECRs. These results can be used to put constraints on the scale of Lorentz violation in the neutrino sector.

The Hubble constant is one of the most important parameters in cosmology. Discrepancies in values of the Hubble constant estimated from various measurements, the so-called Hubble tension, are a serious problem. In this paper, we study the effects of small-scale inhomogeneities of structure formation on the measurement of the Hubble constant using the luminosity distance-redshift relation. By adopting the adhesion model in Newtonian cosmology as the model of structure formation, we investigate whether or not the effects of inhomogeneities can be sufficiently large to affect the current observations of the Hubble constant. We show that inappropriate treatment of the effects of inhomogeneities can cause a large deviation of the measured value of the Hubble constant from the background value, whose magnitude is comparable with the Hubble tension. Our main message is the importance of adopting an appropriate model of structure formation to investigate the effects of inhomogeneities. We also add discussion on the spatial averaging approach used to estimate the measured Hubble constant in the inhomogeneous universe.

We present a phenomenological model to investigate the chiral phase transition characterized by parity doubling in dense, beta equilibrated, cold matter. Our model incorporates effective interactions constrained by SU(3) relations and considers baryonic degrees of freedom. By constraining the model with astrophysical data and nuclear matter properties, we find a first-order phase transition within realistic values of the slope parameter L. The inclusion of the baryon octet and negative parity partners, along with a chiral-invariant mass $m_{0}$, allows for a non-massless chiral symmetric phase. Through exploration of parameter space, we identify parameter sets satisfying mass and radius constraints without requiring a partonic phase. The appearance of the parity partner of the nucleon, the N(1535) resonance, suppresses strangeness, pushing hyperonization to higher densities. We observe a mild first-order phase transition to the chirally restored phase, governed by $m_{0}$. Our calculations of surface tension highlight its strong dependence on $m_{0}$. The existence of mixed phases is ruled out since they become energetically too costly. We compare stars with metastable and stable cores using both branches of the equation of state. Despite limited lifespans due to low surface tension values, phase conversion and star contraction could impact neutron stars with masses around 1.3 solar masses or more. We discuss some applications of this model in its non-zero temperatures generalization and scenarios beyond beta equilibrium that can provide insights into core-collapse supernovae, proto-neutron star evolution, and neutron star mergers. Core-collapse supernovae dynamics, influenced by chiral symmetry restoration and exotic hadronic states, affect explosion mechanisms and nucleosynthesis.