
We derive a description of the family of canonical selfadjoint extensions of
the operator of multiplication in a de Branges space in terms of singular
rankone perturbations using distinguished elements from the set of functions
associated with a de Branges space. The scale of rigged Hilbert spaces
associated with this construction is also studied from the viewpoint of de
Branges's theory.

We present twodimensional particleincell (PIC) simulations of the fully
kinetic collisionless magnetorotational instability (MRI) in weakly magnetized
(high $\beta$) pair plasma. The central result of this numerical analysis is
the emergence of a selfinduced turbulent regime in the saturation state of the
collisionless MRI, which can only be captured for large enough simulation
domains. One of the underlying mechanisms for the development of this turbulent
state is the driftkink instability (DKI) of the current sheets resulting from
the nonlinear evolution of the channel modes. The onset of the DKI can only be
observed for simulation domain sizes exceeding several linear MRI wavelengths.
The DKI, together with ensuing magnetic reconnection, activate the turbulent
motion of the plasma in the late stage of the nonlinear evolution of the MRI.
At steady state, the magnetic energy has an MHDlike spectrum with a slope of
$k^{5/3}$ for $k\rho<1$ and $k^{3}$ for subLarmor scale ($k\rho>1$). We also
examine the role of the collisionless MRI and associated magnetic reconnection
in the development of pressure anisotropy. We study the stability of the system
due to this pressure anisotropy, observing the development of mirror
instability during the earlystage of the MRI. We further discuss the
importance of magnetic reconnection for particle acceleration during the
turbulence regime. In particular, consistent with reconnection studies, we show
that at late times the kinetic energy presents a characteristic slope of
$\epsilon^{2}$ in the highenergy region.

The classical results on oversampling and undersampling (or aliasing), of
functions in PaleyWiener spaces, are generalized to de Branges spaces arising
from regular Schr\"odinger operators with a wide range of potentials.

This work deals with decay bounds for Green matrices and generalized
eigenvectors of block Jacobi matrices when the real part of the spectral
parameter lies in an infinite gap of the operator's essential spectrum. We
consider the cases of commutative and noncommutative matrix entries separately.
An example of a block Jacobi operator with noncommutative entries and
nonnegative essential spectrum is given to illustrate the results.

This work deals with the functional model for extensions of symmetric
operators and its applications to the theory of wave scattering. In terms of
Boris Pavlov's spectral form of this model, we find explicit formulae for the
action of the unitary group of exponentials corresponding to almost solvable
extensions of a given closed symmetric operator with equal deficiency indices.
On the basis of these formulae, we are able to construct wave operators and
derive a new representation for the scattering matrix for pairs of such
extensions in both selfadjoint and nonselfadjoint situations.

On the basis of the explicit formulae for the action of the unitary group of
exponentials corresponding to almost solvable extensions of a given closed
symmetric operator with equal deficiency indices, we derive a new
representation for the scattering matrix for pairs of such extensions. We use
this representation to explicitly recover the coupling constants in the inverse
scattering problem for a finite noncompact quantum graph with $\delta$type
vertex conditions.

This work studies the interplay between Green functions, the index of
determinacy of spectral measures and interior finite rank perturbations of
Jacobi operators. The index of determinacy quantifies the stability of
uniqueness of solutions of the moment problem. We give results on the constancy
of this index in terms of perturbations of the corresponding Jacobi operators.
The permanence of the $N$extremality of a measure is also studied. A measure
$\mu$ is $N$extremal when the polynomials are dense in $L_2(\mathbb{R},\mu)$.
As a byproduct, we give a characterization of the index in terms of cyclic
vectors. We consider a new inverse problem for Jacobi operators in which
information on the place where the interior perturbation occurs is obtained
from the index of determinacy.

In this paper, we investigate the evolution of the energy spread and the
divergence of electron beams while they interact with different laser pulses at
intensities where quantum effects and radiation reaction are of relevance. The
interaction is modeled with a QEDPIC code and the results are compared with
those obtained with a standard PIC code with the addition of a classical
radiation reaction module and with theoretical predictions. While classical
radiation reaction is a continuous process, in QED, radiation emission is
stochastic. The two pictures reconcile in the limit when the emitted photons
energy is small compared to the energy of the emitting electrons. The energy
spread of the electron distribution function always tends to decrease with
classical radiation reaction, whereas the stochastic QED emission can also
enlarge it. These two tendencies compete in the QEDdominated regime. Our
analysis, supported by the QED module, reveals an upper limit to the maximal
attainable energy spread due to stochasticity that depends on laser intensity
and the electron beam average energy. Beyond this limit, the energy spread
decreases. These findings are verified for different laser pulse lengths
ranging from short ~ 30 fs pulses presently available to the long ~ 150 fs
pulses expected in the nearfuture laser facilities, and compared with a
theoretical model. Our results also show that near future experiments will be
able to probe this transition and to demonstrate the competition between
enhanced QED induced energy spread and energy spectrum narrowing from classical
radiation reaction.

We present an analytical and numerical study of multiplelaser QED cascades
induced with linearly polarised laser pulses. We analyse different polarisation
orientations and propose a configuration that maximises the cascade
multiplicity and favours the laser absorption. We generalise the analytical
estimate for the cascade growth rate previously calculated in the field of two
colliding linearly polarised laser pulses and account for multiple laser
interaction. The estimate is verified by a comprehensive numerical study of
fourlaser QED cascades across a range of different laser intensities with QED
PIC module of OSIRIS. We show that by using four linearly polarised 30 fs laser
pulses, one can convert more than 50 % of the total energy to gammarays
already at laser intensity $I\simeq10^{24}\ \mathrm{W/cm^2}$. In this
configuration, the laser conversion efficiency is higher compared with the case
with two colliding lasers.

Under the presence of ultra high intensity lasers or other intense
electromagnetic fields the motion of particles in the ultrarelativistic regime
can be severely affected by radiation reaction. The standard particleincell
(PIC) algorithms do not include radiation reaction effects. Even though this is
a well known mechanism, there is not yet a definite algorithm nor a standard
technique to include radiation reaction in PIC codes. We have compared several
models for the calculation of the radiation reaction force, with the goal of
implementing an algorithm for classical radiation reaction in the Osiris
framework, a stateoftheart PIC code. The results of the different models are
compared with standard analytical results, and the relevance/advantages of each
model are discussed. Numerical issues relevant to PIC codes such as resolution
requirements, application of radiation reaction to macro particles and
computational cost are also addressed. The Landau and Lifshitz reduced model is
chosen for implementation.

This work gives results on the interplay of the spectra of two Jacobi
operators corresponding to an infinite massspring system and a modification of
it obtained by changing one mass and one spring of the system. It is shown that
the system can be recovered from these two spectra. Necessary and sufficient
conditions for two sequences to be the spectra of the massspring system and
the perturbed one are provided.

In this note, we solve an inverse spectral problem for a class of finite band
symmetric matrices. We provide necessary and sufficient conditions for a matrix
valued function to be a spectral function of the operator corresponding to a
matrix in our class and give an algorithm for recovering this matrix from the
spectral function. The reconstructive algorithm is applicable to matrices which
cannot be treated by known inverse block matrix methods. Our approach to the
inverse problem is based on the rational interpolation theory developed in a
previous paper.

Laser wakefield accelerators have great potential as the basis for next
generation compact radiation sources because their accelerating gradients are
three orders of magnitude larger than traditional accelerators. However, Xray
radiation from such devices still lacks of tunability, especially the intensity
and polarization distribution. Here we propose a tunable polarized radiation
source from a helical plasma undulator based on plasma channel guided wakefield
accelerator. When a laser pulse is initially incident with a skew angle
relative to the channel axis, the laser and accelerated electrons experience
collective spiral motions, which leads to elliptically polarized
synchrotronlike radiation with flexible tunability on radiation intensity,
spectra and polarization. We demonstrate that a radiation source with
millimeter size and peak brilliance of $2\times10^{19}
photons/s/mm^{2}/mrad^{2}/0.1%$ bandwidth can be made with moderate laser and
electron beam parameters. This brilliance is comparable with the third
generation of synchrotron radiation facilities running at similar photon
energies, suggesting that laser plasma based radiation sources are promising
for advanced applications.

A model for laser light absorption in electronpositron plasmas
selfconsistently created via QED cascades is described. The laser energy is
mainly absorbed due to hard photon emission via nonlinear Compton scattering.
The degree of absorption depends on the laser intensity and the pulse duration.
The QED cascades are studied with multidimensional particleincell
simulations complemented by a QED module and a macroparticle merging algorithm
that allows to handle the exponential growth of the number of particles.
Results range from moderateintensity regimes ($\sim$ 10 PW) where the laser
absorption is negligible, to extreme intensities (> 100 PW) where the degree of
absorption reaches 80%. Our study demonstrates good agreement between the
analytical model and simulations. The expected properties of the hard photon
emission and the generated pairplasma are investigated, and the experimental
signatures for nearfuture laser facilities are discussed.

This note deals with the direct and inverse spectral analysis for a class of
infinite band symmetric matrices. This class corresponds to operators arising
from difference quations with usual and inner boundary conditions. We give a
characterization of the spectral functions for the operators and provide
necessary and sufficient conditions for a matrixvalued function to be a
spectral function of the operators. Additionally, we give an algorithm for
recovering the matrix from the spectral function. The approach to the inverse
problem is based on the rational interpolation theory.

When modeling laser wakefield acceleration (LWFA) using the particleincell
(PIC) algorithm in a Lorentz boosted frame, the plasma is drifting
relativistically at $\beta_b c$ towards the laser, which can lead to a
computational speedup of $\sim \gamma_b^2=(1\beta_b^2)^{1}$. Meanwhile, when
LWFA is modeled in the quasi3D geometry in which the electromagnetic fields
and current are decomposed into a limited number of azimuthal harmonics,
speedups are achieved by modeling three dimensional problems with the
computation load on the order of two dimensional $rz$ simulations. Here, we
describe how to combine the speed ups from the Lorentz boosted frame and
quasi3D algorithms. The key to the combination is the use of a hybrid YeeFFT
solver in the quasi3D geometry that can be used to effectively eliminate the
Numerical Cerenkov Instability (NCI) that inevitably arises in a Lorentz
boosted frame due to the unphysical coupling of Langmuir modes and EM modes of
the relativistically drifting plasma in these simulations. In addition, based
on the spacetime distribution of the LWFA data in the lab and boosted frame,
we propose to use a moving window to follow the drifting plasma to further
reduce the computational load. We describe the details of how the NCI is
eliminated for the quasi3D geometry, the setups for simulations which combine
the Lorentz boosted frame and quasi3D geometry, the use of a moving window,
and compare the results from these simulations against their corresponding lab
frame cases. Good agreement is obtained, particularly when there is no
selftrapping, which demonstrates it is possible to combine the Lorentz boosted
frame and the quasi3D algorithms when modeling LWFA to achieve unprecedented
speedups.

This paper presents an erroneous proof that if the polynomials are dense in
$L_2(\mathbb{R}, \rho)$, then they are dense in $L_2(\mathbb{R}, \rho+\mu)$
where $\mu$ is a measure supported on a finite set of points.

This work provides a complete characterization of the solutions of a linear
interpolation problem for vector polynomials. The interpolation problem
consists in finding n scalar polynomials such that an equation involving a
linear combination of them is satisfied for each one of the N interpolation
nodes. The results of this work generalize previous results on the socalled
rational interpolation and have applications to direct and inverse spectral
analysis of band matrices.

This work presents a contemporary treatment of Krein's entire operators with
deficiency indices $(1,1)$ and de Branges' Hilbert spaces of entire functions.
Each of these theories played a central role in the research of both renown
mathematicians. Remarkably, entire operators and de Branges spaces are
intimately connected and the interplay between them has had an impact in both
spectral theory and the theory of functions. This work exhibits the
interrelation between Krein's and de Branges' theories by means of a functional
model and discusses recent developments, giving illustrations of the main
objects and applications to the spectral theory of difference and differential
operators.

We analyse the full shock formation process in electronion plasmas in theory
and simulations. It is accepted that electromagnetic shocks in initially
unmagnetised relativistic plasmas are triggered by the filamentation
instability. However, the transition from the first unstable phase to the
quasisteady shock is still missing. We derive a theoretical model for the
shock formation time, taking into account the filament merging in the
nonlinear phase of the filamentation instability. This process is much slower
than in electronpositron pair shocks, so that the shock formation is longer by
a factor proportional to sqrt(m_i/m_e) ln(m_i/m_e).

We investigate the connections between WeylTitchmarshKodaira theory for
onedimensional Schr\"odinger operators and the theory of $n$entire operators.
As our main result we find a necessary and sufficient condition for a
onedimensional Schr\"odinger operator to be $n$entire in terms of square
integrability of derivatives (w.r.t. the spectral parameter) of the Weyl
solution. We also show that this is equivalent to the Weyl function being in a
generalized HerglotzNevanlinna class. As an application we show that perturbed
Bessel operators are $n$entire, improving the previously known conditions on
the perturbation.

A hybrid Maxwell solver for fully relativistic and electromagnetic (EM)
particleincell (PIC) codes is described. In this solver, the EM fields are
solved in $k$ space by performing an FFT in one direction, while using finite
difference operators in the other direction(s). This solver eliminates the
numerical Cerenkov radiation for particles moving in the preferred direction.
Moreover, the numerical Cerenkov instability (NCI) induced by the
relativistically drifting plasma and beam can be eliminated using this hybrid
solver by applying strategies that are similar to those recently developed for
pure FFT solvers. A current correction is applied for the charge conserving
current deposit to correctly account for the EM calculation in hybrid YeeFFT
solver. A theoretical analysis of the dispersion properties in vacuum and in a
drifting plasma for the hybrid solver is presented, and compared with PIC
simulations with good agreement obtained. This hybrid solver is applied to both
2D and 3D Cartesian and quasi3D (in which the fields and current are
decomposed into azimuthal harmonics) geometries. Illustrative results for laser
wakefield accelerator simulation in a Lorentz boosted frame using the hybrid
solver in the 2D Cartesian geometry are presented, and compared against results
from 2D UPICEMMA simulation which uses a pure spectral Maxwell solver, and
from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good
agreement is obtained which demonstrates the feasibility of using the hybrid
solver for high fidelity simulation of relativistically drifting plasma with no
evidence of the numerical Cerenkov instability.

Particleincell merging algorithms aim to resample dynamically the
sixdimensional phase space occupied by particles without distorting
substantially the physical description of the system. Whereas various
approaches have been proposed in previous works, none of them seemed to be able
to conserve fully charge, momentum, energy and their associated distributions.
We describe here an alternative algorithm based on the coalescence of N massive
or massless particles, considered to be close enough in phase space, into two
new macroparticles. The local conservation of charge, momentum and energy are
ensured by the resolution of a system of scalar equations. Various simulation
comparisons have been carried out with and without the merging algorithm, from
classical plasma physics problems to extreme scenarios where quantum
electrodynamics is taken into account, showing in addition to the conservation
of local quantities, the good reproducibility of the particle distributions. In
case where the number of particles ought to increase exponentially in the
simulation box, the dynamical merging permits a considerable speedup, and
significant memory savings that otherwise would make the simulations impossible
to perform.

We study how the spectrum of a Jacobi operator changes when this operator is
modified by a certain finite rank perturbation. The operator corresponds to an
infinite massspring system and the perturbation is obtained by modifying one
interior mass and one spring of this system. In particular, there are detailed
results of what happens in the spectral gaps and which eigenvalues do not move
under the modifications considered. These results were obtained by a new
tecnique of comparative spectral analysis and they generalize and include
previous results for finite and infinite Jacobi matrices.

In this work we consider de Branges spaces where the multiplication operator
by the independent variable is not densely defined. First, we study the
canonical selfadjoint extensions of the multiplication operator as a family of
rankone perturbations from the viewpoint of the theory of de Branges spaces.
Then, on the basis of the obtained results, we provide new necessary and
sufficient conditions for a real, zerofree function to lie in a de Branges
space.