We introduce Astrophysical Hybrid-Kinetic simulations with the flash code ([Formula: see text]) - a new Hybrid particle-in-cell (PIC) code developed within the framework of the multiphysics code flash. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a [Formula: see text] method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfvén and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the turbgen turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the 'hybrid precision' method, which significantly reduces (by a factor [Formula: see text]) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000 cores.

In the context of upcoming large-scale surveys like Euclid, the necessity for the automation of strong lens detection is essential. While existing machine learning pipelines heavily rely on the classification probability (P), this study intends to address the importance of integrating additional metrics, such as Information Content (IC) and the number of pixels above the segmentation threshold ([Formula: see text]), to alleviate the false positive rate in unbalanced data-sets. In this work, we introduce a segmentation algorithm (U-Net) as a supplementary step in the established strong gravitational lens identification pipeline (Denselens), which primarily utilizes [Formula: see text] and [Formula: see text] parameters for the detection and ranking. The results demonstrate that the inclusion of segmentation enables significant reduction of false positives by approximately 25 per cent in the final sample extracted from DenseLens, without compromising the identification of strong lenses. The main objective of this study is to automate the strong lens detection process by integrating these three metrics. To achieve this, a decision tree-based selection process is introduced, applied to the Kilo Degree Survey (KiDS) data. This process involves rank-ordering based on classification scores ([Formula: see text]), filtering based on Information Content ([Formula: see text]), and segmentation score ([Formula: see text]). Additionally, the study presents 14 newly discovered strong lensing candidates identified by the U-Denselens network using the KiDS DR4 data.

Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of time-scales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here, we present a new implicit-explicit method for numerical RHD that has a number of desirable properties that have not previously been combined in a single method. Our scheme is based on moments and allows machine-precision conservation of energy and momentum, making it highly suitable for adaptive mesh refinement applications; it requires no more communication than hydrodynamics and includes no non-local iterative steps, making it highly suitable for massively parallel and Graphics Processing Unit (GPU)-based systems where communication is a bottleneck; and we show that it is asymptotically accurate in the streaming, static diffusion, and dynamic diffusion limits, including in the so-called asymptotic diffusion regime where the computational grid does not resolve the photon mean-free path. We implement our method in the GPU-accelerated RHD code quokka and show that it passes a wide range of numerical tests.

Variations in chemical abundances with evolutionary phase have been identified among stars in globular and open clusters with a wide range of metallicities. In the metal-poor clusters, these variations compare well with predictions from stellar structure and evolution models considering the internal diffusive motions of atoms and ions, collectively known as atomic diffusion, when moderated by an additional mixing process with a fine-tuned efficiency. We present here an investigation of these effects in the Galactic globular cluster NGC 6121 (M4) ([Fe/H] = -1.13) through a detailed chemical abundance analysis of 86 stars using high-resolution ESO Very Large Telescope (VLT) Fibre Large Array Multi Element Spectrograph (FLAMES) spectroscopy. The stars range from the main-sequence turnoff point (TOP) to the red giant branch (RGB) just above the bump. We identify C-N-O and Mg-Al-Si abundance anticorrelations, and confirm the presence of a bimodal population differing by 1 dex in nitrogen abundance. The composition of the second-generation stars imply pollution from both massive (20-40 [Formula: see text]) and asymptotic giant branch stars. We find evolutionary variations in chemical abundances between the TOP and RGB, which are robust to uncertainties in stellar parameters and modelling assumptions. The variations are weak, but match predictions well when employing efficient additional mixing. Without correcting for Galactic production of lithium, we derive an initial lithium abundance 2.63 ± 0.10, which is marginally lower than the predicted primordial big-bang nucleosynthesis value.

To fully take advantage of the data provided by large-scale structure surveys, we need to quantify the potential impact of baryonic effects, such as feedback from active galactic nuclei (AGN) and star formation, on cosmological observables. In simulations, feedback processes originate on scales that remain unresolved. Therefore, they need to be sourced via subgrid models that contain free parameters. We use machine learning to calibrate the AGN and stellar feedback models for the FLAMINGO (Fullhydro Large-scale structure simulations with All-sky Mapping for the Interpretation of Next Generation Observations) cosmological hydrodynamical simulations. Using Gaussian process emulators trained on Latin hypercubes of 32 smaller volume simulations, we model how the galaxy stellar mass function (SMF) and cluster gas fractions change as a function of the subgrid parameters. The emulators are then fit to observational data, allowing for the inclusion of potential observational biases. We apply our method to the three different FLAMINGO resolutions, spanning a factor of 64 in particle mass, recovering the observed relations within the respective resolved mass ranges. We also use the emulators, which link changes in subgrid parameters to changes in observables, to find models that skirt or exceed the observationally allowed range for cluster gas fractions and the SMF. Our method enables us to define model variations in terms of the data that they are calibrated to rather than the values of specific subgrid parameters. This approach is useful, because subgrid parameters are typically not directly linked to particular observables, and predictions for a specific observable are influenced by multiple subgrid parameters.

One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter haloes. There is a large volume of information in strong gravitational lens images - the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter haloes relative to warm dark matter. Since the strong lensing signal receives a contribution from a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter haloes, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with 40 lenses with Bayesian odds of 30:1.

Cosmic rays in star-forming galaxies are a dominant source of both diffuse γ-ray emission and ionization in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for γ-rays and ionization are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, γ-ray luminosities, and ionization rates should all be linked. In this paper, we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate [Formula: see text] and gas depletion time produce a maximum primary ionization rate ζ ≈ 1 × 10( /Gyr) s and a maximum γ-ray luminosity [Formula: see text] erg s in the 0.1-100 GeV band. These budgets imply either that the ionization rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionization in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionization rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of γ-ray luminosities can be used to place constraints on galactic ionization budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration.

At least one in five of all planetary nebulae are the product of a common envelope (CE) interaction, where the companion in-spirals into the envelope of an asymptotic giant branch (AGB) star ejecting the nebula and leaving behind a compact binary. In this work we carry out 3D smoothed particle hydrodynamics simulations of the CE interaction between a 1.7 M AGB star and a 0.6 M companion. We model the AGB structure using a 1D stellar model taken at the seventh thermal pulse. The interaction takes place when the giant is on the expanding phase of the seventh thermal pulse and has a radius of 250 R. The post-CE orbital separations varies between 20 and 31 R, with the inclusion of recombination energy resulting in wider separations. Based on the observed short in-spiral time-scales, we suggest that thermal pulses can trigger CEs, extending the ability of AGB stars to capture companions into CEs, that would lead to the prediction of a larger population of post-AGB, post-CE binaries. Simulations that include a tabulated equation of state unbind a great deal more gas, likely unbinding the entire envelope on short time-scales. The shape of the CE after the in-spiral is more spherical for AGB than red giant branch stars, and even more so if recombination energy is included. We expect that the planetary nebula formed from this CE will have different features from those predicted by Zou et al. 2020.

We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study environmental variation of the initial mass function (IMF). The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of lower mass objects more effectively as the surface density increases, but this only partially compensates for the decreasing Jeans mass in denser environments. Similarly, we find that protostellar outflows are more effective at suppressing the formation of massive stars in higher surface density environments. The combined effect of these two trends is towards an IMF with a lower characteristic mass and a narrower overall mass range in high surface density environments. We discuss the implications for these findings for the interpretation of observational evidence of IMF variation in early type galaxies.

The latticework structure known as the cosmic web provides a valuable insight into the assembly history of large-scale structures. Despite the variety of methods to identify the cosmic web structures, they mostly rely on the assumption that galaxies are embedded in a Euclidean geometric space. Here, we present a novel cosmic web identifier called sconce (pherical and ic osmic wb finder) that inherently considers the 2D (RA, DEC) spherical or the 3D (RA, DEC, ) conic geometry. The proposed algorithms in sconce generalize the well-known subspace constrained mean shift (scms) method and primarily address the predominant filament detection problem. They are intrinsic to the spherical/conic geometry and invariant to data rotations. We further test the efficacy of our method with an artificial cross-shaped filament example and apply it to the SDSS galaxy catalogue, revealing that the 2D spherical version of our algorithms is robust even in regions of high declination. Finally, using -body simulations from Illustris, we show that the 3D conic version of our algorithms is more robust in detecting filaments than the standard scms method under the redshift distortions caused by the peculiar velocities of haloes. Our cosmic web finder is packaged in python as sconce-scms and has been made publicly available.

Gamma-ray bursts (GRBs), observed to redshift = 9.4, are potential probes of the largely unexplored ∼ 2.7-9.4 part of the early Universe. Thus, finding relevant relations among GRB physical properties is crucial. We find that the Platinum GRB data compilation, with 50 long GRBs (with relatively flat plateaus and no flares) in the redshift range 0.553 ≤ ≤ 5.0, and the LGRB95 data compilation, with 95 long GRBs in 0.297 ≤ ≤ 9.4, as well as the 145 GRB combination of the two, strongly favour the 3D Fundamental Plane (Dainotti) correlation (between the peak prompt luminosity, the luminosity at the end of the plateau emission, and its rest-frame duration) over the 2D one (between the luminosity at the end of the plateau emission and its duration). The 3D Dainotti correlations in the three data sets are standardizable. We find that while LGRB95 data have ∼50 per cent larger intrinsic scatter parameter values than the better-quality Platinum data, they provide somewhat tighter constraints on cosmological-model and GRB-correlation parameters, perhaps solely due to the larger number of data points, 95 versus 50. This suggests that when compiling GRB data for the purpose of constraining cosmological parameters, given the quality of current GRB data, intrinsic scatter parameter reduction must be balanced against reduced sample size.

We investigate the environmental properties of the  = 0 radio galaxy population using the simba cosmological hydrodynamic simulation. We identify centrals and satellites from a population of high and low excitation radio galaxies (HERGs and LERGs) in simba, and study their global properties. We find that [Formula: see text] of radio galaxies are satellites, and that there are insignificant differences in the global properties of LERGs based on their central/satellite classification. HERG satellites display lower values of star formation, 1.4 GHz radio luminosity, and Eddington fractions than HERG centrals. We further investigate the environments of radio galaxies and show that HERGs typically live in less dense environments, similar to star-forming galaxies. The environments of high-mass LERGs are similar to non-radio galaxies, but low-mass LERGs live in underdense environments similar to HERGs. LERGs with overmassive black holes reside in the most dense environments, while HERGs with overmassive black holes reside in underdense environments. The richness of a LERG's environment decreases with increasing Eddington fraction, and the environments of all radio galaxies do not depend on radio luminosity for [Formula: see text]. Complementing these results, we find that LERGs cluster on the same scale as the total galaxy population, while multiple HERGs are not found within the same dark matter halo. Finally, we show that high density environments support the growth of HERGs rather than LERGs at  = 2. Simba predicts that with more sensitive surveys, we will find populations of radio galaxies in environments much similar to the total galaxy population.

The hadron-quark phase transition in quantum chromodynamics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of [Formula: see text]. We find explosions only for two low-compactness models (14 and [Formula: see text]) with the DD2F_SF EoS, both with low explosion energies of [Formula: see text]. These weak explosions are characterized by a neutrino signal with several minibursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition-driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as Zr for the [Formula: see text] zero-metallicity model and Zr for the [Formula: see text] solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the protocompact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to 4 /baryon and may have a distinctive gravitational-wave signature, also for a second-order phase transition.

Depending on the stellar densities, protoplanetary discs in stellar clusters undergo: background heating; disc truncation-driven by stellar encounter; and photoevaporation. Disc truncation leads to reduced characteristic sizes and disc masses that eventually halt gas giant planet formation. We investigate how disc truncation impacts planet formation via pebble-based core accretion paradigm, where pebble sizes were derived from the full grain-size distribution within the disc lifetimes. We make the best-case assumption of one embryo and one stellar encounter per disc. Using planet population synthesis techniques, we find that disc truncation shifts the disc mass distributions to the lower margins. This consequently lowered the gas giant occurrence rates. Despite the reduced gas giant formation rates in clustered discs, the encounter models mostly show as in the isolated field; the cold Jupiters are more frequent than the hot Jupiters, consistent with observation. Moreover, the ratio of hot to cold Jupiters depend on the periastron distribution of the perturbers with linear distribution in periastron ratio showing enhanced hot to cold Jupiters ratio in comparison to the remaining models. Our results are valid in the best-case scenario corresponding to our assumptions of: only one disc encounter with a perturber, ambient background heating, and less rampant photoevaporation. It is not known exactly of how much gas giant planet formation would be affected should disc encounter, background heating, and photoevaporation act in a concert. Thus, our study will hopefully serve as motivation for quantitative investigations of the detailed impact of stellar cluster environments on planet formations.

We show that the Platinum gamma-ray burst (GRB) data compilation, probing the redshift range 0.553 ≤ ≤ 5.0, obeys a cosmological-model-independent three-parameter Fundamental Plane (Dainotti) correlation and so is standardizable. While they probe the largely unexplored ∼ 2.3-5 part of cosmological redshift space, the GRB cosmological parameter constraints are consistent with, but less precise than, those from a combination of baryon acoustic oscillation (BAO) and Hubble parameter [()] data. In order to increase the precision of GRB-only cosmological constraints, we exclude common GRBs from the larger Amati-correlated A118 data set composed of 118 GRBs and jointly analyse the remaining 101 Amati-correlated GRBs with the 50 Platinum GRBs. This joint 151 GRB data set probes the largely unexplored ∼ 2.3-8.2 region; the resulting GRB-only cosmological constraints are more restrictive, and consistent with, but less precise than, those from ()  + BAO data.

We use the results of relativistic hydrodynamic simulations of jet-interstellar medium (ISM) interactions in a galaxy with a radio-loud AGN to quantify the extent of ionization in the central few kpcs of the gaseous galactic disc. We perform post-process radiative transfer of AGN radiation through the simulated gaseous jet-perturbed disc to estimate the extent of photo-ionization by the AGN with an incident luminosity of 10 erg s. We also map the gas that is collisionally ionized due to shocks driven by the jet. The analysis was carried out for simulations with similar jet power (10 erg s) but different jet orientations with respect to the gas disc. We find that the shocks from the jets can ionize a significant fraction (up to 33 [Formula: see text]) of dense gas ([Formula: see text]) in the disc, and that the jets clear out the central regions of gas for AGN radiation to penetrate to larger distances in the disc. Jets inclined towards the disc plane couple more strongly with the ISM and ionize a larger fraction of gas in the disc as compared to the vertical jet. However, similar to previous studies, we find that the AGN radiation is quickly absorbed by the outer layers of dense clouds in the disc, and is not able to substantially ionize the disc on a global scale. Thus, compared to jet-ISM interactions, we expect that photo-ionization by the AGN radiation only weakly affects the star-formation activity in the central regions of the galactic disc (≲ 1 kpc), although the jet-induced shocks can spread farther out.

We investigate the spatial, temporal, and spectral properties of 10 microflares from AR12721 on 2018 September 9 and 10 observed in X-rays using the Nuclear Spectroscopic Telescope ARray and the Solar Dynamic Observatory's Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager. We find GOES sub-A class equivalent microflare energies of 10-10 erg reaching temperatures up to 10 MK with consistent quiescent or hot active region (AR) core plasma temperatures of 3-4 MK. One microflare (SOL2018-09-09T10:33), with an equivalent GOES class of A0.1, has non-thermal hard X-ray emission during its impulsive phase (of non-thermal power ~7 × 10 erg s) making it one of the faintest X-ray microflares to have direct evidence for accelerated electrons. In 4 of the 10 microflares, we find that the X-ray time profile matches fainter and more transient sources in the extreme-ultraviolet, highlighting the need for observations sensitive to only the hottest material that reaches temperatures higher than those of the AR core (>5 MK). Evidence for corresponding photospheric magnetic flux cancellation/emergence present at the footpoints of eight microflares is also observed.

WD 0145+234 is a white dwarf that is accreting metals from a circumstellar disc of planetary material. It has exhibited a substantial and sustained increase in 3-5 [Formula: see text]m flux since 2018. Follow-up photometry reveals that emission from the disc had begun to decrease by late 2019. Stochastic brightening events superimposed on the decline in brightness suggest the liberation of dust during collisional evolution of the circumstellar solids. A simple model is used to show that the observations are indeed consistent with ongoing collisions. Rare emission lines from circumstellar gas have been detected at this system, supporting the emerging picture of white dwarf debris discs as sites of collisional gas and dust production.

Accurate line lists are important for the description of the spectroscopic nature of small molecules. While a line list for CN (an important molecule for chemistry and astrophysics) exists, no underlying energy spectroscopic model has been published, which is required to consider the sensitivity of transitions to a variation of the proton-to-electron mass ratio. Here we have developed a Duo energy spectroscopic model as well as a novel hybrid style line list for CN and its isotopologues, combining energy levels that are derived experimentally (Marvel), using the traditional/perturbative approach (Mollist), and the variational approach (from a Duo spectroscopic model using standard ExoMol methodology). The final Trihybrid ExoMol-style line list for CN consists of 28 004 energy levels (6864 experimental, 1574 perturbative, the rest variational) and 2285 103 transitions up to 60 000 cm between the three lowest electronic states (X Σ, A Π, and B Σ). The spectroscopic model created is used to evaluate CN as a molecular probe to constrain the variation of the proton-to-electron mass ratio; no overly promising sensitive transitions for extragalactic study were identified.

We revisit Bondi accretion - steady-state, adiabatic, spherical gas flow on to a Schwarzschild black hole at rest in an asymptotically homogeneous medium - for stiff polytropic equations of state (EOSs) with adiabatic indices Γ > 5/3. A general relativistic treatment is required to determine their accretion rates, for which we provide exact expressions. We discuss several qualitative differences between results for soft and stiff EOSs - including the appearance of a minimum steady-state accretion rate for EOSs with Γ ≥ 5/3 - and explore limiting cases in order to examine these differences. As an example, we highlight results for Γ = 2, which is often used in numerical simulations to model the EOS of neutron stars. We also discuss a special case with this index, the ultrarelativistic 'causal' EOS, = . The latter serves as a useful limit for the still undetermined neutron star EOS above nuclear density. The results are useful, for example, to estimate the accretion rate on to a mini-black hole residing at the centre of a neutron star.

Wavelengths of absorption lines in the spectra of galaxies along the line of sight to distant quasars can be used to probe the variablility of the fine structure constant, , at high redshifts, provided that the laboratory wavelengths are known to better than 6 parts in 10, corresponding to a radial velocity of ≈20 ms. For several lines of Si ii, C ii, Fe i, and Ni ii, previously published wavelengths are inadequate for this purpose. Improved wavelengths for these lines were derived by re-analysing archival Fourier transform (FT) spectra of iron hollow cathode lamps (HCL), a silicon carbide Penning discharge lamp, and with new spectra of nickel HCLs. By reoptimizing the energy levels of Fe i, the absolute uncertainty of 13 resonance lines has been reduced by over a factor of 2. A similar analysis for Si ii gives improved values for 45 lines with wavelength uncertainties over an order of magnitude scer than previous measurements. Improved wavelengths for eight lines of Ni ii were measured and Ritz wavelengths from optimized energy levels determined for an additional three lines at shorter wavelengths. Three lines of C ii near 135 nm were observed using FT spectroscopy and the wavelengths confirm previous measurements.