
We examine the sensitivity of a large scale twophase liquid argon detector
to the directionality of the dark matter signal. This study was performed under
the assumption that, above 50 keV of recoil energy, one can determine (with
some resolution) the direction of the recoil nucleus without headtail
discrimination, as suggested by past studies that proposed to exploit the
dependence of columnar recombination on the angle between the recoil nucleus
direction and the electric field. In this paper we study the differential
interaction recoil rate as a function of the recoil direction angle with
respect to the zenith for a detector located at the Laboratori Nazionali del
Gran Sasso and we determine its diurnal and seasonal modulation. Using a
likelihoodratio based approach we show that, with the angular information
alone, 100 events are enough to reject the isotropic hypothesis at three
standard deviation level. For an exposure of 100 tonne years this would
correspond to a spin independent WIMPnucleon cross section of about 10^46
cm^2 at 200 GeV WIMP mass. The results presented in this paper provide strong
motivation for the experimental determination of directional recoil effects in
twophase liquid argon detectors.

We present the results of a search for dark matter WIMPs in the mass range
below 20 GeV/c^2 using a target of lowradioactivity argon. The data were
obtained using the DarkSide50 apparatus at Laboratori Nazionali del Gran Sasso
(LNGS). The analysis is based on the ionization signal, for which the
DarkSide50 time projection chamber is fully efficient at 0.1 keVee. The
observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and
is almost entirely accounted for by known background sources. We obtain a 90%
C.L. exclusion limit above 1.8 GeV/c^2 for the spinindependent cross section
of dark matter WIMPs on nucleons, extending the exclusion region for dark
matter below previous limits in the range 1.86 GeV/c^2.

We present new constraints on subGeV dark matter particles scattering off
electrons in argon based on an analysis of ionization signal data from the
DarkSide50 detector.

The DarkSide50 directdetection dark matter experiment is a dualphase argon
time projection chamber operating at Laboratori Nazionali del Gran Sasso. This
paper reports on the blind analysis and spinindependent dark matternucleon
coupling results from a 532.4 livedays exposure, using a target of
lowradioactivity argon extracted from underground sources. The backgroundfree
result in the dark matter selection box gave no evidence for dark matter. The
present blind analysis sets a 90% C.L. upper limit on the dark matternucleon
spinindependent cross section of 1.1E44 cm^2 (3.8E44 cm^2, 3.4E43 cm^2) for
a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

We report the measurement of the longitudinal diffusion constant in liquid
argon with the DarkSide50 dualphase time projection chamber. The measurement
is performed at drift electric fields of 100 V/cm, 150 V/cm, and 200 V/cm using
high statistics $^{39}$Ar decays from atmospheric argon. We derive an
expression to describe the pulse shape of the electroluminescence signal (S2)
in dualphase TPCs. The derived S2 pulse shape is fit to events from the
uppermost portion of the TPC in order to characterize the radial dependence of
the signal. The results are provided as inputs to the measurement of the
longitudinal diffusion constant DL, which we find to be (4.12 $\pm$ 0.04)
cm$^2$/s for a selection of 140keV electron recoil events in 200V/cm drift
field and 2.8kV/cm extraction field. To study the systematics of our
measurement we examine datasets of varying event energy, field strength, and
detector volume yielding a weighted average value for the diffusion constant of
(4.09 $\pm$ 0.09) cm$^2$ /s. The measured longitudinal diffusion constant is
observed to have an energy dependence, and within the studied energy range the
result is systematically lower than other results in the literature.

A Geant4based Monte Carlo package named G4DS has been developed to simulate
the response of DarkSide50, an experiment operating since 2013 at LNGS,
designed to detect WIMP interactions in liquid argon. In the process of WIMP
searches, DarkSide50 has achieved two fundamental milestones: the rejection of
electron recoil background with a power of ~10^7, using the pulse shape
discrimination technique, and the measurement of the residual 39Ar
contamination in underground argon, ~3 orders of magnitude lower with respect
to atmospheric argon. These results rely on the accurate simulation of the
detector response to the liquid argon scintillation, its ionization, and
electronion recombination processes. This work provides a complete overview of
the DarkSide Monte Carlo and of its performance, with a particular focus on
PARIS, the custommade liquid argon response model.

We report on the cryogenic characterization of Red Green Blue  High Density
(RGBHD) SiPMs developed at Fondazione Bruno Kessler (FBK) as part of the
DarkSide program of dark matter searches with liquid argon time projection
chambers. A dedicated setup was used to measure the primary dark noise, the
correlated noise, and the gain of the SiPMs at varying temperatures. A
custommade data acquisition system and analysis software were used to
precisely characterize these parameters. We demonstrate that FBK RGBHD SiPMs
with low quenching resistance (RGBHDLR$_q$) can be operated from 40 K to 300
K with gains in the range $10^5$ to $10^6$ and noise rates on the order of a
few Hz/mm$^2$.

Building on the successful experience in operating the DarkSide50 detector,
the DarkSide Collaboration is going to construct DarkSide20k, a direct WIMP
search detector using a twophase Liquid Argon Time Projection Chamber (LArTPC)
with an active (fiducial) mass of 23 t (20 t). The DarkSide20k LArTPC will be
deployed within a shield/veto with a spherical Liquid Scintillator Veto (LSV)
inside a cylindrical Water Cherenkov Veto (WCV). Operation of DarkSide50
demonstrated a major reduction in the dominant $^{39}$Ar background when using
argon extracted from an underground source, before applying pulse shape
analysis. Data from DarkSide50, in combination with MC simulation and
analytical modeling, shows that a rejection factor for discrimination between
electron and nuclear recoils of $\gt3\times10^9$ is achievable. This, along
with the use of the veto system, is the key to unlocking the path to large
LArTPC detector masses, while maintaining an "instrumental backgroundfree"
experiment, an experiment in which less than 0.1 events (other than
$\nu$induced nuclear recoils) is expected to occur within the WIMP search
region during the planned exposure. DarkSide20k will have ultralow
backgrounds than can be measured in situ. This will give sensitivity to
WIMPnucleon cross sections of $1.2\times10^{47}$ cm$^2$ ($1.1\times10^{46}$
cm$^2$) for WIMPs of $1$ TeV$/c^2$ ($10$ TeV$/c^2$) mass, to be achieved during
a 5 yr run producing an exposure of 100 t yr free from any instrumental
background. DarkSide20k could then extend its operation to a decade,
increasing the exposure to 200 t yr, reaching a sensitivity of
$7.4\times10^{48}$ cm$^2$ ($6.9\times10^{47}$ cm$^2$) for WIMPs of $1$
TeV$/c^2$ ($10$ TeV$/c^2$) mass.

Superconducting hybrid junctions are revealing a variety of novel effects.
Some of them are due to the special layout of these devices, which often use a
coplanar configuration with relatively large barrier channels and the
possibility of hosting Pearl vortices. A Josephson junction with a quasi ideal
twodimensional barrier has been realized by growing graphene on SiC with Al
electrodes. Chemical Vapor Deposition offers centimeter size monolayer areas
where it is possible to realize a comparative analysis of different devices
with nominally the same barrier. In samples with a graphene gap below 400 nm,
we have found evidence of Josephson coherence in presence of an incipient
BerezinskiiKosterlitzThouless transition. When the magnetic field is cycled,
a remarkable hysteretic collapse and revival of the Josephson supercurrent
occurs. Similar hysteresis are found in granular systems and are usually
justified within the Bean Critical State model (CSM). We show that the CSM,
with appropriate account for the low dimensional geometry, can partly explain
the odd features measured in these junctions.

The transport properties at finite temperature of crystalline organic
semiconductors are investigated, within the SuSchriefferHeeger model, by
combining exact diagonalization technique, Monte Carlo approaches, and maximum
entropy method. The temperaturedependent mobility data measured in single
crystals of rubrene are successfully reproduced: a crossover from super to
subdiffusive motion occurs in the range $150 \leq T \leq 200$ K, where the
mean free path becomes of the order of the lattice parameter and strong memory
effects start to appear. We provide an effective model which can successfully
explain low frequencies features of the absorption spectra. The observed
response to slowly varying electric field is interpreted by means of a simple
model where the interaction between the charge carrier and lattice polarization
modes is simulated by a harmonic interaction between a fictitious particle and
an electron embedded in a viscous fluid.

We present the first unbiased results for the mobility $\mu$ of
onedimensional Holstein polaron obtained by numerical analytic continuation
combined with diagrammatic and worldline Monte Carlo methods in the
thermodynamic limit. We have identified for the first time, by the
characteristic $\omega$ and $T$ dependence in the wide region of parameters,
several distinct regimes in the $\lambdaT$ plane including band conduction
region, incoherent metallic region, activated hopping region, and high
temperature saturation region. We observe for the first time that although
mobilities and mean free paths at different values of $\lambda$ differ by many
orders of magnitude at small temperatures, their values at $T$ larger than the
bandwidth become very close to each other.

The Kubo formula for the electrical conductivity is rewritten in terms of a
sum of Drudelike contributions associated to the exact eigenstates of the
interacting system, each characterized by its own frequencydependent
relaxation time. The structure of the novel and equivalent formulation,
weighting the contribution from each eigenstate by its Boltzmann occupation
factor, simplifies considerably the access to the static properties (dc
conductivity) and resolves the long standing difficulties to recover the
Boltzmann result for dc conductivity from the Kubo formula. It is shown that
the Boltzmann result, containing the correct transport scattering time instead
of the electron lifetime determined by the Green function, can be recovered in
problems with elastic and inelastic scattering at the lowest order of
interaction.

We study the storage and retrieval of phasecoded patterns as stable
dynamical attractors in recurrent neural networks, for both an analog and a
integrateandfire spiking model. The synaptic strength is determined by a
learning rule based on spiketimedependent plasticity, with an asymmetric time
window depending on the relative timing between pre and postsynaptic
activity. We store multiple patterns and study the network capacity.
For the analog model, we find that the network capacity scales linearly with
the network size, and that both capacity and the oscillation frequency of the
retrieval state depend on the asymmetry of the learning time window. In
addition to fullyconnected networks, we study sparse networks, where each
neuron is connected only to a small number z << N of other neurons. Connections
can be short range, between neighboring neurons placed on a regular lattice, or
long range, between randomly chosen pairs of neurons. We find that a small
fraction of long range connections is able to amplify the capacity of the
network. This imply that a smallworldnetwork topology is optimal, as a
compromise between the cost of long range connections and the capacity
increase.
Also in the spiking integrate and fire model the crucial result of storing
and retrieval of multiple phasecoded patterns is observed. The capacity of the
fullyconnected spiking network is investigated, together with the relation
between oscillation frequency of retrieval state and window asymmetry.

We study the dynamical behavior of a lattice model of glass former on a
random graph, where no corrections to the mean field description are expected.
We find that the behavior of dynamical correlation functions and dynamical
susceptibility are consistent with the quantitative predictions of the Mode
Coupling Theory of the glass transition.

We present a systematic study of dynamical heterogeneity in a model for
permanent gels, upon approaching the gelation threshold. We find that the
fluctuations of the self intermediate scattering function are increasing
functions of time, reaching a plateau whose value, at large length scales,
coincides with the mean cluster size and diverges at the percolation threshold.
Another measure of dynamical heterogeneities, i.e. the fluctuations of the
selfoverlap, displays instead a peak and decays to zero at long times. The
peak, however, also scales as the mean cluster size. Arguments are given for
this difference in the long time behavior. We also find that nonGaussian
parameter reaches a plateau in the long time limit. The value of the plateau of
the nonGaussian parameter, which is connected to the fluctuations of
diffusivity of clusters, increases with the volume fraction and remains finite
at percolation threshold.

We study the formation of a colloidal gel by means of Molecular Dynamics
simulations of a model for colloidal suspensions. A slowing down with gellike
features is observed at low temperatures and low volume fractions, due to the
formation of persistent structures. We show that at low volume fraction the
dynamic susceptibility, which describes dynamic heterogeneities, exhibits a
large plateau, dominated by clusters of long living bonds. At higher volume
fraction, where the effect of the crowding of the particles starts to be
present, it crosses over towards a regime characterized by a peak. We introduce
a suitable mean cluster size of clusters of monomers connected by "persistent"
bonds which well describes the dynamic susceptibility.

We compare the slow dynamics of irreversible gels, colloidal gels, glasses
and spin glasses by analyzing the behavior of the so called nonlinear
dynamical susceptibility, a quantity usually introduced to quantitatively
characterize the dynamical heterogeneities. In glasses this quantity typically
grows with the time, reaches a maximum and then decreases at large time, due to
the transient nature of dynamical heterogeneities and to the absence of a
diverging static correlation length. We have recently shown that in
irreversible gels the dynamical susceptibility is instead an increasing
function of the time, as in the case of spin glasses, and tends asymptotically
to the mean cluster size. On the basis of molecular dynamics simulations, we
here show that in colloidal gelation where clusters are not permanent, at very
low temperature and volume fractions, i.e. when the lifetime of the bonds is
much larger than the structural relaxation time, the nonlinear susceptibility
has a behavior similar to the one of the irreversible gel, followed, at higher
volume fractions, by a crossover towards the behavior of glass forming liquids.

In colloidal suspensions, the competition between attractive and repulsive
interactions gives rise to a rich and complex phenomenology. Here, we study the
equilibrium phase diagram of a model system using a DLVO interaction potential
by means of molecular dynamics simulations and a thermodynamical approach. As a
result, we find tubular and lamellar phases at low volume fraction. Such
phases, extremely relevant for designing new materials, may be not easily
observed in the experiments because of the long relaxation times and the
presence of defects.

We study the structure and the dynamics in the formation of irreversible gels
by means of molecular dynamics simulation of a model system where the gelation
transition is due to the random percolation of permanent bonds between
neighboring particles. We analyze the heterogeneities of the dynamics in terms
of the fluctuations of the intermediate scattering functions: In the sol phase
close to the percolation threshold, we find that this dynamical susceptibility
increases with the time until it reaches a plateau. At the gelation threshold
this plateau scales as a function of the wave vector $k$ as $k^{\eta 2}$, with
$\eta$ being related to the decay of the percolation pair connectedness
function. At the lowest wave vector, approaching the gelation threshold it
diverges with the same exponent $\gamma$ as the mean cluster size. These
findings suggest an alternative way of measuring critical exponents in a system
undergoing chemical gelation.

By means of molecular dynamics, we study a model system for colloidal
suspensions where the interaction is based on a competition between attraction
and repulsion. At low temperatures the relaxation time $\tau$ first increases
as a power law as a function of the volume fraction $\phi$ and then, due to the
finite lifetime of the bonded structures, it deviates from this critical
behavior. We show that colloidal gelation at low temperatures and low volume
fractions is crucially related to the formation of spanning long living
cluster. Besides agreeing with experimental findings in different colloidal
systems, our results shed new light on the different role played by the
formation of long living bonds and the crowding of the particles in colloidal
structural arrest.

In order to study analytically the nature of the jamming transition in
granular material, we have considered a cavity method mean field theory, in the
framework of a statistical mechanics approach, based on Edwards' original idea.
For simplicity we have applied the theory to a lattice model and a transition
with exactly the same nature of the glass transition in mean field models for
usual glass formers is found. The model is also simulated in three dimensions
under tap dynamics and a jamming transition with glassy features is observed.
In particular two step decays appear in the relaxation functions and dynamic
heterogeneities resembling ones usually observed in glassy systems. These
results confirm early speculations about the connection between the jamming
transition in granular media and the glass transition in usual glass formers,
giving moreover a precise interpretation of its nature.

In the framework of schematic hard spheres lattice models for granular media
we investigate the phenomenon of the ``jamming transition''. In particular,
using Edwards' approach, by analytical calculations at a mean field level, we
derive the system phase diagram and show that ``jamming'' corresponds to a
phase transition from a ``fluid'' to a ``glassy'' phase, observed when
crystallization is avoided. Interestingly, the nature of such a ``glassy''
phase turns out to be the same found in mean field models for glass formers.

We study a model for the gel degradation by an enzyme, where the gel is
schematized as a cubic lattice, and the enzyme as a random walker, that cuts
the bonds over which it passes. The model undergoes a (reverse) percolation
transition, which for low density of enzymes falls in a universality class
different from random percolation. In particular we have measured a gel
fraction critical exponent beta=1.0+0.1, in excellent agreement with
experiments made on the real system.

In the framework of schematic hard spheres lattice models we discuss Edwards'
Statistical Mechanics approach to granular media. As this approach appears to
hold here to a very good approximation, by analytical calculations of Edwards'
partition function at a mean field level we derive the system phase diagram and
show that ``jamming'' corresponds to a phase transition from a ``fluid'' to a
``glassy'' phase, observed when crystallization is avoided. The nature of such
a ``glassy'' phase turns out to be the same found in mean field models for
glass formers. In the same context, we also briefly discuss mixing/segregation
phenomena of binary mixtures: the presence of fluidcrystal phase transitions
drives segregation as a form of phase separation and, within a given phase,
gravity can also induce a kind of ``vertical'' segregation, usually not
associated to phase transitions.

Experiments of in vitro formation of blood vessels show that cells randomly
spread on a gel matrix autonomously organize to form a connected vascular
network. We propose a simple model which reproduces many features of the
biological system. We show that both the model and the real system exhibit a
fractal behavior at small scales, due to the process of migration and dynamical
aggregation, followed at large scale by a random percolation behavior due to
the coalescence of aggregates. The results are in good agreement with the
analysis performed on the experimental data.