2024-10-29 | | Total: 114
Galactic-scale simulations rely on sub-grid models to provide prescriptions for the coupling between supernova (SN) feedback and the interstellar medium (ISM). Many of these models are computed in 1-D to allow for an efficient way to account for the variability of properties of their local environment. However, small-scale simulations revealed that the release of energy from SNe within molecular clouds can be highly asymmetrical. This is largely due to the presence of pre-SN feedback, such as ionizing radiation, that are able to carve cavities and channels around the progenitors prior to their detonation. Being partially confined, the SN energy escapes into the outer ISM preferentially through these channels, departing from the spherically symmetric 1-D descriptions. To understand by how much the feedback output could differ, we present a theoretical model for a semi-confined SN. The problem concerns a SN expanding into an evolved HII region, bounded by a molecular cloud with pre-existing vents. With the aid of simple 3-D hydrodynamical simulations, we show that this mode of energy release increases the dynamical impact of the outflows, and extends the timescales over which the SN is energetically coupled to the surrounding matter. We also show that the amount of small-scale solenoidal turbulence driven by semi-confined SNe may be amplified.
The thermal Sunyaev-Zel'dovich effect (tSZ) is a sensitive probe of cosmology, as it traces the abundance of galaxy clusters and groups in the late-time Universe. Upcoming cosmic microwave background experiments such as the Simons Observatory (SO) and CMB-S4 will provide low-noise and high-resolution component-separated tSZ maps covering a large sky fraction. The tSZ signal is highly non-Gaussian; therefore, higher-order statistics are needed to optimally extract information from these maps. In this work, we study the cosmological constraining power of several tSZ statistics -- Minkowski functionals (MFs), peaks, minima, and moments -- that have yielded promising results in capturing non-Gaussian information from other cosmological data. Using a large suite of halo-model-based tSZ simulations with varying $\Omega_{c}$ and $\sigma_{8}$ (154 cosmologies and over $800, 000$ maps, each $10.5\times10.5$ deg$^{2}$), we show that by combining these observables, we can achieve $\approx 29\times$ tighter constraints compared to using the tSZ power spectrum alone in an idealized noiseless case, with the MFs dominating the constraints. We show that much of the MF constraining power arises from halos below the detection threshold of cluster surveys, suggesting promising synergies with cluster-count analyses. Finally, we demonstrate that these statistics have the potential to deliver tight constraints even in the presence of noise. For example, using post-component-separation tSZ noise expected for SO, we obtain $\approx1.6\times$ and $\approx1.8\times$ tighter constraints than the power spectrum with MFs and all statistics combined, respectively. We show that the constraints from MFs approach the noiseless case for white-noise levels $\lesssim 1 \,\, \mu$K-arcmin.
We use a suite of the most recent cosmological observations to test models of dynamical dark energy motivated by quantum gravity. Specifically, we focus on hilltop quintessence scenarios, able to satisfy theoretical constraints from quantum gravity. We discuss their realisation based on axions, their supersymmetric partners, and Higgs-like string constructions. We also examine a specific parameterisation for dynamical dark energy suitable for hilltop quintessence. We then perform an analysis based on Markov Chain Monte-Carlo to assess their predictions against CMB, galaxy surveys, and supernova data. We show to what extent current data can distinguish amongst different hilltop set-ups, providing model parameter constraints that are complementary to and synergetic with theoretical bounds from quantum gravity conjectures, as well as model comparisons across the main dark energy candidates in the literature. However, all these constraints are sensitive to priors based on theoretical assumptions about viable regions of parameter space. Consequently, we discuss theoretical challenges in refining these priors, with the aim of maximizing the informative power of current and forthcoming cosmological datasets for testing dark energy scenarios in quantum gravity.
After reionization, neutral hydrogen (HI) traces the large-scale structure (LSS) of the Universe, enabling HI intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study HI clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of $k$-nearest neighbor ($k$NN) optical galaxies and the HI brightness temperature field smoothed at relevant scales (the $k$NN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate HI data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with $k_\parallel < k_{\rm min,\parallel}$. Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the HI-galaxy cross-correlation signal in our simulations, using the $k$NN-field framework, is detectable at $>30\sigma$ across $r = [3,12] \, h^{-1}$Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp $k_\parallel$ filter can spuriously boost detection to $8\sigma$ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that $k$NN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.
Optical observations of Blue Compact Dwarf galaxies (BCDs) show they typically have high specific star formation rates and low metallicites. A subset of these galaxies (those with the lowest gas phase metallicities) display cometary optical morphologies similar to those found at high redshift. Whether this combination of properties predominantly arises from interactions with neighbours or accretion from the cosmic web or something else remains unclear. We used high resolution HI mapping to gain insights into the processes driving the observed properties of a sample of extremely metal poor (XMP) BCDs. We present Very Large Array B- and C-configuration HI mapping of four BCDs. For three of the targeted BCDs we also detect and map the HI in their nearby companions. In these three cases there is HI evidence for a recent flyby interaction between the BCD and a nearby companion. The HI evidence for recent interactions for these three BCDs is corroborated by our analysis of the tidal forces exerted on the BCDs by companions with available spectroscopic redshifts. For J0204-1009 we had sufficient spatial resolution to determine that it is dark matter dominated and estimate its DM halo mass to be in the range 1.2 x 10^11 to 5.2 x 10^11 solar masses. It is the most isolated BCD in our small sample, J0301-0052, which shows one of the most asymmetric HI morphologies. J0301-0052 has a similar cometary HI morphology to the BCD's optical morphology, although the HI column density maximum is projected at the end of the of the optical tail. Our HI observations suggest J0301-0052 may be undergoing a merger, while the other BCDs show evidence of a recent tidal interaction with a near neighbour. While our selection criteria favoured BCDs with companions our results are consistent with the earlier finding by other authors that most BCDs are associated with either mild tidal interactions or mergers.
The origin of Galactic Cosmic Rays (GCRs) and the potential role of Supernova Remnants (SNRs) as cosmic-ray (CR) accelerators remain subjects of ongoing debate. To shed more light on this topic, we have studied the spectral shape of two SNRs, RX J1713.7-3946 and HAWC J2227+610, performing simulations for the Cherenkov Telescope Array Observatory (CTAO). The previous multi-wavelength (MWL) analysis on these two sources showed an important hadronic contribution at high energies. The interaction of the GCRs accelerated by the SNRs with the medium around the accelerator leads to a process of pion decay (PD) that produces gamma-rays ($\gamma$-rays). These emissions, detectable by CTAO, offer an indirect means of pinpointing the CR source. Two scenarios have been considered: the contribution of heavy CRs and different cut-off sharpnesses ($\beta$) of the particle spectra. The simulations were performed by using different CR composition distributions (protons, CNO, Fe) and different sharpness values ranging from $\beta$=0.5 to $\beta$=1.5. The results show that, in the cases studied here, CTAO will increase the sensitivity to the spectral shape of $\gamma$-rays. This allows us to distinguish protons from heavy CRs and obtain information on $\beta$ values and therefore on different acceleration scenarios.
We present the Variable Star PhasE Curve (VSPEC) Collection, a set of Python packages for simulating combined-light spectroscopic observations of 3-dimensional exoplanet atmospheres in the presence of stellar variability and inhomogeneity. VSPEC uses the Planetary Spectrum Generator's Global Emission Spectra (PSG/GlobES) application along with a custom-built multi-component time-variable stellar model based on a user-defined grid of stellar photosphere models to produce spectroscopic light curves of the planet-host system. VSPEC can be a useful tool for modeling observations of exoplanets in transiting geometries (primary transit, secondary eclipse) as well as orbital phase curve measurements, and is built in a modular and flexible configuration for easy adaptability to new stellar and planetary model inputs. We additionally present a set of codes developed alongside the core VSPEC modules, including the stellar surface model generator vspec-vsm, the stellar spectral grid interpolation code GridPolator, and a Python interface for PSG, libpypsg.
The origin of planetary mass objects (PMOs) wandering in young star clusters remains enigmatic, especially when they come in pairs. They could represent the lowest-mass object formed via molecular cloud collapse or high-mass planets ejected from their host stars. However, neither theory fully accounts for their abundance and multiplicity. Here, we show via hydrodynamic simulations that free-floating PMOs have a unique formation channel via the fragmentation of tidal bridge between encountering circumstellar disks. This process can be highly productive in density clusters like Trapezium forming metal-poor PMOs with disks. Free-floating multiple PMOs also naturally emerge when neighboring PMOs are caught by mutual gravity. PMOs may thus form a distinct population different from stars and planets.
Stellar streams retain a memory of their gravitational interactions with small-scale perturbations. While perturbative models for streams have been formulated in action-angle coordinates, a direct transformation to these coordinates is only available for static and typically axisymmetric models for the galaxy. The real Milky Way potential is in a state of disequilibrium, complicating the application of perturbative methods around an equilibrium system. Here, we utilize a combination of differentiable simulations and Hamiltonian perturbation theory to model the leading-order effect of dark matter subhalos on stream observables. To obtain a perturbative description of streams, we develop a direct and efficient forward mode differentiation of Hamilton's equations of motion. Our model operates in observable coordinates, allowing us to treat the effects of arbitrary subhalo potentials on streams perturbatively, while simultaneously capturing non-linear effects due to other substructures like the infalling LMC or the rotating bar. The model predicts the velocity dispersion of streams as a function of subhalo statistics, allowing us to constrain the low-mass range of subhalos down to $\sim 10^5~M_\odot$. We forecast the velocity dispersion of the GD-1 stream, and find that observations are in agreement with a CDM subhalo population, with a slight preference for more dense subhalos. The method provides a new approach to characterize streams in the presence of substructure, with significantly more modeling flexibility compared to previous works.
The kinematic properties of the Sco-Cen association have been studied using the spatial velocities of young stars. New kinematic age estimates for the three components of the association with the age of UCL and LCC being $17.7\pm2.4$ Myrs and the age of US being $6.4\pm1.7$ Myrs have been obtained. The parameters of the residual velocities US, UCL, and LCC ellipsoid have been estimated.
Non-circular (NC) motions represent the imprints of non-axisymmetric structures in galaxies, providing opportunities to study the physical properties of gas departing from circular rotation. In this work, we have conducted a systematic study of the non-circular motions in a sample of 1624 gas-rich disk galaxies from the MaNGA MPL-11. By using the H$\alpha$ velocity as a tracer of the disk rotation, we find indications that the amplitude of the non-circular motions is related to the stellar mass, with the low mass and late-type galaxies the most affected. In our sample, we find ratios of circular to non-circular rotation ranging from 5% to 20%. By implementing harmonic models to include NC motions associated with spiral arms and stellar bars, we find that the rotational curves traced with H$\alpha$ are barely affected by the NC induced by these structures. Consequently, in our sample, we do not find evidence that NC motions contribute to the scatter of the stellar Tully-Fisher relation. Our results suggest that non-circular motions might have a more localized effect in galaxies rather than a global one.
Demixing properties of planetary major constituents influence the interior structure and evolution of planets. Comparing experimental and computational data on the miscibility of hydrogen and water to adiabatic profiles suggests phase separation between these components occurs in the ice giants Uranus and Neptune. We aim to predict the atmospheric water abundance and transition pressure between the water-poor outer envelope and the water-rich deep interior in Uranus and Neptune. We construct seven H2-H2O phase diagrams from the available experimental and computational data. We compute interior adiabatic structure models and compare these to the phase diagrams to infer whether demixing is occurring. We obtain a strong water depletion in the top layer due to rain-out of water and find upper limits on the atmospheric water mass fraction Z_atm of 0.21 for Uranus and 0.16 for Neptune. The transition from the water-poor to the water-rich layer is sharp and occurs at pressures P_Z between 4 and 11 GPa. Using these constraints on Z_atm and P_Z, we find that the observed gravitational harmonics J2 and J4 can be reproduced if P_Z > 10 GPa in Uranus and > 5 GPa in Neptune, and if the deep interior has a high primordial water mass fraction of 0.8, unless rocks are also present. The agreement with J4 is improved if rocks are confined deeper than P_Z, for instance below a rock cloud level at 2000 K (20-30 GPa). These findings confirm classical few-layer models and suggest that a layered structure may result from a combination of primordial mass accretion and subsequent phase separation. Reduced observational uncertainty in J4 and its dynamic contribution, atmospheric water abundance measurements from an Orbiter with a Probe mission to Uranus (UOP) or Neptune, and better understanding of the mixing behaviour of constituents are needed to constrain the interiors of ice giants.
The advent of large astronomical surveys has made available large and complex data sets. However, the process of discovery and interpretation of each potentially new astronomical source is, many times, still handcrafted. In this context, machine learning algorithms have emerged as a powerful tool to mine large data sets and lower the burden on the domain expert. Active learning strategies are specially good in this task. In this report, we used the PineForest algorithm to search for superluminous supernova (SLSN) candidates in the Zwicky Transient Facility. We showcase how the use of previously confirmed sources can provide important information to boost the convergence of the active learning algorithm. Starting from a data set of $\sim$14 million objects, and using 8 previously confirmed SLSN light curves as priors, we scrutinized 120 candidates and found 8 SLSN candidates, 2 of which have not been reported before (AT 2018moa and AT 2018mob). These results demonstrate how existing spectroscopic samples can be used to improve the efficiency of active learning strategies in searching for rare astronomical sources.
We present an accelerated pipeline, based on high-performance computing techniques and normalizing flows, for joint Bayesian parameter estimation and model selection and demonstrate its efficiency in gravitational wave astrophysics. We integrate the Jim inference toolkit, a normalizing flow-enhanced Markov chain Monte Carlo (MCMC) sampler, with the learned harmonic mean estimator. Our Bayesian evidence estimates run on $1$ GPU are consistent with traditional nested sampling techniques run on $16$ CPU cores, while reducing the computation time by factors of $5\times$ and $15\times$ for $4$-dimensional and $11$-dimensional gravitational wave inference problems, respectively. Our code is available in well-tested and thoroughly documented open-source packages, ensuring accessibility and reproducibility for the wider research community.
We present a novel deep learning framework using Long Short-Term Memory (LSTM) networks to predict the spectra of galactic cosmic rays by leveraging historical solar activity data, addressing the limitations of traditional transport models. By incorporating multiple solar parameters, such as the heliospheric magnetic field, solar wind speed, and sunspot numbers, our model achieves accurate short-term and long-term predictions of cosmic-ray flux. The inclusion of historical cosmic-ray data enhances prediction accuracy, making the model highly effective for space weather forecasting. Moreover, it provides reliable one-day-ahead predictions of full spectra of cosmic rays for different species. Our approach surpasses traditional physics-based methods, providing a scalable, data-driven solution for reliable daily and long-term forecasts. This work paves the way for advanced models that can integrate broader observational data, with significant implications for space weather monitoring and mission planning.
Context. Pulsar wind nebulae (PWNe) are a source of very high energy (VHE) radiation that can reach up to TeV and even PeV energies. Our work uses the pulsar tree, a graph theory tool recently presented to analyze the pulsar population to select the candidates of interest. Aims. We aim to discover possible detectable PWNe. We also aim to test to what extent the pulsar tree groups detectable PWNe despite it considering only pulsars' intrinsic properties. Methods. We select four pulsars as candidates for TeV PWNe based on their positions in the pulsar tree. Using observed and assumed ranges of values for relevant parameters, we anticipate possible spectral energy distributions (SEDs) of the PWNe of four pulsars (PSR J1208-6238, J1341-6220, J1838-0537, and J1844-0346) via a detailed time-dependent, leptonic model that was already found appropriate to describe almost all other detected nebulae. Results. We estimate the likelihood of detection for the four candidates studied by comparing the TeV fluxes predicted by the possible models with the sensitivities of different observatories. In doing so, we provide context for analyzing the advantages and caveats of the pulsar tree position as a marker for properties that go beyond the intrinsic features of pulsars that are considered in producing it.
Human exploration of the moon is expected to resume in the next decade, following the last such activities in the Apollo programme time. One of the major objectives of returning to the Moon is to continue retrieving geological samples, with a focus on collecting high-quality specimens to maximize scientific return. Tools that assist astronauts in making informed decisions about sample collection activities can maximize the scientific value of future lunar missions. A lunar rock classifier is a tool that can potentially provide the necessary information for astronauts to analyze lunar rock samples, allowing them to augment in-situ value identification of samples. Towards demonstrating the value of such a tool, in this paper, we introduce a framework for classifying rock types in thin sections of lunar rocks. We leverage the vast collection of petrographic thin-section images from the Apollo missions, captured under plane-polarized light (PPL), cross-polarised light (XPL), and reflected light at varying magnifications. Advanced machine learning methods, including contrastive learning, are applied to analyze these images and extract meaningful features. The contrastive learning approach fine-tunes a pre-trained Inception-Resnet-v2 network with the SimCLR loss function. The fine-tuned Inception-Resnet-v2 network can then extract essential features effectively from the thin-section images of Apollo rocks. A simple binary classifier is trained using transfer learning from the fine-tuned Inception-ResNet-v2 to 98.44\% ($\pm$1.47) accuracy in separating breccias from basalts.
We assess the impact of QSOs on the high redshift (z > 4) Intergalactic Medium using Monte Carlo realisations of QSO populations and the HeIII regions they generate, applied to the Sherwood-Relics simulations, allowing for uncertainties in the QSO luminosity function, its evolution, and QSO spectra and ages. While QSO luminosity functions based on optical-infrared selection are unable to reproduce the broadening HI Lyman-alpha optical depth distributions at z > 5, much broader distributions are found for the higher numbers of QSOs based on x-ray selection, suggesting a large QSO contribution to the ultra-violet background at z > 5 may offer an alternative to late reionization models to account for the broad HI Lyman-alpha optical depth distributions. Realisations using QSOs based on the higher QSO counts also much better recover the measured pixel flux auto-correlation function at z > 5. The HeIII regions from QSO sources according to both types of luminosity function suppress the pixel flux power spectrum on small scales, k > 0.02 s/km, while enhancing it on larger, both by amounts of up to tens of percent at z > 4, with the magnitude increasing with the intergalactic HeIII filling factor and the boost in temperature within the HeIII regions.
Synchrotron maser emission is a leading candidate to explain the coherent emission from Fast Radio Bursts (FRBs). This mechanism requires a population inversion in order to operate. We show that nonresonant interactions between Alfvén waves and a relativistic plasma result in the formation of population inversions across a wide range of magnetizations, $\sigma\gtrsim10^{-4}$, and temperatures, $10^{-2} \leq k_bT/mc^2 \leq 3$, spanning the parameters expected in FRB environments. We calculate the fraction of energy contained in the inversion across the whole of this parameter space for the first time and we show that energy fractions of $f_{inv}\sim 0.3$ are achieved for high magnetizations $\sigma >1$. The population inversion forms on time-scales compatible with the typical dynamical time-scales of magnetars for all magnetizations. Furthermore, we provide physical explanations for the behaviour of the interaction in different magnetization regimes, and identify the important characteristic values at which this behaviour changes. We also show that the mechanism is capable of producing an FRB signal at GHz frequencies in a relativistic magnetar wind close to the light cylinder and that this signal can escape the magnetar environment without significant damping.
The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-trace tracer of externally driven photoevaporative winds, which agrees with recent excitation models.
This study continues our investigation of early-type binaries using high-precision broad-band polarimetry, focusing on HD 165052, a massive O+O-type binary in the young cluster NGC 6530. Our aim was to monitor linear polarization variations and independently determine the orbital period through polarization data analysis. By examining the phase-locked Stokes parameters, we estimated the orbital inclination at about 55+5/-55 deg, and the orientation at 148+20/-22 deg. The binary's rotation direction is clockwise. Using the Dipol-2 polarimeter with the 60 cm KVA and Tohoku T60 telescopes, we achieved polarization measurement accuracy around 0.01 percent in the B, V, and R passbands. Period analysis identified a strong periodic signal with a cycle of about 1.48 days, corresponding to half of the orbital period of about 2.96 days. The second Fourier harmonics suggest a symmetric scattering geometry, likely due to electron scattering in interacting stellar winds. Our estimates indicate a mass-loss rate of approximately 4.0 x 10^(-7) solar masses per year. Observations of neighboring stars in NGC 6530 revealed complex interstellar polarization behaviors within the cluster, enhancing our understanding of its polarimetric environment.
Tidal forces in close binary systems have diverse impacts on magnetic activity. The synchronicity characteristic of close systems counteracts magnetic braking, thereby sustaining rapid rotation-a key factor in increased levels of magnetic activity. Tidal effects can also work against the slowing down of rotation during stellar evolution, when the star inflates into a red giant. A notable manifestation of the effect of binarity on activity in such systems is the appearance of active longitudes, which are thought to arise from the excitation of non-axisymmetric dynamo modes. Through some recent examples, the dynamo operation in RS CVn and BY Dra type systems is briefly reviewed in terms of spot cycles, active longitudes, flare activity, and differential rotation.
The increasing congestion in the near-Earth space environment has amplified the need for robust and efficient conjunction analysis techniques including the computation of the minimum distance between orbital paths in the presence of perturbations. After showing that classical Minimum Orbit Intersection Distance (MOID) computation schemes are unsuitable to treat Earth orbiting objects, the article presents an analytical approach to provide a more accurate estimate of the true distance between perturbed trajectories by incorporating the effect of zonal harmonics of arbitrary order. Cook's linear secular theory for the motion of the eccentricity vector is extended to include higher order eccentricity effects and applied to the computation of the minimum and maximum radii attained by two orbits at their mutual nodes, which can be employed to estimate the true distance between the two orbital paths and to establish an efficient algorithm for determining or excluding potential conjunctions. Extensive testing and validation are conducted using a high-fidelity propagator and a comprehensive dataset of resident space objects. The results demonstrate an accuracy below the km level for the orbit distance computation in 99\% of cases, which enables high-efficiency conjunction filtering.
Relativistic magnetic reconnection is one of the most fundamental mechanisms considered responsible for the acceleration of relativistic particles in astrophysical jets and magnetospheres of compact objects. Understanding the properties of the dissipation of magnetic fields and the formation of non-ideal electric fields is of paramount importance to quantify the efficiency of reconnection at energizing charged particles. Recent results from particle-in-cell (PIC) simulations suggest that the fundamental properties of how magnetic fields dissipate in a current sheet might be captured by an ``effective resistivity'' formulation, which would locally enhance the amount of magnetic energy dissipated and favor the onset of fast reconnection. Our goal is to assess this ansatz quantitatively by comparing fluid models of magnetic reconnection with a non-constant magnetic diffusivity and fully-kinetic models. We perform 2D resistive relativistic magnetohydrodynamic (ResRMHD) simulations of magnetic reconnection combined to PIC simulations using the same initial conditions (namely a Harris current sheet). We explore the impact of crucial parameters such as the plasma magnetization, its mass density, the grid resolution, and the characteristic plasma skin depth. Our ResRMHD models with effective resistivity can quantitatively reproduce the dynamics of fully-kinetic models of relativistic magnetic reconnection. In particular, they lead to reconnection rates consistent with PIC simulations, while for constant-resistivity fluid models the reconnection dynamics is generally 10 times slower. Even at modest resolutions the adoption of an effective resistivity can qualitatively capture the properties of kinetic reconnection models and produce reconnection rates compatible with collisionless models, i.e. of the order of $\sim10^{-1}$.
The redshift and size distributions of galaxy scale strong lenses depend on the evolution of early-type galaxies (ETGs) in the redshift range 0.2<z<1. We use this dependence to constrain the velocity dispersion function (VDF) evolution from the Strong Lensing Legacy Survey (SL2S) sample of lenses. Our modeling of the lens population includes lens identifiability given survey parameters, and constrains the evolution of the VDF based on the redshift distributions of sources and lenses as well as the distribution of Einstein radii. We consider five different assumptions for the reference VDF at redshift zero and two sets of scaling relations for the VDF. We find that in all cases the observed lens sample favors a slow evolution of both the VDF normalization factor and the VDF characteristic velocity with redshift which is consistent with a VDF that is constant in redshift for z<1.