Papers for Wednesday, Jun 01 2022

From data collection to generating a white dwarf luminosity function, it requires numerous Astrophysical, Mathematical and Computational domain knowledge. The steep learning curve makes it difficult to enter the field and often individuals have to reinvent the wheel to perform identical data reduction and analysis tasks. We have gathered all the publicly available white dwarf cooling models and synthetic photometry to provide a toolkit that allows visualisation of various models, photometric fitting of a white dwarf with or without distance and reddening, and the computing of white dwarf luminosity functions with a choice of initial mass function, main sequence evolution model, star formation history, initial-final mass relation, and white dwarf cooling model. We have recomputed and compared it against the Gaia EDR3 white dwarf catalogue from Gentile Fusillo et al. (2021), we show excellent agreements between the two works.

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

We study star cluster formation at low metallicities of $Z/Z_\odot=10^{-4}$--$10^{-1}$ using three-dimensional hydrodynamics simulations. Particular emphasis is put on how the stellar mass distribution is affected by the cosmic microwave background radiation (CMB), which sets the temperature floor to the gas. Starting from the collapse of a turbulent cloud, we follow the formation of a protostellar system resolving $\sim$au scale. In relatively metal-enriched cases of $Z/Z_\odot \gtrsim 10^{-2}$, where the mass function resembles the present-day one in the absence of the CMB, high temperature CMB suppresses cloud fragmentation and reduces the number of low-mass stars, making the mass function more top-heavy than in the cases without CMB heating at $z\gtrsim10$. In lower-metallicity cases with $Z/Z_\odot \lesssim 10^{-3}$, where the gas temperature is higher than the CMB value due to inefficient cooling, the CMB has only a minor impact on the mass distribution, which is top-heavy regardless of the redshift. In cases either with a low metallicity of $Z/Z_\odot \lesssim 10^{-2}$ or at a high redshift $z\gtrsim10$, the mass spectrum consists of a low-mass Salpeter-like component, peaking at $0.1~M_\odot$, and a top-heavy component with $10$--$50~M_\odot$, with the fraction in the latter increasing with increasing redshift. In galaxies forming at $z\gtrsim10$, the major targets of the future instruments including JWST, CMB heating makes the stellar mass function significantly top-heavy, enhancing the number of supernova explosions by a factor of $1.4$ ($2.8$) at $z=10$ ($20$, respectively) compared to the prediction by Chabrier initial mass function when $Z/Z_\odot=0.1$.

This work is based on high quality integral-field spectroscopic data obtained with the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT). The 21 brightest ($m_B\leq 15$ mag) early-type galaxies (ETGs) inside the virial radius of the Fornax cluster are observed out to distances of $\sim2-3\ R_{\rm e}$. Deep imaging from the VLT Survey Telescope (VST) is also available for the sample ETGs. We investigate the variation of the galaxy structural properties as a function of the total stellar mass and cluster environment. Moreover, we correlate the size scales of the luminous components derived from a multi-component decomposition of the VST surface-brightness radial profiles of the sample ETGs with the MUSE radial profiles of stellar kinematic and population properties. The results are compared with both theoretical predictions and previous observational studies and used to address the assembly history of the massive ETGs of the Fornax cluster. We find that galaxies in the core and north-south clump of the cluster, which have the highest accreted mass fraction, show milder metallicity gradients in their outskirts than the galaxies infalling into the cluster. We also find a segregation in both age and metallicity between the galaxies belonging to the core and north-south clump and the infalling galaxies. The new findings fit well within the general framework for the assembly history of the Fornax cluster.

Measurements of the atomic hydrogen (HI) properties of high-redshift galaxies are critical to understanding the decline in the star-formation rate (SFR) density of the Universe after its peak $\approx8-11$ Gyr ago. Here, we use $\approx510$ hours of observations with the upgraded Giant Metrewave Radio Telescope to measure the dependence of the average HI mass of star-forming galaxies at $z=0.74-1.45$ on their average stellar mass and redshift, by stacking their HI 21 cm emission signals. We divide our sample of 11,419 main-sequence galaxies at $z=0.74-1.45$ into two stellar-mass ($M_*$) subsamples, with $M_*>10^{10} M_\odot$ and $M_*<10^{10} M_\odot$, and obtain clear detections, at $>4.6\sigma$ significance, of the stacked HI 21 cm emission in both subsamples. We find that galaxies with $M_*>10^{10} M_\odot$, which dominate the decline in the cosmic SFR density at $z\lesssim1$, have HI reservoirs that can sustain their SFRs for only a short period, $0.86\pm0.20$ Gyr, unless their HI is replenished via accretion. We also stack the HI 21 cm emission from galaxies in two redshift subsamples, at $z=0.74-1.25$ and $z=1.25-1.45$, again obtaining clear detections of the stacked HI 21 cm emission signals, at $>5.2\sigma$ significance in both subsamples. We find that the average HI mass of galaxies with $\langle M_* \rangle\approx10^{10} M_\odot$ declines steeply over a period of $\approx1$ billion years, from $(33.6\pm6.4) \times 10^9 M_\odot$ at $\langle z\rangle\approx1.3$ to $(10.6\pm1.9)\times10^9 M_\odot$ at $\langle z\rangle\approx1.0$, i.e. by a factor $\gtrsim3$. We thus find direct evidence that accretion of HI onto star-forming galaxies at $z\approx1$ is insufficient to replenish their HI reservoirs and sustain their SFRs, thus resulting in the decline in the cosmic SFR density 8 billion years ago.

The most direct way to confront observed galaxies with those formed in numerical simulations is to forward-model simulated galaxies into synthetic observations. Provided that synthetic galaxy observations include similar constraints and limitations as real observations, they can be used to (1) carry out even-handed comparisons of observation and theory and (2) map the observable characteristics of simulated galaxies to their a priori known origins. In particular, integral field spectroscopy (IFS) expands the scope of such comparisons and mappings to an exceptionally broad set of physical properties. We therefore present RealSim-IFS: a tool for forward-modelling galaxies from hydrodynamical simulations into synthetic IFS observations. The core components of RealSim-IFS model the detailed spatial sampling mechanics of any fibre-bundle, image slicer, or lenslet array IFU and corresponding observing strategy, real or imagined, and support the corresponding propagation of noise adopted by the user. The code is highly generalized and can produce cubes in any light- or mass-weighted quantity (e.g. specific intensity, gas/stellar line-of-sight velocity, stellar age/metallicity, etc.). We show that RealSim-IFS exactly reproduces the spatial reconstruction of specific intensity and variance cubes produced by the MaNGA survey Data Reduction Pipeline using the calibrated fibre spectra as input. We then apply RealSim-IFS by producing a public synthetic MaNGA stellar kinematic survey of 893 galaxies with $\log M_{\star}/M_{\odot}>10$ from the TNG50 cosmological hydrodynamical simulation.

Radiative turbulent mixing layers (TMLs) are ubiquitous in astrophysical environments, e.g., the circumgalactic medium (CGM), and are triggered by the shear velocity at interfaces between different gas phases. To understand the shear velocity dependence of TMLs, we perform a set of 3D hydrodynamic simulations with an emphasis on the TML properties at high Mach numbers $\mathcal{M}$. Since the shear velocity in mixing regions is limited by the local sound speed of mixed gas, high-Mach number TMLs develop into a two-zone structure: a Mach number-independent mixing zone traced by significant cooling and mixing, plus a turbulent zone with large velocity dispersions which expands with greater $\mathcal{M}$. Low-Mach number TMLs do not have distinguishable mixing and turbulent zones. The radiative cooling of TMLs at low and high Mach numbers is predominantly balanced by enthalpy consumption and turbulent dissipation respectively. Both the TML surface brightness and column densities of intermediate-temperature ions (e.g., O VI) scale as $\propto\mathcal{M}^{0.5}$ at $\mathcal{M} \lesssim 1$, but reach saturation ($\propto \mathcal{M}^0$) at $\mathcal{M} \gtrsim 1$. Inflow velocities and hot gas entrainment into TMLs are substantially suppressed at high Mach numbers, and strong turbulent dissipation drives the evaporation of cold gas. This is in contrast to low-Mach number TMLs where the inflow velocities and hot gas entrainment are enhanced with greater $\mathcal{M}$, and cold gas mass increases due to the condensation of entrained hot gas.

Current gravitational-wave observations set the most stringent bounds on the abundance of primordial black holes (PBHs) in the solar mass range. This constraint, however, inherently relies on the merger rate predicted by PBH models. Previous analyses have focused mainly on two binary formation mechanisms: early Universe assembly out of decoupling from the Hubble expansion and dynamical capture in present-day dark matter structures. Using reaction rates of three-body processes studied in the astrophysical context, we show that, under conservative assumptions, three-body interactions in PBH halos efficiently produce binaries and significantly contribute to the overall merger rate, provided PBHs make up a sufficient fraction of the dark matter. Those binaries form at high redshift in Poisson-induced PBH small-scale structures and a fraction is predicted to coalesce and merge within the current age of the Universe, at odds with the dynamical capture scenario where they merge promptly. Our results enable LIGO/Virgo/KAGRA constraints on the PBH abundance to set stronger bounds in this interesting mass range, and have important implications for clustered PBH scenarios that might evade such constraints.

In the last decades, luminous accreting super-massive black holes have been discovered within the first Gyr after the Big Bang, but their origin is still an unsolved mystery. We discuss our state-of-the-art theoretical knowledge of their formation physics and early growth, and describe the results of dedicated observational campaigns in the X-ray band. We also provide an overview of how these systems can be used to derive cosmological parameters. Finally, we point out some open issues, in light of future electro-magnetic and gravitational-wave astronomical facilities.

We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single-binary interactions, that drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step toward the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions.

Recent CHIME/FRB observations of the periodic repeating FRB 180916B have produced a homogeneous sample of 44 bursts. These permit a redetermination of the modulation period and phase window, in agreement with earlier results. If the periodicity results from the precession of an accretion disc, in analogy with those of Her X-1, SS 433, and many other superorbital periods, the width of the observable phase window indicates that the disc axis jitters by an angle of about 0.14 of the inclination angle, similar to the ratio of 0.14 in the well-observed jittering jet source SS 433.

We determine the bright end of the rest-frame UV luminosity function at $z=8-10$ by selecting bright $z\gtrsim 8$ photometric candidates from the largest systematic compilation of HST (pure-)parallel observations to date, the Super-Brightest-of-Reionizing-Galaxies (SuperBoRG) data set. The data set includes $\sim300$ independent sightlines distributed across the sky with WFC3 observations, totalling $800-1300$ arcmin$^2$ (depending on redshift). We identify 31 $z\gtrsim8$ candidates via colour selection and photo-$z$ analysis with observed magnitude ($24.1< H_{160} <26.6$). Following detailed completeness and source recovery simulations in each independent field, as well as modelling of interloper contamination, we derive rest-frame UV luminosity functions at $z=8-10$ down to UV magnitude $\simeq-23$. We find that the bright-end of the galaxy luminosity function can be described both by a Schechter and by a double power-law function, with our space-based large area determination showing some tentative discrepancies with the luminosity functions derived from ground-based observations at the same redshifts. Interestingly, the UV luminosity function we derived at $z=8$ is consistent with no evolution of the bright end from $z=6-7$, suggesting a substantial contribution from AGNs at magnitude $M_{UV}<-22$. On the other hand, our derived UVLFs at all redshifts are still consistent with a scenario where the UV light originates from stars, but there is a vanishing dust content at $z\gtrsim 7.5$. Both scenarios raise interesting prospects to further understand galaxy formation in extreme objects during the epoch of reionization thanks to upcoming James Webb Space Telescope spectroscopic observations.

We generalise existing constraints on primordial black holes to dark objects with extended sizes using the aLIGO design sensitivity. We show that LIGO is sensitive to dark objects with radius $O(10-10^3~{\rm km})$ if they make up more than $\sim O(10^{-2}-10^{-3})$ of dark matter.

Meteor showers occur when streams of meteoroids originating from a common source intersect the Earth. There will be small dissimilarities between the direction of motion of different meteoroids within a stream, and these small differences will act to broaden the radiant, or apparent point of origin, of the shower. This dispersion in meteor radiant can be particularly important when considering the effect of the Earth's gravity on the stream, as it limits the degree of enhancement of the stream's flux due to gravitational focusing. In this paper, we present measurements of the radiant dispersion of twelve showers using observations from the Global Meteor Network. We find that the median offset of individual meteors from the shower radiant ranges from 0.32$^\circ$ for the eta Aquariids to 1.41$^\circ$ for the Southern Taurids. We also find that there is a small but statistically significant drift in Sun-centered ecliptic radiant and/or geocentric speed over time for most showers. Finally, we compare radiant dispersion with shower duration and find that, in contrast with previous results, the two quantities are not correlated in our data.

Jet precessions are widely involved in astrophysical phenomena from galaxies to X-ray binaries and gamma-ray bursts (GRBs). Polarization presents a unique probe of the magnetic fields in GRB jets. The precession of GRBs relativistic jets will change the geometry within the observable emitting region of the jet, which can potentially affect the polarization of the afterglow. In this paper, we take into account jet precession to study the polarization evolution and corresponding light curves in GRB early optical afterglows with ordered and random magnetic field geometries. We find that the jet precession in long-lived engines can significantly reduce the polarization degree (PD) regardless of the magnetic field structure. The strongest PD attenuation is found when the line of sight is aligned with the precession axis. Our results show that jet precession can provide new insight into the low PD measured in the early optical afterglows of GRBs.

Chemical processes on the surface of icy grains play an important role in the chemical evolution in molecular clouds. In particular, reactions involving non-energetic hydrogen atoms accreted from the gaseous phase have been extensively studied. These reactions are believed to effectively proceed only on the surface of the icy grains; thus, molecules embedded in the ice mantle are not considered to react with hydrogen atoms. Recently, Tsuge et al. (2020) suggested that non-energetic hydrogen atoms can react with CO molecules even in ice mantles via diffusive hydrogenation. This investigation was extended to benzene and naphthalene molecules embedded in amorphous solid water (ASW) in the present study, which revealed that a portion of these molecules could be fully hydrogenated in astrophysical environments. The penetration depths of non-energetic hydrogen atoms into porous and non-porous ASW were determined using benzene molecules to be >50 and ~10 monolayers, respectively (1 monolayer ~ 0.3 nm).

This work proposes a Residual Recurrent Neural Network (RRNet) for synthetically extracting spectral information, and estimating stellar atmospheric parameters together with 15 chemical element abundances for medium-resolution spectra from Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). The RRNet consists of two fundamental modules: a residual module and a recurrent module. The residual module extracts spectral features based on the longitudinally driving power from parameters, while the recurrent module recovers spectral information and restrains the negative influences from noises based on Cross-band Belief Enhancement. RRNet is trained by the spectra from common stars between LAMOST DR7 and APOGEE-Payne catalog. The 17 stellar parameters and their uncertainties for 2.37 million medium-resolution spectra from LAMOST DR7 are predicted. For spectra with S/N >= 10, the precision of estimations Teff and log g are 88 K and 0.13 dex respectively, elements C, Mg, Al, Si, Ca, Fe, Ni are 0.05 dex to 0.08 dex, and N, O, S, K, Ti, Cr, Mn are 0.09 dex to 0.14 dex, while that of Cu is 0.19 dex. Compared with StarNet and SPCANet, RRNet shows higher accuracy and robustness. In comparison to Apache Point Observatory Galactic Evolution Experiment and Galactic Archaeology with HERMES surveys, RRNet manifests good consistency within a reasonable range of bias. Finally, this work releases a catalog for 2.37 million medium-resolution spectra from the LAMOST DR7, the source code, the trained model and the experimental data respectively for astronomical science exploration and data processing algorithm research reference.

Using the cross-matched data of Gaia EDR3 and the 2MASS Point Source Catalog, a sample of RC stars with parallax accuracies better than 20% is identified and used to reveal the nearby spiral pattern traced by old stars. As shown in the overdensity distribution of RC stars, there is an arc-like feature extended from $l~\sim$ 90$^\circ$ to $\sim$ 243$^\circ$, which passed close to the Sun. This feature is probably an arm segment traced by old stars, indicating the Galaxy potential in the vicinity of the Sun. By comparing to the spiral arms depicted by young objects, we found that there are considerable offsets between the two different components of Galactic spiral arms. The spiral arm traced by RC stars tends to have a larger pitch angle, hence a more loose wound pattern.

The recent rapid growth of the black hole (BH) catalog from gravitational waves (GWs), has allowed us to study the substructure of black hole mass function (BHMF) beyond the simplest Power-Law distribution. However, the BH masses inferred from binary BH merger events, may be systematically 'brightened' or 'dimmed' by gravitational lensing effect. In this work, we investigate the impact of gravitational lensing on the BHMF inference considering the detection of the third-generation GW detector -- the Einstein Telescope (ET). We focus on high redshift, $z=10$, in order to obtain the upper-limits of this effect. We use Monte Carlo(MC) method to simulate the data adopting 3 original BHMFs under Un-Lensed and Lensed scenarios, then recovery the parameters of BHMFs from the mock data, and compare the difference of results, respectively. We found that all the parameters are well recovered within one standard deviation(std., 1$\sigma$), and all 3 BHMF models are reconstructed within 68\% credible interval, suggesting that lensing would not change the main structure drastically, even at very high redshifts and with high precision of ET. And the modest influence beyond $50M_{\odot}$, depends on the modeling of the high mass tail or substructure of BHMF. We conclude that the impact of lensing on BHMF inference with ET can be safely ignored in the foreseeable future. Careful handling of lensing effects is required only when focus on an accurate estimation of the high mass end of BHMF at high redshifts.

We use 38 C IV quasar (QSO) reverberation-measured observations, which span eight orders of magnitude in luminosity and the redshift range $0.001064 \leq z \leq 3.368$, to simultaneously constrain cosmological-model and QSO radius-luminosity ($R-L$) relation parameters in six cosmological models, using an improved technique that more correctly accounts for the asymmetric errors bars of the time-lag measurements. We find that $R-L$ relation parameters are independent of the cosmological models used in the analysis and so the $R-L$ relation can be used to standardize the C IV QSOs. The C IV QSO cosmological constraints are consistent with those from Mg II QSOs, allowing us to derive joint C IV + Mg II QSO cosmological constraints which are consistent with currently accelerated cosmological expansion, as well as consistent with cosmological constraints derived using better-established baryon acoustic oscillation (BAO) and Hubble parameter [$H(z)$] measurements. When jointly analyzed with $H(z)$ + BAO data, current C IV + Mg II QSO data mildly tighten current $H(z)$ + BAO data cosmological constraints.

We present the system equivalent flux density (SEFD) expressions for all four Stokes parameters: I, Q, U, V. The expressions were derived based on our derivation of SEFD I (for Stokes I) and subsequent extensions of that work to phased array and multipole interferometers. The key to the derivation of the SEFD Q, U, V expressions is to recognize that the noisy estimates of Q, U, V can be written as the trace of a matrix product. This shows that the SEFD I is a special case, where the general case involves a diagonal or anti-diagonal 2x2 matrix interposed in the matrix multiplication. Following this step, the relation between the SEFD for I and Q, U, V becomes immediately evident. We present example calculations for a crossed dipole based on the formulas derived and the comparison between simulation and observation using the Murchison Widefield Array (MWA).

OJ 287 is a BL Lacertae type quasar in which the active galactic nucleus (AGN) outshines the host galaxy by an order of magnitude. The only exception to this may be at minimum light when the AGN activity is so low that the host galaxy may make quite a considerable contribution to the photometric intensity of the source. Such a dip or a fade in the intensity of OJ 287 occurred in November 2017, when its brightness was about 1.75 magnitudes lower than the recent mean level. We compare the observations of this fade with similar fades in OJ 287 observed earlier in 1989, 1999, and 2010. It appears that there is a relatively strong reddening of the B$-$V colours of OJ 287 when its V-band brightness drops below magnitude 17. Similar changes are also seen V$-$R, V$-$I, and R$-$I colours during these deep fades. These data support the conclusion that the total magnitude of the host galaxy is $V=18.0 \pm 0.3$, corresponding to $M_{K}=-26.5 \pm 0.3$ in the K-band. This is in agreement with the results, obtained using the integrated surface brightness method, from recent surface photometry of the host. These results should encourage us to use the colour separation method also in other host galaxies with strongly variable AGN nuclei. In the case of OJ 287, both the host galaxy and its central black hole are among the biggest known, and its position in the black hole mass-galaxy mass diagram lies close to the mean correlation.

Early-time radiative signals from type Ia supernovae (SNe Ia) can provide important constraints on the explosion mechanism and the progenitor system. We present observations and analysis of SN 2019np, a nearby SN Ia discovered within 1-2 days after the explosion. Follow-up observations were conducted in optical, ultraviolet, and near-infrared bands, covering the phases from $\sim-$16.7 days to $\sim$+367.8 days relative to its $B-$band peak luminosity. The photometric and spectral evolutions of SN 2019np resembles the average behavior of normal SNe Ia. The absolute B-band peak magnitude and the post-peak decline rate are $M_{\rm max}(B)=-19.52 \pm 0.47$mag and $\Delta m_{\rm15}(B) =1.04 \pm 0.04$mag, respectively. No Hydrogen line has been detected in the near-infrared and nebular-phase spectra of SN 2019np. Assuming that the $^{56}$Ni powering the light curve is centrally located, we find that the bolometric light curve of SN 2019np shows a flux excess up to 5.0% in the early phase compared to the radiative diffusion model. Such an extra radiation perhaps suggests the presence of an additional energy source beyond the radioactive decay of central nickel. Comparing the observed color evolution with that predicted by different models such as interactions of SN ejecta with circumstellar matter (CSM)/companion star, a double-detonation explosion from a sub-Chandrasekhar mass white dwarf (WD), and surface $^{56}$Ni mixing, the latter one is favored.

The proper motions (PMs) of OB stars in Cygnus have recently been found to exhibit two large-scale kinematic patterns suggestive of expansion. We perform a 3D traceback on these OB stars, the newly-identified OB associations and related open clusters in the region. We find that there are two groups of stars, associations and clusters and that they were each more compact in the past, reaching their closest approach $7.9^{+3.0}_{-1.8}$ and $8.5^{+0.8}_{-2.8}$ Myr ago. We consider two main scenarios for the driver of these large-scale expansion patterns: feedback-driven expansion from a previous generation of massive stars, and expansion as a result of the turbulent velocity field in the primordial molecular cloud. While it is tempting to attribute such large-scale expansion patterns to feedback processes, we find that the observed kinematics are fully consistent with the turbulent origin, and therefore that the injection of further energy or momentum from feedback is not required. Similar conclusions may be drawn for other star forming regions with large-scale expansion patterns.

Insights on stellar surface large-scale magnetic field topologies are usually drawn by applying Zeeman-Doppler-Imaging (ZDI) to the observed spectropolarimetric time series. However, ZDI requires experience for reliable results to be reached and is based on a number of prior assumptions that may not be valid, e.g., when the magnetic topology is evolving on timescales comparable to or shorter than the time span over which observations are collected. In this paper, we present a method based on Principal Component Analysis (PCA) applied to circularly polarised (Stokes~$V$) line profiles of magnetic stars to retrieve the main characteristics of the parent large-scale magnetic topologies, like for instance, the relative strength of the poloidal and toroidal components, and the degree of axisymmetry of the dominant field component and its complexity (dipolar or more complex). We show that this method can also be used to diagnose the temporal variability of the large-scale magnetic field. Performing best for stars with moderate projected equatorial velocities hosting relatively simple magnetic field topologies, this new method is simpler than ZDI, making it convenient to rapidly diagnose the main characteristics of the large-scale fields of non-degenerate stars and to provide insights into the temporal evolution of the field topology.

We examine the solar-cycle variation of the power in the low-degree helioseismic modes by looking at binned power spectra from 45 years of observations with the Birmingham Solar Oscillations Network, which provides a more robust estimate of the mode power than that obtained by peak fitting. The solar-cycle variation of acoustic mode power in the five-minute band is clearly seen. Unusually, even though Cycle 24 was substantially weaker in terms of surface magnetic activity than Cycle 23, the reduction in mode power at solar maximum is very similar for the two cycles, suggesting that the relationship between mode power and magnetic activity is more complex than has previously been thought. This is in contrast to the mode frequencies, which show a strong correlation with activity with only subtle differences between in the response across different solar cycles.

We present the Tenerife Inversion Code (TIC), which has been developed to infer the magnetic and plasma properties of the solar chromosphere and transition region via full-Stokes inversion of polarized spectral lines. The code is based on the HanleRT forward engine, which takes into account many of the physical mechanisms that are critical for a proper modeling of the Stokes profiles of spectral lines originating in the tenuous and highly dynamic plasmas of the chromosphere and transition region: quantum level population imbalance and interference (atomic polarization), frequency coherence effects in polarized resonance scattering (partial frequency redistribution), and the impact of arbitrary magnetic fields on the atomic polarization and the radiation field. We present first results of atmospheric and magnetic inversions, and discuss future developments for the project.

We discuss the statistical distribution of galaxy shapes and viewing angles under the assumption of triaxiality by deprojecting observed Surface Brightness (SB) profiles of 56 Brightest Cluster Galaxies coming from a recently published large deep-photometry sample. For the first time, we address this issue by directly measuring axis ratio profiles without limiting ourselves to a statistical analysis of average ellipticities. We show that these objects are strongly triaxial, with triaxiality parameters 0.39 $ \leq T \leq $ 0.72, have on average axis ratios $< p(r) > = $ 0.84 and $< q(r) > =$ 0.68, and are more spherical in the central regions but flatten out at large radii. Measured shapes in the outskirts agree well with the shapes found for simulated massive galaxies and their dark matter halos from both the IllustrisTNG and the Magneticum simulations, possibly probing the nature of dark matter. In contrast, both simulations fail to reproduce the observed inner regions of BCGs, producing too flattened objects.

The dynamics of ionized gas around the W33 Main ultracompact HII region is studied using observations of hydrogen radio recombination lines and a detailed multiwavelength characterization of the massive star-forming region W33 Main is performed. We used the Giant Meterwave Radio Telescope (GMRT) to observe the H167$\alpha$ recombination line at 1.4 GHz at an angular resolution of 10 arcsec, and Karl. G. Jansky Very Large Array (VLA) data acquired in the GLOSTAR survey to study the dynamics of ionized gas. We also observed the radio continuum at 1.4 GHz and 610 MHz with the GMRT and used GLOSTAR 4$-$8 GHz continuum data to characterize the nature of the radio emission. In addition, archival data from submillimeter to near-infrared wavelengths were used to study the dust emission and identify YSOs in the W33 Main star-forming region. The radio recombination lines were detected at good signal to noise in the GLOSTAR data, while the H167$\alpha$ radio recombination line was marginally detected with the GMRT. The spectral index of radio emission in the region determined from GMRT and GLOSTAR shows the emission to be thermal in the entire region. Along with W33 Main, an arc-shaped diffuse continuum source, G12.81$-$0.22, was detected with the GMRT data. The GLOSTAR recombination line data reveal a velocity gradient across W33 Main and G12.81$-$0.22. The electron temperature is found to be 6343 K and 4843 K in W33 Main and G12.81$-$0.22, respectively. The physical properties of the W33 Main molecular clump were derived by modeling the dust emission using data from the ATLASGAL and Hi-GAL surveys and they are consistent with the region being a relatively evolved site of massive star formation. The gas dynamics and physical properties of G12.81$-$0.22 are consistent with the HII region being in an evolved phase and its expansion on account of the pressure difference is slowing down.

We study the explicitly time-dependent response of a razor-thin axisymmetric disc to externally imposed perturbations by recasting the linearized Collisionless Boltzmann equation as an integral equation and applying Kalnajs' matrix method. As an application we consider the idealized problem of calculating the dynamical friction torque on a steadily rotating, two-dimensional bar. We consider two choices of basis functions in the matrix method, showing that both lead to comparable results. The torques from our linearised calculation are in excellent agreement with those measured from $N$-body simulation, as long as the bar perturbation does not resonate with a significant fraction of the disc's stars.

Bayesian parameter inference is an essential tool in modern cosmology, and typically requires the calculation of $10^5$--$10^6$ theoretical models for each inference of model parameters for a given dataset combination. Computing these models by solving the linearised Einstein-Boltzmann system usually takes tens of CPU core-seconds per model, making the entire process very computationally expensive. In this paper we present \textsc{connect}, a neural network framework emulating \textsc{class} computations as an easy-to-use plug-in for the popular sampler \textsc{MontePython}. \textsc{connect} uses an iteratively trained neural network which emulates the observables usually computed by \textsc{class}. The training data is generated using \textsc{class}, but using a novel algorithm for generating favourable points in parameter space for training data, the required number of \textsc{class}-evaluations can be reduced by two orders of magnitude compared to a traditional inference run. Once \textsc{connect} has been trained for a given model, no additional training is required for different dataset combinations, making \textsc{connect} many orders of magnitude faster than \textsc{class} (and making the inference process entirely dominated by the speed of the likelihood calculation). For the models investigated in this paper we find that cosmological parameter inference run with \textsc{connect} produces posteriors which differ from the posteriors derived using \textsc{class} by typically less than $0.01$--$0.1$ standard deviations for all parameters. We also stress that the training data can be produced in parallel, making efficient use of all available compute resources. The \textsc{connect} code is publicly available for download at \url{https://github.com/AarhusCosmology}.

Over the south pole of Enceladus, an icy moon of Saturn, geysers eject water into space in a striped pattern, making Enceladus one of the most attractive destinations in the search for extraterrestrial life. We explore the ocean dynamics and tracer/heat transport associated with geysers as a function of the assumed salinity of the ocean and various core-shell heat partitions and bottom heating patterns. We find that, even if heating is concentrated into a narrow band on the seafloor directly beneath the south pole, the warm fluid becomes quickly mixed with its surroundings due to baroclinic instability. The warming signal beneath the ice is diffuse and insufficient to prevent the geyser from freezing over. Instead, if heating is assumed to be local to the geyser, emanating from tidal dissipation in the ice itself, the geyser can be sustained. In this case, the upper ocean beneath the ice becomes stably stratified and thus a barrier to vertical communication, leading to transit timescales from the core to the ice shell of hundreds of years.

According to adaptive-optics observations by Ferrais et al., (22) Kalliope is a 150-km, dense and differentiated body. Here, we interpret (22) Kalliope in the context of bodies in its surroundings. While there is a known moon Linus, with a 5:1 size ratio, no family has been reported in the literature, which is in contradiction with the existence of the moon. Using the hierarchical clustering method (HCM) along with physical data, we identified the Kalliope family. Previously, it was associated to (7481) San Marcello. We then used various models (N-body, Monte-Carlo, SPH) of its orbital and collisional evolution, including the break-up of the parent body, to estimate the dynamical age of the family and address its link to Linus. The best-fit age is (900+-100) My according to our collisional model, in agreement with the position of (22) Kalliope, which was modified by chaotic diffusion due to 4-1-1 three-body resonance with Jupiter and Saturn. It seems possible to create Linus and the Kalliope family at the same time, although our SPH simulations show a variety of outcomes, for both satellite size and the family size-frequency distribution. The shape of (22) Kalliope itself was most likely affected by gravitational reaccumulation of `streams', which creates characteristic hills observed on the surface. If the body was differentiated, its internal structure is surely asymmetric.

This paper presents an analytical model for collision probability assessments between de-orbiting or injecting space objects and satellite constellations. Considering the first to be subjected to a continuous tangential acceleration, its spiraling motion would result in a series of close approaches in the proximity of a constellation. The proposed methodology involves the integration of the collision probability density function on the encounter plane, from which two analytical formulas, one for the number of close approaches and one for their respective average collision probability, are obtained. The mathematical description of the crossing dynamics relies on the assumption of circular orbits and independent collision probabilities, but does not require to propagate the satellites' orbit. A comparison with a conventional propagation method has been performed for validation purposes, proving its accuracy also in case of elliptical crossing orbits. The model developed has been used to assess the risk connected to constellation's satellites replacement, once they have reached their programmed End-of-Life. The environmental impact of the full replacement of 12 approved constellations is analysed by means of average collision probability. In particular, it is shown that the key features for space exploitation sustainability are the maximum propulsion available from the thruster, the selection of an optimal crossing orbit and the true anomaly phases between constellations' and crossing satellites. The consequences of an in-orbit collision are also investigated by assessing the collision risk generated by the formation of a debris cloud. The results corroborate the need for international standards for space traffic management as an exponentially increasing satellites population could trigger a chain reaction of collisions, making LEO inaccessible for decades.

The kinetic Sunyaev Zel'dovich (kSZ) effect, cosmic microwave background (CMB) temperature anisotropies induced by the scattering of CMB photons from free electrons, will be measured by near-term CMB experiments at high significance. By combining CMB temperature anisotropies with a tracer of structure, such as a galaxy redshift survey, previous literature introduced a number of techniques to reconstruct the radial velocity field. This reconstructed radial velocity field encapsulates the majority of the cosmological information contained in the kSZ temperature anisotropies, and can provide powerful new tests of the standard cosmological model and theories beyond it. In this paper, we introduce a new estimator for the radial velocity field based on a coarse-grained maximum likelihood fit for the kSZ component of the temperature anisotropies, given a tracer of the optical depth. We demonstrate that this maximum likelihood estimator yields a higher fidelity reconstruction than existing quadratic estimators in the low-noise and high-resolution regime targeted by upcoming CMB experiments. We describe implementations of the maximum likelihood estimator in both harmonic space and map space, using either direct measurements of the optical depth field or a galaxy survey as a tracer. We comment briefly on the impact of biases introduced by imperfect reconstruction of the optical depth field.

The imprints of stellar nucleosynthesis and chemical evolution of the galaxy can be seen in different stellar populations, with older generation stars showing higher $\alpha$-element abundances while the later generations becoming enriched with iron-peak elements. The evolutionary connections and chemical characteristics of circumstellar disks, stars, and their planetary companions can be inferred by studying the interdependence of planetary and host star properties. Numerous studies in the past have confirmed that high-mass giant planets are commonly found around metal-rich stars, while the stellar hosts of low-mass planets have a wide range of metallicity. In this work, we analyzed the detailed chemical abundances for a sample of $>900$ exoplanet hosting stars drawn from different radial velocity and transit surveys. We correlate the stellar abundance trends for $\alpha$ and iron-peak elements with the planets' mass. We find the planet mass-abundance correlation to be primarily negative for $\alpha$-elements and marginally positive or zero for the iron-peak elements, indicating that stars hosting giant planets are relatively younger. This is further validated by the age of the host stars obtained from isochrone fitting. The later enrichment of protoplanetary material with iron and iron-peak elements is also consistent with the formation of the giant planets via the core accretion process. A higher metal fraction in the protoplanetary disk is conducive to rapid core growth, thus providing a plausible route for the formation of giant planets. This study, therefore, indicates the observed trends in stellar abundances and planet mass are most likely a natural consequence of Galactic chemical evolution.

Exoplanetary properties depend on stellar properties: to know the planet with accuracy and precision it is necessary to know the star as accurately and precisely as possible. Our immediate aim is to characterize in a homogeneous and accurate way a sample of 27 transiting planet-hosting stars observed within the GAPS program. We determined stellar parameters (effective temperature, surface gravity, rotational velocity) and abundances of 26 elements (Li,C,N,O,Na,Mg,Al,Si,S,Ca,Sc,Ti,V,Cr,Fe,Mn,Co,Ni,Cu,Zn,Y,Zr,Ba,La,Nd,Eu). Our study is based on high-resolution HARPS-N@TNG and FEROS@ESO spectra and uniform techniques. We derived kinematic properties from Gaia data and estimated for the first time in exoplanet host stars ages using elemental ratios as chemical clocks. Teff of our stars is of 4400-6700 K, while [Fe/H] is within -0.3 and 0.4 dex. Lithium is present in 7 stars. [X/H] and [X/Fe] abundances vs [Fe/H] are consistent with the Galactic Chemical Evolution. The dependence of [X/Fe] with the condensation temperature is critically analyzed with respect to stellar and kinematic properties. All targets with measured C and O abundances show C/O<0.8, compatible with Si present in rock-forming minerals. Most of targets show 1.0<Mg/Si<1.5, compatible with Mg distributed between olivine and pyroxene. HAT-P-26, the target hosting the lowest-mass planet, shows the highest Mg/Si ratio. From our chemo-dinamical analysis we find agreement between ages and position within the Galactic disk. We note a tendency for higher density planets to be around metal-rich stars and hints of higher stellar abundances of some volatiles for lower mass planets. We cannot exclude that part of our results could be also related to the location of the stars within the Galactic disk. We trace the planetary migration scenario from the composition of the planets related to the chemical composition of the hosting stars

The mass (M) of a star can be evaluated from its spectroscopically determined effective temperature (T _eff) and metallicity ([Fe/H]) along with the luminosity (L; derived from parallax), while comparing them with grids of theoretical evolutionary tracks. It has been argued, however, that such a track-based mass (M_trk) may tend to be overestimated for the case of red giants. Meanwhile, there is an alternative approach of evaluating mass (M_gLT) directly from surface gravity (g), L, and T_eff. The practical reliability of M_gLT was examined for ~100 benchmark giants in the Kepler field, for which atmospheric parameters are already determined and the reliable mass (M_seis) along with the evolutionary status are known from asteroseismology. In addition, similar check was also made for the accuracy of M_trk for comparison. It turned out that, while a reasonable correlation is seen between M_gLT and M_seis almost irrespective of the stellar property, its precision is rather insufficient because log(M_gLT/M_seis) distributes rather widely within ~+/-0.2--0.3dex. In contrast, the reliability of M_trk was found to depend on the evolutionary status. Although M_trk and M_seis are satisfactorily consistent with each other (typical dispersion of log(M_trk}/M_seis) is within ~+/-0.1dex) for H-burning red giants as well as He-burning 2nd clump giants of higher mass, M_trk tends to be considerably overestimated as compared to M_seis by up to ~<0.4~dex for He-burning 1st clump giants of lower mass. Accordingly, M_gLT and M_trk are complementary with each other in terms of their characteristic merit and demerit.

Photospheric abundances of C, N, O, and Na were determined by applying the synthetic spectrum-fitting technique to 34 snap-shot high-dispersion spectra of 22 RR Lyr stars covering a metallicity range of -1.8 <[Fe/H] < 0.0, with an aim of investigating the mixing mechanism in the interior of low-mass giant stars by examining the abundance anomalies of these elements possibly affected by the evolution-induced dredge-up of nuclear burning products. Special attention was paid to check the recent theoretical stellar evolution simulations indicating the importance of thermohaline mixing in low-mass stars (M <~1 M_sun), which is expected to be more significant as the metallicity is lowered. By inspecting the resulting abundances in comparison with those of unevolved metal-poor dwarfs at the same metallicity, the deficiency in C as well as enrichment in N was confirmed (while O is almost unchanged), the extent of peculiarities tending to increase with a decrease in [Fe/H]. Accordingly, the [C/N] ratio turned out to progressively decrease towards lower metallicity from ~0 (Fe/H]~0) to ~-1 ([Fe/H]~-1.5), which is reasonably consistent with the theoretical prediction in the presence of thermohaline mixing. However, these RR Lyr stars do not show any apparent Na anomaly (i.e., essentially the same [Na/Fe] vs. [Fe/H] trends as those of dwarfs), despite that metallicity-dependent overabundance in Na is theoretically expected for the case of non-canonical mixing. This inconsistency between C/N and Na may suggest a necessity of further improvement in the current theory.

We report on a study that compares energetic particle fluxes in corotating interaction regions (CIRs) associated with type III radio storm with those in nonstorm CIRs. In a case study, we compare the CIR particle events on 2010October 21 and 2005 November 2. The two events have similar solar and solar wind circumstances, except that the former is associated with a type III radio storm and has a higher CIR particle flux and fluence. We also perform a statistical study, which shows that the proton and electron fluences are higher in the storm associated CIRs by factor of about 6 and 8, respectively than those in the storm-free CIRs.

Ultra-high-energy (UHE) neutrinos ($>10^{16}$ eV) can be measured cost-effectively using in-ice radio detection, which has been explored successfully in pilot arrays. A large radio detector is currently being constructed in Greenland with the potential to measure the first UHE neutrino, and an order-of-magnitude more sensitive detector is being planned with IceCube-Gen2. For such shallow radio detector stations, we present an end-to-end reconstruction of the neutrino energy and direction using deep neural networks (DNNs). The DNN determines the energy with a standard deviation of a factor of two around the true energy, which meets the science requirements of UHE neutrino detectors. For the first time, we are able to predict the neutrino direction well for all event topologies including the complicated electron neutrino charged-current ($\nu_e$-CC) interactions, a significant improvement compared to previous approaches. The obtained angular resolution follows a Gaussian distribution with $\sigma \approx 0.6^\circ (0.8^\circ)$ with extended tails that push the 68\% quantile to $4^\circ (5^\circ)$ for non-$\nu_e$-CC and $\nu_e$-CC interactions, respectively. This highlights the advantages of DNNs for modeling the complex correlations in radio detector data, thereby enabling measurement of neutrino energy and direction.

In this work, we employ two publicly available analysis tools to study four hydrogen(H)--stripped core--collapse supernovae (CCSNe) namely, SN 2009jf, iPTF13bvn, SN 2015ap, and SN 2016bau. We use the Modular Open-Source Fitter for Transients ({\tt MOSFiT}) to model the multi band light curves. {\tt MOSFiT} analyses show ejecta masses (log M$_{ej}$) of $0.80_{-0.13}^{+0.18}$ M$_{\odot}$, $0.15_{-0.09}^{+0.13}$ M$_{\odot}$, $0.19_{-0.03}^{+0.03}$ M$_{\odot}$, and $0.19_{+0.02}^{-0.01}$ M$_{\odot}$ for SN 2009jf, iPTF13vn, SN 2015ap, and SN 2016au, respectively. Later, Modules for Experiments in Stellar Astrophysics ({\tt MESA}), is used to construct models of stars from pre-main sequence upto core collapse which serve as the possible progenitors of these H-stripped CCSNe. Based on literature, we model a 12 M$_{\odot}$ ZAMS star as the possible progenitor for iPTF13vn, SN 2015ap, and SN 2016bau while a 20 M$_{\odot}$ ZAMS star is modeled as the possible progenitor for SN 2009jf. Glimpses of stellar engineering and the physical properties of models at various stages of their lifetime have been presented to demonstrate the usefulness of these analysis threads to understand the observed properties of several classes of transients in detail.

In this work, we study the synthetic explosions of a massive star. We take a 100 M$_{\odot}$ zero--age main--sequence (ZAMS) star and evolve it until the onset of core-collapse using {\tt MESA}. Then, the resulting star model is exploded using the publicly available stellar explosion code, {\tt STELLA}. The outputs of {\tt STELLA} calculations provide us the bolometric light curve and photospheric velocity evolution along with other physical properties of the underlying supernova. In this paper, the effects of having large Hydrogen-envelope on the supernova light curve have been explored. We also explore the effects of the presence of different amounts of nickel mass and the effect of changing the explosion energy of the resulting supernovae from such heavy progenitors, on their bolometric light curves and photospheric velocities.

Recently, a single-line spectroscopic binary, LTD064402+245919, has been discovered by Yang et al. Using data from LAMOST and ZTF, the unseen companion is estimated to have a mass of 1-3 $M_{\odot}$, orbiting a subgiant with orbital period of 14.50 days, making it a good compact binary candidate without X-ray emission. However, new light curves from ZTF and ASAS-SN, have shown the depth of one dip increases towards a bluer wavelength, indicating LTD064402+245919 is more likely to be a subgiant with a red star. Using both Wilson-Devinney code and Phoebe, the derived $T_{eff}$ of secondary is about 3400 K, corresponding to a red M2/3 star. Additionally, the 20% error of parallax from Gaia is large. The mass of subgiant will be 1.28 $M_{\odot}$ instead of 2.77 $M_{\odot}$, if the refined distance of 5.0kpc is used. Nevertheless, new multi-colour photometry are warranted for the final confirmation of binary properties.

The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (<1 MK) loops seen in the core of an active region to address the role of field-line braiding in energising these structures. We report observations from the High-resolution Coronal imager (Hi-C) that have been acquired in a coordinated campaign with the Interface Region Imaging Spectrograph (IRIS). In the core of the active region, the 172 A band of Hi-C and the 1400 A channel of IRIS show plasma loops at different temperatures that run in parallel. There is a small but detectable spatial offset of less than 1 arcsec between the loops seen in the two bands. Most importantly, we do not see observational signatures that these loops might be twisted around each other. Considering the scenario of magnetic braiding, our observations of parallel loops imply that the stresses put into the magnetic field have to relax while the braiding is applied: the magnetic field never reaches a highly braided state on these length-scales comparable to the separation of the loops. This supports recent numerical 3D models of loop braiding in which the effective dissipation is sufficiently large that it keeps the magnetic field from getting highly twisted within a loop.

Chaotic three-body interactions may lead to the formation of gravitational-wave sources. Here, by modelling the encounter as a series of close, non-hierarchical, triple approaches, interspersed with hierarchical phases, in which the system consists of an inner binary and a star that orbits it, we compute the pericentre probability distribution, and thereby the in-spiral probability in any given binary-single encounter. We then consider the indirect influence of binary-single encounters on the population of gravitational-wave sources, by changing the eccentricity distribution of hard binaries in clusters; we calculate this distribution analytically, by requiring that it be invariant under interactions with single stars.

We use 10 years of publicly available IceCube data to investigate the correlations between hight-energy neutrinos and various Fermi-LAT gamma-ray samples. This work considers the following gamma-ray samples:the third Fermi-LAT catalog of high-energy sources(3FHL), >100GeV Fermi-LAT events, LAT 12-year source catalog(4FGL), the fourth catalog of activate galactic nuclei(4LAC) and subsets of these samples. For each sample, both a single-source analysis and a joint likelihood analysis are performed. We find no indication that the sources in these samples produce significant high-energy neutrinos .From the null search result, we infer that each source population can produce no more than ~0.3%-27% (at the 95% confidence level) of the IceCube's diffuse neutrino flux. Since we are using a larger(10 years) dataset of IceCube neutrinos , the constriants are improved by a factor of ~2 compared to those based on 3 years of data.

Primordial non-Gaussianity can source $\mu$-distortion anisotropies that are correlated with the large-scale temperature and polarization signals of the cosmic microwave background (CMB). A measurement of $\mu T$ and $\mu E$ correlations can therefore be used to constrain it on wavelengths of perturbations not directly probed by the standard CMB anisotropies. In this work, we carry out a first rigorous search for $\mu$-type spectral distortion anisotropies with \Planck data, applying the well-tested constrained ILC component-separation method combined with the needlet framework. We reconstruct a $\mu$ map from \Planck data, which we then correlate with the CMB anisotropies to derive constraints on the amplitude $\fNL$ of the local form bispectrum, specifically on the highly squeezed configurations with effective wavenumbers $k_s \simeq \SI{740}{Mpc^{-1}}$ and $k_L \simeq \SI{0.05}{Mpc^{-1}}$. We improve previously estimated constraints by more than an order of magnitude. This enhancement is owing to the fact that for the first time we are able to use the full multipole information by carefully controlling biases and systematic effects in the final analysis. We also for the first time incorporate constraints from measurements of $\mu E$ correlations, which further tighten the limits. A combination of the derived \Planck $\mu T$ and $\mu E$ power spectra yields $|\fNL| \lesssim 6800$ (95\% c.l.) on this highly squeezed bispectrum. This is only $\simeq 3$ times weaker than the anticipated constraint from \LiteBIRD alone. We show that a combination of \LiteBIRD with \Planck will improve the expected future constraint by $\simeq 20\%$ over \LiteBIRD alone. These limits can be used to constrain multi-field inflation models and primordial black hole formation scenarios, thus providing a promising novel avenue forward in CMB cosmology.

Modern astronomical potentials modeling galaxies or stellar systems can be rather involved, and deriving their first derivatives (accelerations) and second derivatives (variational equations) in order to compute orbits and their chaoticity may be a formidable task. We present here a fully automated routine, dubbed Smart, with which the accelerations and the variational equations of an arbitrary potential that has been written in the Fortran 77 language can be computed. Almost any Fortran 77 statement is admitted in the potential, and the output are standard Fortran 77 routines ready to use. We validate our algorithm with a set of potentials including time-dependent, velocity-dependent and very complex potentials that even involve auxiliary routines. We also describe with some detail a realistic seven-component Galactic potential, MilkyWayHydra, which yields very involved derivatives, thus being a good test bed for Smart.

Upcoming neutrino telescopes may discover ultra-high-energy (UHE) cosmic neutrinos, with energies beyond 100~PeV, in the next 10--20 years. Finding their sources would expose the long-sought origin of UHE cosmic rays. We search for sources by looking for multiplets of UHE neutrinos arriving from similar directions. Our forecasts are state-of-the-art, geared at neutrino radio-detection in IceCube-Gen2. They account for detector energy and angular response, and for critical, but uncertain backgrounds. We report powerful insight. Sources at declination of $-45^\circ$ to $0^\circ$ will be easiest to discover. Discovering even one steady-state source in 10~years would disfavor most known steady-state source classes as dominant. Discovering no transient source would disfavor most known transient source classes as dominant. Our results aim to inform the design of upcoming detectors.

A key ingredient for semi-analytic models (SAMs) of galaxy formation is the mass assembly history of haloes, encoded in a tree structure. The most commonly used method to construct halo merger histories is based on the outcomes of high-resolution, computationally intensive N-body simulations. We show that machine learning (ML) techniques, in particular Generative Adversarial Networks (GANs), are a promising new tool to tackle this problem with a modest computational cost and retaining the best features of merger trees from simulations. We train our GAN model with a limited sample of merger trees from the EAGLE simulation suite, constructed using two halo finders-tree builder algorithms: SUBFIND-D-TREES and ROCKSTAR-ConsistentTrees. Our GAN model successfully learns to generate well-constructed merger tree structures with high temporal resolution, and to reproduce the statistical features of the sample of merger trees used for training, when considering up to three variables in the training process. These inputs, whose representations are also learned by our GAN model, are mass of the halo progenitors and the final descendant, progenitor type (main halo or satellite) and distance of a progenitor to that in the main branch. The inclusion of the latter two inputs greatly improves the final learned representation of the halo mass growth history, especially for SUBFIND-like ML trees. When comparing equally sized samples of ML merger trees with those of the EAGLE simulation, we find better agreement for SUBFIND-like ML trees. Finally, our GAN-based framework can be utilised to construct merger histories of low and intermediate mass haloes, the most abundant in cosmological simulations.

The sky observed by space telescopes in Low Earth Orbit (LEO) can be dominated by stray light from multiple sources including the Earth, Sun and Moon. This stray light presents a significant challenge to missions that aim to make a secure measurement of the Extragalactic Background Light (EBL). In this work we quantify the impact of stray light on sky observations made by the Hubble Space Telescope (HST) Advanced Camera for Surveys. By selecting on orbital parameters we successfully isolate images with sky that contain minimal and high levels of Earthshine. In addition, we find weather observations from CERES satellites correlates with the observed HST sky surface brightness indicating the value of incorporating such data to characterise the sky. Finally we present a machine learning model of the sky trained on the data used in this work to predict the total observed sky surface brightness. We demonstrate that our initial model is able to predict the total sky brightness under a range of conditions to within 3.9% of the true measured sky. Moreover, we find that the model matches the stray light-free observations better than current physical Zodiacal light models.

Experimental and theoretical results are presented for double, triple, and quadruple photoionization of Si$^+$ and Si$^{2+}$ ions and for double photoionization of Si$^{3+}$ ions by a single photon. The experiments employed the photon-ion merged-beams technique at a synchrotron light source. The experimental photon-energy range 1835--1900 eV comprises resonances associated with the excitation of a $1s$ electron to higher subshells and subsequent autoionization. Energies, widths, and strengths of these resonances are extracted from high-resolution photoionization measurements, and the core-hole lifetime of K-shell ionized neutral silicon is inferred. In addition, theoretical cross sections for photoabsorption and multiple photoionization were obtained from large-scale Multi-Configuration Dirac-Hartree-Fock (MCDHF) calculations. The present calculations agree with the experiment much better than previously published theoretical results. The importance of an accurate energy calibration of laboratory data is pointed out. The present benchmark results are particularly useful for discriminating between silicon absorption in the gaseous and in the solid component (dust grains) of the interstellar medium.

It has been recently pointed out that in certain axion models it is possible to suppress simultaneously both the axion couplings to nucleons and electrons, realising the so-called astrophobic axion scenarios, wherein the tight bounds from SN1987A and from stellar evolution of red giants and white dwarfs are greatly relaxed. So far, however, the conditions for realising astrophobia have only been set out in tree-level analyses. Here we study whether these conditions can still be consistently implemented once renormalization group effects are included in the running of axion couplings. We find that axion astrophobia keeps holding, albeit within fairly different parameter space regions, and we provide transparent insights into this result. Given that astrophobic axion models generally feature flavour violating axion couplings, we also assess the impact of renormalization group effects on axion-mediated flavour violating observables.

Complex metrics are a double-edged sword: they allow one to replace singular spacetimes, such as those containing a big bang, with regular metrics, yet they can also describe unphysical solutions in which quantum transitions may be more probable than ordinary classical evolution. In the cosmological context, we investigate a criterion proposed by Witten (based on works of Kontsevich & Segal and of Louko & Sorkin) to decide whether a complex metric is allowable or not. Because of the freedom to deform complex metrics using Cauchy's theorem, deciding whether a metric is allowable in general requires solving a complicated optimisation problem. We describe a method that allows one to quickly determine the allowability of minisuperspace metrics. This enables us to study the off-shell structure of minisuperspace path integrals, which we investigate for various boundary conditions. Classical transitions always reside on the boundary of the domain of allowable metrics, and care must be taken in defining appropriate integration contours for the corresponding gravitational path integral. Perhaps more surprisingly, we find that proposed quantum (`tunnelling') transitions from a contracting to an expanding universe violate the allowability criterion and may thus be unphysical. No-boundary solutions, by contrast, are found to be allowable, and moreover we demonstrate that with an initial momentum condition an integration contour over allowable metrics may be explicitly described in arbitrary spacetime dimensions.

A sub-fraction of dark matter or new particles trapped inside celestial objects can significantly alter their macroscopic properties. We investigate the new physics imprint on celestial objects by using a generic framework to solve the Tolman-Oppenheimer-Volkoff (TOV) equations for up to two fluids. We test the impact of populations of new particles on celestial objects, including the sensitivity to self-interaction sizes, new particle mass, and net population mass. Applying our setup to neutron stars and boson stars, we find rich phenomenology for a range of these parameters, including the creation of extended atmospheres. These atmospheres are detectable by their impact on the tidal love number, which can be measured at upcoming gravitational wave experiments such as Advanced LIGO, the Einstein Telescope, and LISA. We release our calculation framework as a publicly available code, allowing the TOV equations to be generically solved for arbitrary new physics models in novel and admixed celestial objects.

In this work we shall consider the effects of a geometrically generated post-inflationary era on the energy spectrum of the primordial gravitational waves. Specifically, we shall consider a post-inflationary constant equation of state era, generated by the synergistic effect of $f(R)$ gravity and of radiation and matter perfect fluids. Two cases of interest shall be studied, one with equation of state parameter $w=-1/3$, in which case the Universe neither accelerates nor decelerates, and one with $w=0$ so an early matter domination era. For the evaluation of the inflationary observational indices which is relevant for the calculation of the gravitational waves energy spectrum, we also took into account the effects of the constant equation of state parameter era, on the $e$-foldings number. In both the $w=-1/3$ and $w=0$ cases, the energy spectrum of the primordial gravitational waves is amplified, but for the $w=0$ case, the effect is stronger.

Water has been stored in the Martian mantle since its formation, primarily in nominally anhydrous minerals. The short-lived early hydrosphere and intermittently flowing water on the Martian surface may have been supplied and replenished by magmatic degassing of water from the mantle. Estimating the water storage capacity of the solid Martian mantle places important constraints on its water inventory and helps elucidate the sources, sinks, and temporal variations of water on Mars. In this study, we applied a bootstrap aggregation method to investigate the effects of iron on water storage capacities in olivine, wadsleyite, and ringwoodite, based on high-pressure experimental data compiled from the literature, and we provide a quantitative estimate of the upper bound of the bulk water storage capacity in the FeO-rich solid Martian mantle. Along a series of areotherms at different mantle potential temperatures ($T_{p}$), we estimated a water storage capacity equal to $9.0_{-2.2} ^{+2.8}$ km Global Equivalent Layer (GEL) for the present-day Martian mantle at $T_{p}$ = 1600 K and $4.9_{-1.5}^{+1.7}$ km GEL for the initial Martian mantle at $T_{p}$ = 1900 K. The water storage capacity of the Martian mantle increases with secular cooling through time, but due to the lack of an efficient water recycling mechanism on Mars, its actual mantle water content may be significantly lower than its water storage capacity today.

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

Io's movement relative to the plasma in Jupiter's magnetosphere creates Alfv\'en waves propagating along the magnetic field lines which are partially reflected along their path. These waves are the root cause for auroral emission, which is subdivided into the Io Footprint (IFP), its tail and leading spot. New observations of the Juno spacecraft by Mura et al. (2018) have shown puzzling substructure of the footprint and its tail. In these observations, the symmetry between the poleward and equatorward part of the footprint tail is broken and the tail spots are alternatingly displaced. We show that the location of these bright spots in the tail are consistent with Alfv\'en waves reflected at the boundary of the Io torus and Jupiter's ionosphere. Then, we investigate three different mechanisms to explain this phenomenon: (1) The Hall effect in Io's ionosphere, (2) travel time differences of Alfv\'en waves between Io's Jupiter facing and its opposing side and (3) asymmetries in Io's atmosphere. For that, we use magnetohydrodynamic simulations within an idealized geometry of the system. We use the Poynting flux near the Jovian ionosphere as a proxy for the morphology of the generated footprint and its tail. We find that the Hall effect is the most important mechanism under consideration to break the symmetry causing the "Alternating Alfv\'en spot street". The travel time differences contributes to enhance this effect. We find no evidence that the inhomogeneities in Io's atmosphere contribute significantly to the location or shape of the tail spots.

The recently observed accelerated expansion of the universe has put a challenge for its theoretical understanding. As a possible explanation of this, it is considered that the most part of the present universe is filled with a form of energy that exerts a negative pressure called dark energy, which drives the acceleration. In the present work, we assume a dynamical dark energy model, where dark energy interacts with matter and grows at the expense of the latter. Using this model, we discuss the evolution of the universe within the context of loop quantum cosmology. Our work successfully explains the presently observed accelerated expansion of the universe, by predicting that the present universe is phantom dominated. We also found that in the past, the expansion of the universe was decelerated one and transition from deceleration to acceleration would occur at $t_{q=0}=0.688t_0$, where $t_0$ is the present age of the universe. Again, our analysis predicted that at the transition time, the universe would be dominated with quintessence type dark energy.

We study the generation and propagation of gravitational waves in scalar-tensor gravity using numerical relativity simulations of scalar field collapses beyond spherical symmetry. This allows us to compare the tensor and additional massive scalar waves that are excited. As shown in previous work in spherical symmetry, massive propagating scalar waves decay faster than 1/r and disperse, resulting in an inverse chirp. These effects obscure the ringdown in any extracted signal by mixing it with the transient responses of the collapse during propagation. In this paper we present a simple method to rewind the extracted signals to horizon formation, which allows us to clearly identify the ringdown phase and extract the amplitudes of the scalar quasinormal modes, quantifying their excitation in strong gravity events and verifying the frequencies to perturbative calculations. The effects studied are relevant to any theories in which the propagating waves have a dispersion relation, including the tensor case.

We present an improved measurement of the CNO solar neutrino interaction rate at Earth obtained with the complete Borexino Phase-III dataset. The measured rate R$_{\rm CNO}$ = $6.7^{+2.0}_{-0.8}$ counts/(day$ \cdot$ 100 tonnes), allows us to exclude the absence of the CNO signal with about 7$\sigma$ C.L. The correspondent CNO neutrino flux is $6.6^{+2.0}_{-0.9} \times 10^8$ cm$^{-2}$ s$^{-1}$, taking into account the neutrino flavor conversion. We use the new CNO measurement to evaluate the C and N abundances in the Sun with respect to the H abundance for the first time with solar neutrinos. Our result of $N_{\rm CN}$ = $(5.78^{+1.86}_{-1.00})\times10^{-4}$ displays a $\sim$2$\sigma$ tension with the "low metallicity" spectroscopic photospheric measurements. On the other hand, our result used together with the $^7$Be and $^8$B solar neutrino fluxes, also measured by Borexino, permits to disfavour at 3.1$\sigma$ C.L. the "low metallicity" SSM B16-AGSS09met as an alternative to the "high metallicity" SSM B16-GS98.