• 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 rank-one 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 two-dimensional particle-in-cell (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 self-induced 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 drift-kink 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 MHD-like spectrum with a slope of $k^{-5/3}$ for $k\rho<1$ and $k^{-3}$ for sub-Larmor 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 early-stage 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 high-energy region.
  • The classical results on oversampling and undersampling (or aliasing), of functions in Paley-Wiener 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 self-adjoint and non-self-adjoint 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 non-compact 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 by-product, 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 QED-PIC 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 QED-dominated 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 near-future 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 multiple-laser 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 four-laser 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 gamma-rays 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 particle-in-cell (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 state-of-the-art 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 mass-spring 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 mass-spring 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, X-ray 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 synchrotron-like 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 electron-positron plasmas self-consistently 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 multi-dimensional particle-in-cell simulations complemented by a QED module and a macro-particle merging algorithm that allows to handle the exponential growth of the number of particles. Results range from moderate-intensity 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 pair-plasma are investigated, and the experimental signatures for near-future 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 matrix-valued 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 particle-in-cell (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 quasi-3D 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 $r-z$ simulations. Here, we describe how to combine the speed ups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D 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 space-time 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 quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame and quasi-3D 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 self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D 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 so-called 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 electron-ion 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 quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking into account the filament merging in the non-linear phase of the filamentation instability. This process is much slower than in electron-positron 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 Weyl-Titchmarsh-Kodaira theory for one-dimensional Schr\"odinger operators and the theory of $n$-entire operators. As our main result we find a necessary and sufficient condition for a one-dimensional 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 Herglotz-Nevanlinna 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) particle-in-cell (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 Yee-FFT 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 quasi-3D (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 UPIC-EMMA 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.
  • Particle-in-cell merging algorithms aim to resample dynamically the six-dimensional 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 macro-particles. 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 mass-spring 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 rank-one 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, zero-free function to lie in a de Branges space.