
We investigate the lowenergy scaling behavior of an interacting 3D Weyl
semimetal in the presence of disorder. In order to achieve a renormalization
group analysis of the theory, we focus on the effects of a
shortrangedcorrelated disorder potential, checking nevertheless that this
choice is not essential to locate the different phases of the Weyl semimetal.
We show that there is a line of fixedpoints in the renormalization group flow
of the interacting theory, corresponding to the disorderdriven transition to a
diffusive metal phase. Along that boundary, the critical disorder strength
undergoes a strong increase with respect to the noninteracting theory, as a
consequence of the unconventional screening of the Coulomb and disorderinduced
interactions. A complementary resolution of the SchwingerDyson equations
allows us to determine the full phase diagram of the system, showing the
prevalence of a renormalized semimetallic phase in the regime of intermediate
interaction strength, and adjacent to the nonFermi liquid phase characteristic
of the strong interaction regime of 3D Weyl semimetals.

We investigate the topological protection of surface states in Weyl and
nodalline semimetals by characterizing them as evanescent states when the band
structure is extended to complex momenta. We find in this way a sequence of
exceptional points that is, branch points with zero energy in the complex
spectrum allowing us to identify the set of surface states with complex
momentum signaling the decay into the 3D semimetal. From this point of view,
Weyl and nodalline semimetals can be classified in two types depending on the
way surface states decay. Type A semimetals have surface states with smaller
penetration length and oscillating decay while type B semimetals have longer
simple exponential decays. The difference between both types reflects in the
way the branch cuts in the spectrum accommodate in the complex plane. The
stability of the surface states stems in this approach from the complex
structure that develops around the exceptional points, with a topological
protection which is based on the fact that the branch cuts cannot be closed by
small perturbations. We check this property when nodalline semimetals are
placed under circularly polarized light, where we observe that the exceptional
points survive the effect of such a perturbation, though appropriate boundary
conditions for zeroenergy surface states cannot be satisfied in general due to
the breakdown of timereversal invariance by the radiation field.

The study of cosmology, galaxy formation and exoplanets has now advanced to a
stage where a cosmic inventory of terrestrial planets may be attempted. By
coupling semianalytic models of galaxy formation to a recipe that relates the
occurrence of planets to the mass and metallicity of their host stars, we trace
the population of terrestrial planets around both solarmass (FGK type) and
lowermass (M dwarf) stars throughout all of cosmic history. We find that the
mean age of terrestrial planets in the local Universe is $7\pm{}1$ Gyr for FGK
hosts and $8\pm{}1$ Gyr for M dwarfs. We estimate that hot Jupiters have
depleted the population of terrestrial planets around FGK stars by no more than
$\approx 10\%$, and that only $\approx 10\%$ of the terrestrial planets at the
current epoch are orbiting stars in a metallicity range for which such planets
have yet to be confirmed. The typical terrestrial planet in the local Universe
is located in a spheroiddominated galaxy with a total stellar mass comparable
to that of the Milky Way. When looking at the inventory of planets throughout
the whole observable Universe, we argue for a total of $\approx 1\times
10^{19}$ and $\approx 5\times 10^{20}$ terrestrial planets around FGK and M
stars, respectively. Due to light travel time effects, the terrestrial planets
on our past light cone exhibit a mean age of just $1.7\pm 0.2$ Gyr. These
results are discussed in the context of cosmic habitability, the Copernican
principle and searches for extraterrestrial intelligence at cosmological
distances.

We investigate the effect of shear and strain in graphene bilayers, under
conditions where the distortion of the lattice gives rise to a smooth
onedimensional modulation in the stacking sequence of the bilayer. We show
that strain and shear produce characteristic Moir\'e patterns which can have
the same visual appearance on a large scale, but representing graphene bilayers
with quite different electronic properties. The different features in the
lowenergy electronic bands can be ascribed to the effect of a fictitious
nonAbelian gauge field mimicking the smooth modulation of the stacking order.
Strained and sheared bilayers show a complementary behavior, which can be
understood from the fact that the nonAbelian gauge field acts as a repulsive
interaction in the former, expelling the electron density away from the
stacking domain walls, while behaving as a confining interaction leading to
localization of the electronic states in the sheared bilayers. In this latter
case, the presence of the effective gauge field explains the development of
almost flat lowenergy bands, resembling the form of the zeroth Landau level
characteristic of a Dirac fermion field. The estimate of the gauge field
strength in those systems gives a magnitude of the order of several tens of
Tesla, implying a robust phenomenology that should be susceptible of being
observed in suitably distorted bilayer samples.

We give a short compendium of the main ongoing detectors and concepts capable
of performing accurate sub100 ps timing at high particle fluxes and on large
areas, through technologies based on gaseous media. We briefly discuss the
stateoftheart, technological limitations and prospects, and a new bizarre
idea.

A three species onedimensional kinetic model is presented for a spatially
homogeneous weakly ionized plasma subjected to the action of a time varying
electric field. Planar geometry is assumed, which means that the plasma
dynamics evolves in the privileged direction of the field. The energy
transmitted to the charges is be channelized to the neutrals thanks to
collisions and impacting the plasma dynamics. Chargecharge interactions have
been designed as a one dimensional collision term equivalent to the Landau
operator used for fully ionized plasmas. Chargeneutral collisions are modelled
by a conservative drift diffusion operator in the Doughertys form. The
resulting set of coupled drift diffusion equations is solved with the stable
and robust Propagator Integral Method. This method feasibility accounts for
nonlinear effects without appealing to linearisation or simplifications,
providing conservative physically meaningful solutions. It is found that charge
neutral collisions exert a significant effect since a quite different plasma
dynamics arises if compared to the collisionless limit. In addition,
substantial differences in the system evolution are found for constant and
temperature dependent collision frequencies cases.

We report results of a wideband search for periodic gravitational waves from
isolated neutron stars within the Orion spur towards both the inner and outer
regions of our Galaxy. As gravitational waves interact very weakly with matter,
the search is unimpeded by dust and concentrations of stars. One search disk
(A) is $6.87^\circ$ in diameter and centered on
$20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other
(B) is $7.45^\circ$ in diameter and centered on
$8^\textrm{h}35^\textrm{m}20.61^\textrm{s}46^\circ49'25.151"$. We explored the
frequency range of 501500 Hz and frequency derivative from $0$ to $5\times
10^{9}$ Hz/s. A multistage, loosely coherent search program allowed probing
more deeply than before in these two regions, while increasing coherence length
with every stage.
Rigorous followup parameters have winnowed initial coincidence set to only 70
candidates, to be examined manually. None of those 70 candidates proved to be
consistent with an isolated gravitational wave emitter, and 95% confidence
level upper limits were placed on continuouswave strain amplitudes. Near $169$
Hz we achieve our lowest 95% CL upper limit on worstcase linearly polarized
strain amplitude $h_0$ of $6.3\times 10^{25}$, while at the high end of our
frequency range we achieve a worstcase upper limit of $3.4\times 10^{24}$ for
all polarizations and sky locations.

In this paper we present the results of the first low frequency allsky
search of continuous gravitational wave signals conducted on Virgo VSR2 and
VSR4 data. The search covered the full sky, a frequency range between 20 Hz and
128 Hz with a range of spindown between $1.0 \times 10^{10}$ Hz/s and $+1.5
\times 10^{11}$ Hz/s, and was based on a hierarchical approach. The starting
point was a set of short Fast Fourier Transforms (FFT), of length 8192 seconds,
built from the calibrated strain data. Aggressive data cleaning, both in the
time and frequency domains, has been done in order to remove, as much as
possible, the effect of disturbances of instrumental origin. On each dataset a
number of candidates has been selected, using the FrequencyHough transform in
an incoherent step. Only coincident candidates among VSR2 and VSR4 have been
examined in order to strongly reduce the false alarm probability, and the most
significant candidates have been selected. Selected candidates have been
subject to a followup by constructing a new set of longer FFTs followed by a
further incoherent analysis, still based on the FrequencyHough transform. No
evidence for continuous gravitational wave signals was found, therefore we have
set a populationbased joint VSR2VSR4 90$\%$ confidence level upper limit on
the dimensionless gravitational wave strain in the frequency range between 20
Hz and 128 Hz. This is the first allsky search for continuous gravitational
waves conducted, on data of groundbased interferometric detectors, at
frequencies below 50 Hz. We set upper limits in the range between about
$10^{24}$ and $2\times 10^{23}$ at most frequencies. Our upper limits on
signal strain show an improvement of up to a factor of $\sim$2 with respect to
the results of previous allsky searches at frequencies below $80~\mathrm{Hz}$.

We present our extensive observational campaign on the Swiftdiscovered
GRB141121A, al most ten years after its launch. Our observations covers radio
through Xrays, and extends for more than 30 days after discovery. The prompt
phase of GRB 141121A lasted 1410 s and, at the derived redshift of z = 1.469,
the isotropic energy is E{\gamma},iso = 8.0x10^52 erg. Due to the long prompt
duration, GRB141121A falls into the recently discovered class of ULGRBs.
Peculiar features of this burst are a flat earlytime optical light curve and a
radiotoXray rebrightening around 3 days after the burst. The latter is
followed by a steep opticaltoXray decay and a much shallower radio fading.
We analyze GRB 141121A in the context of the standard forwardreverse shock
(FS,RS) scenario and we disentangle the FS and RS contributions. Finally, we
comment on the puzzling earlytime (t ~3 d) behavior of GRB 141121A, and
suggest that its interpretation may require a twocomponent jet model. Overall,
our analysis confirms that the class of ULGRBs represents our best opportunity
to firmly establish the prominent emission mechanisms in action during powerful
GRB explosions, and future missions (like SVOM, XTiDE, or ISSLobster) will
provide many more of such objects.

Gammaray bursts (GRBs) are the most luminous explosions in the universe, yet
the nature and physical properties of their energy sources are far from
understood. Very important clues, however, can be inferred by studying the
afterglows of these events. We present optical and Xray observations of GRB
130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and
GTC. This burst shows a steep drop in the Xray lightcurve at $\simeq 10^5$ s
after the trigger, with a powerlaw decay index of $\alpha \sim 6$. Such a rare
behaviour cannot be explained by the standard forward shock (FS) model and
indicates that the emission, up to the fast decay at $10^5$ s, must be of
"internal origin", produced by a dissipation process within an
ultrarelativistic outflow. We propose that the source of such an outflow, which
must produce the Xray flux for $\simeq 1$ day in the cosmological rest frame,
is a newly born magnetar or black hole. After the drop, the faint Xray
afterglow continues with a much shallower decay. The optical emission, on the
other hand, shows no break across the Xray steep decrease, and the latetime
decays of both the Xray and optical are consistent. Using both the Xray and
optical data, we show that the emission after $\simeq 10^5$ s can be explained
well by the FS model. We model our data to derive the kinetic energy of the
ejecta and thus measure the efficiency of the central engine of a GRB with
emission of internal origin visible for a long time. Furthermore, we break down
the energy budget of this GRB into the prompt emission, the late internal
dissipation, the kinetic energy of the relativistic ejecta, and compare it with
the energy of the associated supernova, SN 2013fu.

We study the Quantum Electrodynamics of 2D and 3D Dirac semimetals by means
of a selfconsistent resolution of the SchwingerDyson equations, aiming to
obtain the respective phase diagrams in terms of the relative strength of the
Coulomb interaction and the number N of Dirac fermions. In this framework, 2D
Dirac semimetals have just a strongcoupling instability characterized by
exciton condensation (and dynamical generation of mass) that we find at a
critical coupling well above previous theoretical estimates, thus explaining
the absence of that instability in freestanding graphene samples. On the other
hand, we show that 3D Dirac semimetals have a richer phase diagram, with a
strongcoupling instability leading to dynamical mass generation up to N = 4
and a line of critical points for larger values of N characterized by the
vanishing of the electron quasiparticle weight in the lowenergy limit. Such a
critical behavior signals the transition to a strongly correlated liquid,
characterized by noninteger scaling dimensions that imply the absence of a pole
in the electron propagator and are the signature of nonFermi liquid behavior
with no stable electron quasiparticles.

We investigate the strongcoupling phases that may arise in 3D Dirac and Weyl
semimetals under the effect of the longrange Coulomb interaction, considering
the manybody theory of these electron systems as a variant of the conventional
fully relativistic Quantum Electrodynamics (QED). For this purpose, we apply
two different nonperturbative approaches, consisting in the sum of ladder
diagrams and taking the limit of a large number N of fermion flavors. We
benefit from the renormalizability that the theory shows in both cases to
compute the anomalous scaling dimensions of different operators exclusively in
terms of the renormalized coupling constant, allowing us to determine the
precise location of the singularities signaling the onset of the
strongcoupling phases. We show then that the QED of 3D Dirac semimetals has
two competing effects at strong coupling. One of them is the tendency to chiral
symmetry breaking and dynamical mass generation, which are analogous to the
same phenomena arising in the conventional QED at strong coupling. This trend
is however outweighed by the strong suppression of electron quasiparticles that
takes place at large N, leading to a different type of critical point at
sufficiently large interaction strength, shared also by the 3D Weyl semimetals.
Overall, the phase diagram of the 3D Dirac semimetals turns out to be richer
than that of their 2D counterparts, displaying a transition to a phase with
nonFermi liquid behavior which may be observed in materials hosting a
sufficiently large number of Dirac or Weyl points.

The third catalog of active galactic nuclei (AGNs) detected by the FermiLAT
(3LAC) is presented. It is based on the third FermiLAT catalog (3FGL) of
sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater
than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs
located at high Galactic latitudes (b>10{\deg}), a 71% increase over the
second catalog based on 2 years of data. There are 28 duplicate associations,
thus 1563 of the 2192 highlatitude gammaray sources of the 3FGL catalog are
AGNs. Most of them (98%) are blazars. About half of the newly detected blazars
are of unknown type, i.e., they lack spectroscopic information of sufficient
quality to determine the strength of their emission lines. Based on their
gammaray spectral properties, these sources are evenly split between
flatspectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected
BL~Lacs are of the highsynchrotronpeaked (HSP) type. About 50% of the BL~Lacs
have no measured redshifts. A few new rare outliers (HSPFSRQs and
highluminosity HSP BL~Lacs) are reported. The general properties of the 3LAC
sample confirm previous findings from earlier catalogs. The fraction of 3LAC
blazars in the total population of blazars listed in BZCAT remains
nonnegligible even at the faint ends of the BZCATblazar radio, optical and
Xray flux distributions, which is a clue that even the faintest known blazars
could eventually shine in gamma rays at LATdetection levels. The energyflux
distributions of the different blazar populations are in good agreement with
extrapolation from earlier catalogs.

The microwave behavior of polymer metacomposites containing parallel Febased
and continuous/shortcut Cobased microwire arrays has been investigated. A
magnetic fieldtunable metacomposite feature has been identified in the dense
continuous hybrid composite confirmed by the transmission windows in the
frequency band of 1 to 3.5 GHz. The complex magnetically tuned
redshiftblueshift evolution of the transmission window is reasoned to result
from the competition between the dynamic wirewire interaction and the
ferromagnetic resonance of Febased wires. Increasing Cobased interwire
spacing to 10 mm in the continuous hybrid composites, a remarkable dualband
transmission window in the 1.53.5 GHz and 917 GHz is respectively induced by
the ferromagnetic resonance of Febased wires and the magnetic resonance
arising between FeCo wire couples. The hybridization of parallel Febased and
shortcut Cobased wires in the polymer composite leads to a significant
enhancement of the transmission window in the frequency band of 1 to 6 GHz due
to the bandstop nature of Cobased wires. The advanced hybridized microwire
metacomposites are arguably demonstrated to be particularly attractive for
microwave cloaking and radio frequency barcoding applications.

We present optical and nearinfrared (NIR) photometry of 28 gammaray bursts
(GRBs) detected by the \textit{Swift} satellite and rapidly observed by the
Reionization and Transients Infrared/Optical (RATIR) camera. We compare the
optical flux at fiducial times of 5.5 and 11 hours after the highenergy
trigger to that in the Xray regime to quantify optical darkness. 46$\pm$9 per
cent (13/28) of all bursts in our sample and 55$\pm$10 per cent (13/26) of long
GRBs are optically dark, which is statistically consistently with previous
studies. Fitting RATIR optical and NIR spectral energy distributions (SEDs) of
19 GRBs, most (6/7) optically dark GRBs either occur at highredshift ($z>4.5$)
or have a high dust content in their host galaxies ($A_{\rm V} > 0.3$).
Performing KS tests, we compare the RATIR sample to those previously presented
in the literature, finding our distributions of redshift, optical darkness,
host dust extinction and Xray derived column density to be consistent. The one
reported discrepancy is with host galaxy dust content in the BAT6 sample, which
appears inconsistent with our sample and other previous literature. Comparing
Xray derived host galaxy hydrogen column densities to host galaxy dust
extinction, we find that GRBs tend to occur in host galaxies with a higher
metaltodust ratio than our own Galaxy, more akin to the Large and Small
Magellanic Clouds. Finally, to mitigate time evolution of optical darkness, we
measure $\beta_{\rm OX,rest}$ at a fixed rest frame time, $t_{\rm rest}=1.5$
hours and fixed rest frame energies in the Xray and optical regimes. Choosing
to evaluate optical flux at $\lambda_{\rm rest}=0.25~\mu$m, we remove
highredshift as a source of optical darkness, demonstrating that optical
darkness must result from either highredshift, dust content in the host galaxy
along the GRB sight line, or a combination of the two.

The study of the Gmode pressure coefficients of carbon nanotubes, reflecting
the stiff sp2 bond pressure dependence, is essential to the understanding of
their extraordinary mechanical properties as well as fundamental mechanics.
However, it is hindered by the availability of carbon nanotubes samples only as
bundles or isolated with surfactants. Octadecylamine functionalized carbon
nanotubes are mostly of a single diameter and can be stably dispersed in 1,
2dichloroethane and chloroform without surfactants. Here we perform high
pressure Raman spectroscopy on these tubes and obtain their experimental Gmode
pressure coefficients for individual tubes and bundles. The Gplus pressure
coefficient for bundles is only about half of that for individual tubes in 1,
2dichloroethane and is about twothirds in chloroform. The Gminus pressure
coefficient for bundles is about onethird of Gplus in 1, 2dichloroethane and
about the same in chloroform. These results for the first time provide
unambiguous experimental evidence of the significant effect of bundling on
carbon nanotubes' Gmode pressure coefficients, identifying it as one of the
major reasons for the lack of consensus on what the values should be in the
literature.

We measured the sizes and morphological parameters of LARS galaxies in the
continuum, Lya, and Ha images. We studied morphology by using the Gini
coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated
LARS galaxies at z~2 and 5.7, performing the same morphological measurements.
We also investigated the detectability of LARS galaxies in current deep field
observations. The subsample of LAEs within LARS (LARSLAEs) was stacked to
provide a comparison to stacking studies performed at high redshift. LARS
galaxies have continuum size, stellar mass, and restframe absolute magnitude
typical of Lyman break analogues in the local Universe and also similar to
2<z<3 starforming galaxies and massive LAEs. LARS optical morphology is
consistent with the one of merging systems, and irregular or starburst
galaxies. For the first time we quantify the morphology in Lya images: even if
a variety of intrinsic conditions of the interstellar medium can favour the
escape of Lya photons, LARSLAEs appear small in the continuum, and their Lya
is compact. LARS galaxies tend to be more extended in Lya than in the
restframe UV. It means that Lya photons escape by forming haloes around HII
regions of LARS galaxies. The stack of LARSLAE Lya images is peaked in the
centre, indicating that the conditions, which make a galaxy an LAE, tend to
produce a concentrated surface brightness profile. On the other hand, the stack
of all LARS galaxies is shallower and more extended. This can be caused by the
variety of dust and HI amount and distribution, which produces a more complex,
patchy, and extended profile, like the one observed for Lyman break galaxies
that can contribute to the stack. We cannot identify a single morphological
property that controls whether a galaxy emits a net positive Lya flux. However,
the LARSLAEs have continuum properties consistent with merging systems.

We study the different phases in the Quantum Electrodynamics of 3D Dirac
semimetals depending on the number $N$ of Dirac fermions, using renormalization
group methods and the selfconsistent resolution of the SchwingerDyson
equation. We find that, for $N < 4$, a phase with dynamical generation of mass
prevails at sufficiently strong coupling, sharing the same physics of the
excitonic instability in 2D Dirac semimetals. For $N \geq 4$, we show that the
phase diagram has instead a line of critical points characterized by the
suppression of the quasiparticle weight at low energies, making the system to
fall into the class of marginal Fermi liquids. Such a boundary marks the
transition to a kind of strange metal which can be still defined in terms of
electron quasiparticles, but with parameters that have large imaginary parts
implying an increasing deviation from the conventional Fermi liquid picture.

A quantum field theory approach is applied to investigate the dynamics of
flexural phonons in a metallic membrane like graphene, looking for the effects
deriving from the strong interaction between the electronic excitations and
elastic deformations. Relying on a selfconsistent screening approximation to
the phonon selfenergy, we show that the theory has a critical point
characterized by the vanishing of the effective bending rigidity of the
membrane at low momentum. We also check that the instability in the sector of
flexural phonons takes place without the development of an inplane static
distortion, which is avoided due to the significant reduction of the
electronphonon couplings for inplane phonons at large momenta. Furthermore,
we analyze the scaling properties of the manybody theory to identify the order
parameter that opens up at the point of the transition. We find that the
vanishing of the effective bending rigidity and the onset of a nonzero
expectation value of the mean curvature are concurrent manifestations of the
critical behavior. The results presented here imply that, even in the absence
of tension, the theory has a critical point at which the flat geometry becomes
an unstable configuration of the metallic membrane, with a condensation of the
mean curvature field that may well reproduce the smooth distribution of ripples
observed in freestanding graphene.

We present a template fitting algorithm for determining photometric
redshifts, $z_{\rm phot}$, of candidate highredshift gammaray bursts (GRBs).
Using afterglow photometry, obtained by the Reionization And Transients
InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB
afterglow spectral energy distribution (SED), host dust extinction and the
effect of neutral hydrogen (local and cosmological) along the line of sight. We
present the results obtained by this algorithm and RATIR photometry of GRB
130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for
models of several host dust extinction laws (none, MW, LMC and SMC), consistent
with spectroscopic measurements of the redshift of this GRB. Using simulated
RATIR photometry, we find our algorithm provides precise measures of $z_{\rm
phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$
and can robustly determine when $z_{\rm phot}>4$. Further testing highlights
the required caution in cases of highly dust extincted host galaxies. These
tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$
when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to
rapidly identify all potential highredshift events.

We investigated the microwave properties of polymer based metacomposites
containing hybridized parallel Fe and Cobased microwire arrays. A dualband
lefthanded feature was observed in the frequency bands of 1.5 to 5.5 GHz and 9
to 17 GHz, indicated by two transmission windows associated with ferromagnetic
resonance of Febased microwires and long range dipolar resonance between the
wire arrays. The plasma frequency after hybridization is significantly
increased due to the enhanced effective diameter through the wirewire
interactions between the Fe and Co microwire couples. These results offer
essential perspectives in designing the multiband metamaterial for microwave
applications such as sensors and cloaking devices.

The microwave properties of glassfibers reinforced polymer composite
embedded with an orthogonal array of Fe77Si10B10C3 microwires have been
investigated. The composites containing orthogonal wire arrays display a
remarkable transmission window in the frequency band of 1 to 6 GHz under zero
external magnetic field indicating an intrinsic doublenegativeindex
characteristic. The polymer matrices have proved to exert a synergistic effect
on the microwave properties, which is responsible for the disappearance of the
transmission windows when Ek is perpendicular to the glass fiber direction. The
plasma frequency of the orthogonal microwire array composite is higher than
that of the parallel microwire array with identical wire spacing; this could be
attributed to the enhanced microwirewave interaction induced by the axial
electrical components on the additional layer of perpendicular wires. All these
features make this new kind of orthogonal microwire composites promising for
potential cloaking and sensing applications.

We investigate the manybody instabilities of electrons interacting near Van
Hove singularities arising in monolayer and twisted bilayer graphene. We show
that a pairing instability must be dominant over the tendency to magnetic order
as the Fermi level is tuned to the Van Hove singularity in the conduction band
of graphene. As a result of the extended character of the saddle points in the
dispersion, we find that the pairing of the electrons takes place
preferentially in a channel of fwave symmetry, with an order parameter
vanishing at the position of the saddle points along the Fermi line. In the
case of the twisted bilayers, the dispersion has instead its symmetry reduced
down to the C_{3v} group and, most importantly, it leads to susceptibilities
that diverge at the saddle points but are integrable along the Fermi line. This
implies that a ferromagnetic instability becomes dominant in the twisted
graphene bilayers near the Van Hove singularity, with a strength which is
amplified as the lowest subband of the electron system becomes flatter for
decreasing twist angle.

We show that, at sufficiently large strength of the longrange Coulomb
interaction, a mass term breaking parity (socalled Haldane mass) is
dynamically generated in the manybody theory of Dirac fermions describing the
graphene layer. While the tendency towards chiral symmetry breaking is stronger
than for the dynamical breakdown of parity at spatial dimension D > 2, we find
that the situation is reversed at D = 2. The need to regularize the manybody
theory in a gaugeinvariant manner (taking the limit D = 2  \epsilon) is what
leads to the dominance of the paritybreaking pattern in graphene. We compute
the critical coupling for the generation of a paritybreaking mass from the
finite radius of convergence of the ladder series supplemented with electron
selfenergy corrections, finding a value quite close to the effective
interaction strength for graphene in vacuum after including Fermi velocity
renormalization and static RPA screening of the Coulomb interaction.

We study the manybody theory of graphene Dirac quasiparticles interacting
via the longrange Coulomb potential, taking as a starting point the ladder
approximation to different vertex functions. We test in this way the lowenergy
behavior of the electron system beyond the simple logarithmic dependence of
electronic correlators on the highenergy cutoff, which is characteristic of
the largeN approximation. We show that the graphene manybody theory is
perfectly renormalizable in the ladder approximation, as all higher powers in
the cutoff dependence can be absorbed into the redefinition of a finite number
of parameters (namely, the Fermi velocity and the weight of the fields) that
remain free of infrared divergences even at the charge neutrality point. We
illustrate this fact in the case of the vertex for the current density, where a
complete cancellation between the cutoff dependences of vertex and electron
selfenergy corrections becomes crucial for the preservation of the gauge
invariance of the theory. The other potentially divergent vertex corresponds to
the staggered (sublattice odd) charge density, which is made cutoff independent
by a redefinition in the scale of the density operator. This allows to compute
a welldefined, scale invariant anomalous dimension to all orders in the ladder
series, which becomes singular at a value of the interaction strength marking
the onset of chiral symmetry breaking (and gap opening) in the Dirac field
theory. The critical coupling we obtain in this way matches with great accuracy
the value found with a quite different method, based on the resolution of the
gap equation, thus reassuring the predictability of our renormalization
approach.