title: A novel least squares method for Helmholtz equations with large wave numbers

abstract: In this paper we are concerned with numerical methods for nonhomogeneous Helmholtz equations in inhomogeneous media. We design a least squares method for discretization of the considered Helmholtz equations. In this method, an auxiliary unknown is introduced on the common interface of any two neighboring elements and a quadratic subject functional is defined by the jumps of the traces of the solutions of local Helmholtz equations across all the common interfaces, where the local Helmholtz equations are defined on elements and are imposed Robin-type boundary conditions given by the auxiliary unknowns. A minimization problem with the subject functional is proposed to determine the auxiliary unknowns. The resulting discrete system of the auxiliary unknowns is Hermitian positive definite and so it can be solved by the PCG method. Under some assumptions we show that the generated approximate solutions possess almost the optimal error estimates with little "wave number pollution". Moreover, we construct a substructuring preconditioner for the discrete system of the auxiliary unknowns. Numerical experiments show that the proposed methods are very effective for the tested Helmholtz equations with large wave numbers.


title: Stationary solutions to a chemotaxis-consumption model with realistic boundary conditions

abstract: Previous studies of chemotaxis models with consumption of the chemoattractant (with or without fluid) have not been successful in explaining pattern formation even in the simplest form of concentration near the boundary, which had been experimentally observed. Following the suggestions that the main reason for that is usage of inappropriate boundary conditions, in this article we study solutions to the stationary chemotaxis system \[ \begin{cases} 0 = \Delta n - \nabla\cdot(n\nabla c) \\ 0 = \Delta c - nc \end{cases} \] in bounded domains $\Omega\subset\mathbb{R}^N$, $N\ge 1$, under no-flux boundary conditions for $n$ and the physically meaningful condition \[ \partial_{\nu} c = (\gamma-c)g \] on $c$, with given parameter $\gamma>0$ and $g\in C^{1+\beta}(\Omega)$ satisfying $g\ge 0$, $g \not\equiv 0$ on $\partial \Omega$. We prove existence and uniqueness of solutions for any given mass $\int_\Omega n > 0$. These solutions are non-constant.


title: Distributed-memory parallelization of the aggregated unfitted finite element method

abstract: The aggregated unfitted finite element method (AgFEM) is a methodology recently introduced in order to address conditioning and stability problems associated with embedded, unfitted, or extended finite element methods. The method is based on removal of basis functions associated with badly cut cells by introducing carefully designed constraints, which results in well-posed systems of linear algebraic equations, while preserving the optimal approximation order of the underlying finite element spaces. The specific goal of this work is to present the implementation and performance of the method on distributed-memory platforms aiming at the efficient solution of large-scale problems. In particular, we show that, by considering AgFEM, the resulting systems of linear algebraic equations can be effectively solved using standard algebraic multigrid preconditioners. This is in contrast with previous works that consider highly customized preconditioners in order to allow one the usage of iterative solvers in combination with unfitted techniques. Another novelty with respect to the methods available in the literature is the problem sizes that can be handled with the proposed approach. While most of previous references discussing linear solvers for unfitted methods are based on serial non-scalable algorithms, we propose a parallel distributed-memory method able to efficiently solve problems at large scales. This is demonstrated by means of a weak scaling test defined on complex 3D domains up to 300M degrees of freedom and one billion cells on 16K CPU cores in the Marenostrum-IV platform. The parallel implementation of the AgFEM method is available in the large-scale finite element package FEMPAR.


title: Extension of B-spline Material Point Method for unstructured triangular grids using Powell-Sabin splines

abstract: The Material Point Method (MPM) is a numerical technique that combines a fixed Eulerian background grid and Lagrangian point masses to simulate materials which undergo large deformations. Within the original MPM, discontinuous gradients of the piecewise-linear basis functions lead to so-called `grid-crossing errors' when particles cross element boundaries. Previous research has shown that B-spline MPM (BSMPM) is a viable alternative not only to MPM, but also to more advanced versions of the method that are designed to reduce the grid-crossing errors. In contrast to many other MPM-related methods, BSMPM has been used exclusively on structured rectangular domains, considerably limiting its range of applicability. In this paper, we present an extension of BSMPM to unstructured triangulations. The proposed approach combines MPM with $C^1$-continuous high-order Powell-Sabin (PS) spline basis functions. Numerical results demonstrate the potential of these basis functions within MPM in terms of grid-crossing-error elimination and higher-order convergence.


title: UVES Spectroscopy of T Chamaeleontis: Line Variability, Mass Accretion Rate and Spectro-astrometric Analysis

abstract: Although advances in exoplanet detection techniques have seen an increase in discoveries, observing a planet in the earliest stages of formation still remains a difficult task. Here four epochs of spectra of the transitional disk object T Cha are analysed to determine whether spectro-astrometry can be used to detect a signal from its proposed protoplanet, T Cha b. The unique properties of T Cha are also further constrained. H{\alpha} and [O I]{\lambda} 6300, the most prominent lines, were analysed using spectro-astrometry. H{\alpha} being a direct accretion tracer is the target for the T Cha b detection while [O I]{\lambda} 6300 is considered to be an indirect tracer of accretion. [O I]{\lambda} 6300 is classified as a broad low velocity component (BC LVC). The mass accretion rate was derived for all epochs using new [O I]{\lambda} 6300 LVC relationships and the H{\alpha} line luminosity. It is shown that a wind is the likely origin of the [O I]{\lambda} 6300 line and that the [O I]{\lambda} 6300 line serves as a better accretion tracer than H{\alpha} in this case. From the comparison between M acc([OI]) and M acc(H{\alpha}) it is concluded that T Cha is not an intrinsically weak accretor but rather that a significant proportion of the H{\alpha} emission tracing accretion is obscured. T Cha b is not detected in the spectro-astrometric analysis yet a detection limit of 0.5 mas is derived. The analysis in this case was hampered by spectro-astrometric artefacts and by the unique properties of T Cha. While it seems that spectro-astrometry as a means of detecting exoplanets in TDs can be challenging it can be used to put an limit on the strength of the H{\alpha} emission from accreting planetary companions and thus can have an important input into the planning of high angular resolution observations.


title: Random graphs and their subgraphs

abstract: Random graphs are more and more used for modeling real world networks such as evolutionary networks of proteins. For this purpose we look at two different models and analyze how properties like connectedness and degree distributions are inherited by differently constructed subgraphs. We also give a formula for the variance of the degrees of fixed nodes in the preferential attachment model and additionally draw a connection between weighted graphs and electrical networks.


title: A Reeb flow on the three-sphere without a disk-like global surface of section

abstract: We show that there are Reeb flows on the standard, tight three-sphere that do not admit global surfaces of section with one boundary component. In particular, the Reeb flows that we construct do not admit disk-like global surfaces of section. These Reeb flows are constructed using integrable systems, and a connected sum construction that extends the integrable system.


title: Potential Vorticity of Saturn's Polar Regions: Seasonality and Instabilities

abstract: We analyse the potential vorticity of Saturn's polar regions, as it is a fundamental dynamical tracer that enables us to improve our understanding of the dynamics of these regions and their seasonal variability. In particular, we present zonally averaged quasi-geostrophic potential vorticity maps between 68{\deg} planetographic latitude and the poles at altitudes between 500 mbar and 1mbar for three different epochs: (i) June 2013 (early northern summer) for the north polar region, (ii) December 2008 (late northern winter) for both polar regions and (iii) October 2006 (southern summer) for the south, computed using temperature profiles retrieved from Cassini Composite Infrared Spectrometer (CIRS) data and wind profiles obtained from Cassini's Imaging Science Subsystem (ISS). The results show that quasi-geostrophic potential vorticity maps are very similar at all the studied epochs, showing positive vorticities at the north and negative at the south, indicative of the dominance of the Coriolis parameter 2${\Omega}sin{\phi}$ at all latitudes, except near the pole. The meridional gradients of the quasi-geostrophic potential vorticity show that dynamical instabilities, mainly due to the barotropic term, could develop at the flanks of the Hexagon at 78{\deg}N, the jet at 73.9{\deg}S and on the equatorward flank of both polar jets. There are no differences in potential vorticity gradients between the two hemispheres that could explain why a hexagon forms in the north and not in the south. No seasonal variability of the potential vorticity and its meridional gradient has been found, despite significant changes in the atmospheric temperatures over time.


title: Existence of solution for a system involving fractional Laplacians and a Radon measure

abstract: An existence of a nontrivial solution in some `weaker' sense of the following system of equations \begin{align*} (-\Delta)^{s}u+l(x)\phi u+w(x)|u|^{k-1}u&=\mu~\text{in}~\Omega\nonumber\\ (-\Delta)^{s}\phi&= l(x)u^2~\text{in}~\Omega\nonumber\\ u=\phi&=0 ~\text{in}~\mathbb{R}^N\setminus\Omega \end{align*} has been proved. Here $s \in (0,1)$, $l,w$ are bounded nonnegative functions in $\Omega$, $\mu$ is a Radon measure and $k > 1$ belongs to a certain range.


title: Multi-wavelength torus-jet model for Sgr~A*

abstract: Context. The properties of the accretion flow surrounding the supermassive central black hole of the Galaxy, Sgr A*, will be scrutinized by the new-generation instrument GRAVITY and the Event Horizon Telescope (EHT). Developing fast, robust, and simple models of such flows is thus important and very timely. Aims. We want to model the quiescent emission of Sgr A* from radio to mid-infrared by considering a magnetized compact torus and an extended jet. Results. We find perfect spectral fit both for face-on and edge-on views. These best fits give parameters values very close to that found by the most recent numerical simulations, which are much more complex than our model. The intrinsic radio size of Sgr A* is found to be in reasonable agreement with the centimetric observed constraints. Our best-fit infrared spectral index is in perfect agreement with the latest constraints. Our emission region at 1.3 mm, although larger than the Doeleman et al. (2008) Gaussian best-fit, does contain bright features at the 40 microarcsec scale. EHT-reconstructed images show that torus/jet-specific features persist after the reconstruction procedure, and that these features are sensitive to inclination. Conclusions. The main interest of our model is to give a simple and fast model of the quiescent state of Sgr A*, which gives extremely similar results as compared to state-of-the-art numerical simulations. Our model is easy to use and we publish all the material necessary to reproduce our spectra and images, so that anyone interested can use our results rather straightforwardly. We hope that such a public tool can be useful in the context of the recent and near-future GRAVITY and EHT results. Our model can in particular be easily used to test alternative compact objects models, or alternative gravity theories. The main limitation of our model is that we do not yet treat the X-ray emission.


title: Position Representation of Effective Electron-Electron Interactions in Solids

abstract: An essential ingredient in many model Hamiltonians, such as the Hubbard model, is the effective electron-electron interaction $U$, which enters as matrix elements in some localized basis. These matrix elements provide the necessary information in the model, but the localized basis is incomplete for describing $U$. We present a systematic scheme for computing the manifestly basis-independent dynamical interaction in position representation, $U({\bf r},{\bf r}';\omega)$, and its Fourier transform to time domain, $U({\bf r},{\bf r}';\tau)$. These functions can serve as an unbiased tool for the construction of model Hamiltonians. For illustration we apply the scheme within the constrained random-phase approximation to the cuprate parent compounds La$_2$CuO$_4$ and HgBa$_2$CuO$_4$ within the commonly used 1- and 3-band models, and to non-superconducting SrVO$_{3}$ within the $t_{2g}$ model. Our method is used to investigate the shape and strength of screening channels in the compounds. We show that the O 2$p_{x,y}-$Cu 3$d_{x^2-y^2}$ screening gives rise to regions with strong attractive static interaction in the minimal (1-band) model in both cuprates. On the other hand, in the minimal ($t_{2g}$) model of SrVO$_3$ only regions with a minute attractive interaction are found. The temporal interaction exhibits generic damped oscillations in all compounds, and its time-integral is shown to be the potential caused by inserting a frozen point charge at $\tau=0$. When studying the latter within the three-band model for the cuprates, short time intervals are found to produce a negative potential.


title: Unsupervised Clinical Language Translation

abstract: As patients' access to their doctors' clinical notes becomes common, translating professional, clinical jargon to layperson-understandable language is essential to improve patient-clinician communication. Such translation yields better clinical outcomes by enhancing patients' understanding of their own health conditions, and thus improving patients' involvement in their own care. Existing research has used dictionary-based word replacement or definition insertion to approach the need. However, these methods are limited by expert curation, which is hard to scale and has trouble generalizing to unseen datasets that do not share an overlapping vocabulary. In contrast, we approach the clinical word and sentence translation problem in a completely unsupervised manner. We show that a framework using representation learning, bilingual dictionary induction and statistical machine translation yields the best precision at 10 of 0.827 on professional-to-consumer word translation, and mean opinion scores of 4.10 and 4.28 out of 5 for clinical correctness and layperson readability, respectively, on sentence translation. Our fully-unsupervised strategy overcomes the curation problem, and the clinically meaningful evaluation reduces biases from inappropriate evaluators, which are critical in clinical machine learning.


title: Improved Decoding of Staircase Codes: The Soft-aided Bit-marking (SABM) Algorithm

abstract: Staircase codes (SCCs) are typically decoded using iterative bounded-distance decoding (BDD) and hard decisions. In this paper, a novel decoding algorithm is proposed, which partially uses soft information from the channel. The proposed algorithm is based on marking certain number of highly reliable and highly unreliable bits. These marked bits are used to improve the miscorrection-detection capability of the SCC decoder and the error-correcting capability of BDD. For SCCs with $2$-error-correcting Bose-Chaudhuri-Hocquenghem component codes, our algorithm improves upon standard SCC decoding by up to $0.30$~dB at a bit-error rate (BER) of $10^{-7}$. The proposed algorithm is shown to achieve almost half of the gain achievable by an idealized decoder with this structure. A complexity analysis based on the number of additional calls to the component BDD decoder shows that the relative complexity increase is only around $4\%$ at a BER of $10^{-4}$. This additional complexity is shown to decrease as the channel quality improves. Our algorithm is also extended (with minor modifications) to product codes. The simulation results show that in this case, the algorithm offers gains of up to $0.44$~dB at a BER of $10^{-8}$.


title: Quantum paramagnetism and helimagnetic orders in the Heisenberg model on the body centered cubic lattice

abstract: We investigate the spin $S=1/2$ Heisenberg model on the body centered cubic lattice in the presence of ferromagnetic and antiferromagnetic nearest-neighbor $J_{1}$, second-neighbor $J_{2}$, and third-neighbor $J_{3}$ exchange interactions. The classical ground state phase diagram obtained by a Luttinger-Tisza analysis is shown to host six different (noncollinear) helimagnetic orders in addition to ferromagnetic, N\'eel, stripe and planar antiferromagnetic orders. Employing the pseudofermion functional renormalization group (PFFRG) method for quantum spins ($S=1/2$) we find an extended nonmagnetic region, and significant shifts to the classical phase boundaries and helimagnetic pitch vectors caused by quantum fluctuations while no new long-range dipolar magnetic orders are stabilized. The nonmagnetic phase is found to disappear for $S=1$. We calculate the magnetic ordering temperatures from PFFRG and quantum Monte Carlo methods, and make comparisons to available data


title: Single-photon light emitting diodes based on pre-selected quantum dots using a deterministic lithography technique

abstract: In the present study, we developed a fabrication process of an electrically driven single-photon LED based on InP QDs emitting in the red spectral range, the wavelength of interest coinciding with the high efficiency window of Si APDs. A deterministic lithography technique allowed for the pre-selection of a suitable QD, here exclusively operated under electrical carrier injection. The final device was characterized under micro-electroluminescence in direct current, as well as in pulsed excitation mode. In particular, under pulsed excitation of one device, single-photon emission of a spectral line, identified as an exciton, has been observed with $g^{(2)}_\mathrm{raw}(0)=0.42\pm0.02$, where the non-zero $g^{(2)}$-value is mainly caused by background contribution in the spectrum and re-excitation processes due to the electrical pulse length. The obtained results constitute an important step forward in the fabrication of electrically driven single-photon sources, where deterministic lithography techniques can be used to sensibly improve the device performances. In principle, the developed process can be extended to any desired emitter wavelength above $600\,\mathrm{nm}$ up to the telecom bands.


title: On the structure of defects in the Fe7Mo6 $\mu$-Phase

abstract: Topologically close packed phases, among them the $\mu$-phase studied here, are commonly considered as being hard and brittle due to their close packed and complex structure. Nanoindentation enables plastic deformation and therefore investigation of the structure of mobile defects in the $\mu$-phase, which, in contrast to grown-in defects, has not been examined yet. High resolution transmission electron microscopy (HR-TEM) performed on samples deformed by nanoindentation revealed stacking faults which are likely induced by plastic deformation. These defects were compared to theoretically possible stacking faults within the $\mu$-phase building blocks, and in particular Laves phase layers. The experimentally observed stacking faults were found resulting from synchroshear assumed to be associated with deformation in the Laves-phase building blocks.


title: Constructing the Matrix Multilayer Perceptron and its Application to the VAE

abstract: Like most learning algorithms, the multilayer perceptrons (MLP) is designed to learn a vector of parameters from data. However, in certain scenarios we are interested in learning structured parameters (predictions) in the form of symmetric positive definite matrices. Here, we introduce a variant of the MLP, referred to as the matrix MLP, that is specialized at learning symmetric positive definite matrices. We also present an application of the model within the context of the variational autoencoder (VAE). Our formulation of the VAE extends the vanilla formulation to the cases where the recognition and the generative networks can be from the parametric family of distributions with dense covariance matrices. Two specific examples are discussed in more detail: the dense covariance Gaussian and its generalization, the power exponential distribution. Our new developments are illustrated using both synthetic and real data.


title: ALMA FITS header keywords: a study from the archive User perspective

abstract: ALMA products are stored in the Science Archive in the form of FITS images. It is a common idea that the FITS image headers should collect in their keywords all the information that an archive User might want to search for in order to quickly select, compare, or discard datasets. With this perspective in mind, we first present a short description of the current status of the ALMA FITS archive and images. We realized that at the moment most of the parameters that could be useful for a general User are still missing in the archived data. We then provide a CASA task generating the image header keywords that we suggest to be relevant for the scientific exploitation of the ALMA archival data. The proposed tool could be also applied to several types of interferometer data and part of it is implemented in a web interface. An example of the scientific application of the keywords is also discussed.


title: Performance Analysis of Joint Radar and Communication using OFDM and OTFS

abstract: We consider a joint radar estimation and communication system using orthogonal frequency division multiplexing (OFDM) and orthogonal time frequency space (OTFS) modulations. The scenario is motivated by vehicular applications where a vehicle equipped with a mono-static radar wishes to communicate data to its target receiver, while estimating parameters of interest related to this receiver. By focusing on the case of a single target, we derive the maximum likelihood (ML) estimator and the Cram\'er-Rao lower bound on joint velocity and range estimation. Numerical examples demonstrate that both digital modulation formats can achieve as accurate range/velocity estimation as state-of-the-art radar waveforms such as frequency modulated continuous wave (FMCW) while sending digital information at their full achievable rate. We conclude that it is possible to obtain significant data transmission rate without compromising the radar estimation capabilities of the system.


title: Sex-ratio bias induced by mutation

abstract: A question in evolutionary biology is why the number of males is approximately equal to that of females in many species, and Fisher's theory of equal investment answers that it is the evolutionarily stable state. The Fisherian mechanism can be given a concrete form by a genetic model based on the following assumptions: (1) Males and females mate at random. (2) An allele acts on the father to determine the expected progeny sex ratio. (3) The offspring inherits the allele from either side of the parents with equal probability. The model is known to achieve the 1:1 sex ratio due to the invasion of mutant alleles with different progeny sex ratios. In this study, however, we argue that mutation plays a more subtle role in that fluctuations caused by mutation renormalize the sex ratio and thereby keep it away from 1:1 in general. This finding shows how the sex ratio is affected by mutation in a systematic way, whereby the effective mutation rate can be estimated from an observed sex ratio.


title: A double EP-based proposal for turbo equalization

abstract: This letter deals with the application of the expectation propagation (EP) algorithm to turbo equalization. The EP has been successfully applied to obtain either a better approximation at the output of the equalizer or at the output of the channel decoder to better initialize the Gaussian prior used by the equalizer. In this letter we combine both trends to propose a novel double EP-based equalizer that is able to decrease the number of iterations needed, reducing the computational complexity to twice that of the linear MMSE. This novel equalizer is developed in three different implementations: a block design that exploits the whole vector of observations, a Wiener filter-type approach that just uses the observations within a predefined window and a Kalman smoothing filter-type approach that emulates the BCJR behavior. Finally, we include some experimental results to compare the three different implementations and to detail their improvements with respect to other EP-based proposals in the literature.


title: The class of supernova progenitors that result from fatal common envelope evolution

abstract: I construct the class of supernovae and supernova progenitors that result from fatal common envelope evolution (CEE). The fatal CEE progenitors are stellar binary systems where a companion spirals-in inside the envelope of a giant star and merges with the core. The companion can be a neutron star (NS; or a black hole) that destroys the core and by that forms a common envelope jets supernova (CEJSN), a white dwarf (WD) that merges with the core to form a massive WD that later might explode as a Type Ia supernova (the core degenerate scenario), or a main sequence companion. In the latter case the outcome might be a core collapse supernova (CCSN) of a blue giant, a CCSN of type IIb or of type Ib. In another member of this class two giant stars merge and the two cores spiral-in toward each other to form a massive core that later explodes as a CCSN with a massive circumstellar matter (CSM). I discuss the members of this class, their characteristics, and their common properties. I find that fatal CEE events account for $\approx 6-10 \%$ of all CCSNe, and raise the possibility that a large fraction of peculiar and rare supernovae result from the fatal CEE. The study of these supernova progenitors as a class will bring insights on other types of supernova progenitors, as well as on the outcome of the CEE.


title: Proper time in atom interferometers: Diffractive versus specular mirrors

abstract: We compare a conventional Mach-Zehnder light-pulse atom interferometer based on diffractive mirrors with one that uses specular reflection. In contrast to diffractive mirrors that generate a symmetric configuration, specular mirrors realized, for example, by evanescent fields lead under the influence of gravity to an asymmetric geometry. In such an arrangement the interferometer phase contains nonrelativistic signatures of proper time.


title: Improved bound for the dimension of posets of treewidth two

abstract: Joret et al. proved that posets with cover graphs of treewidth at most 2 have dimension at most 1276. Their proof is long and very complex. We give a short and much simpler proof that the dimension of such posets is at most 12.


title: Dynamics on immediate basins for parabolic Newton maps

abstract: Dynamics on parabolic immediate basins for rational Newton maps of entire functions have been studied. It is proved that every parabolic immediate basin contains invariant accesses to the parabolic fixed point at infinity. Moreover, among these accesses there exists a unique dynamically defined access where dynamics are attracted towards the parabolic fixed point, whereas for other accesses, if there is any, the dynamics are repelled.


title: Magnetic absorption of VHE photons in the magnetosphere of the Crab pulsar

abstract: The detection of the pulsed $\sim 1 $~TeV gamma-ray emission from the Crab pulsar reported by MAGIC and VERITAS collaborations demands a substantial revision of existing models of particle acceleration in the pulsar magnetosphere. In this regard model independent restrictions on the possible production site of the VHE photons become an important issue. In this paper, we consider limitations imposed by the process of conversion of VHE gamma rays into $e^{\pm}$ pairs in the magnetic field of the pulsar magnetosphere. Photons with energies exceeding 1~TeV are effectively absorbed even at large distances from the surface of the neutron star. Our calculations of magnetic absorption in the force-free magnetosphere show that the twisting of the magnetic field due to the pulsar rotation makes the magnetosphere more transparent compared to the dipole magnetosphere. The gamma-ray absorption appears stronger for photons emitted in the direction of rotation than in the opposite direction. There is a small angular cone inside which the magnetosphere is relatively transparent and photons with energy $1.5$~TeV can escape from distances beyond $0.1$~light cylinder radius ($R_{\rm{lc}}$). The emission surface from where photons can be emitted in the observer's direction further restricts the sites of VHE gamma-ray production. For the observation angle $57^{\circ}$ relative to the Crab pulsar axis of rotation and the orthogonal rotation, the emission surface in the open field line region is located as close as $0.4\,R_{\rm{lc}}$ from the stellar surface for a dipole magnetic field, and $0.1\,R_{\rm{lc}}$ for a force-free magnetic field.


title: Advances in the Treatment of Trimmed CAD Models due to Isogeometric Analysis

abstract: Trimming is a core technique in geometric modeling. Unfortunately, the resulting objects do not take the requirements of numerical simulations into account and yield various problems. This paper outlines principal issues of trimmed models and highlights different analysis-suitable strategies to address them. It is discussed that these concepts not only provide important computational tools for isogeometric analysis, but can also improve the treatment of trimmed models in a design context.


title: Solving Nurse Scheduling Problem Using Constraint Programming Technique

abstract: Staff scheduling is a universal problem that can be encountered in many organizations, such as call centers, educational institution, industry, hospital, and any other public services. It is one of the most important aspects of workforce management strategy and the one that is most prone to errors or issues as there are many entities should be considered, such as the staff turnover, employee availability, time between rotations, unusual periods of activity, and even the last-minute shift changes. The nurse scheduling problem is a variant of staff scheduling problems which appoints nurses to shifts as well as rooms per day taking both hard constraints, i.e., hospital requirements, and soft constraints, i.e., nurse preferences, into account. Most algorithms used for scheduling problems fall short when it comes to the number of inputs they can handle. In this paper, constraint programming was developed to solve the nurse scheduling problem. The developed constraint programming model was then implemented using python programming language.


title: Deep One-Class Classification Using Intra-Class Splitting

abstract: This paper introduces a generic method which enables to use conventional deep neural networks as end-to-end one-class classifiers. The method is based on splitting given data from one class into two subsets. In one-class classification, only samples of one normal class are available for training. During inference, a closed and tight decision boundary around the training samples is sought which conventional binary or multi-class neural networks are not able to provide. By splitting data into typical and atypical normal subsets, the proposed method can use a binary loss and defines an auxiliary subnetwork for distance constraints in the latent space. Various experiments on three well-known image datasets showed the effectiveness of the proposed method which outperformed seven baselines and had a better or comparable performance to the state-of-the-art.


title: Approximation of solutions of the stochastic wave equation by using the Fourier series

abstract: A one-dimensional stochastic wave equation driven by a general stochastic measure is studied in this paper. The Fourier series expansion of stochastic measures is considered. It is proved that changing the integrator by the corresponding partial sums or by Fej\`{e}r sums we obtain the approximations of mild solution of the equation.


title: Application of image processing in optical method, Moire deflectometry for investigating the optical properties of zinc oxide nanoparticle

abstract: Nowadays for measurement of refractive index of nanomaterials usually spectro-photometric and mechanical methods are used which are expensive and indirect. In this paper for measuring these parameters of zinc oxide nanomaterial with two different stabilizers, a simple optical method, Moire deflectometry, which is based on wave front analysis and geometric optics is used. In the Moire deflectometry method, the beam of light from the laser diode passes through the sample. As a result of that, a change in the sample environment is observed as deflections of the fringes. By recording of these deflections using CCD and image processing with MATLAB, the nanomaterials refractive indices can be calculated. Due to the high accuracy of this method and improved the image processing code in Matlab, the smallest changes of the refractive index in the sample can be measured. Digital Image processing is used for processing images in a way that features can be selected and being shown. The results obtained in this method show a good improvement over the other used methods.


title: Late time kilonova light curves and implications to GW 170817

abstract: We discuss the late time (tens of days) emission from the radioactive ejecta of mergers involving neutron stars, when the ionization energy loss time of beta-decay electrons and positrons exceeds the expansion time. We show that if the e$^\pm$ are confined to the plasma (by magnetic fields), then the time dependence of the plasma heating rate, $\dot{\varepsilon}_d$, and hence of the bolometric luminosity $L=\dot{\varepsilon}_d$, are given by $d\log L/d\log t\simeq-2.8$, nearly independent of the composition and of the instantaneous radioactive energy release rate, $\dot{\varepsilon}$. This universality of the late time behavior is due to the weak dependence of the ionization loss rate on composition and on e$^\pm$ energy. The late time IR and optical measurements of GW 170817 are consistent with this expected behavior provided that the ionization loss time exceeds the expansion time at $t>t_\varepsilon\approx 7$~d, as predicted based on the early (few day) electromagnetic emission.


title: DBSCAN for nonlinear equalization in high-capacity multi-carrier optical communications

abstract: Coherent optical multi-carrier communications have recently dominated metro-regional and long-haul optical communications. However, the major obstacle of networks involving coherent multi-carrier signals such as coherent optical orthogonal frequency-division multiplexing (CO-OFDM) is the fiber-induced nonlinearity and the parametric noise amplification from cascaded optical amplifiers which results in significant nonlinear distortion among subcarriers. Here, we present the first nonlinear equalizer in optical communications using the traditional Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm and a novel modified version of DBSCAN which combines K-means clustering on the noisy un-clustered symbols. For a 24.72 Gbit/sec differential quaternary phase-shift keying (DQPSK) CO-OFDM system, the modified DBSCAN can increase the signal quality-factor by up to 2.158 dB compared to linear equalization at 500 km of transmission. The modified DBSCAN slightly outperforms the traditional DBSCAN, fuzzy-logic C-means, hierarchical and conventional K-means clustering at high launched optical powers.


title: Numerical Reconstruction in Magnetic Particle Imaging

abstract: Magnetic particle imaging (MPI) is a medical imaging modality of recent origin, and it exploits the nonlinear magnetization phenomenon to recover the spatially dependent concentration of the nanoparticles. Currently, image reconstruction in MPI is frequently carried out by standard Tikhonov regularization with nonnegativity constraint, which is then minimized by a Kaczmarz type method. In this work, we revisit several issues in the numerical reconstruction in MPI from the perspective of modern inverse theory, i.e., the choice of data fidelity, and choosing a suitable regularization parameter and accelerating Kaczmarz iteration via randomized singular value decomposition. These algorithmic tricks are straightforward to implement and easy to incorporate in existing reconstruction algorithms. Their significant potentials are illustrated by extensive numerical experiments on a publicly available dataset.


title: Non-Fickian macroscopic model of axial diffusion of granular materials in a long cylindrical tumbler

abstract: We provide new concepts for understanding transport phenomena in flows of granular materials by using a non-Fickian macroscopic model of axial diffusion of a granular material in a finite cylindrical tumbler. The model accounts for diffusion induced by particle collisions only in a thin surface flowing layer due to localization of shear within the cross-section of the drum. All model parameters are related to measurable quantities in a granular flow. It is shown that the proposed model is a member of the general class of linear constitutive relations with memory. An exact solution for the spreading of a finite-width pulse initial condition under the proposed non-Fickian model is derived and compared to the solution of a Fickian model (i.e., the ``classical'' diffusion equation).


title: Spectral Efficiency Analysis in Downlink Massive MIMO System for Perfect CSI with Precoding

abstract: In this paper we first derived the mathematical expression for lower bound spectral efficiency (SE) calculation for zero-force (ZF), and minimum mean square error (MMSE). Secondly, we calculated the simulation SE with three algorithms for ZF and MMSE precoding. We compared the simulation and theoretical results and found that the theoretical results are 1 to 1.5 bits less than the simulation values which implied that the theoretical lower bounds are actually the lower bounds. To achieve the maximum spectral efficiency in downlink massive MIMO systems we assumed perfect CSI, ZF and MMSE precoding in this paper. We also considered that the channel has the characteristics of small and large scale fading (SSF and LSF) as the model is like a practical. We investigated the effect of different SNR, base station (M) and radius (R) of the cell on spectral efficiency for simulation and theoretical results. We also evaluated the performance of SE of each algorithms and precoding schemes for different configurations. From the results we have observed that algorithms-1 and ZF outperform other algorithms and MMSE. From our investigation we noticed that the LSF parameters are the most dominated factor in SE in massive MIMO systems.


title: Probabilistic Power Flow Calculation using Non-intrusive Low-rank Approximation Method

abstract: In this paper, a novel non-intrusive probabilistic power flow (PPF) analysis method based on the low-rank approximation (LRA) is proposed, which can accurately and efficiently estimate the probabilistic characteristics (e.g., mean, variance, probability density function) of the PPF solutions. This method aims at building up a statistically-equivalent surrogate for the PPF solutions through a small number of power flow evaluations. By exploiting the retained tensor-product form of the univariate polynomial basis, a sequential correction-updating scheme is applied, making the total number of unknowns to be linear rather than exponential to the number of random inputs. Consequently, the LRA method is particularly promising for dealing with high-dimensional problems with a large number of random inputs. Numerical studies on the IEEE 39-bus, 118-bus, and 1354-bus systems show that the proposed method can achieve accurate probabilistic characteristics of the PPF solutions with much less computational effort compared to the Monte Carlo simulations. Even compared to the polynomial chaos expansion method, the LRA method can achieve comparable accuracy, while the LRA method is more capable of handling higher-dimensional problems. Moreover, numerical results reveal that the randomness brought about by the renewable energy resources and loads may inevitably affect the feasibility of dispatch/planning schemes.


title: Determination of hybrid and direct laser acceleration dominated regimes in a 55fs laser driven plasma accelerator with ionization induced injection

abstract: An experimental study on 55fs laser driven plasma accelerator using mixed gas-jet target with varying plasma density is used to identify the role of different acceleration mechanisms, viz. Direct Laser Acceleration (DLA) and wakefield. At lower electron density electron acceleration could be attributed mainly to DLA with ionization induced injection. With increase in density, increasing role of wakefield was observed leading to hybrid regime, and at densities higher than self-injection threshold, observed experimentally for He target contribution of DLA and wakefield was found to be comparable. Dominant DLA mechanism was also observed in case of pure nitrogen target. 2D PIC simulations performed using the EPOCH code corroborate the above scenario, and also showed generation of surface waves, considered as a potential mechanism of pre-acceleration to DLA.


title: A general model for dynamic contact angle over full speed regime

abstract: The prevailing models for advancing dynamic contact angle are under intensive debates, and the fitting performances are far from satisfying in practice. The present study proposes a model based on the recent understanding of the multi-scale structure and the local friction at the contact line. The model has unprecedented fitting performance for dynamic contact angle over wide spreading speed regime across more than five orders of magnitude. The model also well applies for non-monotonous angle variations which have long been considered as abnormal. The model has three fitting parameters, and two of them are nearly predictable at this stage. The one denoting the multi-scale ratio was nearly constant and the one representing the primary frictional coefficient, has a simple correlation with the liquid bulk viscosity.


title: Taylor States in Stellarators: A Fast High-order Boundary Integral Solver

abstract: We present a boundary integral equation solver for computing Taylor relaxed states in non-axisymmetric solid and shell-like toroidal geometries. The computation of Taylor states in these geometries is a key element for the calculation of stepped pressure stellarator equilibria. The integral representation of the magnetic field in this work is based on the generalized Debye source formulation, and results in a well-conditioned second-kind boundary integral equation. The integral equation solver is based on a spectral discretization of the geometry and unknowns, and the computation of the associated weakly-singular integrals is performed with high-order quadrature based on a partition of unity. The resulting scheme for applying the integral operator is then coupled with an iterative solver and suitable preconditioners. Several numerical examples are provided to demonstrate the accuracy and efficiency of our method, and a direct comparison with the leading code in the field is reported.


title: Recycling Solutions for Vertex Coloring Heuristics

abstract: The vertex coloring problem is a well-known NP-hard problem and has many applications in operations research and in scheduling. A conventional approach to the problem solves the k-colorability problem iteratively, decreasing k one by one. Whether a heuristic algorithm finds a legal k-coloring quickly or not is largely affected by an initial solution. We highlight a simple initial solution generator, which we call the recycle method, which makes use of the legal (k+1)-coloring that has been found. An initial solution generated by the method is expected to guide a heuristic algorithm to find a legal k-coloring more quickly than conventional methods, as demonstrated by experimental studies. The results suggest that the recycle method should be used as the standard initial solution generator for both local search algorithms and modern hybrid methods.


title: Islanding, growth mode and ordering in Si heteroepitaxy on Ge(001) substrates structured by thermal annealing

abstract: Si/Ge heteroepitaxial dots under tensile strain are grown on nanostructured Ge substrates produced by high-temperature flash heating exploiting the spontaneous faceting of the Ge(001) surface close to the onset of surface melting. A very diverse growth mode is obtained depending on the specific atomic structure and step density of nearby surface domains with different vicinal crystallographic orientations. On highly-miscut areas of the Ge(001) substrate, the critical thickness for islanding is lowered to about 5 ML, in contrast to the 11 ML reported for the flat Ge(001) surface, while on unreconstructed (1x1) domains the growth is Volmer-Weber driven. An explanation is proposed considering the diverse relative contributions of step and surface energies on misoriented substrates. In addition, we show that the bottom-up pattern of the substrate naturally formed by thermal annealing determines a spatial correlation for the dot sites.


title: A New Family of Neural Networks Provably Resistant to Adversarial Attacks

abstract: Adversarial attacks add perturbations to the input features with the intent of changing the classification produced by a machine learning system. Small perturbations can yield adversarial examples which are misclassified despite being virtually indistinguishable from the unperturbed input. Classifiers trained with standard neural network techniques are highly susceptible to adversarial examples, allowing an adversary to create misclassifications of their choice. We introduce a new type of network unit, called MWD (max of weighed distance) units that have a built-in resistant to adversarial attacks. These units are highly non-linear, and we develop the techniques needed to effectively train them. We show that simple interval techniques for propagating perturbation effects through the network enables the efficient computation of robustness (i.e., accuracy guarantees) for MWD networks under any perturbations, including adversarial attacks. MWD networks are significantly more robust to input perturbations than ReLU networks. On permutation invariant MNIST, when test examples can be perturbed by 20% of the input range, MWD networks provably retain accuracy above 83%, while the accuracy of ReLU networks drops below 5%. The provable accuracy of MWD networks is superior even to the observed accuracy of ReLU networks trained with the help of adversarial examples. In the absence of adversarial attacks, MWD networks match the performance of sigmoid networks, and have accuracy only slightly below that of ReLU networks.


title: Minimal-excitation single-particle emitters: A comparison of charge and energy transport properties

abstract: We investigate different types of time-dependently driven single-particle sources whose common feature is that they produce pulses of integer charge and minimally excite the Fermi sea. These sources are: a slowly driven mesoscopic capacitor, a Lorentzian-shaped time-dependent bias voltage, and a local gate-voltage modulation of a quantum Hall edge state. They differ by their specific driving protocols, e.g., they have a pure ac driving or a driving with a dc component. In addition, only in the first of these setups, strong confinement leading to a discrete energy spectrum of the conductor, is exploited for the single-particle emission. Here, we study if and how these basic differences impact transport properties. Specifically, we address time- and energy-resolved charge and energy currents, as well as their zero-frequency correlators (charge-, energy- and mixed noise), as they are frequently used to characterize experiments in quantum optics with electrons. Beyond disparities due to a different number and polarity of particles emitted per period, we in particular identify differences in the impact, which temperature has on the observables for sources with and without energy-dependent scattering properties. We trace back these characteristics to a small set of relevant parameter ratios.


title: Epileptiform spikes in specific left temporal and mesial temporal structures disrupt verbal episodic memory encoding

abstract: Patients diagnosed with epilepsy experience cognitive dysfunction that may be due to a transient cognitive/memory impairment (TCI/TMI) caused by spontaneous epileptiform spikes. We asked in a cohort of 166 adult patients with medically refractory focal epilepsy if spikes in specific neuroanatomical regions during verbal episodic memory encoding would significantly decrease the probability of recall. We found using a na\"ive Bayesian machine learning model that the probability of correct word recall decreased significantly by 11.9% when spikes occurred in left Brodmann area 21 (BA)21, (p<0.001), 49.7% in left BA38 (p=0.01), and 32.2% in right BA38 (p<0.001), and 21.4% in left BA36 (p<0.01). We also examined the influence of the seizure-onset zone and the language dominant hemisphere on this effect. Our results demonstrate that spontaneous epileptiform spikes produce a large effect TCI/TMI in brain regions known to be important in semantic processing and episodic memory. Thus memory impairment in patients with epilepsy may be attributable to cellular events associated with abnormal inter-ictal electrical events.


title: A next-generation LHC heavy-ion experiment

abstract: The present document discusses plans for a compact, next-generation multi-purpose detector at the LHC as a follow-up to the present ALICE experiment. The aim is to build a nearly massless barrel detector consisting of truly cylindrical layers based on curved wafer-scale ultra-thin silicon sensors with MAPS technology, featuring an unprecedented low material budget of 0.05% X$_0$ per layer, with the innermost layers possibly positioned inside the beam pipe. In addition to superior tracking and vertexing capabilities over a wide momentum range down to a few tens of MeV/$c$, the detector will provide particle identification via time-of-flight determination with about 20~ps resolution. In addition, electron and photon identification will be performed in a separate shower detector. The proposed detector is conceived for studies of pp, pA and AA collisions at luminosities a factor of 20 to 50 times higher than possible with the upgraded ALICE detector, enabling a rich physics program ranging from measurements with electromagnetic probes at ultra-low transverse momenta to precision physics in the charm and beauty sector.


title: Model-Based Detector for SSDs in the Presence of Inter-cell Interference

abstract: In this paper, we consider the problem of reducing the bit error rate of flash-based solid state drives (SSDs) when cells are subject to inter-cell interference (ICI). By observing that the outputs of adjacent victim cells can be correlated due to common aggressors, we propose a novel channel model to accurately represent the true flash channel. This model, equivalent to a finite-state Markov channel model, allows the use of the sum-product algorithm to calculate more accurate posterior distributions of individual cell inputs given the joint outputs of victim cells. These posteriors can be easily mapped to the log-likelihood ratios that are passed as inputs to the soft LDPC decoder. When the output is available with high precision, our simulation showed that a significant reduction in the bit-error rate can be obtained, reaching $99.99\%$ reduction compared to current methods, when the diagonal coupling is very strong. In the realistic case of low-precision output, our scheme provides less impressive improvements due to information loss in the process of quantization. To improve the performance of the new detector in the quantized case, we propose a new iterative scheme that alternates multiple times between the detector and the decoder. Our simulations showed that the iterative scheme can significantly improve the bit error rate even in the quantized case.


title: Phase-Resolved Spectroscopy of Gaia14aae: Line Emission From Near the White Dwarf Surface

abstract: AM CVn binaries are a class of ultracompact, hydrogen-deficient binaries, each consisting of a white dwarf accreting helium-dominated material from a degenerate or semi-degenerate donor star. Of the 56 known systems, only Gaia14aae undergoes complete eclipses of its central white dwarf, allowing the parameters of its stellar components to be tightly constrained. Here, we present phase-resolved optical spectroscopy of Gaia14aae. We use the spectra to test the assumption that the narrow emission feature known as the `central spike' traces the motion of the central white dwarf. We measure a central spike velocity amplitude of $13.8 \pm 3.2$ km/s, which agrees at the 1 $\sigma$ level with the predicted value of $17.6 \pm 1.0$ km/s based on eclipse-derived system parameters. The orbital phase offset of the central spike from its expected position is $4 \pm 15$ $^\circ$, consistent with 0 $^\circ$. Doppler maps of the He I lines in Gaia14aae show two accretion disc bright spots, as seen in many AM CVn systems. The formation mechanism for the second spot remains unclear. We detect no hydrogen in the system, but we estimate a 3 $\sigma$ limit on H$\alpha$ emission with an equivalent width of -1.14 \AA. Our detection of nitrogen and oxygen with no corresponding detection of carbon, in conjunction with evidence from recent studies, mildly favours a formation channel in which Gaia14aae is descended from a cataclysmic variable with a significantly evolved donor.


title: Unconventional Bloch-Gr\"{u}neisen scattering in hybrid Bose-Fermi systems

abstract: We report on the novel mechanism of electron scattering in hybrid Bose-Fermi systems consisting of a two-dimensional electron gas in the vicinity of an exciton condensate: We show that a pair-of-bogolons--mediated scattering proves to be dominating over the conventional acoustic phonon channel and over the single-bogolon scattering, even if the screening is taken into account. We develop a microscopic theory of this effect, focusing on GaAs and MoS$_2$ materials, and find the principal temperature dependence of resistivity, distinct from the conventional phonon--mediated processes. Further, we scrutinize parameters and suggest a way to design composite samples with predefined electron mobilities and propose a mechanism of electron pairing for superconductivity.


title: New Risk Bounds for 2D Total Variation Denoising

abstract: 2D Total Variation Denoising (TVD) is a widely used technique for image denoising. It is also an important nonparametric regression method for estimating functions with heterogenous smoothness. Recent results have shown the TVD estimator to be nearly minimax rate optimal for the class of functions with bounded variation. In this paper, we complement these worst case guarantees by investigating the adaptivity of the TVD estimator to functions which are piecewise constant on axis aligned rectangles. We rigorously show that, when the truth is piecewise constant, the ideally tuned TVD estimator performs better than in the worst case. We also study the issue of choosing the tuning parameter. In particular, we propose a fully data driven version of the TVD estimator which enjoys similar worst case risk guarantees as the ideally tuned TVD estimator.


title: Laser refrigeration using exciplex resonances in gas filled hollow-core fibres

abstract: We theoretically study prospects and limitations of a new route towards macroscopic scale laser refrigeration based on exciplex-mediated frequency up-conversion in gas filled hollow-core fibres. Using proven quantum optical rate equations we model the dynamics of a dopant-buffer gas mixture filling an optically pumped waveguide. In the particular example of alkali-noble gas mixtures, recent high pressure gas cell setup experiments have shown that efficient kinetic energy extraction cycles appear via the creation of transient exciplex excited electronic bound states. The cooling cycle consists of absorption of lower energy laser photons during collisions followed by blue-shifted spontaneous emission on the atomic line of the alkali atoms. For any arbitrary dopant-buffer gas mixture, we derive scaling laws for cooling power, cooling rates and temperature drops with varying input laser power, dopant and buffer gas concentration, fibre geometry and particularities of the exciplex ground and excited state potential landscapes.


title: Fermion Dark Matter with Scalar Triplet at Direct and Collider Searches

abstract: Fermion dark matter (DM) as an admixture of additional singlet and doublet vector like fermions provides an attractive and allowed framework by relic density and direct search constraints within TeV scale, although limited by its discovery potential at the Large Hadron Collider (LHC). An extension of the model with scalar triplet can yield neutrino masses and provide some cushion to the direct search constraint of the DM through pseudo-Dirac mass splitting. This in turn, allow the model to live in a larger region of the parameter space and open the door for detection at LHC, even if slightly. The model however can see an early discovery at International Linear Collider (ILC) without too much of fine-tuning. The complementarity of LHC, ILC and direct search prospect of this framework is studied in this paper.


title: First-order continuous- and discontinuous-Galerkin moment models for a linear kinetic equation: model derivation and realizability theory

abstract: We provide two new classes of moment models for linear kinetic equations in slab and three-dimensional geometry. They are based on classical finite elements and low-order discontinuous-Galerkin approximations on the unit sphere. We investigate their realizability conditions and other basic properties. Numerical tests show that these models are more efficient than classical full-moment models in a space-homogeneous test, when the analytical solution is not smooth.


title: Local minimax rates for closeness testing of discrete distributions

abstract: We consider the closeness testing problem for discrete distributions. The goal is to distinguish whether two samples are drawn from the same unspecified distribution, or whether their respective distributions are separated in $L_1$-norm. In this paper, we focus on adapting the rate to the shape of the underlying distributions, i.e. we consider \textit{a local minimax setting}. We provide, to the best of our knowledge, the first local minimax rate for the separation distance up to logarithmic factors, together with a test that achieves it. In view of the rate, closeness testing turns out to be substantially harder than the related one-sample testing problem over a wide range of cases.

categories: math.ST cs.IT cs.LG math.IT stat.ML stat.TH


title: Adaptive Gradient for Adversarial Perturbations Generation

abstract: Deep Neural Networks have achieved remarkable success in computer vision, natural language processing, and audio tasks.


title: Precision and informational limits in inelastic optical spectroscopy

abstract: Using Fisher information and the Cram\'er-Rao lower bound, we analyse fundamental precision limits in the determination of spectral parameters in inelastic optical scattering. General analytic formulae are derived which account for the instrument response functions of the dispersive element and relay optics found in practical Raman and Brillouin spectrometers. Limiting cases of dispersion and diffraction limited spectrometers, corresponding to measurement of Lorentzian and Voigt lineshapes respectively, are discussed in detail allowing optimal configurations to be identified. Effects of defocus, spherical aberration, detector pixellation and a finite detector size are also considered.


title: Preliminary model of the outer disk of RU Lup presently showing only four dark gaps

abstract: Recent ALMA/DSHARP observations reveal that the disk of RU Lup shows presently only four dark gaps at large radii in which planets are presumably already orbiting. This gap arrangement is quite similar to the current structure of the outer planets in our solar system. Although there may be more undetected planets forming in the inner disk of RU Lup, the information obtained from the current observations prompted us to fit preliminary models with an isothermal oscillatory density profile to the outer disk of RU Lup. The best-fit model confirms that several physical properties of the disk of RU Lup are comparable to those of our solar nebula to within factors of 1.0-2.3; but there are some differences as well: central density ($\rho_0$), centrifugal support ($\beta_0$), and core angular velocity ($\Omega_0$) in RU Lup are smaller by factors of 5.5, 8.0, and 18, respectively, whereas the inner core radius ($R_1$) of RU Lup is 25 times larger. The smaller values found for RU Lup are linked to the (apparently gap-empty) inner 10 AU of its disk, a region that is presently inadequately resolved.


title: Compositional Abstraction of Large-Scale Stochastic Systems: A Relaxed Dissipativity Approach

abstract: In this paper, we propose a compositional approach for the construction of finite abstractions (a.k.a. finite Markov decision processes (MDPs)) for networks of discrete-time stochastic control subsystems that are not necessarily stabilizable. The proposed approach leverages the interconnection topology and a notion of finite-step stochastic storage functions, that describes joint dissipativity-type properties of subsystems and their abstractions, and establishes a finite-step stochastic simulation function as a relation between the network and its abstraction. To this end, we first develop a new type of compositionality conditions which is less conservative than the existing ones. In particular, using a relaxation via a finite-step stochastic simulation function, it is possible to construct finite abstractions such that stabilizability of each subsystem is not necessarily required. We then propose an approach to construct finite MDPs together with their corresponding finite-step storage functions for general discrete-time stochastic control systems satisfying an incremental passivablity property. We also construct finite MDPs for a particular class of nonlinear stochastic control systems. To demonstrate the effectiveness of the proposed results, we apply our results on three different case studies.


title: Estimating the Mixing Time of Ergodic Markov Chains

abstract: We address the problem of estimating the mixing time $t_{\mathsf{mix}}$ of an arbitrary ergodic finite-state Markov chain from a single trajectory of length $m$. The reversible case was addressed by Hsu et al. [2019], who left the general case as an open problem. In the reversible case, the analysis is greatly facilitated by the fact that the Markov operator is self-adjoint, and Weyl's inequality allows for a dimension-free perturbation analysis of the empirical eigenvalues. As Hsu et al. point out, in the absence of reversibility (which induces asymmetric pair probabilities matrices), the existing perturbation analysis has a worst-case exponential dependence on the number of states $d$. Furthermore, even if an eigenvalue perturbation analysis with better dependence on $d$ were available, in the non-reversible case the connection between the spectral gap and the mixing time is not nearly as straightforward as in the reversible case. Our key insight is to estimate the pseudo-spectral gap $\gamma_{\mathsf{ps}}$ instead, which allows us to overcome the loss of symmetry and to achieve a polynomial dependence on the minimal stationary probability $\pi_\star$ and $\gamma_{\mathsf{ps}}$. Additionally, in the reversible case, we obtain simultaneous nearly (up to logarithmic factors) minimax rates in $t_{\mathsf{mix}}$ and precision $\varepsilon$, closing a gap in Hsu et al., who treated $\varepsilon$ as constant in the lower bounds. Finally, we construct fully empirical confidence intervals for $\gamma_{\mathsf{ps}}$, which shrink to zero at a rate of roughly $1/\sqrt{m}$, and improve the state of the art in even the reversible case.


title: Semantic Factorization and Descent

abstract: Let $\mathbb{A}$ be a $2$-category with suitable opcomma objects and pushouts. We give a direct proof that, provided that the codensity monad of a morphism $p$ exists and is preserved by a suitable morphism, the factorization given by the lax descent object of the higher cokernel of $p$ is up to isomorphism the same as the semantic factorization of $p$, either one existing if the other does. The result can be seen as a counterpart account to the celebrated B\'{e}nabou-Roubaud theorem. This leads in particular to a monadicity theorem, since it characterizes monadicity via descent. It should be noted that all the conditions on the codensity monad of $p$ trivially hold whenever $p$ has a left adjoint and, hence, in this case, we find monadicity to be a $2$-dimensional exact condition on $p$, namely, to be an effective faithful morphism of the $2$-category $\mathbb{A} $.


title: Analysis and Application of Optimal Transport For Challenging Seismic Inverse Problems

abstract: In seismic exploration, sources and measurements of seismic waves on the surface are used to determine model parameters representing geophysical properties of the earth. Full-waveform inversion (FWI) is a nonlinear seismic inverse technique that inverts the model parameters by minimizing the difference between the synthetic data from the forward wave propagation and the observed true data in PDE-constrained optimization. The traditional least-squares method of measuring this difference suffers from three main drawbacks including local minima trapping, sensitivity to noise, and difficulties in reconstruction below reflecting layers. Unlike the local amplitude comparison in the least-squares method, the quadratic Wasserstein distance from the optimal transport theory considers both the amplitude differences and the phase mismatches when measuring data misfit. We will briefly review our earlier development and analysis of optimal transport-based inversion and include improvements, for example, a stronger convexity proof. The main focus will be on the third "challenge" with new results on sub-reflection recovery.


title: Stochastic models associated to a Nonlocal Porous Medium Equation

abstract: The nonlocal porous medium equation considered in this paper is a degenerate nonlinear evolution equation involving a space pseudo-differential operator of fractional order. This space-fractional equation admits an explicit, nonnegative, compactly supported weak solution representing a probability density function. In this paper we analyze the link between isotropic transport processes, or random flights, and the nonlocal porous medium equation. In particular, we focus our attention on the interpretation of the weak solution of the nonlinear diffusion equation by means of random flights.


title: Stabilization Time in Weighted Minority Processes

abstract: A minority process in a weighted graph is a dynamically changing coloring. Each node repeatedly changes its color in order to minimize the sum of weighted conflicts with its neighbors. We study the number of steps until such a process stabilizes. Our main contribution is an exponential lower bound on stabilization time. We first present a construction showing this bound in the adversarial sequential model, and then we show how to extend the construction to establish the same bound in the benevolent sequential model, as well as in any reasonable concurrent model. Furthermore, we show that the stabilization time of our construction remains exponential even for very strict switching conditions, namely, if a node only changes color when almost all (i.e., any specific fraction) of its neighbors have the same color. Our lower bound works in a wide range of settings, both for node-weighted and edge-weighted graphs, or if we restrict minority processes to the class of sparse graphs.


title: The boundary of the Milnor fibre of certain non-isolated singularities

abstract: Let $ \Phi: ({\mathbb C}^2, 0) \to ( {\mathbb C}^3, 0) $ be a finitely determined complex analytic germ and let $(\{f=0\},0)$ be the reduced equation of its image, a non-isolated hypersurface singularity. We provide the plumbing graph of the boundary of the Milnor fibre of $f$ from the double-point-geometry of $\Phi$.


title: On generalized stochastic fractional integrals and related inequalities

abstract: The generalized mean-square fractional integrals $\mathcal{J}_{\rho,\lambda,u+;\omega}^{\sigma}$ and $\mathcal{J}_{\rho,\lambda,v-;\omega}^{\sigma}$ of the stochastic process $X$ are introduced. Then, for Jensen-convex and strongly convex stochastic proceses, the generalized fractional Hermite--Hadamard inequality is establish via generalized stochastic fractional integrals.


title: Huygens' Dipole for Polarization-Controlled Nanoscale Light Routing

abstract: Structured illumination allows for satisfying the first Kerker condition of in-phase perpendicular electric and magnetic dipole moments in any isotropic scatterer that supports electric and magnetic dipolar resonances. The induced Huygens' dipole may be utilized for unidirectional coupling to waveguide modes that propagate transverse to the excitation beam. We study two configurations of a Huygens' dipole -- longitudinal electric and transverse magnetic dipole moments or vice versa. We experimentally show that only the radially polarized emission of the first and azimuthally polarized emission of the second configuration are directional in the far-field. This polarization selectivity implies that directional excitation of either TM or TE waveguide modes is possible. Applying this concept to a single nanoantenna excited with structured light, we are able to experimentally achieve scattering directivities of around 23 dB and 18 dB in TM and TE modes, respectively. This strong directivity paves the way for tunable polarization-controlled nanoscale light routing and applications in optical metrology, localization microscopy and on-chip optical devices.


title: Can Oscillating Neutrino States Be Formulated Universally?

abstract: A standing problem in neutrino physics is the consistent and universal definition of oscillating neutrino states as coherent superpositions of massive neutrino states. This problem is solved in a quantum field theoretical framework of neutrino mixing developed in analogy with the Nambu--Jona-Lasinio model for the dynamical generation of nucleon masses. The massive neutrino states are Bogoliubov quasiparticles and their vacuum is a condensate of "Cooper pairs" of massless flavour neutrinos. Their superpositions as oscillating neutrino states have intrinsic quantum coherence by construction. In this quantization framework, the standard phenomenological flavour neutrino states and oscillation probability formula are validated in the ultrarelativistic approximation.


title: Improved Visualisation of Patient-Specific Heart Structure Using Three-Dimensional Printing Coupled with Image-Processing Techniques Inspired by Astrophysical Methods

abstract: The aim of this study is to use image-processing techniques developed in the field of astrophysics as inspiration for a novel approach to the three-dimensional (3D) imaging of periprocedural medical data, with the intention of providing improved visualisation of patient-specific heart structure and thereby allowing for an improved quality of procedural planning with regards to individualized cardiovascular healthcare. Using anonymized patient DICOM data for a cardiac computed tomography (CT) angiography, two-dimensional slices of the patient's heart were processed using a series of software packages in order to produce an accurate 3D representation of the patient's heart tissue as a computer-generated stereolithography (STL) file, followed by the creation of a tactile 3D printout. We find that the models produced provide clear definition of heart structure, in particular in the left atrium, left ventricle and aorta. This level of clarity also allows for the aortic valve to be observed and 3D printed. This study provides a step-by-step blueprint of how this can be achieved using open source software, specifically Slicer 4.8.1, MeshLab and AutoDesk Netfabb. In addition, the implementation of astrophysical image-processing techniques shows an improvement in modelling of the heart based on the CT data, in particular in the case of small-scale features where echocardiography has previously been required for more reliable results.


title: EBIT Observation of Ar Dielectronic Recombination Lines Near the Unknown Faint X-Ray Feature Found in the Stacked Spectrum of Galaxy Clusters

abstract: Motivated by possible atomic origins of the unidentified emission line detected at 3.55 keV to 3.57 keV in a stacked spectrum of galaxy clusters (Bulbul et al. 2014), an electron beam ion trap (EBIT) was used to investigate the resonant dielectronic recombination (DR) process in highly-charged argon ions as a possible contributor to the emission feature. The He-like Ar DR-induced transition 1s$^2$2l - 1s2l3l$^\prime$ was suggested to produce a 3.62 keV photon (Bulbul et al. 2014) near the unidentified line at 3.57 keV and was the starting point of our investigation. The collisional-radiative model NOMAD was used to create synthetic spectra for comparison with both our EBIT measurements and with spectra produced with the AtomDB database/Astrophysical Plasma Emission Code (APEC) used in the Bulbul et al. (2014) work. Excellent agreement was found between the NOMAD and EBIT spectra, providing a high level of confidence in the atomic data used. Comparison of the NOMAD and APEC spectra revealed a number of missing features in the AtomDB database near the unidentified line. At an electron temperature of $T_e$ = 1.72 keV, the inclusion of the missing lines in AtomDB increases the total flux in the 3.5 keV to 3.66 keV energy band by a factor of 2. While important, this extra emission is not enough to explain the unidentified line found in the galaxy cluster spectra.


title: Robustness Certificates Against Adversarial Examples for ReLU Networks

abstract: While neural networks have achieved high performance in different learning tasks, their accuracy drops significantly in the presence of small adversarial perturbations to inputs. Defenses based on regularization and adversarial training are often followed by new attacks to defeat them. In this paper, we propose attack-agnostic robustness certificates for a multi-label classification problem using a deep ReLU network. Although computing the exact distance of a given input sample to the classification decision boundary requires solving a non-convex optimization, we characterize two lower bounds for such distances, namely the simplex certificate and the decision boundary certificate. These robustness certificates leverage the piece-wise linear structure of ReLU networks and use the fact that in a polyhedron around a given sample, the prediction function is linear. In particular, the proposed simplex certificate has a closed-form, is differentiable and is an order of magnitude faster to compute than the existing methods even for deep networks. In addition to theoretical bounds, we provide numerical results for our certificates over MNIST and compare them with some existing upper bounds.


title: Monolithic integration of broadband optical isolators for polarization-diverse silicon photonics

abstract: Integrated optical isolators have been a longstanding challenge for photonic integrated circuits (PIC). An ideal integrated optical isolator for PIC should be made by a monolithic process, have a small footprint, exhibit broadband and polarization-diverse operation, and be compatible with multiple materials platforms. Despite significant progress, the optical isolators reported so far do not meet all these requirements. In this article we present monolithically integrated broadband magneto-optical isolators on silicon and silicon nitride (SiN) platforms operating for both TE and TM modes with record high performances, fulfilling all the essential characteristics for PIC applications. In particular, we demonstrate fully-TE broadband isolators by depositing high quality magneto-optical garnet thin films on the sidewalls of Si and SiN waveguides, a critical result for applications in TE-polarized on-chip lasers and amplifiers. This work demonstrates monolithic integration of high performance optical isolators on chip for polarization-diverse silicon photonic systems, enabling new pathways to impart nonreciprocal photonic functionality to a variety of integrated photonic devices.


title: Ordinal Patterns in Clusters of Subsequent Extremes of Regularly Varying Time Series

abstract: In this paper, we investigate temporal clusters of extremes defined as subsequent exceedances of high thresholds in a stationary time series. Two meaningful features of these clusters are the probability distribution of the cluster size and the ordinal patterns within a cluster. Since these patterns take only the ordinal structure of consecutive data points into account the method is robust under monotone transformations and measurement errors. We verify the existence of the corresponding limit distributions in the framework of regularly varying time series, develop non-parametric estimators and show their asymptotic normality under appropriate mixing conditions. The performance of the estimators is demonstrated in a simulated example and a real data application to discharge data of the river Rhine.


title: Dynamic robust stabilization of fractional-order linear systems with nonlinear uncertain parameters: An LMI approach

abstract: This paper considers the problem of robust stability and stabilization for linear fractional-order system with nonlinear uncertain parameters, with fractional order 0<a<2. A dynamic output feedback controller, with predetermined order, for asymptotically stabilizing such uncertain fractional-order systems is designed. The derived stabilization conditions are in LMI form. Simulation results of two numerical examples illustrate that the proposed sufficient theoretical results are applicable and effective for tackling robust stabilization problems. Keywords: Fractional-order system, nonlinear uncertain parameters, linear matrix inequality (LMI), robust stabilization, dynamic output feedback.


title: A Practical Algorithm for Multiplayer Bandits when Arm Means Vary Among Players

abstract: We study a multiplayer stochastic multi-armed bandit problem in which players cannot communicate, and if two or more players pull the same arm, a collision occurs and the involved players receive zero reward. We consider the challenging heterogeneous setting, in which different arms may have different means for different players, and propose a new and efficient algorithm that combines the idea of leveraging forced collisions for implicit communication and that of performing matching eliminations. We present a finite-time analysis of our algorithm, giving the first sublinear minimax regret bound for this problem, and prove that if the optimal assignment of players to arms is unique, our algorithm attains the optimal $O(\ln(T))$ regret, solving an open question raised at NeurIPS 2018.


title: PIPPS: Flexible Model-Based Policy Search Robust to the Curse of Chaos

abstract: Previously, the exploding gradient problem has been explained to be central in deep learning and model-based reinforcement learning, because it causes numerical issues and instability in optimization. Our experiments in model-based reinforcement learning imply that the problem is not just a numerical issue, but it may be caused by a fundamental chaos-like nature of long chains of nonlinear computations. Not only do the magnitudes of the gradients become large, the direction of the gradients becomes essentially random. We show that reparameterization gradients suffer from the problem, while likelihood ratio gradients are robust. Using our insights, we develop a model-based policy search framework, Probabilistic Inference for Particle-Based Policy Search (PIPPS), which is easily extensible, and allows for almost arbitrary models and policies, while simultaneously matching the performance of previous data-efficient learning algorithms. Finally, we invent the total propagation algorithm, which efficiently computes a union over all pathwise derivative depths during a single backwards pass, automatically giving greater weight to estimators with lower variance, sometimes improving over reparameterization gradients by $10^6$ times.


title: Equivariant Kazhdan-Lusztig polynomials of thagomizer matroids

abstract: The equivariant Kazhdan-Lusztig polynomial of a matroid was introduced by Gedeon, Proudfoot, and Young. Gedeon conjectured an explicit formula for the equivariant Kazhdan-Lusztig polynomials of thagomizer matroids with an action of symmetric groups. In this paper, we discover a new formula for these polynomials which is related to the equivariant Kazhdan-Lusztig polynomials of uniform matroids. Based on our new formula, we confirm Gedeon's conjecture by the Pieri rule.


title: The electromagnetic form factors of $\Lambda$ hyperon in the vector meson dominance model

abstract: We perform an analysis on the electromagnetic form factors of the $\Lambda$ hyperon in the time-like reaction $e^+e^-\rightarrow \Lambda\bar\Lambda$ by using a modified vector meson dominance model. We consider both the intrinsic structure components and the meson clouds components. For the latter one, we not only include the contributions from the $\phi$ and $\omega$ mesons, but also take into account the contributions from the resonance states $\omega(1420)$, $\omega(1650)$, $\phi(1680)$ and $\phi(2170)$. We extract the model parameters by combined fit to the time-like effective form factor $|G_{\rm{eff}}|$, the electromagnetic form factor ratio $|G_E/G_M|$ and the relative phase $\Delta\Phi$ of the $\Lambda$ hyperon from the BaBar and BESIII Collaborations. We find that the vector meson dominance model can simultaneously describe these observables. Particularly, the inclusion of the resonance states in the model is necessary for explaining the ratio $|G_E/G_M|$ in a wide range of $\sqrt{s}$ as well as the large phase angle. With the fitted parameters, we predict the single and double polarization observables, which could be measured in polarized annihilation reactions. Moreover, we analytically continue the expression of the form factors to space-like region and estimate the space-like form factors of $\Lambda$ hyperon.


title: Individually addressable double quantum dots formed with nanowire polytypes and identified by epitaxial markers

abstract: Double quantum dots (DQDs) hold great promise as building blocks for quantum technology as they allow for two electronic states to coherently couple. Defining QDs with materials rather than using electrostatic gating allows for QDs with a hard-wall confinement potential and more robust charge and spin states. An unresolved problem is how to individually address these quantum dots, which is necessary for controlling quantum states. We here report the fabrication of double quantum dot devices defined by the conduction band edge offset at the interface of the wurtzite and zinc blende crystal phases of InAs in nanowires. By using sacrifical epitaxial GaSb markers selectively forming on one crystal phase, we are able to precisely align gate electrodes allowing us to probe and control each QD independently. We hence observe textbook-like charge stability diagrams, a discrete energy spectrum and electron numbers consistent with theoretical estimates and investigate the tunability of the devices, finding that changing the electron number can be used to tune the tunnel barrier as expected by simple band diagram arguments.


title: Martingale-like sequences in Banach lattices

abstract: Martingale-like sequences in vector lattice and Banach lattice frameworks are defined in the same way as martingales are defined in [Positivity 9 (2005), 437--456]. In these frameworks, a collection of bounded $X$-martingales is shown to be a Banach space under the supremum norm, and under some conditions it is also a Banach lattice with coordinate-wise order. Moreover, a necessary and sufficient condition is presented for the collection of $\mathcal{E}$-martingales to be a vector lattice with coordinate-wise order. It is also shown that the collection of bounded $\mathcal{E}$-martingales is a normed lattice but not necessarily a Banach space under the supremum norm.


title: Wavenumber Selection in Pattern Forming Systems

abstract: Wavenumber selection in pattern forming systems remains a long standing puzzle in physics. Previous studies have shown that external noise is a possible mechanism for wavenumber selection. We conduct an extensive numerical study of the noisy stabilized Kuramoto Sivashinsky equation. We use a fast spectral method of integration, which enables us to investigate long time behavior for large system sizes that could not be investigated by earlier work. We find that a state with a unique wavenumber has the highest probability of occurring at very long times. We also find that this state is independent of the strength of the noise and initial conditions, thus making a convincing case for the role of noise as a mechanism of state selection.


title: On the fundamentals of Richtmyer-Meshkov dynamics with variable acceleration

abstract: Richtmyer-Meshkov instability (RMI) plays important role in nature and technology, from supernovae and fusion to scramjets and nano-fabrication. Canonical Richtmyer-Meshkov instability is induced by a steady shock and impulsive acceleration, whereas in realistic environments the acceleration is usually variable. This work focuses on RMI induced by acceleration with a power-law time-dependence, and applies group theory to solve the classical problem. For early-time dynamics, we find the dependence of RMI growth-rate on the initial conditions and show it is free from the acceleration parameters. For late time dynamics, we find a continuous family of regular asymptotic solutions, including their curvature, velocity, Fourier amplitudes, and interfacial shear, and we study the solutions stability. For each of the solutions, the interface dynamics is directly linked to the interfacial shear, and the non-equilibrium velocity field has intense fluid motion near the interface and effectively no motion in the bulk. The quasi-invariance of the fastest stable solution suggests that nonlinear coherent dynamics in RMI is characterized by two macroscopic length-scales - the wavelength and the amplitude, in excellent agreement with observations. We elaborate new theory benchmarks for experiments and simulations, and put forward a hypothesis on the role of viscous effects in interfacial nonlinear RMI.


title: The effect of time-correlated noise on the Kuramoto model studied via the unified colored noise approximation

abstract: Many natural and social phenomena are characterized by synchronization. The Kuramoto model, taking into account the basic ingredients for observing synchronized states, allows to study mathematically synchronization in a simplified but nontrivial picture. Here we study how a noise that is correlated on a finite time-scale $\tau$ impacts the ability of the Kuramoto model to achieve synchronization. We develop an approximated theory that allows to compute the critical coupling constant $k_c$ as a function of the correlation time $\tau$. We obtain that that $k_c(\tau)$ decreases as $\tau$ increases indicating that time-correlated noise promotes synchronization. Moreover, we show that theory describes qualitatively well the degree of synchronization near $k_c$ obtained numerically. Finally, we show that, independently on the value of $\tau$, the curves of the order parameter versus $k$ scale on the same master curve even at values of $k$ very far from $k_c$.


title: The fate of the Konishi multiplet in the \beta-deformed Quantum Spectral Curve

abstract: We investigate the solution space of the $\beta$-deformed Quantum Spectral Curve by studying a sample of solutions corresponding to single-trace operators that in the undeformed theory belong to the Konishi multiplet. We discuss how to set the precise boundary conditions for the leading Q-system for a given state, how to solve it, and how to build perturbative corrections to the $\mathbf{P}\mu$-system. We confirm and add several loop orders to known results in the literature.


title: A local contact systolic inequality in dimension three

abstract: Let $\alpha$ be a contact form on a connected closed three-manifold $\Sigma$. The systolic ratio of $\alpha$ is defined as $\rho_{\mathrm{sys}}(\alpha):=\tfrac{1}{\mathrm{Vol}(\alpha)}T_{\min}(\alpha)^2$, where $T_{\min}(\alpha)$ and $\mathrm{Vol}(\alpha)$ denote the minimal period of periodic Reeb orbits and the contact volume. The form $\alpha$ is said to be Zoll if its Reeb flow generates a free $S^1$-action on $\Sigma$. We prove that the set of Zoll contact forms on $\Sigma$ locally maximises the systolic ratio in the $C^3$-topology. More precisely, we show that every Zoll form $\alpha_*$ admits a $C^3$-neighbourhood $\mathcal U$ in the space of contact forms such that, for every $\alpha\in\mathcal U$, there holds $\rho_{\mathrm{sys}}(\alpha)\leq \rho_{\mathrm{sys}}(\alpha_*)$ with equality if and only if $\alpha$ is Zoll.


title: Dissipation and plastic deformation in collisions between metallic nanoparticles

abstract: Collisions between amorphous Fe nanoparticles were studied using molecular-dynamics simulation. For head-on collisions of nanoparticles with radii $R =$ 1.4 nm, $R =$ 5.2 nm, and $R =$ 11 nm, sticking was observed at all simulated velocities. The results were compared to the description provided by the JKR model. It was found that strong disagreement exists between the predictions of JKR and the results of the molecular-dynamics simulation due to the presence of additional dissipative processes which strengthen sticking behavior. First, it is demonstrated that very strong dissipation into atomic vibrations occurs during the collision. The dissipation is strong enough to prevent significant rebound of the nanoparticles. Additionally, the morphology of the adhered nanoparticles includes a ``neck'' that increases in radius with increasing collision velocity which results in amplified irreversibility and adhesion. Approximate calculation of the stress during the collision indicates that stress levels are well above typical yield stress values even for low velocity collisions, consistent with the observation of plastic deformation. Furthermore, it is shown that for nanoparticles with $R \leq$ 11 nm, the dominance of surface attraction results in large effective collision velocities and plastic deformation. By obtaining scaling relations for computed quantities, predictions are made for larger nanoparticles up to $R$ $\sim$ 1 $\mu$m. This work provides a new perspective on collisional dissipation and adhesion with an important connection to the modern understanding of tribology and friction.


title: Cosmological Shapes of Higher-Spin Gravity

abstract: We explore non-Gaussian features of a massless spin-two field in the Vasiliev theory of higher-spin gravity. The theory contains an infinite tower of interacting gauge fields with increasing spin, and admits four-dimensional asymptotically de Sitter configurations. Using a recent proposal for calculating late-time quantum correlations in Vasiliev theory, we provide an exact formula for the tensor non-Gaussianities of the massless spin-two graviton field. By general symmetry considerations, we relate our result to that produced by a tree-level calculation in a gravitational theory containing an Einstein term and a term cubic in the Weyl tensor. The relative coefficient between the two terms is calculated explicitly, exhibiting a significant contribution from the Weyl cubed term. We discuss potential cosmological implications of our results.


title: Equivalent definitions for (degree one) Cameron-Liebler classes of generators in finite classical polar spaces

abstract: In this article, we study degree one Cameron-Liebler sets of generators in all finite classical polar spaces, which is a particular type of a Cameron-Liebler set of generators in this polar space, [9]. These degree one Cameron-Liebler sets are defined similar to the Boolean degree one functions, [15]. We summarize the equivalent definitions for these sets and give a classification result for the degree one Cameron-Liebler sets in the polar spaces W(5,q) and Q(6,q).


title: Semidefinite programming bounds for error-correcting codes

abstract: This chapter is written for the forthcoming book "A Concise Encyclopedia of Coding Theory" (CRC press), edited by W. Cary Huffman, Jon-Lark Kim, and Patrick Sol\'e. This book will collect short but foundational articles, emphasizing definitions, examples, exhaustive references, and basic facts on the model of the Handbook of Finite Fields. The target audience of the Encyclopedia is upper level undergraduates and graduate students.


title: Gd(III)-Gd(III) RIDME for In-Cell EPR Distance Determination

abstract: In-cell distance determination by EPR reveals essential structural information about biomacromolecules under native conditions. We demonstrate that the pulsed EPR technique RIDME (relaxation induced dipolar modulation enhancement) can be utilized for such distance determination. The performance of in-cell RIDME has been assessed at Q band using stiff molecular rulers labelled with Gd(III)-PyMTA tags and microinjected into X. laevis oocytes. The overtone coefficients are determined to be the same for protonated aqueous solutions and inside cells. As compared to in-cell DEER (double electron-electron resonance, also abbreviated as PELDOR), in-cell RIDME features approximately 5 times larger modulation depth and does not show artificial broadening in the distance distributions due to the effect of pseudo-secular terms.


title: Central Limit Theorems for Moving Average Random Fields with Non-Random and Random Sampling On Lattices

abstract: For a L\'evy basis $L$ on $\mathbb{R}^d$ and a suitable kernel function $f:\mathbb{R}^d \to \mathbb{R}$, consider the continuous spatial moving average field $X=(X_t)_{t\in \mathbb{R}^d}$ defined by $X_t = \int_{\mathbb{R}^d} f(t-s) \, dL(s)$. Based on observations on finite subsets $\Gamma_n$ of $\mathbb{Z}^d$, we obtain central limit theorems for the sample mean and the sample autocovariance function of this process. We allow sequences $(\Gamma_n)$ of deterministic subsets of $\mathbb{Z}^d$ and of random subsets of $\mathbb{Z}^d$. The results generalise existing results for time indexed stochastic processes (i.e. $d=1$) to random fields with arbitrary spatial dimension $d$, and additionally allow for random sampling. The results are applied to obtain a consistent and asymptotically normal estimator of $\mu>0$ in the stochastic partial differential equation $(\mu - \Delta) X = dL$ in dimension 3, where $L$ is L\'evy noise.


title: Nuclear collective dynamics in the lattice Hamiltonian Vlasov method

abstract: The lattice Hamiltonian method is developed for solving the Vlasov equation with nuclear mean-field based on the Skyrme pseudopotential up to next-to-next-to-next-to leading order. The ground states of nuclei are obtained through varying the total energy with respect to the density distribution of nucleons. Owing to the self-consistent treatment of initial nuclear ground state and the exact energy conservation in the lattice Hamiltonian method, the present framework of solving the Vlasov equation exhibits very stable nuclear ground state evolution. As a first application of the new lattice Hamiltonian Vlasov method, we explore the iso-scalar giant monopole and iso-vector giant dipole modes of finite nuclei. The obtained results are shown to be comparable to that from random-phase approximation and consistent with the experimental data, indicating the capability of the present method in dealing with the long-time near-equilibrium nuclear dynamics.


title: Cortical thickness and functional networks modules by cortical lobes

abstract: This study aims to investigate topological organization of cortical thickness and functional networks by cortical lobes. First, I demonstrated modular organization of these networks by the cortical surface frontal, temporal, parietal and occipital divisions. Secondly, I mapped the overlapping edges of cortical thickness and functional networks for positive and negative correlations. Finally, I showed that overlapping positive edges map onto within-lobe cortical interactions and negative onto between-lobes interactions.


title: Analysis of fine-tuning measures in models with extended Higgs sectors

abstract: In the literature measures of fine-tuning have been discussed as one of the tools to assess the feasibility of beyond the Standard Model theories. In this paper we focus on two specific measures and investigate what kind of fine-tuning they actually quantify. First we apply both measures to the two Higgs doublet model, for which one can analyze the numerical results in terms of available analytic expressions. After drawing various conclusions about the fine-tuning measures, we investigate a particular left-right symmetric model for which it has been claimed that already at tree-level it suffers from a high amount of fine-tuning. We will reach a different conclusion, although left-right symmetric models may require a modest amount of fine-tuning if phenomenological constraints are imposed. Our analysis shows that the two considered measures can probe different aspects of fine-tuning and are both useful if applied and interpreted in the appropriate way.