title: Tunable room-temperature ferromagnetism in the SiC monolayer

abstract: It is essential to explore two-dimensional (2D) material with magnetic ordering in new generation spintronic devices. Particularly, the seeking of room-temperature 2D ferromagnetic (FM) materials is a hot topic of current research. Here, we study magnetism of the Mn-doped and electron-doped SiC monolayer using first-principle calculations. For the Mn-doped SiC monolayer, we find that either electron or hole could mediate the ferromagnetism in the system and the Curie temperature ($T_C$) can be improved by appropriate carrier doping. The codoping strategy is also discussed on improving $T_C$. The transition between antiferromagnetic and FM phase can be found by strain engineering. The $T_C$ is improved above room temperature (RT) under the strain larger than $0.06$. Moreover, the Mn-doped SiC monolayer develops half-metal at the strain range of $0.05-0.1$. On the other hand, the direct electron doping can induce ferromagnetism due to the van Hove singularity in density of states of the conduction band edge of the SiC monolayer. The $T_C$ is found to be around RT. These fascinating controllable electronic and magnetic properties are desired for spintronic applications.


title: Existence of solutions for fractional p-Kirchhoff type equations with a generalized Choquard nonlinearities

abstract: In this article, we establish the existence of solutions to the fractional $p-$Kirchhoff type equations with a generalized Choquard nonlinearities without assuming the Ambrosetti-Rabinowitz condition.


title: Non-asymptotic bounds for percentiles of independent non-identical random variables

abstract: This note displays an interesting phenomenon for percentiles of independent but non-identical random variables. Let $X_1,\cdots,X_n$ be independent random variables obeying non-identical continuous distributions and $X^{(1)}\geq \cdots\geq X^{(n)}$ be the corresponding order statistics. For any $p\in(0,1)$, we investigate the $100(1-p)$%-th percentile $X^{(pn)}$ and prove non-asymptotic bounds for $X^{(pn)}$. In particular, for a wide class of distributions, we discover an intriguing connection between their median and the harmonic mean of the associated standard deviations. For example, if $X_k\sim\mathcal{N}(0,\sigma_k^2)$ for $k=1,\cdots,n$ and $p=\frac{1}{2}$, we show that its median $\big|{\rm Med}\big(X_1,\cdots,X_n\big)\big|= O_P\Big(n^{1/2}\cdot\big(\sum_{k=1}^n\sigma_k^{-1}\big)^{-1}\Big)$ as long as $\{\sigma_k\}_{k=1}^n$ satisfy certain mild non-dispersion property.


title: Online Static Security Assessment of Power Systems Based on Lasso Algorithm

abstract: As one important means of ensuring secure operation in a power system, the contingency selection and ranking methods need to be more rapid and accurate. A novel method-based least absolute shrinkage and selection operator (Lasso) algorithm is proposed in this paper to apply to online static security assessment (OSSA). The assessment is based on a security index, which is applied to select and screen contingencies. Firstly, the multi-step adaptive Lasso (MSA-Lasso) regression algorithm is introduced based on the regression algorithm, whose predictive performance has an advantage. Then, an OSSA module is proposed to evaluate and select contingencies in different load conditions. In addition, the Lasso algorithm is employed to predict the security index of each power system operation state with the consideration of bus voltages and power flows, according to Newton-Raphson load flow (NRLF) analysis in post-contingency states. Finally, the numerical results of applying the proposed approach to the IEEE 14-bus, 118-bus, and 300-bus test systems demonstrate the accuracy and rapidity of OSSA.


title: Features of word similarity

abstract: In this theoretical note we compare different types of computational models of word similarity and association in their ability to predict a set of about 900 rating data. Using regression and predictive modeling tools (neural net, decision tree) the performance of a total of 28 models using different combinations of both surface and semantic word features is evaluated. The results present evidence for the hypothesis that word similarity ratings are based on more than only semantic relatedness. The limited cross-validated performance of the models asks for the development of psychological process models of the word similarity rating task.


title: Transverse spin transfer to $\Lambda$ and $\bar{\Lambda}$ hyperons in polarized proton-proton collisions at $\sqrt{s}=200\,\mathrm{GeV}$

abstract: The transverse spin transfer from polarized protons to $\Lambda$ and $\bar{\Lambda}$ hyperons is expected to provide sensitivity to the transversity distribution of the nucleon and to the transversely polarized fragmentation functions. We report the first measurement of the transverse spin transfer to $\Lambda$ and $\bar{\Lambda}$ along the polarization direction of the fragmenting quark, $D_\mathrm{TT}$, in transversely polarized proton-proton collisions at $\sqrt{s}=200\,\mathrm{GeV}$ with the STAR detector at RHIC. The data correspond to an integrated luminosity of $18\,\mathrm{pb}^{-1}$ and cover the pseudorapidity range $\left|\eta\right| < 1.2$ and transverse momenta $p_{\mathrm{T}}$ up to $8\,\mathrm{GeV}/c$. The dependence on $p_\mathrm{T}$ and $\eta$ are presented. The $D_\mathrm{TT}$ results are found to be comparable with a model prediction, and are also consistent with zero within uncertainties.


title: Non-Cartan Mordell-Weil lattices of rational elliptic surfaces and heterotic/F-theory compactifications

abstract: The Mordell-Weil lattices (MW lattices) associated to rational elliptic surfaces are classified into 74 types. Among them, there are cases in which the MW lattice is none of the weight lattices of simple Lie algebras or direct sums thereof. We study how such "non-Cartan MW lattices" are realized in the six-dimensional heterotic/F-theory compactifications. In this paper, we focus on non-Cartan MW lattices that are torsion free and whose associated singularity lattices are sublattices of $A_7$. For the heterotic string compactification, a non-Cartan MW lattice yields an instanton gauge group $H$ with one or more $U(1)$ group(s). We give a method for computing massless spectra via the index theorem and show that the $U(1)$ instanton number is limited to be a multiple of some particular non-one integer. On the F-theory side, we examine whether we can construct the corresponding threefold geometries, i.e., rational elliptic surface fibrations over $P^1$. Except for some cases, we obtain such geometries for specific distributions of instantons. All the spectrum derived from those geometries completely match with the heterotic results.


title: Polynomial Chaos-Based Adaptive Control for Nonlinear Systems

abstract: In this work we present a nonlinear adaptive suboptimal control strategy for uncertain nonlinear systems. Stochastic parametric uncertainty is dealt with by employing spectral decomposition of the random variables by means of the generalized polynomial chaos expansion. The projection of uncertainty onto the orthogonal polynomial basis, prescribed by the Wiener-Askey scheme, provides a deterministic model from which the control laws are designed. We apply the nonlinear feedback control law to an automatic pilot system for recovering an aircraft with uncertain aerodynamic data from stall while providing acceptable dynamic response. The control law is obtained by approximating the Hamilton-Jacobi-Bellman equation by a pertubational procedure and asymptotic weak stability in probability of the controlled nonlinear system is verified in a deterministic way.


title: Approximate Distribution Matching for Sequence-to-Sequence Learning

abstract: Sequence-to-Sequence models were introduced to tackle many real-life problems like machine translation, summarization, image captioning, etc. The standard optimization algorithms are mainly based on example-to-example matching like maximum likelihood estimation, which is known to suffer from data sparsity problem. Here we present an alternate view to explain sequence-to-sequence learning as a distribution matching problem, where each source or target example is viewed to represent a local latent distribution in the source or target domain. Then, we interpret sequence-to-sequence learning as learning a transductive model to transform the source local latent distributions to match their corresponding target distributions. In our framework, we approximate both the source and target latent distributions with recurrent neural networks (augmenter). During training, the parallel augmenters learn to better approximate the local latent distributions, while the sequence prediction model learns to minimize the KL-divergence of the transformed source distributions and the approximated target distributions. This algorithm can alleviate the data sparsity issues in sequence learning by locally augmenting more unseen data pairs and increasing the model's robustness. Experiments conducted on machine translation and image captioning consistently demonstrate the superiority of our proposed algorithm over the other competing algorithms.


title: Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

abstract: Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.


title: Photoluminescence Blinking beyond Quantum-Confinement: Spatiotemporally Correlated Intermittency over Entire Micron Sized Perovskite Polycrystalline Disks

abstract: Abrupt fluorescence intermittency or blinking is long recognized to be characteristic of single nano-emitters. Extended quantum-confined nanostructures also undergo spatially heterogeneous blinking, however, there is no such precedence in dimensionally unconfined (bulk) materials. Here, we report multi-level blinking of entire individual organo-lead bromide perovskite micro-crystals (volume 0.1-3 micron-cuble) under ambient conditions. Extremely high spatiotemporal correlation (>0.9) in intra-crystal emission intensity fluctuations signifies effective communication amongst photogenerated carriers at distal locations (up to ~4 microns) within each crystal. Fused polycrystalline grains also exhibit this intriguing phenomenon, which is rationalized by correlated and efficient migration of carriers to a few transient non-radiative traps, the nature and population of which determine blinking propensity. Observation of spatiotemporally correlated emission intermittency in bulk semiconductor crystals opens up the possibility to design novel devices involving long range (mesoscopic) electronic communication.


title: Energy-Efficient Massive IoT Shared Spectrum Access over UAV-enabled Cellular Networks

abstract: Data aggregation has become an emerging paradigm to support massive Internet-of-things (IoT), a new and critical use case for fifth-generation new radio (5G-NR). Indeed, data aggregators can complement cellular base stations and process IoT traffic to reduce network congestion. In this paper, we consider using mobile data aggregators, e.g., drones, that collect IoT traffic and aggregate them to the network. Specifically, we first discuss how the spectrum can be shared between cellular users (UEs) and IoT devices in the presence of drones, proposing a time-division duplexing protocol. We use stochastic geometry to analyze this protocol, comparing it to the standard spectrum sharing and orthogonal allocation protocols. We then formulate a stochastic optimization problem to optimize the nominal IoT transmit power, maximizing the average energy-efficiency (EE) of the IoT device subject to interference constraints to protect UEs. Simulations are presented to validate the theoretical insights and the effectiveness of the proposed protocol. It is shown that using drones, to aggregate IoT traffic, improves the EE of IoT devices, yet the EE degrades as their altitudes increases. Equally important, optimizing the transmit power is critical to further improve the EE, while ensuring fair coexistence with UEs.


title: Remarks on the higher dimensional Suita conjecture

abstract: To study the analog of Suita's conjecture for domains $D \subset \mathbb{C}^n$, $n \ge 2$, B\l ocki introduced the invariant $F^k_D(z)=K_D(z)\lambda\big(I^k_D(z)\big)$, where $K_D(z)$ is the Bergman kernel of $D$ along the diagonal and $\lambda\big(I^k_D(z)\big)$ is the Lebesgue measure of the Kobayashi indicatrix at the point $z$. In this note, we study the behaviour of $F^k_D(z)$ (and other similar invariants using different metrics) on strongly pseudconvex domains and also compute its limiting behaviour explicitly at certain points of decoupled egg domains in $\mathbb{C}^2$.


title: Microscopic theory of electric polarization induced by skyrmionic order in GaV$_{4}$S$_{8}$

abstract: The lacunar spinel GaV$_{4}$S$_{8}$ was recently suggested to be a prototype multiferroic material hosting skyrmion lattice states with a sizeable polarization $\boldsymbol{P}$ coupled to magnetic order. We explain this phenomenon on the microscopic level. On the basis of density functional theory, we construct an effective model describing the behavior of magnetically active electrons in a weakly coupled lattice formed by molecular orbitals of the (V$_{4}$S$_{4}$)$^{5+}$ clusters. By applying superexchange theory combined with the Berry-phase theory for $\boldsymbol{P}$, we derive a compass model relating the energy and polarization change with the directions of spins $\boldsymbol{e}_{i}$ in magnetic bonds. We argue that, although each skyrmion layer is mainly formed by superexchange interactions in the same plane, the spin-dependence of $\boldsymbol{P}$ arises from the stacking misalignment of such planes in the perpendicular direction, which is inherent to the lacunar spinel structure. We predict a strong competition of isotropic, $\sim \boldsymbol{e}_{i}\boldsymbol{e}_{j}$, and antisymmetric, $\sim \boldsymbol{e}_{i} \times \boldsymbol{e}_{j}$, contributions to $\boldsymbol{P}$ that explains the experimentally observed effect.


title: On determinants identity minus Hankel matrix

abstract: In this note, we study the asymptotics of the determinant $\det(I_N - \beta H_N)$ for $N$ large, where $H_N$ is the $N\times N$ restriction of a Hankel matrix $H$ with finitely many jump discontinuities in its symbol satisfying $\|H\|\leq 1$. Moreover, we assume $\beta\in\mathbb C$ with $|\beta|<1$ and $I_N$ denotes the identity matrix. We determine the first order asymtoptics as $N\to\infty$ of such determinants and show that they exhibit power-like asymptotic behaviour, with exponent depending on the height of the jumps. For example, for the $N \times N$ truncation of the Hilbert matrix $\mathbf{H}$ with matrix elements $\pi^{-1}(j+k+1)^{-1}$, where $j,k\in \mathbb Z_+$ we obtain $$ \log \det(I_N - \beta \mathbf{H}_N) = -\frac{\log N}{2\pi^2} \big(\pi\arcsin(\beta)+\arcsin^2(\beta)+o(1)\big),\qquad N\to\infty. $$


title: Novel electronic states realized by the interplay between Li diffusion and Co$^{\rm 3+}$$/$Co$^{\rm 4+}$ charge ordering in Li$_{\rm x}$CoO$_2$

abstract: Measurements of dc magnetization ($M$) and electrical resistivity ($\rho$) have been carried out as a function of temperature ($T$) for layered oxide Li$_{\rm x}$CoO$_2$ (0.51$\leq$$x$$\leq$1.0) using single crystal specimens. After slow cool of the specimens down to 10 K, both of the $M$($T$) and $\rho$($T$) curves are found to exhibit a clear anomaly due to the occurrence of Co$^{\rm 3+}$$/$Co$^{\rm 4+}$ charge ordering (CO) at $T_{\rm S}$$\sim$155 K for 0.6$\leq$$x$$\leq$0.98 (at $T_{\rm S}$$\sim$180$-$190 K for 0.5$\leq$$x$$\leq$0.55). After rapid cool of the specimens, additional anomalies are observed related to the onset of Li diffusion at $T_{\rm F1}$$\sim$370 K and$/$or $T_{\rm F2}$=120$-$130 K. Due to phase mixing with compositions of nearly LiCoO$_2$ and Li$_{\rm 2/3}$CoO$_2$, the specimens with 0.7$\lesssim$$x$$\lesssim$0.9 show anomlies both at $T_{\rm F1}$ and $T_{\rm F2}$. For 0.6$\lesssim$$x$$\lesssim$0.9, the resistivity measured after rapid cool is found to be fairly larger than that measured after slow cool below $T_{\rm S}$. The enhanced resistivity can be explained by the scenario that disordered Co$^{\rm 3+}$$/$Co$^{\rm 4+}$ arrangements, which have been observed and revealed to have an insulating electronic structure contrasting to the regular CO state in the previous scanning tunneling microscopy measurements [Phys. Rev. Lett. {\bf 111}, 126104 (2013)], are realized due to the formation of an amorphous-like structure of Li ions after rapid cool via the interlayer Coulomb coupling. An electronic phase diagram for 0.5$\leq$$x$$\leq$1.0 is proposed.


title: A new approach for rain gush formation associated with ionic wind

abstract: Based upon experimental observation in the laboratory, we propose that ionic wind from corona discharge inside a thundercloud would play an important role in producing a rain gush. A cyclic chain of events inside a super-saturated environment in a thundercloud is proposed, each event enhancing the successive ones until lightning occurs. These successive events are collision between snowflakes and rimers, charge separation, corona discharge, avalanche ionization, ionic wind originating from the positively and negatively charged masses of cloud, vortex motion and turbulence when mixed with the updraft, more collision, more charge separation, stronger corona discharge, and so on. Meanwhile, avalanche ionization would produce more CCN (cloud condensation nuclei) resulting in more precipitation and hence rimers formation in the super-saturated environment. More collision in the buoyant turbulence would lead to more fusion of droplets and the formation of larger rimers. The cyclic processes would repeat themselves until the electric field between the two oppositely charged masses of cloud was strong enough to induce a breakdown. The latter would create a sudden short circuit between the two charged masses of cloud neutralizing the charges. There would be no more ionic wind, hence, much less buoyant turbulence. The updraft alone would not be sufficiently strong to support larger rimers which would fall down 'suddenly' to the earth surface as a rain gush.


title: The quantum Zeno and anti-Zeno effects with non-selective measurements

abstract: In studies of the quantum Zeno and anti-Zeno effects, it is usual to consider rapid projective measurements with equal time intervals being performed on the system to check whether or not the system is in the initial state. These projective measurements are selective measurements in the sense that the measurement results are read out and only the case where all the measurement results correspond to the initial state is considered in the analysis of the effect of the measurements. In this paper, we extend such a treatment to consider the effect of repeated non-selective measurements - only the final measurement is required to correspond to the initial state, while we do not know the results of the intermediate measurements. We present a general formalism to derive the effective decay rate of the initial quantum state with such non-selective measurements. Importantly, we show that there is a difference between using non-selective measurements and the usual approach of considering only selective measurements only if we go beyond the weak system-environment coupling regime in models other than the usual population decay models. As such, we then apply our formalism to investigate the quantum Zeno and anti-Zeno effects for three exactly solvable system-environment models: a single two-level system undergoing dephasing, a single two-level system interacting with an environment of two-level systems, and a large spin undergoing dephasing. Our results show that the quantum Zeno and anti-Zeno effects in the presence of non-selective measurements can differ very significantly as compared to the repeated selective measurement scenario.


title: Reinforcement Learning for Relation Classification from Noisy Data

abstract: Existing relation classification methods that rely on distant supervision assume that a bag of sentences mentioning an entity pair are all describing a relation for the entity pair. Such methods, performing classification at the bag level, cannot identify the mapping between a relation and a sentence, and largely suffers from the noisy labeling problem. In this paper, we propose a novel model for relation classification at the sentence level from noisy data. The model has two modules: an instance selector and a relation classifier. The instance selector chooses high-quality sentences with reinforcement learning and feeds the selected sentences into the relation classifier, and the relation classifier makes sentence level prediction and provides rewards to the instance selector. The two modules are trained jointly to optimize the instance selection and relation classification processes. Experiment results show that our model can deal with the noise of data effectively and obtains better performance for relation classification at the sentence level.


title: Topological axion states in magnetic insulator MnBi$_2$Te$_4$ with the quantized magnetoelectric effect

abstract: Topological states of quantum matter have attracted great attention in condensed matter physics and materials science. The study of time-reversal-invariant (TRI) topological states in quantum materials has made tremendous progress in both theories and experiments. As a great success, thousands of TRI topological materials are predicted through sweeping search. Richer exotic phenomena are expected to appear in magnetic topological materials because of varied magnetic configurations, but this study falls much behind due to the complex magnetic structures and transitions. Here, we predict the tetradymite-type compound MnBi$_2$Te$_4$ and its related materials host interesting magnetic topological states. The magnetic ground state of MnBi$_2$Te$_4$ is an antiferromagnetic phase which leads to an antiferromagetic topological insulator state with a large topologically non-trivial energy gap ($\sim$0.2~eV). It is the parent state for the axion state, which has gapped bulk and surface states, and quantized topological magnetoelectric effect. The ferromagnetic phase of MnBi$_2$Te$_4$ leads to an ideal minimal type-II Weyl semimetal with two Weyl points accompanied by one hole-type and one electron-type Fermi pocket at the Fermi level, which has never been discovered elsewhere. We further present a simple and unified continuum model to capture the salient topological features of this kind of materials.


title: An Enhanced SCMA Detector Enabled by Deep Neural Network

abstract: In this paper, we propose a learning approach for sparse code multiple access (SCMA) signal detection by using a deep neural network via unfolding the procedure of message passing algorithm (MPA). The MPA can be converted to a sparsely connected neural network if we treat the weights as the parameters of a neural network. The neural network can be trained off-line and then deployed for online detection. By further refining the network weights corresponding to the edges of a factor graph, the proposed method achieves a better performance. Moreover, the deep neural network based detection is a computationally efficient since highly paralleled computations in the network are enabled in emerging Artificial Intelligence (AI) chips.


title: Eddy magnetization from the chiral Barnett effect

abstract: We discuss the spin, the angular momentum, and the magnetic moment of rotating chiral fermions using a kinetic theory. We find that, in addition to the chiral vortical contribution along the rotation axis, finite circular spin polarization is induced by the spin-momentum correlation of chiral fermions, which is canceled by a change in the orbital angular momentum. We point out that the eddy magnetic moment is nonvanishing due to the $g$-factors, exhibiting the chiral Barnett effect.


title: De-distorting ionospheric effects in the image plane

abstract: The Earth's ionosphere refracts radio waves incident on an interferometer, resulting in shifts to the measured positions of radio sources. We present a method to smoothly remove these shifts and restore sources to their reference positions, in both the catalogue and image domains. The method is applicable to instruments and ionospheric weather such that all antennas see the same ionosphere. The method is generalisable to repairing any sparsely-sampled vector field distortion to some input data. The code is available under the Academic Free License (https://opensource.org/licenses/AFL-3.0) from https://github.com/nhurleywalker/fits_warp


title: Magnetic domain walls in strain-patterned ultrathin films

abstract: We present a comparison of the characteristics of the magnetic domain walls in an atomic monolayer of Co on Pt(111) and a Ni/Fe atomic bilayer on Ir(111), based on spin-polarized scanning tunneling microscopy measurements. In both cases, the films exhibit a roughly triangular dislocation line pattern created by epitaxial strain relief, as well as out-of-plane ferromagnetic order. Domains with opposite magnetization are separated by domain walls with a unique rotational sense, demonstrating the important role of the Dzyaloshinskii-Moriya interaction induced by the Co/Pt and the Fe/Ir interfaces. The domain walls in Co/Pt(111) are straight and usually found in geometrical constrictions of the film, where they can minimize their length. In contrast, the domain walls in Ni/Fe/Ir(111) follow complicated paths, which can be correlated to the structural triangular pattern. The comparison between the two systems shows that the structural patterns, despite their similarity, have a different impact on the domain walls. In the Co/Pt(111) case, the magnetic state is not influenced by the dislocation line network, in contrast to the Ni/Fe/Ir(111) system in which the formation of the walls is favored at specific positions of the structural pattern.


title: Linear complexity of generalized cyclotomic sequences of period $2p^{m}$

abstract: In this paper, we construct two generalized cyclotomic binary sequences of period $2p^{m}$ based on the generalized cyclotomy and compute their linear complexity, showing that they are of high linear complexity when $m\geq 2$.


title: The Operadic Nerve, Relative Nerve, and the Grothendieck Construction

abstract: We relate the relative nerve $\mathrm{N}_f(\mathcal{D})$ of a diagram of simplicial sets $f \colon \mathcal{D} \to \mathsf{sSet}$ with the Grothendieck construction $\mathsf{Gr} F$ of a simplicial functor $F \colon \mathcal{D} \to \mathsf{sCat}$ in the case where $f = \mathrm{N} F$. We further show that any strict monoidal simplicial category $\mathcal{C}$ gives rise to a functor $\mathcal{C}^\bullet \colon \Delta^\mathrm{op} \to \mathsf{sCat}$, and that the relative nerve of $\mathrm{N} \mathcal{C}^\bullet$ is the operadic nerve $\mathrm{N}^\otimes(\mathcal{C})$. Finally, we show that all the above constructions commute with appropriately defined opposite functors.


title: Deep demosaicking for multispectral filter arrays

abstract: We propose a novel demosaicking method for multispectral filter arrays based on a deep convolutional neural network. The proposed method first interpolates mosaicked multispectral images utilizing a bilinear approach, then applies a residual network to initial demosaicked images. The residual network consists of various three-dimensional convolutional layers and a rectified linear unit for describing the features of a multispectral data cube. Experimental results reveal that the proposed method outperforms conventional demosaicking methods.


title: Semi-infinite highest weight categories

abstract: We develop axiomatics of highest weight categories and quasi-hereditary algebras in order to incorporate two semi-infinite situations which are in Ringel duality with each other; the underlying posets are either upper finite or lower finite. We also consider various more general sorts of stratified categories. In the upper finite cases, we give an alternative characterization of these categories in terms of based quasi-hereditary algebras and based stratified algebras, which are certain locally unital algebras possessing triangular bases.


title: A Jointly Learned Context-Aware Place of Interest Embedding for Trip Recommendations

abstract: Trip recommendation is an important location-based service that helps relieve users from the time and efforts for trip planning. It aims to recommend a sequence of places of interest (POIs) for a user to visit that maximizes the user's satisfaction. When adding a POI to a recommended trip, it is essential to understand the context of the recommendation, including the POI popularity, other POIs co-occurring in the trip, and the preferences of the user. These contextual factors are learned separately in existing studies, while in reality, they impact jointly on a user's choice of a POI to visit. In this study, we propose a POI embedding model to jointly learn the impact of these contextual factors. We call the learned POI embedding a context-aware POI embedding. To showcase the effectiveness of this embedding, we apply it to generate trip recommendations given a user and a time budget. We propose two trip recommendation algorithms based on our context-aware POI embedding. The first algorithm finds the exact optimal trip by transforming and solving the trip recommendation problem as an integer linear programming problem. To achieve a high computation efficiency, the second algorithm finds a heuristically optimal trip based on adaptive large neighborhood search. We perform extensive experiments on real datasets. The results show that our proposed algorithms consistently outperform state-of-the-art algorithms in trip recommendation quality, with an advantage of up to 43% in F1-score.


title: Decision fusion with multiple spatial supports by conditional random fields

abstract: Classification of remotely sensed images into land cover or land use is highly dependent on geographical information at least at two levels. First, land cover classes are observed in a spatially smooth domain separated by sharp region boundaries. Second, land classes and observation scale are also tightly intertwined: they tend to be consistent within areas of homogeneous appearance, or regions, in the sense that all pixels within a roof should be classified as roof, independently on the spatial support used for the classification. In this paper, we follow these two observations and encode them as priors in an energy minimization framework based on conditional random fields (CRFs), where classification results obtained at pixel and region levels are probabilistically fused. The aim is to enforce the final maps to be consistent not only in their own spatial supports (pixel and region) but also across supports, i.e., by getting the predictions on the pixel lattice and on the set of regions to agree. To this end, we define an energy function with three terms: 1) a data term for the individual elements in each support (support-specific nodes); 2) spatial regularization terms in a neighborhood for each of the supports (support-specific edges); and 3) a regularization term between individual pixels and the region containing each of them (intersupports edges). We utilize these priors in a unified energy minimization problem that can be optimized by standard solvers. The proposed 2LCRF model consists of a CRF defined over a bipartite graph, i.e., two interconnected layers within a single graph accounting for interlattice connections.


title: Three-dimensional finite temperature Z$_2$ gauge theory with tensor network scheme

abstract: We apply a tensor network scheme to finite temperature Z$_2$ gauge theory in 2+1 dimensions. Finite size scaling analysis with the spatial extension up to $N_{\sigma}=4096$ at the temporal extension of $N_\tau=2,3,5$ allows us to determine the transition temperature and the critical exponent $\nu$ at high level of precision, which shows the consistency with the Svetitsky-Yaffe conjecture.


title: A Numerical Exploration of the Spherically Symmetric SU(2) Einstein-Yang-Mills Equations

abstract: The Einstein-Yang-Mills equations are the source of many interesting solutions within general relativity, including families of particle-like and black hole solutions, and critical phenomena of more than one type. These solutions, discovered in the last thirty years, all assume a restricted form for the Yang-Mills gauge potential known as the "magnetic" ansatz. In this thesis we relax that assumption and investigate the most general solutions of the Einstein-Yang-Mills system assuming spherically symmetry, a Yang-Mills gauge group of SU(2), and zero cosmological constant. We proceed primarily by numerically integrating the equations and find new static solutions, for both regular and black hole boundary conditions, which are not asymptotically flat, and attempt to classify the possible static behaviours. We develop a code to solve the dynamic equations that uses a novel adaptive mesh refinement algorithm making full use of double-null coordinates. We find that the "type II" critical behaviour in the general case exhibits non-universal critical solutions, in contrast to the magnetic case and to all previously observed type II critical behaviour.


title: Measuring Motional Dynamics of (CH$_3$)$_2$ NH$_2^+$ in the Perovskite-Like Metal--Organic Framework [(CH$_3$)$_2$ NH$_2$][Zn(HCOO)$_3$]: The Value of Low-Frequency Electron Paramagnetic Resonance

abstract: Dimethylammonium zinc formate (DMAZnF) is the precursor for a large family of multiferroics, materials which display co-existing magnetic and dielectric ordering. However, the mechanism underlying these orderings remains unclear. While it is generally believed that the dielectric transition is related to the freezing of the order-disorder dynamics of the dimethylammonium (DMA+) cation, no quantitative data on this motion are available. We surmise that this is due to the fact that the timescale of this cationic motion is on the borderline of the timescales of experimental techniques used in earlier reports. Using multifrequency EPR, we find that the timescale of this motion is ~ 5 x 10 -9 s. Thus, S-band (4 GHz) EPR spectroscopy is presented as the technique of choice for studying these motional dynamics. This work highlights the value of the lower-frequency end of EPR spectroscopy. The data are interpreted using DFT calculations and provide direct evidence for the motional freezing model of the ferroelectric transition in these metal-organic frameworks with the ABX3 perovskite-like architecture.


title: The XDEM Multi-physics and Multi-scale Simulation Technology: Review on DEM-CFD Coupling, Methodology and Engineering Applications

abstract: The XDEM multi-physics and multi-scale simulation platform roots in the Ex- tended Discrete Element Method (XDEM) and is being developed at the In- stitute of Computational Engineering at the University of Luxembourg. The platform is an advanced multi- physics simulation technology that combines flexibility and versatility to establish the next generation of multi-physics and multi-scale simulation tools. For this purpose the simulation framework relies on coupling various predictive tools based on both an Eulerian and Lagrangian approach. Eulerian approaches represent the wide field of continuum models while the Lagrange approach is perfectly suited to characterise discrete phases. Thus, continuum models include classical simulation tools such as Computa- tional Fluid Dynamics (CFD) or Finite Element Analysis (FEA) while an ex- tended configuration of the classical Discrete Element Method (DEM) addresses the discrete e.g. particulate phase. Apart from predicting the trajectories of individual particles, XDEM extends the application to estimating the thermo- dynamic state of each particle by advanced and optimised algorithms. The thermodynamic state may include temperature and species distributions due to chemical reaction and external heat sources. Hence, coupling these extended features with either CFD or FEA opens up a wide range of applications as diverse as pharmaceutical industry e.g. drug production, agriculture food and processing industry, mining, construction and agricultural machinery, metals manufacturing, energy production and systems biology.


title: Residual spin susceptibility in the spin-triplet, orbital-singlet model

abstract: Nuclear magnetic resonance (NMR) and Knight shift measurements are critical tools in the identification of spin-triplet superconductors. We discuss the effects of spin orbit coupling on the Knight shift and susceptibilities for a variety of spin triplet multi-orbital gap functions with orbital-singlet character and compare their responses to "traditional" single band spin-triplet ($p_x+ip_y$) superconductors. We observe a non-negligible residual spin-susceptibility at low temperature.


title: Note about Stiefel-Whitney classes on real Bott manifolds

abstract: Real Bott manifolds is a class of flat manifolds with holonomy group $\mathbb Z_2^k$ of diagonal type. In this paper we want to show how we can compute even Stiefel - Whitney classes on real Bott manifolds. This paper is an answer to the question of professor Masuda if is it possible to extend A. G\k{a}sior "Spin-structures on real Bott manifolds" (J. Korean Math. Soc. {\bf 54}, (2017), no. 2, 507 - 516) and compute any Stiefel-Whitney classes for real Bott manifolds. It also extends results of A. G\k{a}sior, A. Szczepa\'nski "Flat manifolds with holonomy group $Z_2^k$ of diagonal type" (Osaka J. Math. {\bf 51} (2014), 1015 - 1025).


title: CRIRES high-resolution infrared spectroscopy of the long-period Cepheid l Car

abstract: The dynamical structure of the atmosphere of Cepheids has been well studied in the optical. Several authors have found very interesting spectral features in the J band, but little data have been secured beyond 1.6 um. However, such observations can probe different radial velocities and line asymmetry regimes, and are able to provide crucial insights into stellar physics. Our goal was to investigate the infrared line-forming region in the K-band domain, and its impact on the projection factor and the k-term of Cepheids. We secured CRIRES observations for the long-period Cepheid l Car, with a focus on the unblended spectral line NaI2208.969 nm. We measured the corresponding radial velocities (by using the first moment method) and the line asymmetries (by using the bi-Gaussian method). These quantities are compared to the HARPS visible spectra we previously obtained on l Car. The optical and infrared radial velocity curves show the same amplitude (only about 3% of difference), with a slight radial velocity shift of about 0.5 +/- 0.3 km s^-1 between the two curves. Around the minimum radius (phase ~ 0.9) the visible radial velocity curve is found in advance compared to the infrared one (phase lag), which is consistent with an infrared line forming higher in the atmosphere (compared to the visible line) and with a compression wave moving from the bottom to the top of the atmosphere during maximum outward velocity. The asymmetry of the K-band line is also found to be significantly different from that of the optical line.


title: Structure of the $\Xi_b(6227)^-$ Resonance

abstract: We explore the recently observed $\Xi_b(6227)^{-}$ resonance to fix its quantum numbers. To this end, we consider various possible scenarios: It can be considered as either 1P/2S excitations of the $\Xi_b^-$ and $\Xi_b'(5935)^{-}$ ground state baryons with spin-$\frac{1}{2}$ or 1P/2S excitations of the ground state $\Xi_b(5955)^{-}$ with spin-$\frac{3}{2}$. We calculate the masses of the possible angular-orbital $1P$ and $2S$ excited states corresponding to each channel employing the QCD sum rule technique. It is seen that all the obtained masses are in agreement with the experimentally observed value, implying that the mass calculations are not enough to distinguish the quantum numbers of the state under question. Therefore, we extend the analysis to investigate the possible decays of the same excited states into $\Lambda_b^0 K^-$ and $\Xi_b^-\pi$. Using the light cone QCD sum rule method, we calculate the corresponding strong coupling constants, which are used to extract the decay widths of the modes under consideration. Our results on decay widths indicate that the $\Xi_b(6227)^{-}$ is $1P$ angular-orbital excited state of the $\Xi_b(5955)^{-}$ baryon with quantum numbers $J^P=\frac{3}{2}^{-}$.


title: Circular spectropolarimetric sensing of higher plant and algal chloroplast structural variations

abstract: Photosynthetic eukaryotes show a remarkable variability in photosynthesis, including large differences in light harvesting proteins and pigment composition. In vivo circular spectropolarimetry enables us to probe the molecular architecture of photosynthesis in a non-invasive and non-destructive way and, as such, can offer a wealth of physiological and structural information. In the present study we have measured the circular polarizance of several multicellular green, red and brown algae and higher plants, which show large variations in circular spectropolarimetric signals with differences in both spectral shape and magnitude. Many of the algae display spectral characteristics not previously reported, indicating a larger variation in molecular organization than previously assumed. As the strengths of these signals vary by three orders of magnitude, these results also have important implications in terms of detectability for the use of circular polarization as a signature of life.


title: Finite-dimensional complex manifolds on commutative Banach algebras and continuous families of compact complex manifolds

abstract: An n-dimensional complex manifold is a manifold by biholomorphic mappings between open sets of the finite direct product of the complex number field. On the other hand, when A is a commutative Banach algebra, Lorch gave a definition that an A-valued function on an open set of A is holomorphic. The definition of a holomorphic function by Lorch can be straightforwardly generalized to an A-valued function on an open set of the finite direct product of A. Therefore, a manifold modeled on the finite direct product of A (an n-dimensional A-manifold) is easily defined. However, in my opinion, it seems that so many nontrivial examples were not known (including the case of n=1, that is, Riemann surfaces). By the way, if X is a compact Hausdorff space, then the algebra C(X) of all complex valued continuous functions on X is the most basic example of a commutative Banach algebra (furthermore, a commutative C*-algebra). In this note, we see that if the set of all continuous cross sections of a continuous family M of compact complex manifolds (a topological deformation M of compact complex analytic structures) on X is denoted by G(M), then the structure of a C(X)-manifold modeled on the C(X)-modules of all continuous cross sections of complex vector bundles on X is introduced into G(M). Therefore, especially, if X is contractible, then G(M) is a finite-dimensional C(X)-manifold.


title: A beam of ultrahigh energy particles of the same magnetic rigidity

abstract: Two EAS arrays during a day have recorded 3 particles with energies above 30 EeV arriving from the same sky region. Two events were registered by the Yakutsk array and one - by the Telescope Array. Two Yakutsk events were estimated to be same energy and the Telescope Array shower's energy was almost 2 times higher. This indicates the same magnetic rigidity of all three particles, if the charge differs by a factor of 2. The relatively short sequence of all three events and their "monochromaticity" in rigidity can be due to the magnetic separation of particles in the acceleration process and propagation of the beam.


title: Guaranteed and computable error bounds for approximations constructed by an iterative decoupling of the Biot problem

abstract: The paper is concerned with guaranteed a posteriori error estimates for a class of evolutionary problems related to poroelastic media governed by the quasi-static linear Biot equations. The system is decoupled employing the fixed-stress split scheme, which leads to a semi-discrete system solved iteratively. The error bounds are derived by combining a posteriori estimates for contractive mappings with those of the functional type for elliptic partial differential equations. The estimates are applicable for any approximation in the admissible functional space and are independent of the discretization method. They are fully computable, do not contain mesh dependent constants, and provide reliable global estimates of the error measured in the energy norm. Moreover, they suggest efficient error indicators for the distribution of local errors, which can be used in adaptive procedures.


title: Analogy between thermal emission of nano objects and Hawking's radiation

abstract: We analyze in this work some analogies between thermal emission of nano objects and Hawking's radiation. We first focus on the famous expression of the black hole radiating temperature derived by Hawking in 1974 and consider the case of thermal emission of a small aperture made into a cavity (Ideal Blackbody). We show that an expression very similar to Hawking's temperature determines a temperature below which an aperture in a cavity cannot be considered as standard blackbody radiating like T^4. Hawking's radiation therefore appear as a radiation at a typical wavelength which is of the size of the horizon radius. In a second part, we make the analogy between the emission of particle-anti particle pairs near the black hole horizon and the scattering and coupling of thermally populated evanescent waves by a nano objects. We show here again that a temperature similar to the Hawking temperature determines the regimes where the scattering occur or where it is negligible.


title: Measurement of the $\gamma^{\star}\gamma^{\star} \to \eta'$ transition form factor

abstract: We study the process $e^+e^-\to e^+e^- \eta'$ in the double-tag mode and measure for the first time the $\gamma^{\star}\gamma^{\star} \to \eta'$ transition form factor $F_{\eta'}(Q_1^2, Q_2^2)$ in the momentum-transfer range \mbox{2 $<Q_1^2, Q_2^2<$ 60 GeV$^2$}. The analysis is based on a data sample corresponding to an integrated luminosity of around 469 fb$^{-1}$ collected at the PEP-II $e^+e^-$ collider with the BaBar detector at center-of-mass energies near 10.6 GeV.


title: Guiding-center motion for electrons in strong laser fields

abstract: We consider the dynamics of electrons in combined strong laser and Coulomb fields. Under a timescale separation condition, we reduce this dynamics to a guiding-center framework. More precisely, we derive a hierarchy of models for the guiding-center dynamics based on averaging over the fast motion of the electron using Lie transforms. The reduced models we obtain describe well the different ionization channels, in particular, the conditions under which an electron is rescattered by the ionic core or is directly ionized. The comparison between these models highlights the models which are best suited for a qualitative and quantitative agreement with the parent dynamics.


title: Hybrid Job-driven Scheduling for Virtual MapReduce Clusters

abstract: It is cost-efficient for a tenant with a limited budget to establish a virtual MapReduce cluster by renting multiple virtual private servers (VPSs) from a VPS provider. To provide an appropriate scheduling scheme for this type of computing environment, we propose in this paper a hybrid job-driven scheduling scheme (JoSS for short) from a tenant's perspective. JoSS provides not only job level scheduling, but also map-task level scheduling and reduce-task level scheduling. JoSS classifies MapReduce jobs based on job scale and job type and designs an appropriate scheduling policy to schedule each class of jobs. The goal is to improve data locality for both map tasks and reduce tasks, avoid job starvation, and improve job execution performance. Two variations of JoSS are further introduced to separately achieve a better map-data locality and a faster task assignment. We conduct extensive experiments to evaluate and compare the two variations with current scheduling algorithms supported by Hadoop. The results show that the two variations outperform the other tested algorithms in terms of map-data locality, reduce-data locality, and network overhead without incurring significant overhead. In addition, the two variations are separately suitable for different MapReduce-workload scenarios and provide the best job performance among all tested algorithms.


title: Characterization of a novel pixelated Silicon Drift Detector (PixDD) for high-throughput X-ray astrophysics

abstract: Multi-pixel fast silicon detectors represent the enabling technology for the next generation of space-borne experiments devoted to high-resolution spectral-timing studies of low-flux compact cosmic sources. Several imaging detectors based on frame-integration have been developed as focal plane devices for X-ray space-borne missions but, when coupled to large-area concentrator X-ray optics, these detectors are affected by strong pile-up and dead-time effects, thus limiting the time and energy resolution as well as the overall system sensitivity. The current technological gap in the capability to realize pixelated silicon detectors for soft X-rays with fast, photon-by-photon response and nearly Fano-limited energy resolution therefore translates into the unavailability of sparse read-out sensors suitable for high throughput X-ray astronomy applications. In the framework of the ReDSoX Italian collaboration, we developed a new, sparse read-out, pixelated silicon drift detector which operates in the energy range 0.5-15 keV with nearly Fano-limited energy resolution ($\leq$150 eV FWHM @ 6 keV) at room temperature or with moderate cooling ($\sim$0 {\deg}C to +20 {\deg}C). In this paper, we present the design and the laboratory characterization of the first 16-pixel (4$\times$4) drift detector prototype (PixDD), read-out by individual ultra low-noise charge sensitive preamplifiers (SIRIO) and we discuss the future PixDD prototype developments.


title: Using SWISH to realise interactive web based tutorials for logic based languages

abstract: Programming environments have evolved from purely text based to using graphical user interfaces, and now we see a move towards web based interfaces, such as Jupyter. Web based interfaces allow for the creation of interactive documents that consist of text and programs, as well as their output. The output can be rendered using web technology as, e.g., text, tables, charts or graphs. This approach is particularly suitable for capturing data analysis workflows and creating interactive educational material. This article describes SWISH, a web front-end for Prolog that consists of a web server implemented in SWI-Prolog and a client web application written in JavaScript. SWISH provides a web server where multiple users can manipulate and run the same material, and it can be adapted to support Prolog extensions. In this paper we describe the architecture of SWISH, and describe two case studies of extensions of Prolog, namely Probabilistic Logic Programming (PLP) and Logic Production System (LPS), which have used SWISH to provide tutorial sites.


title: On the control of Rydberg state population with realistic femtosecond laser pulses

abstract: We investigate computationally a method for ultrafast preparation of alkali metal atoms in their Rydberg states using a three-dimensional model potential in the single active electron approximation. By optimizing laser pulse shapes that can be generated with modern waveform synthesizers, we propose pulses for controlling the population transfer from the ground state to a preselected set of Rydberg states. Dynamical processes under the optimized pulses are shown to be much more complicated than in the traditional optical two-photon preparation of Rydberg states.


title: Hunting for Statistical Anisotropy in Tensor Modes with B-mode Observations

abstract: We investigate a possibility of constraining statistical anisotropies of the primordial tensor perturbations by using future observations for the Cosmic Microwave Background (CMB) B-mode polarization. By parameterizing a statistically-anisotropic tensor power spectrum as $P_h ({\boldsymbol{k}}) = P_h (k) \sum_n g_n \cos^n \theta_{\boldsymbol{k}}$, where $\theta_{\boldsymbol{k}}$ is an angle of the direction of $\hat{k}={\boldsymbol{k}}/k$ from a preferred direction, we find that it would be possible for future B-mode observations such as CMB-S4 to detect the tensor statistical anisotropy at the level of $g_n \sim {\mathcal O} (0.1)$.


title: On the ascent-descent spectrum

abstract: We establish the various properties as well as diverse relations of the ascent and descent spectra for bounded linear operators. We specially focus on the theory of subspectrum. Furthermore, we construct a new concept of convergence for such spectra.


title: Collective Modes in Excitonic Magnets: Dynamical Mean-Field Study

abstract: We present a dynamical mean-field study of dynamical susceptibilities in two-band Hubbard model. Varying the model parameters we analyze the two-particle excitations in the normal as well as in the ordered phase, an excitonic condensate. The two-particle DMFT spectra in the ordered phase reveal the gapless Goldstone modes arising from spontaneous breaking of continuous symmetries. We also observe gapped Higgs mode, characterized by vanishing of the gap at the phase boundary. Qualitative changes observed in the spin susceptibility can be used as an experimental probe to identify the excitonic condensation.


title: Role Semantics for Better Models of Implicit Discourse Relations

abstract: Predicting the structure of a discourse is challenging because relations between discourse segments are often implicit and thus hard to distinguish computationally. I extend previous work to classify implicit discourse relations by introducing a novel set of features on the level of semantic roles. My results demonstrate that such features are helpful, yielding results competitive with other feature-rich approaches on the PDTB. My main contribution is an analysis of improvements that can be traced back to role-based features, providing insights into why and when role semantics is helpful.


title: Premetric teleparallel theory of gravity and its local and linear constitutive law

abstract: We continue to investigate the premetric teleparallel theory of gravity (TG) with the coframe (tetrad) as gravitational potential. We start from the field equations and a local and linear constitutive law. We create a Tonti diagram of TG in order to disclose the structure of TG. Subsequently we irreducibly decompose the 6th order constitutive tensor under the linear group. Moreover, we construct the most general constitutive tensors from the metric and the totally antisymmetric Levi-Civita symbol, and we demonstrate that they encompass nontrivial axion and skewon type pieces. Using these tools, we derive for TG in the geometric-optics approximation propagating massless spin 0, 1, and 2 waves, including the special case of Einstein's general relativity.


title: Mollweide's formula in teaching trigonometry

abstract: Trigonometry is one of the topics in mathematics that the students in both high school and pre-undergraduate levels need to learn. Generally, the topic covers trigonometric functions, trigonometric equations, trigonometric identities and solving oblique triangles using the Laws of Sines and Cosines. However, when solving the oblique triangles, Mollweide's formula is most likely to be omitted from the discussion. Mollweide's formula--which exhibits a cyclical nature--is particularly useful in checking one's result after solving an oblique triangle since all six components of the triangle are involved. It is interesting to note that proving Mollweide's formula can be performed without words. Furthermore, the Law of Tangents can be derived straightforwardly from this equation. In this article, we revisit Mollweide's formula and provide classroom examples where this equation comes into alive. In addition, we suggest that this seemingly less-known equation is to be included in the mathematics syllabus on the topic of Trigonometry.


title: Multiple multivariate subdivision schemes: matrix and operator approaches

abstract: This paper extends the matrix based approach to the setting of multiple subdivision schemes studied in [Sauer 2012]. Multiple subdivision schemes, in contrast to stationary and non-stationary schemes, allow for level dependent subdivision weights and for level dependent choice of the dilation matrices. The latter property of multiple subdivision makes the standard definition of the transition matrices, crucial ingredient of the matrix approach in the stationary and non-stationary settings, inapplicable. We show how to avoid this obstacle and characterize the convergence of multiple subdivision schemes in terms of the joint spectral radius of certain square matrices derived from subdivision weights. We illustrate our results with several examples.


title: Integration with an Adaptive Harmonic Mean Algorithm

abstract: Numerically estimating the integral of functions in high dimensional spaces is a non-trivial task. A oft-encountered example is the calculation of the marginal likelihood in Bayesian inference, in a context where a sampling algorithm such as a Markov Chain Monte Carlo provides samples of the function. We present an Adaptive Harmonic Mean Integration (AHMI) algorithm. Given samples drawn according to a probability distribution proportional to the function, the algorithm will estimate the integral of the function and the uncertainty of the estimate by applying a harmonic mean estimator to adaptively chosen regions of the parameter space. We describe the algorithm and its mathematical properties, and report the results using it on multiple test cases.


title: A New Probability-one Homotopy Method for Solving Minimum-Time Low-Thrust Orbital Transfer Problems

abstract: Homotopy methods have been widely utilized to solve low-thrust orbital transfer problems, however, it is not guaranteed that the optimal solution can be obtained by the existing homotopy methods. In this paper, a new homotopy method is presented, by which the optimal solution can be found with probability one. Generalized sufficient conditions, which are derived from the parametrized Sard's theorem, are first developed. A new type of probability-one homotopy formulation, which is custom-designed for solving minimum-time low-thrust trajectory optimization problems and satisfies all these sufficient conditions, is then constructed. By tracking the continuous zero curve initiated by an initial problem with known solution, the optimal solution of the original problem is guaranteed to be solved with probability one. Numerical demonstrations in a three-dimensional time-optimal low-thrust orbital transfer problem with 43 revolutions is presented to illustrate the applications of the method.


title: Explicit rates of convergence in the multivariate CLT for nonlinear statistics

abstract: We investigate the multivariate central limit theorem for nonlinear statistics by means of Stein's method and Slepian's smart path interpolation method. Based on certain difference operators in theory of concentration inequalities, we obtain two explicit bounds for the rate of convergence. Applications to Rademacher functionals, the runs and quadratic forms are provided as well.


title: The Co-Terminal Swap Market Model with Bergomi Stochastic Volatility

abstract: In this article, we apply the forward variance modeling approach by L.Bergomi to the co-terminal swap market model. We build an interest rate model for which all the market price changes of hedging instruments, interest rate swaps and European swaptions, are interpreted as the state variable variations, and no diffusion parameter calibration procedure is required. The model provides quite simple profit and loss (PnL) formula, with which we can easily understand where a material PnL trend comes from when it appears, and consider how we should modify the model parameters. The model has high flexibility to control the model dynamics because parameter calibration is unnecessary and the model parameters can be used solely for the purpose of the model dynamics control. With the model, the position management of the exotic interest rate products, e.g. Bermudan swaptions, can be carried out in a more sophisticated and systematic manner. A numerical experiment is performed to show the effectiveness of the approach for a Canary swaption, which is a special form of a Bermudan swaption.


title: Jupiter radio emission induced by Ganymede and consequences for the radio detection of exoplanets

abstract: By analysing a database of 26 years of observations of Jupiter with the Nancay Decameter Array, we unambiguously identify the radio emissions caused by the Ganymede-Jupiter interaction. We study the energetics of these emissions via the distributions of their intensities, duration, and power, and compare them to the energetics of the Io-Jupiter radio emissions. This allows us to demonstrate that the average emitted radio power is proportional to the Poynting flux from the rotating Jupiter's magnetosphere intercepted by the obstacle. We then generalize this result to the radio-magnetic scaling law that appears to apply to all plasma interactions between a magnetized flow and an obstacle, magnetized or not. Extrapolating this scaling law to the parameter range corresponding to hot Jupiters, we predict large radio powers emitted by these objects, that should result in detectable radio flux with new-generation radiotelescopes. Comparing the distributions of the durations of Ganymede-Jupiter and Io-Jupiter emission events also suggests that while the latter results from quasi-permanent Alfven wave excitation by Io, the former likely results from sporadic reconnection between magnetic fields Ganymede and Jupiter, controlled by Jupiter's magnetic field geometry and modulated by its rotation.


title: Independent Low-Rank Matrix Analysis Based on Time-Variant Sub-Gaussian Source Model

abstract: Independent low-rank matrix analysis (ILRMA) is a fast and stable method for blind audio source separation. Conventional ILRMAs assume time-variant (super-)Gaussian source models, which can only represent signals that follow a super-Gaussian distribution. In this paper, we focus on ILRMA based on a generalized Gaussian distribution (GGD-ILRMA) and propose a new type of GGD-ILRMA that adopts a time-variant sub-Gaussian distribution for the source model. By using a new update scheme called generalized iterative projection for homogeneous source models, we obtain a convergence-guaranteed update rule for demixing spatial parameters. In the experimental evaluation, we show the versatility of the proposed method, i.e., the proposed time-variant sub-Gaussian source model can be applied to various types of source signal.


title: A non-local approach to waves of maximal height for the Degasperis-Procesi equation

abstract: We consider the non-local formulation of the Degasperis-Procesi equation $u_t+uu_x+L(\frac{3}{2}u^2)_x=0$, where $L$ is the non-local Fourier multiplier operator with symbol $m(\xi)=(1+\xi^2)^{-1}$. We show that all $L^\infty$, pointwise travelling-wave solutions are bounded above by the wave-speed and that if the maximal height is achieved they are peaked at those points, otherwise they are smooth. For sufficiently small periods we find the highest, peaked, travelling-wave solution as the limiting case at the end of the main bifurcation curve of $P$-periodic solutions. The results imply that the Degasperis-Procesi equation does not admit cuspon solutions.


title: The Slow Invariant Manifold of the Lorenz--Krishnamurthy Model

abstract: During this last decades, several attempts to construct slow invariant manifold of the Lorenz-Krishnamurthy five-mode model of slow-fast interactions in the atmosphere have been made by various authors. Unfortunately, as in the case of many two-time scales singularly perturbed dynamical systems the various asymptotic procedures involved for such a construction diverge. So, it seems that till now only the first-order and third-order approximations of this slow manifold have been analytically obtained. While using the Flow Curvature Method we show in this work that one can provide the eighteenth-order approximation of the slow manifold of the generalized Lorenz-Krishnamurthy model and the thirteenth-order approximation of the "conservative" Lorenz-Krishnamurthy model. The invariance of each slow manifold is then established according to Darboux invariance theorem.


title: Capture into first-order resonances and long-term stability of pairs of equal-mass planets

abstract: Massive planets form within the lifetime of protoplanetary disks and undergo orbital migration due to planet-disk interactions. When the first planet reaches the inner edge of the disk its migration stops and the second planet is locked in resonance. We detail the resonant trapping comparing semi-analytical formulae and numerical simulations in the case of two equal-mass coplanar planets trapped in first order resonances. We describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using expanded and non-expanded Hamiltonians. We show that during convergent migration the planets evolve along these families of equilibria. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that the resonant configurations achieved by migration can become unstable after the removal of the disk. We probe the stability of the resonances as a function of planetary mass. To do this, we fictitiously increase the masses of the planets, adiabatically maintaining the low-amplitude libration regime until instability. We discuss two hypotheses for instability: low-order secondary resonance of the libration frequency with a fast synodic frequency, and minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low-amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter still play the decisive role in the destabilisation of resonant pairs. We show that, as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.


title: Global N-Body Simulation of Galactic Spiral Arms

abstract: The origin of galactic spiral arms is one of fundamental problems in astrophysics. Based on the local analysis Toomre (1981) proposed the swing amplification mechanism in which the self-gravity forms spiral arms as leading waves of stars rotate to trailing ones due to galactic shear. The structure of spiral arms is characterized by their number and pitch angle. We perform global $N$-body simulations of spiral galaxies to investigate the dependence of the spiral structure on disk parameters and compare the simulation results with the swing amplification model. We find that the spiral structure in the $N$-body simulations agrees well with that predicted by the swing amplification for the wide range of parameters. The pitch angle decreases with increasing the shear rate and is independent of the disk mass fraction. The number of spiral arms decreases with both increasing the shear rate and the disk mass fraction. If the disk mass fraction is fixed, the pitch angle increases with the number of spiral arms.


title: Bloch-like energy oscillations

abstract: We identify a new type of periodic evolution that appears in driven quantum systems. Provided that the instantaneous (adiabatic) energies are equidistant we show how such systems can be mapped to (time-dependent) tilted single-band lattice models. Having established this mapping, the dynamics can be understood in terms of Bloch oscillations in the instantaneous energy basis. In our lattice model the site-localized states are the adiabatic ones, and the Bloch oscillations manifest as a periodic repopulation among these states, or equivalently a periodic change in the system's instantaneous energy. Our predictions are confirmed by considering two different models: a driven harmonic oscillator and a Landau-Zener grid model. Both models indeed show convincing, or even perfect, oscillations. To strengthen the link between our energy Bloch oscillations and the original spatial Bloch oscillations we add a random disorder that breaks the translational invariance of the spectrum. This verifies that the oscillating evolution breaks down and instead turns into a ballistic spreading.


title: A note on the complex SYK model and warped CFTs

abstract: We discuss the connections between the complex SYK model at the conformal limit and warped conformal field theories. Both theories have an $SL(2,R) \times U(1)$ global symmetry. We present comparisons on symmetries, correlation functions, the effective action and the entropy formula. We also use modular covariance to reinterpret results in the complex SYK model.


title: A Generalized Open Quantum System Approach for the Electron Paramagnetic Resonance of Magnetic Atoms

abstract: A recent experimental breakthrough allowed to probe electronic parametric resonance of a single magnetic atom in an STM setup. The results present intriguing features, such as an asymmetric lineshape and unusually large ratio of the decoherence and decay rates, which defy standard approaches using the conventional Bloch equations. To address these issues we employ novel generalized Bloch equations, together with proper microscopic modeling of the magnetic adatom, and show how all the experimental features can naturally be accounted for. The proposed approach may also be useful in treating any future similar experiments, as well as next generation hybrid quantum devices.


title: Conformally symmetric triangular lattices and discrete $\vartheta$-conformal maps

abstract: Two immersed triangulations in the plane with the same combinatorics are considered as preimage and image of a discrete immersion $F$. We compare the cross-ratios $Q$ and $q$ of corresponding pairs of adjacent triangles in the two triangulations. If for every pair the arguments of these cross-ratios (i.e. intersection angles of circumcircles) agree, $F$ is a discrete conformal map based on circle patterns. Similarly, if for every pair the absolute values of the corresponding cross-ratios $Q$ and $q$ (i.e. length cross-ratios) agree, the two triangulations are discrete conformally equivalent. We introduce a new notion, discrete $\vartheta$-conformal maps, which interpolates between these two known definitions of discrete conformality for planar triangulations. We prove that there exists an associated variational principle. In particular, discrete $\vartheta$-conformal maps are unique minimizers of a locally defined convex functional ${\cal F}_\vartheta$ in suitable variables. Furthermore, we study conformally symmetric triangular lattices which contain examples of discrete $\vartheta$-conformal maps.


title: Towards Machine Learning-Based Optimal HAS

abstract: Mobile video consumption is increasing and sophisticated video quality adaptation strategies are required to deal with mobile throughput fluctuations. These adaptation strategies have to keep the switching frequency low, the average quality high and prevent stalling occurrences to ensure customer satisfaction. This paper proposes a novel methodology for the design of machine learning-based adaptation logics named HASBRAIN. Furthermore, the performance of a trained neural network against two algorithms from the literature is evaluated. We first use a modified existing optimization formulation to calculate optimal adaptation paths with a minimum number of quality switches for a wide range of videos and for challenging mobile throughput patterns. Afterwards we use the resulting optimal adaptation paths to train and compare different machine learning models. The evaluation shows that an artificial neural network-based model can reach a high average quality with a low number of switches in the mobile scenario. The proposed methodology is general enough to be extended for further designs of machine learning-based algorithms and the provided model can be deployed in on-demand streaming scenarios or be further refined using reward-based mechanisms such as reinforcement learning. All tools, models and datasets created during the work are provided as open-source software.


title: Game Theory Meets Network Security: A Tutorial at ACM CCS

abstract: The increasingly pervasive connectivity of today's information systems brings up new challenges to security. Traditional security has accomplished a long way toward protecting well-defined goals such as confidentiality, integrity, availability, and authenticity. However, with the growing sophistication of the attacks and the complexity of the system, the protection using traditional methods could be cost-prohibitive. A new perspective and a new theoretical foundation are needed to understand security from a strategic and decision-making perspective. Game theory provides a natural framework to capture the adversarial and defensive interactions between an attacker and a defender. It provides a quantitative assessment of security, prediction of security outcomes, and a mechanism design tool that can enable security-by-design and reverse the attacker's advantage. This tutorial provides an overview of diverse methodologies from game theory that includes games of incomplete information, dynamic games, mechanism design theory to offer a modern theoretic underpinning of a science of cybersecurity. The tutorial will also discuss open problems and research challenges that the CCS community can address and contribute with an objective to build a multidisciplinary bridge between cybersecurity, economics, game and decision theory.


title: Minimal covers of hypergraphs

abstract: For a hypergraph $H=(V,\mathcal E)$, a subfamily $\mathcal C\subseteq \mathcal E$ is called a cover of the hypergraph if $\bigcup\mathcal C=\bigcup\mathcal E$. A cover $\mathcal C$ is called minimal if each cover $\mathcal D\subseteq\mathcal C$ of the hypergraph $H$ coincides with $\mathcal C$. We prove that for a hypergraph $H$ the following conditions are equivalent: (i) each countable subhypergraph of $H$ has a minimal cover; (ii) each non-empty subhypergraph of $H$ has a maximal edge; (iii) $H$ contains no isomorphic copy of the hypergraph $(\omega,\omega)$. This characterization implies that a countable hypergraph $(V,\mathcal E)$ has a minimal cover if every infinite set $I\subseteq V$ contains a finite subset $F\subseteq I$ such that the family of edges $\mathcal E_F:=\{E\in\mathcal E:F\subseteq E\}$ is finite. Also we prove that a hypergraph $(V,\mathcal E)$ has a minimal cover if $\sup\{|E|:E\in\mathcal E\}<\omega$ or for every $v\in V$ the family $\mathcal E_v:=\{E\in\mathcal E:v\in E\}$ is finite.


title: Self-Paced Multi-Task Clustering

abstract: Multi-task clustering (MTC) has attracted a lot of research attentions in machine learning due to its ability in utilizing the relationship among different tasks. Despite the success of traditional MTC models, they are either easy to stuck into local optima, or sensitive to outliers and noisy data. To alleviate these problems, we propose a novel self-paced multi-task clustering (SPMTC) paradigm. In detail, SPMTC progressively selects data examples to train a series of MTC models with increasing complexity, thus highly decreases the risk of trapping into poor local optima. Furthermore, to reduce the negative influence of outliers and noisy data, we design a soft version of SPMTC to further improve the clustering performance. The corresponding SPMTC framework can be easily solved by an alternating optimization method. The proposed model is guaranteed to converge and experiments on real data sets have demonstrated its promising results compared with state-of-the-art multi-task clustering methods.


title: Topological quantum quench dynamics carrying arbitrary Hopf and second-Chern numbers

abstract: A quantum quench is a nonequilibrium dynamics governed by the unitary evolution. We propose a two-band model whose quench dynamics is characterized by an arbitrary Hopf number belonging to the homotopy group $\pi _{3}(S^{2})=\mathbb{Z}$. When we quench a system from an insulator with the Chern number $C_{i}\in \pi _{2}(S^{2})=\mathbb{Z}$ to another insulator with the Chern number $C_{f} $, the preimage of the Hamiltonian vector forms links having the Hopf number $C_{f}-C_{i}$. We also investigate a quantum-quench dynamics for a four-band model carrying an arbitrary second-Chern number $N\in \pi _{4}(S^{4})=\mathbb{Z}$, which can be realized by quenching a three-dimensional topological insulator having the three-dimensional winding number $N\in \pi _{3}(S^{3})=\mathbb{Z}$.


title: The Open Energy Modelling Framework (oemof) - A new approach to facilitate open science in energy system modelling

abstract: Energy system models have become indispensable tools for planning future energy systems by providing insights into different development trajectories. However, sustainable systems with high shares of renewable energy are characterized by growing cross-sectoral interdependencies and decentralized structures. To capture important properties of increasingly complex energy systems, sophisticated and flexible modelling tools are needed. At the same time, open science is becoming increasingly important in energy system modelling. This paper presents the Open Energy Modelling Framework (oemof) as a novel approach to energy system modelling, representation and analysis. The framework provides a toolbox to construct comprehensive energy system models and has been published open source under a free licence. Through collaborative development based on open processes, the framework supports a maximum level of participation, transparency and open science principles in energy system modelling. Based on a generic graph-based description of energy systems, it is well-suited to flexibly model complex cross-sectoral systems and incorporate various modelling approaches. This makes the framework a multi-purpose modelling environment for modelling and analyzing different systems at scales ranging from urban to transnational.


title: Proceedings Combined 25th International Workshop on Expressiveness in Concurrency and 15th Workshop on Structural Operational Semantics

abstract: This volume contains the proceedings of the Combined 25th International Workshop on Expressiveness in Concurrency and the 15th Workshop on Structural Operational Semantics (EXPRESS/SOS 2018), which was held on September 3, 2018, in Beijing, China, as an affiliated workshop of CONCUR 2018, the 29th International Conference on Concurrency Theory. The EXPRESS workshops aim at bringing together researchers interested in the expressiveness of various formal systems and semantic notions, particularly in the field of concurrency. Their focus has traditionally been on the comparison between programming concepts (such as concurrent, functional, imperative, logic and object-oriented programming) and between mathematical models of computation (such as process algebras, Petri nets, event structures, modal logics, and rewrite systems) on the basis of their relative expressive power. The SOS workshops aim at being a forum for researchers, students and practitioners interested in new developments, and directions for future investigation, in the field of structural operational semantics. One of the specific goals of the SOS workshop series is to establish synergies between the concurrency and programming language communities working on the theory and practice of SOS. Since 2012, the EXPRESS and SOS communities have organized an annual combined EXPRESS/SOS workshop on the expressiveness of mathematical models of computation and the formal semantics of systems and programming concepts.


title: Control of a quantum emitter's bandwidth by managing its reactive power

abstract: Reactive power plays a crucial role in the design of small antenna systems, but its impact on the bandwidth of quantum emitters is typically disregarded. Here, we theoretically demonstrate that there is an intermediate domain between the usual weak and strong coupling regimes where the bandwidth of a quantum emitter is directly related to the dispersion properties of the reactive power. This result emphasizes that reactive power must be understood as an additional degree of freedom in engineering the bandwidth of quantum emitters. We illustrate the applicability of this concept by revisiting typical configurations of quantum emitters coupled to resonant cavities and waveguides. Analysis of the reactive power in these system unveils new functionalities, including the design of effcient but narrowband photon sources, as well as quantum emitters exhibiting a bandwidth narrower than its nonradiative linewidth.


title: Homotopy classes of proper maps out of vector bundles

abstract: In this paper we classify the homotopy classes of proper maps $E\rightarrow \mathbb R^k$, where $E$ is a vector bundle over a compact Hausdorff space. As a corollary we compute the homotopy classes of proper maps $\mathbb R^n\rightarrow \mathbb R^k$. We find a stability range of such maps. We conclude with some remarks on framed submanifolds of non-compact manifolds, the relationship with proper homotopy classes of maps and the Pontryagin-Thom construction.


title: Large deviations and one-sided scaling limit of randomized multicolor box-ball system

abstract: The basic $\kappa$-color box-ball (BBS) system is an integrable cellular automaton on one dimensional lattice whose local states take $\{0,1,\cdots,\kappa \}$ with $0$ regarded as an empty box. The time evolution is defined by a combinatorial rule of quantum group theoretical origin, and the complete set of conserved quantities is given by a $\kappa$-tuple of Young diagrams. In the randomized BBS, a probability distribution on $\{0,1,\cdots,\kappa \}$ to independently fill the consecutive $n$ sites in the initial state induces a highly nontrivial probability measure on the $\kappa$-tuple of those invariant Young diagrams. In a recent work \cite{kuniba2018randomized}, their large $n$ `equilibrium shape' has been determined in terms of Schur polynomials by a Markov chain method and also by a very different approach of Thermodynamic Bethe Ansatz (TBA). In this paper, we establish a large deviations principle for the row lengths of the invariant Young diagrams. As a corollary, they are shown to converge almost surely to the equilibrium shape at an exponential rate. We also refine the TBA analysis and obtain the exact scaling form of the vacancy, the row length and the column multiplicity, which exhibit nontrivial factorization in a one-parameter specialization.


title: Chaotic dynamics of fractional Vallis system for El-Nino

abstract: Vallis proposed a simple model for El-Nino weather phenomenon (referred as Vallis system) by adding an additional parameter p to the Lorenz system. He showed that the chaotic behavior of the Vallis system is related to the El-Nino effect. In the present article we study fractional version of Vallis system in depth. We investigate bifurcations and chaos present in the fractional Vallis system along with the effect of variation of system parameter p. It is observed that the range of values of parameter p for which the Vallis system is chaotic, reduces with the reduction of the fractional order. Further we analyze the incommensurate fractional Vallis system and find the critical value below which the system loses chaos. We also synchronize Vallis system with Bhalekar-Gejji system.


title: A continuity principle equivalent to the monotone $\Pi^{0}_{1}$ fan theorem

abstract: The strong continuity principle reads "every pointwise continuous function from a complete separable metric space to a metric space is uniformly continuous near each compact image." We show that this principle is equivalent to the fan theorem for monotone $\Pi^{0}_{1}$ bars. We work in the context of constructive reverse mathematics.


title: Adiabatic approach to large-amplitude collective motion with the higher-order collective-coordinate operator

abstract: We propose a new set of equations to determine the collective Hamiltonian including the second-order collective-coordinate operator on the basis of the adiabatic self-consistent collective-coordinate (ASCC) theory. We illustrate, with the two-level Lipkin model, that the collective operators including the second-order one are self-consistently determined. We compare the results of the calculations with and without the second-order operator and show that, without the second-order operator, the agreement with the exact solution becomes worse as the excitation energy increases, but that, with the second-order operator included, the exact solution is well reproduced even for highly excited states. We also reconsider which equations one should adopt as the basic equations in the case where only the first-order operator is taken into account, and suggest an alternative set of fundamental equations instead of the conventional ASCC equations. Moreover, we briefly discuss the gauge symmetry of the new basic equations we propose in this paper.


title: On the Moser-Trudinger inequality in complex space

abstract: In this paper we prove the pluricomplex counterpart of the Moser-Trudinger and Sobolev inequalities in complex space. We consider these inequalities for plurisubharmonic functions with finite pluricomplex energy, and we estimate the concerned constants.


title: Under the Hood: Using Diagnostic Classifiers to Investigate and Improve how Language Models Track Agreement Information

abstract: How do neural language models keep track of number agreement between subject and verb? We show that `diagnostic classifiers', trained to predict number from the internal states of a language model, provide a detailed understanding of how, when, and where this information is represented. Moreover, they give us insight into when and where number information is corrupted in cases where the language model ends up making agreement errors. To demonstrate the causal role played by the representations we find, we then use agreement information to influence the course of the LSTM during the processing of difficult sentences. Results from such an intervention reveal a large increase in the language model's accuracy. Together, these results show that diagnostic classifiers give us an unrivalled detailed look into the representation of linguistic information in neural models, and demonstrate that this knowledge can be used to improve their performance.


title: On the Finite Horizon Optimal Switching Problem with Random Lag

abstract: We consider an optimal switching problem with random lag and possibility of component failure. The random lag is modeled by letting the operation mode follow a regime switching Markov-model with transition intensities that depend on the switching mode. The possibility of failures is modeled by having absorbing components. We show existence of an optimal control for the problem by applying a probabilistic technique based on the concept of Snell envelopes.


title: Looking for a Black Cat in a Dark Room: Security Visualization for Cyber-Physical System Design and Analysis

abstract: Today, there is a plethora of software security tools employing visualizations that enable the creation of useful and effective interactive security analyst dashboards. Such dashboards can assist the analyst to understand the data at hand and, consequently, to conceive more targeted preemption and mitigation security strategies. Despite the recent advances, model-based security analysis is lacking tools that employ effective dashboards---to manage potential attack vectors, system components, and requirements. This problem is further exacerbated because model-based security analysis produces significantly larger result spaces than security analysis applied to realized systems---where platform specific information, software versions, and system element dependencies are known. Therefore, there is a need to manage the analysis complexity in model-based security through better visualization techniques. Towards that goal, we propose an interactive security analysis dashboard that provides different views largely centered around the system, its requirements, and its associated attack vector space. This tool makes it possible to start analysis earlier in the system lifecycle. We apply this tool in a significant area of engineering design---the design of cyber-physical systems---where security violations can lead to safety hazards.


title: New constraints on inelastic dark matter from IceCube

abstract: We study the capture and subsequent annihilation of inelastic dark matter (DM) in the Sun, placing constraints on the DM-nucleon scattering cross section from the null result of the IceCube neutrino telescope. We then compare such constraints with exclusion limits on the same cross section that we derive from XENON1T, PICO and CRESST results. We calculate the cross section for inelastic DM-nucleon scattering within an extension of the effective theory of DM-nucleon interactions which applies to the case of inelastic DM models characterised by a mass splitting between the incoming and outgoing DM particle. We find that for values of the mass splitting parameter larger than about 200 keV, neutrino telescopes place limits on the DM-nucleon scattering cross section which are stronger than the ones from current DM direct detection experiments. The exact mass splitting value for which this occurs depends on whether DM thermalises in the Sun or not. This result applies to all DM-nucleon interactions that generate DM-nucleus scattering cross sections which are independent of the nuclear spin, including the "canonical" spin-independent interaction. We explicitly perform our calculations for a DM candidate with mass of 1 TeV, but our conclusions qualitatively also apply to different masses. Furthermore, we find that exclusion limits from IceCube on the coupling constants of this family of spin-independent interactions are more stringent than the ones from a (hypothetical) reanalysis of XENON1T data based on an extended signal region in nuclear recoil energy. Our results should be taken into account in global analyses of inelastic DM models.


title: Automated shape differentiation in the Unified Form Language

abstract: We discuss automating the calculation of weak shape derivatives in the Unified Form Language (Aln{\ae}s et al., ACM Trans. Math. Softw., 2014) by introducing an appropriate additional step in the pullback from physical to reference space that computes G\^ateaux derivatives with respect to the coordinate field. We illustrate the ease of use with several examples.


title: The Forward-Backward-Forward Method from continuous and discrete perspective for pseudo-monotone variational inequalities in Hilbert spaces

abstract: Tseng's forward-backward-forward algorithm is a valuable alternative for Korpelevich's extragradient method when solving variational inequalities over a convex and closed set governed by monotone and Lipschitz continuous operators, as it requires in every step only one projection operation. However, it is well-known that Korpelevich's method converges and can therefore be used also for solving variational inequalities governed by pseudo-monotone and Lipschitz continuous operators. In this paper, we first associate to a pseudo-monotone variational inequality a forward-backward-forward dynamical system and carry out an asymptotic analysis for the generated trajectories. The explicit time discretization of this system results into Tseng's forward-backward-forward algorithm with relaxation parameters, which we prove to converge also when it is applied to pseudo-monotone variational inequalities. In addition, we show that linear convergence is guaranteed under strong pseudo-monotonicity. Numerical experiments are carried out for pseudo-monotone variational inequalities over polyhedral sets and fractional programming problems.


title: A dynamical approach to privacy preserving average consensus

abstract: In this paper we propose a novel method for achieving average consensus in a continuous-time multiagent network while avoiding to disclose the initial states of the individual agents. In order to achieve privacy protection of the state variables, we introduce maps, called output masks, which alter the value of the states before transmitting them. These output masks are local (i.e., implemented independently by each agent), deterministic, time-varying and converging asymptotically to the true state. The resulting masked system is also time-varying and has the original (unmasked) system as its limit system. It is shown in the paper that the masked system has the original average consensus value as a global attractor. However, in order to preserve privacy, it cannot share an equilibrium point with the unmasked system, meaning that in the masked system the global attractor cannot be also stable.


title: Bayesian Multi--Dipole Modeling in the Frequency Domain

abstract: Background: Magneto- and Electro-encephalography record the electromagnetic field generated by neural currents with high temporal frequency and good spatial resolution, and are therefore well suited for source localization in the time and in the frequency domain. In particular, localization of the generators of neural oscillations is very important in the study of cognitive processes in the healthy and in the pathological brain. New method: We introduce the use of a Bayesian multi-dipole localization method in the frequency domain. Given the Fourier Transform of the data at one or multiple frequencies and/or trials, the algorithm approximates numerically the posterior distribution with Monte Carlo techniques. Results: We use synthetic data to show that the proposed method behaves well under a wide range of experimental conditions, including low signal-to-noise ratios and correlated sources. We use dipole clusters to mimic the effect of extended sources. In addition, we test the algorithm on real MEG data to confirm its feasibility. Comparison with existing method(s): Throughout the whole study, DICS (Dynamic Imaging of Coherent Sources) is used systematically as a benchmark. The two methods provide similar general pictures; the posterior distributions of the Bayesian approach contain much richer information at the price of a higher computational cost. Conclusions: The Bayesian method described in this paper represents a reliable approach for localization of multiple dipoles in the frequency domain.


title: Tidal stripping as a mechanism for placing globular clusters on wide orbits: the case of MGC1 in M31

abstract: The globular clusters of large spiral galaxies can be divided into two populations: one which formed in-situ and one which comprises clusters tidally stripped away from other galaxies. In this paper we investigate the contribution to the outer globular cluster population in the M31 galaxy through donation of clusters from dwarf galaxies. We test this numerically by comparing the contribution of globular clusters from simulated encounters to the observed M31 globular cluster population. To constrain our simulations, we specifically investigate the outermost globular cluster in the M31 system, MGC1. The remote location of MGC1 favours the idea of it being captured, however, the cluster is devoid of features associated with tidal interactions. Hence we separate simulations where tidal features are present and where they are hidden. We find that our simulated encounters can place clusters on MGC1-like orbits. In addition, we find that tidal stripping of clusters from dwarf galaxies leaves them on orbits having a range of separations, broadly matching those observed in M31. We find that the specific energies of globular clusters captured by M31 closely matches those of the incoming host dwarf galaxies. Furthermore, in our simulations we find an equal number of accreted clusters on co-rotating and counter-rotating orbits within M31 and use this to infer the fraction of clusters that has been accreted. We find that even close in roughly 50% of the clusters are accreted, whilst this figure increases to over 80% further out.