Search

We report the first measurements of deep inelastic scattering spin-dependent azimuthal asymmetries in back-to-back dihadron electroproduction, where two hadrons are produced in opposite hemispheres along the z-axis in the center-of-mass frame, with the first hadron produced in the current-fragmentation region and the second in the target-fragmentation region. The data were taken with longitudinally polarized electron beams of 10.2 and 10.6 GeV incident on an unpolarized liquid-hydrogen target using the CLAS12 spectrometer at Jefferson Lab. Observed non-zero $\sin\Delta\phi$ modulations in $ep \rightarrow e'p\pi^+X$ events, where $\Delta\phi$ is the difference of the azimuthal angles of the proton and pion in the virtual photon and target nucleon center-of-mass frame, indicate that correlations between the spin and transverse momenta of hadrons produced in the target- and current-fragmentation regions may be significant. The measured beam-spin asymmetries provide a first access in dihadron production to a previously unobserved leading-twist spin- and transverse-momentum-dependent fracture function. The fracture functions describe the hadronization of the target remnant after the hard scattering of a virtual photon off a quark in the target particle and provide a new avenue for studying nucleonic structure and hadronization., Comment: 7 pages, 6 figures

Several factors can contribute to the difficulty of aligning the sensors of tracking detectors, including a large number of modules, multiple types of detector technologies, and non-linear strip patterns on the sensors. All three of these factors apply to the CLAS12 CVT, which is a hybrid detector consisting of planar silicon sensors with non-parallel strips, and cylindrical micromegas sensors with longitudinal and arc-shaped strips located within a 5~T superconducting solenoid. To align this detector, we used the Kalman Alignment Algorithm, which accounts for correlations between the alignment parameters without requiring the time-consuming inversion of large matrices. This is the first time that this algorithm has been adapted for use with hybrid technologies, non-parallel strips, and curved sensors. We present the results for the first alignment of the CLAS12 CVT using straight tracks from cosmic rays and from a target with the magnetic field turned off. After running this procedure, we achieved alignment at the level of 10~$\mu$m, and the widths of the residual spectra were greatly reduced. These results attest to the flexibility of this algorithm and its applicability to future use in the CLAS12 CVT and other hybrid or curved trackers, such as those proposed for the future Electron-Ion Collider.

The study of reaction-diffusion systems on networks is of paramount relevance for the understanding of nonlinear processes in systems where the topology is intrinsically discrete, such as the brain. Until now reaction-diffusion systems have been studied only when species are defined on the nodes of a network. However, in a number of real systems including, e.g., the brain and the climate, dynamical variables are not only defined on nodes but also on links, faces and higher-dimensional cells of simplicial or cell complexes, leading to topological signals. In this work we study reaction-diffusion processes of topological signals coupled through the Dirac operator. The Dirac operator allows topological signals of different dimension to interact or cross-diffuse as it projects the topological signals defined on simplices or cells of a given dimension to simplices or cells of one dimension up or one dimension down. By focusing on the framework involving nodes and links we establish the conditions for the emergence of Turing patterns and we show that the latter are never localized only on nodes or only on links of the network. Moreover when the topological signals display Turing pattern their projection does as well. We validate the theory hereby developed on a benchmark network model and on square lattices with periodic boundary conditions.

We present the first measurement of di-hadron angular correlations in electron-nucleus scattering. The data were taken with the CLAS detector and a 5.0 GeV electron beam incident on deuterium, carbon, iron, and lead targets. Relative to deuterium, the nuclear yields of charged-pion pairs show a strong suppression for azimuthally opposite pairs, no suppression for azimuthally nearby pairs, and an enhancement of pairs with large invariant mass. These effects grow with increased nuclear size. The data are qualitatively described by the GiBUU model, which suggests that hadrons form near the nuclear surface and undergo multiple-scattering in nuclei. These results show that angular correlation studies can open a new way to elucidate how hadrons form and interact inside nuclei, Comment: 6 pages, 4 figures

Hypergraphs and simplical complexes both capture the higher-order interactions of complex systems, ranging from higher-order collaboration networks to brain networks. One open problem in the field is what should drive the choice of the adopted mathematical framework to describe higher-order networks starting from data of higher-order interactions. Unweighted simplicial complexes typically involve a loss of information of the data, though having the benefit to capture the higher-order topology of the data. In this work we show that weighted simplicial complexes allow to circumvent all the limitations of unweighted simplicial complexes to represent higher-order interactions. In particular, weighted simplicial complexes can represent higher-order networks without loss of information, allowing at the same time to capture the weighted topology of the data. The higher-order topology is probed by studying the spectral properties of suitably defined weighted Hodge Laplacians displaying a normalized spectrum. The higher-order spectrum of (weighted) normalized Hodge Laplacians is here studied combining cohomology theory with information theory. In the proposed framework, we quantify and compare the information content of higher-order spectra of different dimension using higher-order spectral entropies and spectral relative entropies. The proposed methodology is tested on real higher-order collaboration networks and on the weighted version of the simplicial complex model "Network Geometry with Flavor"., Comment: 20 pages, 11 figures

This chapter provides a comprehensive and self-contained discussion of the most recent developments of information theory of networks. Maximum entropy models of networks are the least biased ensembles enforcing a set of constraints and are used in a number of application to produce null model of networks. Here maximum entropy ensembles of networks are introduced by distinguishing between microcanonical and canonical ensembles revealing the the non-equivalence of these two classes of ensembles in the case in which an extensive number of constraints is imposed. It is very common that network data includes also meta-data describing feature of the nodes such as their position in a real or in an abstract space. The features of the nodes can be treated as latent variables that determine the cost associated to each link. Maximum entropy network ensembles with latent variables include spatial networks and their generalisation. In this chapter we cover the case of transportation networks including airport and rail networks. Maximum entropy network ensemble satisfy a given set of constraints. However traditional approaches do not provide any insight on the origin of such constraints. We use information theory principles to find the optimal distribution of latent variables in the framework of the classical information theory of networks. This theory explains the "blessing of non-uniformity" of data guaranteeing the efficiency of inference algorithms., Comment: 51 pages, 9 figures

Multiplex networks describe systems whose interactions can be of different nature, and are fundamental to understand complexity of networks beyond the framework of simple graphs. Recently it has been pointed out that restricting the attention to pairwise interactions is also a limitation, as the vast majority of complex systems include higher-order interactions that strongly affect their dynamics. Here, we propose hyper-diffusion on multiplex networks, a dynamical process in which diffusion on each single layer is coupled with the diffusion in other layers thanks to the presence of higher-order interactions occurring when there exists link overlap. We show that hyper-diffusion on a duplex network can be described by the Hyper-Laplacian in which the strength of four-body interactions among every set of four replica nodes connected by a multilink $(1,1)$ can be tuned by a parameter $\delta_{11}\ge 0$. The Hyper-Laplacian reduces to the standard lower Laplacian, capturing pairwise interactions at the two layers, when $\delta_{11}=0$. By combining tools of spectral graph theory, applied topology and network science we provide a general understanding of hyper-diffusion on duplex networks when $\delta_{11}>0$, including theoretical bounds on the Fiedler and the largest eigenvalue of Hyper-Laplacians and the asymptotic expansion of their spectrum for $\delta_{11}\ll1$ and $\delta_{11}\gg1$. Although hyper-diffusion on multiplex networks does not imply a direct ``transfer of mass" among the layers (i.e. the average state of replica nodes in each layer is a conserved quantity of the dynamics), we find that the dynamics of the two layers is coupled as the relaxation to the steady state becomes synchronous when higher-order interactions are taken into account and the Fiedler eigenvalue of the Hyper-Laplacian is not localized in a single layer of the duplex network., Comment: 27 pages, 5 figures, 8 appendices

Cross-referencing a watchlist of galaxy groups and clusters with transient detections from real-time streams of wide-field survey data is a promising method for discovering gravitationally lensed explosive transients including supernovae, kilonovae, gravitational waves and gamma-ray bursts in the next ten years. However, currently there exists no catalogue of objects with both sufficient angular extent and depth to adequately perform such a search. In this study, we develop a cluster-finding method capable of creating an all-sky list of galaxy group- and cluster-scale objects out to $z\simeq1$ based on their lens-plane properties and using only existing data from wide-field infrared surveys such as VHS and UHS, and all-sky \textit{WISE} data. In testing this method, we recover 91 per cent of a sample containing known and candidate lensing objects with Einstein radii of $\theta_E \geq 5\arcsec$. We also search the surrounding regions of this test sample for other groups and clusters using our method and verify the existence of any significant findings by visual inspection, deriving estimates of the false positive rate that are as low as 6 per cent. The method is also tested on simulated Rubin data from their DP0 programme, which yields complementary results of a good recovery rate of $\gtrsim 80$ per cent for $M_{200}\geq7\times10^{13}$M$_\odot$ clusters and with no false positives produced in our test region. Importantly, our method is positioned to create a watchlist in advance of Rubin's LSST, as it utilises only existing data, therefore enabling the discovery of lensed transients early within the survey's lifetime., Comment: 11 pages, 9 figures. Published, MNRAS

We search for gravitational wave (GW) events from LIGO-Virgo's third run that may have been affected by gravitational lensing. Gravitational lensing delays the arrival of GWs, and alters their amplitude -- thus biasing the inferred progenitor masses. This would provide a physically well-understood interpretation of GW detections in the ''mass gap'' between neutron stars and black holes, as gravitationally lensed binary neutron star (BNS) mergers. We selected three GW detections in LIGO-Virgo's third run for which the probability of at least one of the constituent compact objects being in the mass gap was reported as high with low latency -- i.e. candidate lensed BNS mergers. Our observations of powerful strong lensing clusters located adjacent to the peak of their sky localisation error maps reached a sensitivity $\rm AB\simeq25.5$ in the $z'$-band with the GMOS instruments on the Gemini telescopes, and detected no candidate lensed optical counterparts. We combine recent kilonova lightcurve models with recent predictions of the lensed BNS population and the properties of the objects that we followed up to show that realistic optical counterparts were detectable in our observations. Further detailed analysis of two of the candidates suggests that they are a plausible pair of images of the same low-mass binary black hole merger, lensed by a local galaxy or small group of galaxies. This further underlines that access to accurate mass information with low latency would improve the efficiency of candidate lensed BNS selection., Comment: Accepted for publication on MNRAS

Secure confirmation that a gravitational wave (GW) has been gravitationally lensed would bring together these two pillars of General Relativity for the first time. This breakthrough is challenging for many reasons, including: GW sky localization uncertainties dwarf the angular scale of gravitational lensing, the mass and structure of gravitational lenses is diverse, the mass function of stellar remnant compact objects is not yet well constrained, and GW detectors do not operate continuously. We introduce a new approach that is agnostic to the mass and structure of the lenses, compare the efficiency of different methods for lensed GW discovery, and explore detection of lensed kilonova counterparts as a direct method for localising candidates. Our main conclusions are: (1) lensed neutron star mergers (NS-NS) are magnified into the "mass gap" between NS and black holes, therefore selecting candidates from public GW alerts with high mass gap probability is efficient, (2) the rate of detectable lensed NS-NS will approach one per year in the mid-2020s, (3) the arrival time difference between lensed NS-NS images is $1\,\rm sec\lesssim\Delta t\lesssim1\,year$, and thus well-matched to the operations of GW detectors and optical telescopes, (4) lensed kilonova counterparts are faint at peak (e.g.\ $r_{\rm AB}\simeq24-26$ in the mid-2020s), fade quickly ($d<2\,\rm days$), and are detectable with target of opportunity observations with large wide-field telescopes. For example, just $\lesssim0.25$ per cent of Vera C.\ Rubin Observatory's observing time will be sufficient to follow up one well-localized candidate per year. Our predictions also provide a physically well-defined basis for exploring electromagnetically the exciting new "mass gap" discovery space., Comment: Accepted by MNRAS. 20 pages, 10 figures

Percolation establishes the connectivity of complex networks and is one of the most fundamental critical phenomena for the study of complex systems. On simple networks, percolation displays a second-order phase transition; on multiplex networks, the percolation transition can become discontinuous. However, little is known about percolation in networks with higher-order interactions. Here, we show that percolation can be turned into a fully-fledged dynamical process when higher-order interactions are taken into account. By introducing signed triadic interactions, in which a node can regulate the interactions between two other nodes, we define triadic percolation. We uncover that in this paradigmatic model the connectivity of the network changes in time and that the order parameter undergoes a period-doubling and a route to chaos. We provide a general theory for triadic percolation which accurately predicts the full phase diagram on random graphs as confirmed by extensive numerical simulations. We find that triadic percolation on real network topologies reveals a similar phenomenology. These results radically change our understanding of percolation and may be used to study complex systems in which the functional connectivity is changing in time dynamically and in a non-trivial way, such as in neural and climate networks., Comment: 62 pages, 20 figures

Maximum entropy network ensembles have been very successful in modelling sparse network topologies and in solving challenging inference problems. However the sparse maximum entropy network models proposed so far have fixed number of nodes and are typically not exchangeable. Here we consider hierarchical models for exchangeable networks in the sparse limit, i.e. with the total number of links scaling linearly with the total number of nodes. The approach is grand canonical, i.e. the number of nodes of the network is not fixed a priori: it is finite but can be arbitrarily large. In this way the grand canonical network ensembles circumvent the difficulties in treating infinite sparse exchangeable networks which according to the Aldous-Hoover theorem must vanish. The approach can treat networks with given degree distribution or networks with given distribution of latent variables. When only a subgraph induced by a subset of nodes is known, this model allows a Bayesian estimation of the network size and the degree sequence (or the sequence of latent variables) of the entire network which can be used for network reconstruction., Comment: (30 pages, 4 figures)

Data on hyperon form factors recently obtained from annihilation reactions are compared each other, and to proton and neutron form factors, in terms of two kinematical variables: the transferred momentum square, $q^2$, and the modulus $P$ of the relative momentum between the outgoing baryons. They are critically discussed in terms of possible correlated structures. The present status of the time-like form factor data for baryons is described and suggestions are given on the reactions and the kinematical range where data are desirable in order to clarify the arisen questions., Comment: 20 pages, 10 figures

New results for the exclusive and quasi-free cross sections off neutrons bound in deuterium ${\gamma}_vn(p) \rightarrow p{\pi}^{-} (p)$ are presented over a wide final state hadron angle range with a kinematic coverage of the invariant mass ($W$) up to 1.825 GeV and the virtual photon four-momentum transfer squared ($Q^{2}$) from 0.4 to 1.0 GeV$^2$. The exclusive structure functions were extracted and their Legendre moments were obtained. Final-state-interaction contributions have been kinematically separated from the extracted quasi-free cross sections off bound neutrons solely based on the analysis of the experimental data. These new results will serve as long-awaited input for phenomenological analyses to extract the $Q^{2}$ evolution of previously unavailable $n \to N^{*}$ electroexcitation amplitudes and to improve state-of-the-art models of neutrino scattering off nuclei by augmenting the already available results from free protons., Comment: The author list has been updated

While stripe phases with broken rotational symmetry of charge density appear in many complex correlated systems, the heterogeneity of spatial ordering and dynamics remains elusive. This missing info is at the heart of understanding the structure and function relation in quantum complex materials. We focus here on the spatial heterogeneity of the motion of charge density wave (CDW) at nanoscale in the archetypal case of La$_{2-x}$Sr$_x$NiO$_{4+y}$ perovskite at low temperature. We report compelling evidence that the unconventional increasing motion of CDW at T < 50K is related with the decreasing of its correlation length using resonant soft X-ray photon correlation spectroscopy (XPCS). The key result of this work is the direct visualization of nanoscale spatial inhomogeneity of CDW relaxation dynamics by scanning micro X-ray diffraction (SmXRD) showing a nanoscale landscape of percolating short range dynamic CDW puddles competing with large quasi-static CDW puddles giving rise to a novel form of nanoscale phase separation of the incommensurate stripes order landscape., Comment: 15 pages, 3 figures

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2)., Comment: 29 pages, 25 figures; accepted for publication in MNRAS following minor corrections

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE's 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters., Comment: 20 pages, 12 figures, 3 tables (including 1 in an appendix). Accepted in MNRAS. Data Release products available at https://astro.dur.ac.uk/~hbpn39/pilot-wings.html

Beam-recoil transferred polarizations for the exclusive electroproduction of $K^+\Lambda$ and $K^+\Sigma^0$ final states from an unpolarized proton target have been measured using the CLAS12 spectrometer at Jefferson Laboratory. The measurements at beam energies of 6.535~GeV and 7.546~GeV span the range of four-momentum transfer $Q^2$ from 0.3 to 4.5~GeV$^2$ and invariant energy $W$ from 1.6 to 2.4~GeV, while covering the full center-of-mass angular range of the $K^+$. These new data extend the existing hyperon polarization data from CLAS in a similar kinematic range but from a significantly larger dataset. They represent an important addition to the world data, allowing for better exploration of the reaction mechanism in strangeness production processes, for further understanding of the spectrum and structure of excited nucleon states, and for improved insight into the strong interaction in the regime of non-perturbative dynamics., Comment: 25 pages, 22 figures, 4 tables

While it is known that the double-layer Bi$_2$Sr$_2$CaCu$_2$O$_{8+y}$ (BSCCO) cuprate superconductor exhibits a one-dimensional (1D) incommensurate superlattice (IS), the effect of IS on the electronic structure remains elusive. Following the recent shift of the interest from underdoped to optimum and overdoped phase in BSCCO by increasing the hole doping $x$, controlled by the oxygen interstitials concentration $y$, here we focus on the multiple splitting of the density of states (DOS) peaks and emergence of higher order van Hove singularities (VHS) due to the 1D incommensurate superlattice. It is known that 1D incommensurate wave vector $\textbf{q} = \epsilon \textbf{b}$ (where $\textbf{b}$ is the reciprocal lattice vector of the orthorhombic lattice) is controlled by the misfit strain between different atomic layers in the range $0.209$ - $0.215$ in BSCCO and in the range $0.209$ - $0.25$ in Bi$_2$Sr$_2$Ca$_{1-x}$Y$_x$Cu$_2$O$_{8+y}$ (BSCYCO). This work reports the theoretical calculation of a complex pattern of VHS due to the 1D incommensurate superlattice with large 1D quasi-commensurate supercells with the wave vector $\epsilon =9/\eta$ in the range $36 > \eta > 43$. The similarity of the complex VHS splitting and appearing of higher order VHS in a mismatched CuO$_2$ bilayer with VHS due to the moir{\'e} lattice in strained twisted bilayer graphene is discussed. This makes mismatched CuO$_2$ bilayer quite promising for constructing quantum devices with tuned physical characteristics., Comment: 11 pages, 7 figures

During the COVID pandemic, periods of exponential growth of the disease have been mitigated by containment measures that in different occasions have resulted in a power-law growth of the number of cases. The first observation of such behaviour has been obtained from 2020 late spring data coming from China by Ziff and Ziff in Ref. [1]. After this important observation the power-law scaling (albeit with different exponents) has also been observed in other countries during periods of containment of the spread. Early interpretations of these results suggest that this phenomenon might be due to spatial effects of the spread. Here we show that temporal modulations of infectivity of individuals due to containment measures can also cause power-law growth of the number of cases over time. To this end we propose a stochastic well-mixed Susceptible-Infected-Removed (SIR) model of epidemic spreading in presence of containment measures resulting in time-dependent infectivity and we explore the statistical properties of the resulting branching process at criticality. We show that at criticality it is possible to observe power-law growth of the number of cases with exponents ranging between one and two. Our asymptotic analytical results are confirmed by extensive Monte Carlo simulations. Although these results do not exclude that spatial effects might be important in modulating the power-law growth of the number of cases at criticality, this work shows that even well-mixed populations may already feature non-trivial power-law exponents at criticality., Comment: 24 pages, 6 figures

While understanding of periodic recurrent waves of Covid-19 epidemics would aid to combat the pandemics, quantitative analysis of data over a two years period from the outbreak, is lacking. The complexity of Covid-19 recurrent waves is related with the concurrent role of i) the containment measures enforced to mitigate the epidemics spreading ii) the rate of viral gene mutations, and iii) the variable immune response of the host implemented by vaccination. This work focuses on the effect of massive vaccination and gene variants on the recurrent waves in a representative case of countries enforcing mitigation and vaccination strategy. The spreading rate is measured by the ratio between the reproductive number Rt(t) and the doubling time Td(t) called RIC-index and the daily fatalities number. The dynamics of the Covid-19 epidemics has been studied by wavelet analysis and represented by a non-linear helicoid vortex in a 3D space where both RIC-index and fatalities change with time. The onset of periodic recurrent waves has been identified by the transition from convergent to divergent trajectories on the helicoid vortex. We report a main period of recurrent waves of 120 days and the elongation of this period after the vaccination campaign, Comment: 15 pages, 3 figures

The first results on the $\sigma_{LT'}$ structure function in exclusive $\pi^0p$ electroproduction at invariant masses of the final state of 1.5 GeV $<$ $W$ $<$ 1.8 GeV and in the range of photon virtualities 0.4 GeV$^2 < Q^2 < 1.0$ GeV$^2$ were obtained from data on beam spin asymmetries and differential cross sections measured with the CLAS detector at Jefferson Lab. The Legendre moments determined from the $\sigma_{LT'}$ structure function have demonstrated sensitivity to the contributions from the nucleon resonances in the second and third resonance regions. These new data on the beam spin asymmetries in $\pi^0p$ electroproduction extend the opportunities for the extraction of the nucleon resonance electroexcitation amplitudes in the mass range above 1.6 GeV., Comment: 7 pages, 3 figures

Nanoscale phase separation (NPS) characterized by particular types of correlated disorder plays an important role in the functionality of high temperature superconductors (HTS). Our results show that multiscale heterogeneity is an essential ingredient of quantum functionality in complex materials. Here, the interactions developing between different structural units cause dynamical spatiotemporal conformations with correlated disorder. Thus visualizing conformations landscape is needed for understanding the physical properties of complex matter and requires advanced methodologies based on high precision X ray measurements. We discuss the dynamical correlated disorder at nanoscale and the related functionality in oxygen doped perovskite superconducting materials, Comment: 11 pages, 5 figures

Exchangeability is a desired statistical property of network ensembles requiring their invariance upon relabelling of the nodes. However combining sparsity of network ensembles with exchangeability is challenging. Here we propose a statistical physics framework and a Metropolis-Hastings algorithm defining exchangeable sparse network ensembles. The model generates networks with heterogeneous degree distributions by enforcing only global constraints while existing (non exchangeable) exponential random graphs enforce an extensive number of local constraints. This very general theoretical framework to describe exchangeable networks is here first formulated for uncorrelated simple networks and then it is extended to treat simple networks with degree correlations, directed networks, bipartite networks and generalized network structures including multiplex networks and simplicial complexes. In particular here we formulate and treat both uncorrelated and correlated exchangeable ensembles of simplicial complexes using statistical mechanics approaches., Comment: (25 pages, 1 figure)

Complex networks have become the main paradigm for modelling the dynamics of interacting systems. However, networks are intrinsically limited to describing pairwise interactions, whereas real-world systems are often characterized by higher-order interactions involving groups of three or more units. Higher-order structures, such as hypergraphs and simplicial complexes, are therefore a better tool to map the real organization of many social, biological and man-made systems. Here, we highlight recent evidence of collective behaviours induced by higher-order interactions, and we outline three key challenges for the physics of higher-order systems., Comment: pre-peer-reviewed version of the Nature Physics perspective, 7 pages, 4 figures

Background: Energetic quarks in nuclear DIS propagate through the nuclear medium. Processes that are believed to occur inside nuclei include quark energy loss through medium-stimulated gluon bremsstrahlung and intra-nuclear interactions of forming hadrons. More data are required to gain a more complete understanding of these effects. Purpose: To test the theoretical models of parton transport and hadron formation, we compared their predictions for the nuclear and kinematic dependence of pion production in nuclei. Methods: We have measured charged-pion production in semi-inclusive DIS off D, C, Fe, and Pb using the CLAS detector and the CEBAF 5.014 GeV electron beam. We report results on the nuclear-to-deuterium multiplicity ratio for $\pi^{+}$ and $\pi^{-}$ as a function of energy transfer, four-momentum transfer, and pion energy fraction or transverse momentum - the first three-dimensional study of its kind. Results: The $\pi^{+}$ multiplicity ratio is found to depend strongly on the pion fractional energy $z$, and reaches minimum values of $0.67\pm0.03$, $0.43\pm0.02$, and $0.27\pm0.01$ for the C, Fe, and Pb targets, respectively. The $z$ dependences of the multiplicity ratios for $\pi^{+}$ and $\pi^{-}$ are equal within uncertainties for C and Fe targets but show differences at the level of 10$\%$ for the Pb-target data. The results are qualitatively described by the GiBUU transport model, as well as with a model based on hadron absorption, but are in tension with calculations based on nuclear fragmentation functions. Conclusions: These precise results will strongly constrain the kinematic and flavor dependence of nuclear effects in hadron production, probing an unexplored kinematic region. They will help to reveal how the nucleus reacts to a fast quark, thereby shedding light on its color structure, transport properties, and on the mechanisms of the hadronization process., Comment: 22 pages, 5 figures, 12 tables

We present a Spitzer/IRAC survey of H-faint ($H_{160} \gtrsim 26.4$, $<5\sigma$) sources in 101 lensing cluster fields. Across a CANDELS/Wide-like survey area of $\sim$648 arcmin$^2$ (effectively $\sim$221 arcmin$^2$ in the source plane), we have securely discovered 53 sources in the IRAC Channel-2 band (CH2, 4.5 $\mathrm{\mu m}$; median CH2$=22.46\pm0.11$ AB mag) that lack robust HST/WFC3-IR F160W counterparts. The most remarkable source in our sample, namely ES-009 in the field of Abell 2813, is the brightest H-faint galaxy at 4.5 $\mathrm{\mu m}$ known so far ($\mathrm{CH2}=20.48\pm0.03$ AB mag). We show that the H-faint sources in our sample are massive (median $M_\mathrm{star} = 10^{10.3\pm 0.3}$ $M_\odot$), star-forming (median star formation rate $=100_{-40}^{+60}$ $M_\odot$yr$^{-1}$) and dust-obscured ($A_V=2.6\pm0.3$) galaxies around a median photometric redshift of $z=3.9\pm0.4$. The stellar continua of 14 H-faint galaxies can be resolved in the CH2 band, suggesting a median circularized effective radius ($R_\mathrm{e,circ}$; lensing corrected) of $1.9\pm0.2$ kpc and $<1.5$ kpc for the resolved and whole samples, respectively. This is consistent with the sizes of massive unobscured galaxies at $z\sim4$, indicating that H-faint galaxies represent the dusty tail of the distribution of a wider galaxy population. Comparing with the ALMA dust continuum sizes of similar galaxies reported previously, we conclude that the heavy dust obscuration in H-faint galaxies is related to the compactness of both stellar and dust continua ($R_\mathrm{e,circ}\sim 1$ kpc). These H-faint galaxies make up $16_{-7}^{+13}$% of the galaxies in the stellar mass range of $10^{10}-10^{11.2}$ $M_\odot$ at $z=3\sim5$, contributing to $8_{-4}^{+8}$% of the cosmic star formation rate density in this epoch and likely tracing the early phase of massive galaxy formation., Comment: 26 pages, 10 figures, 3 tables. Accepted for publication in ApJ