-
The COSINE-100 dark matter search experiment is an array of NaI(Tl) crystal
detectors located in the Yangyang Underground Laboratory (Y2L). To understand
measured backgrounds in the NaI(Tl) crystals we have performed Monte Carlo
simulations using the Geant4 toolkit and developed background models for each
crystal that consider contributions from both internal and external sources,
including cosmogenic nuclides. The background models are based on comparisons
of measurement data with Monte Carlo simulations that are guided by a campaign
of material assays and are used to evaluate backgrounds and identify their
sources. The average background level for the six crystals (70 kg total mass)
that are studied is 3.5 counts/day/keV/kg in the (2-6) keV energy interval. The
dominant contributors in this energy region are found to be $^{210}$Pb and
$^3$H.
-
This paper reports measurements of final-state proton multiplicity, muon and
proton kinematics, and their correlations in charged-current pionless neutrino
interactions, measured by the T2K ND280 near detector in its plastic
scintillator (C$_8$H$_8$) target. The data were taken between years 2010 and
2013, corresponding to approximately 6$\times10^{20}$ protons on target. Thanks
to their exploration of the proton kinematics and of kinematic imbalances
between the proton and muon kinematics, the results offer a novel probe of the
nuclear-medium effects most pertinent to the (sub-)GeV neutrino-nucleus
interactions that are used in accelerator-based long-baseline neutrino
oscillation measurements. These results are compared to many neutrino-nucleus
interaction models which all fail to describe at least part of the observed
phase space. In case of events without a proton above a detection threshold in
the final state, a fully consistent implementation of the local Fermi gas model
with multinucleon interactions gives the best description of the data. In the
case of at least one proton in the final state the spectral function model
agrees well with the data, most notably when measuring the kinematic imbalance
between the muon and the proton in the plane transverse to the incoming
neutrino. A clear indication of existence of multinucleon interactions is
observed. The effect of final-state interactions is also discussed.
-
This paper reports the first differential measurement of the charged-current
interaction cross section of $\nu_{\mu}$ on water with no pions in the final
state. This flux-averaged measurement has been made using the T2K experiment's
off-axis near detector, and is reported in doubly-differential bins of muon
momentum and angle. The flux-averaged total cross section in a restricted
region of phase space was found to be $ \sigma= (0.95 \pm 0.08 (\mbox{stat})
\pm 0.06 (\mbox{det. syst.}) \pm 0.04(\mbox{model syst.}) \pm{}
0.08(\mbox{flux}) ) \times 10^{-38} \mbox{cm}^2$ per n.
-
We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628
(M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II).
Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from
the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with
archival HST/ACS data, we correct thousands of individual stars for extinction
across these five galaxies using an isochrone-matching (reddening-free Q)
method. We generate extinction maps for each galaxy from the individual stellar
extinctions using both adaptive and fixed resolution techniques, and correlate
these maps with neutral HI and CO gas maps from literature, including The HI
Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey
(HERACLES). We calculate dust-to-gas ratios and investigate variations in the
dust-to-gas ratio with galaxy metallicity. We find a power law relationship
between dust-to-gas ratio and metallicity, consistent with other studies of
dust-to-gas ratio compared to metallicity. We find a change in the relation
when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in
low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$
could have produced too low a slope in the derived relationship between
dust-to-gas ratio and metallicity. We also compare our extinctions to those
derived from fitting the spectral energy distribution (SED) using the Bayesian
Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower
extinctions from SED-fitting as compared to isochrone matching.
-
COSINE is a dark matter search experiment based on an array of low background
NaI(Tl) crystals located at the Yangyang underground laboratory. The assembly
of COSINE-100 was completed in the summer of 2016 and the detector is currently
collecting physics quality data aimed at reproducing the DAMA/LIBRA experiment
that reported an annual modulation signal. Stable operation has been achieved
and will continue for at least two years. Here, we describe the design of
COSINE-100, including the shielding arrangement, the configuration of the
NaI(Tl) crystal detection elements, the veto systems, and the associated
operational systems, and we show the current performance of the experiment.
-
We report a measurement of the flux-integrated cross section for inclusive
muon neutrino charged-current interactions on carbon. The double differential
measurements are given as function of the muon momentum and angle. Relative to
our previous publication on this topic, these results have an increased angular
acceptance and higher statistics. The data sample presented here corresponds to
$5.7 \times 10^{20}$ protons-on-target. The total flux-integrated cross section
is measured to be $(6.950 \pm 0.662) \times 10^{-39}$ cm$^2$nucleon$^{-1}$ and
is consistent with our simulation.
-
The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21
Treasury program on the Hubble Space Telescope. It studied 50 nearby
star-forming galaxies in five bands from the near UV to the I-band, combining
new Wide Field Camera 3 observations with archival Advanced Camera for Surveys
data. LEGUS was designed to investigate how star formation occurs and develops
on both small and large scales, and how it relates to the galactic
environments. In this paper we present the photometric catalogs for all the
apparently single stars identified in the 50 LEGUS galaxies. Photometric
catalogs and mosaicked images for all filters are available for download.
We present optical and near UV color-magnitude diagrams for all the galaxies.
For each galaxy we derived the distance from the tip of the red giant branch.
We then used the NUV color-magnitude diagrams to identify stars more massive
than 14 Mo, and compared their number with the number of massive stars expected
from the GALEX FUV luminosity. Our analysis shows that the fraction of massive
stars forming in star clusters and stellar associations is about constant with
the star formation rate. This lack of a relation suggests that the time scale
for evaporation of unbound structures is comparable or longer than 10 Myr. At
low star formation rates this translates to an excess of mass in clustered
environments as compared to model predictions of cluster evolution, suggesting
that a significant fraction of stars form in unbound systems.
-
The COSINE-100 dark matter search experiment has started taking physics data
with the goal of performing an independent measurement of the annual modulation
signal observed by DAMA/LIBRA. A muon detector was constructed by using plastic
scintillator panels in the outermost layer of the shield surrounding the
COSINE-100 detector. It is used to detect cosmic ray muons in order to
understand the impact of the muon annual modulation on dark matter analysis.
Assembly and initial performance test of each module have been performed at a
ground laboratory. The installation of the detector in Yangyang Underground
Laboratory (Y2L) was completed in the summer of 2016. Using three months of
data, the muon underground flux was measured to be 328 $\pm$ 1(stat.)$\pm$
10(syst.) muons/m$^2$/day. In this report, the assembly of the muon detector
and the results from the analysis are presented.
-
The single $\pi^0$ production rate in neutral current neutrino interactions
on water in a neutrino beam with a peak neutrino energy of 0.6 GeV has been
measured using the P{\O}D, one of the subdetectors of the T2K near detector.
The production rate was measured for data taking periods when the P{\O}D
contained water ($2.64\times{}10^{20}$ protons-on-target) and also periods
without water ($3.49 \times 10^{20}$ protons-on-target). A measurement of the
neutral current single $\pi^0$ production rate on water is made using
appropriate subtraction of the production rate with water in from the rate with
water out of the target region. The subtraction analysis yields 106 $\pm$ 41
(stat.) $\pm$ 69 (sys.) signal events, which is consistent with the prediction
of 157 events from the nominal simulation. The measured to expected ratio is
0.68 $\pm$ 0.26 (stat.) $\pm$ 0.44 (sys.) $\pm$ 0.12 (flux). The nominal
simulation uses a flux integrated cross section of $7.63\times{}10^{-39}$
cm${}^2$ per nucleon with an average neutrino interaction energy of 1.3 GeV.
-
Merging neutron stars offer an exquisite laboratory for simultaneously
studying strong-field gravity and matter in extreme environments. We establish
the physical association of an electromagnetic counterpart EM170817 to
gravitational waves (GW170817) detected from merging neutron stars. By
synthesizing a panchromatic dataset, we demonstrate that merging neutron stars
are a long-sought production site forging heavy elements by r-process
nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to
classical short gamma-ray bursts with ultra-relativistic jets. Instead, we
suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the
jet elegantly explains the low-luminosity gamma-rays, the high-luminosity
ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit
that all merging neutron stars may lead to a wide-angle cocoon breakout;
sometimes accompanied by a successful jet and sometimes a choked jet.
-
We generalize the notions of elliptic pseudoprimes and elliptic Carmichael
numbers introduced by Silverman to analogues of Euler-Jacobi and strong
pseudoprimes. We investigate the relationships among Euler Elliptic Carmichael
numbers , strong elliptic Carmichael numbers, products of anomalous primes and
elliptic Korselt numbers of Type I: The former two of these are introduced in
this paper, and the latter two of these were introduced by Mazur (1973) and
Silverman (2012) respectively. In particular, we expand upon a previous work of
Babinkostova et al. by proving a conjecture about the density of certain
elliptic Korselt numbers of Type I that are products of anomalous primes.
-
KAGRA is a 3-km interferometric gravitational wave telescope located in the
Kamioka mine in Japan. It is the first km-class gravitational wave telescope
constructed underground to reduce seismic noise, and the first km-class
telescope to use cryogenic cooling of test masses to reduce thermal noise. The
construction of the infrastructure to house the interferometer in the tunnel,
and the initial phase operation of the interferometer with a simple 3-km
Michelson configuration have been completed. The first cryogenic operation is
expected in 2018, and the observing runs with a full interferometer are
expected in 2020s. The basic interferometer configuration and the current
status of KAGRA are described.
-
Advanced satellite-based frequency transfers by TWCP and IPPP have been
performed between NICT and KRISS. We confirm that the disagreement between them
is less than 1x10^{-16} at an averaging time of several days. Additionally, an
intercontinental frequency ratio measurement of Sr and Yb optical lattice
clocks was directly performed by TWCP. We achieved an uncertainty at the
mid-10^{-16} level after a total measurement time of 12 hours. The frequency
ratio was consistent with the recently reported values within the uncertainty.
-
Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy
Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used
to study the young star cluster (YSC) population of this interacting system.
Our newly extracted source catalogue contains 2834 cluster candidates,
morphologically classified to be compact and uniform in colour, for which ages,
masses and extinction are derived. In this first work we study the main
properties of the YSC population of the whole galaxy, considering a
mass-limited sample. Both luminosity and mass functions follow a power law
shape with slope -2, but at high luminosities and masses a dearth of sources is
observed. The analysis of the mass function suggests that it is best fitted by
a Schechter function with slope -2 and a truncation mass at
$1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm
this result and link the shape of the luminosity function to the presence of a
truncation in the mass function. A mass limited age function analysis, between
10 and 200 Myr, suggests that the cluster population is undergoing only
moderate disruption. We observe little variation in the shape of the mass
function at masses above $1\times10^4\ M_\odot$, over this age range. The
fraction of star formation happening in the form of bound clusters in M51 is
$\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over
the whole range from 1 to 200 Myr.
-
Subgroups $H_1$ and $H_2$ of a group $G$ are said to be locally conjugate if
there is a bijection $f: H_1 \rightarrow H_2$ such that $h$ and $f(h)$ are
conjugate in $G$ for every $h \in H_1$. This paper studies local conjugacy
among subgroups of $\text{GL}_2(\mathbb{Z}/p^2\mathbb{Z})$, where $p$ is an odd
prime, building on Sutherland's categorizations of subgroups of
$\text{GL}_2(\mathbb{Z}/p\mathbb{Z})$ and local conjugacy among them. There are
two conditions that locally conjugate subgroups $H_1$ and $H_2$ of
$\text{GL}_2(\mathbb{Z}/p^2\mathbb{Z})$ must satisfy: letting $\varphi:
\text{GL}_2(\mathbb{Z}/p^2\mathbb{Z}) \rightarrow
\text{GL}_2(\mathbb{Z}/p\mathbb{Z})$ be the natural homomorphism, $H_1 \cap
\ker \varphi$ and $H_2 \cap \ker \varphi$ must be locally conjugate in
$\text{GL}_2(\mathbb{Z}/p^2\mathbb{Z})$ and $\varphi(H_1)$ and $\varphi(H_2)$
must be locally conjugate in $\text{GL}_2(\mathbb{Z}/p\mathbb{Z})$. To identify
$H_1$ and $H_2$ up to conjugation, we choose $\varphi(H_1)$ and $\varphi(H_2)$
to be similar to each other, then understand the possibilities for $H_1 \cap
\ker \varphi$ and $H_2 \cap \ker \varphi$. This study fully categorizes local
conjugacy in $\text{GL}_2(\mathbb{Z}/p^2\mathbb{Z})$ through such casework.
-
The T2K experiment reports an updated analysis of neutrino and antineutrino
oscillations in appearance and disappearance channels. A sample of electron
neutrino candidates at Super-Kamiokande in which a pion decay has been tagged
is added to the four single-ring samples used in previous T2K oscillation
analyses. Through combined analyses of these five samples, simultaneous
measurements of four oscillation parameters, $|\Delta m^2_{32}|$,
$\sin^2(\theta_{23})$, $\sin^2(\theta_{13})$, and $\delta_{CP}$ and of the mass
ordering are made. A set of studies of simulated data indicates that the
sensitivity to the oscillation parameters is not limited by neutrino
interaction model uncertainty. Multiple oscillation analyses are performed, and
frequentist and Bayesian intervals are presented for combinations of the
oscillation parameters with and without the inclusion of reactor constraints on
$\sin^2(\theta_{13})$. When combined with reactor measurements, the hypothesis
of CP conservation ($\delta_{CP}$$=0$ or $\pi$) is excluded at 90% confidence
level. The 90% confidence region for $\delta_{CP}$ is [-2.95,-0.44] ([-1.47,
-1.27]) for normal (inverted) ordering. The central values and 68\% confidence
intervals for the other oscillation parameters for normal (inverted) ordering
are $\Delta m^{2}_{32}=2.54\pm0.08$ ($2.51\pm0.08$) $\times 10^{-3}$ eV$^2 /
c^4$ and $\sin^2(\theta_{23}) = 0.55^{+0.05}_{-0.09}$ ($0.55^{+0.05}_{-0.08}$),
compatible with maximal mixing. In the Bayesian analysis, the data weakly
prefer normal ordering (Bayes factor 3.7) and the upper octant for
$\sin^2(\theta_{23})$ (Bayes factor 2.4).
-
We report measurements by the T2K experiment of the parameters $\theta_{23}$
and $\Delta m^{2}_{32}$ governing the disappearance of muon neutrinos and
antineutrinos in the three flavor neutrino oscillation model. Utilizing the
ability of the experiment to run with either a mainly neutrino or a mainly
antineutrino beam, the parameters are measured separately for neutrinos and
antineutrinos. Using $7.482 \times 10^{20}$ POT in neutrino running mode and
$7.471 \times 10^{20}$ POT in antineutrino mode, T2K obtained,
$\sin^{2}(\theta_{23})=0.51^{+0.08}_{-0.07}$ and $\Delta m^{2}_{32} =
2.53^{+0.15}_{-0.13} \times 10^{-3}$eV$^{2}$/c$^{4}$ for neutrinos, and
$\sin^{2}({\overline{\theta}}_{23})=0.42^{+0.25}_{-0.07}$ and
${\Delta\overline{m}^2}_{32} = 2.55^{+0.33}_{-0.27} \times
10^{-3}$eV$^{2}$/c$^{4}$ for antineutrinos (assuming normal mass ordering). No
significant differences between the values of the parameters describing the
disappearance of muon neutrinos and antineutrinos were observed.
-
We report a measurement of cross section $\sigma(\nu_{\mu}+{\rm
nucleus}\rightarrow\mu^{-}+X)$ and the first measurements of the cross section
$\sigma(\bar{\nu}_{\mu}+{\rm nucleus}\rightarrow\mu^{+}+X)$ and their ratio
$R(\frac{\sigma(\bar \nu)}{\sigma(\nu)})$ at (anti-)neutrino energies below 1.5
GeV. We determine the single momentum bin cross section measurements, averaged
over the T2K $\bar{\nu}/\nu$-flux, for the detector target material (mainly
Carbon, Oxygen, Hydrogen and Copper) with phase space restricted laboratory
frame kinematics of $\theta_{\mu}<32^\circ$ and $p_{\mu}>$500 MeV/c. The
results are $\sigma(\bar{\nu})=\left( 0.900\pm0.029{\rm (stat.)}\pm0.088{\rm
(syst.)}\right)\times10^{-39}$ and $\sigma(\nu)=\left( 2.41\
\pm0.022{\rm{(stat.)}}\pm0.231{\rm (syst.)}\ \right)\times10^{-39}$ in units of
cm$^{2}$/nucleon and $R\left(\frac{\sigma(\bar{\nu})}{\sigma(\nu)}\right)=
0.373\pm0.012{\rm (stat.)}\pm0.015{\rm (syst.)}$.
-
A class of extensions of the Standard Model allows Lorentz and CPT
violations, which can be identified by the observation of sidereal modulations
in the neutrino interaction rate. A search for such modulations was performed
using the T2K on-axis near detector. Two complementary methods were used in
this study, both of which resulted in no evidence of a signal. Limits on
associated Lorentz and CPT violating terms from the Standard Model Extension
have been derived taking into account their correlations in this model for the
first time. These results imply such symmetry violations are suppressed by a
factor of more than $10^{20}$ at the GeV scale.
-
We present an analysis of the positions and ages of young star clusters in
eight local galaxies to investigate the connection between the age difference
and separation of cluster pairs. We find that star clusters do not form
uniformly but instead are distributed such that the age difference increases
with the cluster pair separation to the 0.25-0.6 power, and that the maximum
size over which star formation is physically correlated ranges from ~200 pc to
~1 kpc. The observed trends between age difference and separation suggest that
cluster formation is hierarchical both in space and time: clusters that are
close to each other are more similar in age than clusters born further apart.
The temporal correlations between stellar aggregates have slopes that are
consistent with turbulence acting as the primary driver of star formation. The
velocity associated with the maximum size is proportional to the galaxy's
shear, suggesting that the galactic environment influences the maximum size of
the star-forming structures.
-
Precise neutrino--nucleus interaction measurements in the sub-multi GeV
region are important to reduce the systematic uncertainty in future neutrino
oscillation experiments. Furthermore, the excess of ${\nu_e}$ interactions, as
a possible interpretation of the existence of a sterile neutrino has been
observed in such an energy region. The nuclear emulsion technique can measure
all the final state particles with low energy threshold for a variety of
targets (Fe, C, H${_2}$O, and so on). Its sub-$\mu$m position resolution allows
measurements of the ${\nu_e}$ cross-section with good electron/gamma separation
capability. We started a new experiment at J-PARC to study sub-multi GeV
neutrino interactions by introducing the nuclear emulsion technique. The J-PARC
T60 experiment has been implemented as a first step of such a project.
Systematic neutrino event analysis with full scanning data in the nuclear
emulsion detector was performed for the first time. The first neutrino event
detection and its analysis is described in this paper.
-
We present a study of the effective (half-light) radii and other structural
properties of a systematically selected sample of young, massive star clusters
(YMCs, $\geq$$5\times10^3$ M$_{\odot}$ and $\leq$200 Myr) in two nearby spiral
galaxies, NGC 628 and NGC 1313. We use Hubble Space Telescope WFC3/UVIS and
archival ACS/WFC data obtained by the Legacy Extragalactic UV Survey (LEGUS),
an HST Treasury Program. We measure effective radii with GALFIT, a
two-dimensional image-fitting package, and with a new technique to estimate
effective radii from the concentration index (CI) of observed clusters. The
distribution of effective radii from both techniques spans $\sim$0.5-10 pc and
peaks at 2-3 pc for both galaxies. We find slight positive correlations between
effective radius and cluster age in both galaxies, but no significant
relationship between effective radius and galactocentric distance. Clusters in
NGC 1313 display a mild increase in effective radius with cluster mass, but the
trend disappears when the sample is divided into age bins. We show that the
vast majority of the clusters in both galaxies are much older than their
dynamical times, suggesting they are gravitationally bound objects. We find
that about half of the clusters in NGC 628 are underfilling their Roche lobes,
based on their Jacobi radii. Our results suggest that the young, massive
clusters in NGC 628 and NGC 1313 are expanding due to stellar mass loss or
two-body relaxation and are not significantly influenced by the tidal fields of
their host galaxies.
-
We report the large effort which is producing comprehensive high-level young
star cluster (YSC) catalogues for a significant fraction of galaxies observed
with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We
present the methodology developed to extract cluster positions, verify their
genuine nature, produce multiband photometry (from NUV to NIR), and derive
their physical properties via spectral energy distribution fitting analyses. We
use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact
that LEGUS will have on our understanding of the formation and evolution of
YSCs and compact stellar associations within their host galaxy. Our analysis of
the cluster luminosity function from the UV to the NIR finds a steepening at
the bright end and at all wavelengths suggesting a dearth of luminous clusters.
The cluster mass function of NGC628 is consistent with a power-law distribution
of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After
their formation YSCs and compact associations follow different evolutionary
paths. YSCs survive for a longer timeframe, confirming their being potentially
bound systems. Associations disappear on time scales comparable to
hierarchically organized star-forming regions, suggesting that they are
expanding systems. We find mass-independent cluster disruption in the inner
region of NGC628, while in the outer part of the galaxy there is little or no
disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$
M$_\odot$) clusters suggesting that a mass-dependent component is necessary to
fully describe the YSC disruption process in NGC628.
-
SS~433 is an X-ray binary and the source of sub-relativistic, precessing,
baryonic jets. We present high-resolution spectrograms of SS 433 in the
infrared H and K bands. The spectrum is dominated by hydrogen and helium
emission lines. The precession phase of the emission lines from the jet
continues to be described by a constant period, P_jet= 162.375 d. The limit on
any secularly changing period is $|\dot P| \lesssim 10^{-5}$. The He I 2.0587
micron line has complex and variable P Cygni absorption features produced by an
inhomogeneous wind with a maximum outflow velocity near 900 km/s. The He II
emission lines in the spectrum also arise in this wind. The higher members of
the hydrogen Brackett lines show a double-peaked profile with symmetric wings
extending more than +/-1500 km/s from the line center. The lines display radial
velocity variations in phase with the radial velocity variation expected of the
compact star, and they show a distortion during disk eclipse that we interpret
as a rotational distortion. We fit the line profiles with a model in which the
emission comes from the surface of a symmetric, Keplerian accretion disk around
the compact object. The outer edge of the disk has velocities that vary from
110 to 190 km/s. These comparatively low velocities place an important
constraint on the mass of the compact star: Its mass must be less than 2.2
M_solar and is probably less than 1.6 M_solar.
-
We present a study of the hierarchical clustering of the young stellar
clusters in six local (3--15 Mpc) star-forming galaxies using Hubble Space
Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program
LEGUS (Legacy ExtraGalactic UV Survey). We have identified 3685 likely clusters
and associations, each visually classified by their morphology, and we use the
angular two-point correlation function to study the clustering of these stellar
systems. We find that the spatial distribution of the young clusters and
associations are clustered with respect to each other, forming large, unbound
hierarchical star-forming complexes that are in general very young. The
strength of the clustering decreases with increasing age of the star clusters
and stellar associations, becoming more homogeneously distributed after ~40--60
Myr and on scales larger than a few hundred parsecs. In all galaxies, the
associations exhibit a global behavior that is distinct and more strongly
correlated from compact clusters. Thus, populations of clusters are more
evolved than associations in terms of their spatial distribution, traveling
significantly from their birth site within a few tens of Myr whereas
associations show evidence of disruption occurring very quickly after their
formation. The clustering of the stellar systems resembles that of a turbulent
interstellar medium that drives the star formation process, correlating the
components in unbound star-forming complexes in a hierarchical manner,
dispersing shortly after formation, suggestive of a single, continuous mode of
star formation across all galaxies.