title: Transient effects at resonant light scattering by particles: Anapole as a storage for "frozen light"?

abstract: Transient effects at excitation of high-Q resonant modes at a leading front of a laser pulse and their decomposition at its trailing edge during the pulse scattering by a particle are discussed. The main attention is paid to the nonradiating (anapole) mode excited in a high-index cylinder. The problem is studied both analytically and numerically. It is shown that the anapole is a steady but essentially nonequilibrium mode, which can exist if and only if the host particle is irradiated by a continuous incident wave. As soon as the incident wave is switched off the anapole becomes a radiating mode and collapses very fast owing to extensive radiative losses. It ruins the hopes to employ the anapole-like excitations as storages for the "frozen light."


title: Elastic scattering of virtual photons via quark loop in Double-Logarithmic Approximation

abstract: We calculate the amplitude of elastic photon-photon scattering via a single quark loop in the Double-Logarithmic Approximation, presuming all external photons to be off-shell and unpolarized.. At the same time we account for the running coupling effects. We consider this process in the forward kinematics at arbitrary relations between t and the external photon virtualities. We obtain explicit expressions for the photon-photon scattering amplitudes in all double-logarithmic kinematic regions. Then we calculate the small-x asymptotics of the obtained amplitudes and compare them with the parent amplitudes, thereby fixing the applicability regions of the asymptorics, i.e. fixing the applicability region for the non-vacuum Reggeons. We find that these Reggeons should be used at x < 10^{-8} only.


title: Majorana Neutrino as Bogoliubov Quasiparticle

abstract: We suggest that the Majorana neutrino should be regarded as a Bogoliubov quasiparticle that is consistently understood only by use of a relativistic analogue of the Bogoliubov transformation. The unitary charge conjugation condition ${\cal C}\psi{\cal C}^{\dagger}=\psi$ is not maintained in the definition of a quantum Majorana fermion from a Weyl fermion. This is remedied by the Bogoliubov transformation accompanying a redefinition of the charge conjugation properties of vacuum, such that a C-noninvariant fermion number violating term (condensate) is converted to a Dirac mass. We also comment on the chiral symmetry of a Majorana fermion; a massless Majorana fermion is invariant under a global chiral transformation $\psi\rightarrow \exp[i\alpha\gamma_{5}]\psi$ and different Majorana fermions are distinguished by different chiral $U(1)$ charge assignments. The reversed process, namely, the definition of a Weyl fermion from a well-defined massless Majorana fermion is also briefly discussed.


title: Quantum transport and the Wigner distribution function for Bloch electrons in spatially homogeneous electric and magnetic fields

abstract: The theory of Bloch electron dynamics for carriers in homogeneous electric and magnetic fields of arbitrary time dependence is developed in the framework of the Liouville equation. The Wigner distribution function (WDF) is determined from the single particle density matrix in the ballistic regime, i.e., collision effects are excluded. The single particle transport equation is established with the electric field described in the vector potential gauge, and the magnetic field is treated in the symmetric gauge. The general approach is to employ the accelerated Bloch state representation (ABR) as a basis so that the dependence upon the electric field, including multiband Zener tunneling, is treated exactly. In the formulation of the WDF, we transform to a new set of variables so that the final WDF is gauge invariant and is expressed explicitly in terms of the position, kinetic momentum, and time. The methodology for developing the WDF is illustrated by deriving the exact WDF equation for free electrons in homogeneous electric and magnetic fields. The methodology is then extended to the case of electrons described by an effective Hamiltonian corresponding to an arbitrary energy band function. In treating the problem of Bloch electrons in a periodic potential, the methodology for deriving the WDF reveals a multiband character due to the inherent nature of the Bloch states. In examining the single-band WDF, it is found that the collisionless WDF equation matches the equivalent Boltzmann transport equation to first order in the magnetic field. These results are necessarily extended to second order in the magnetic field by employing a unitary transformation that diagonalizes the Hamiltonian using the ABR to second order. The work includes a discussion of the multiband WDF transport analysis and the identification of the combined Zener-magnetic field induced tunneling.


title: Comparison of Distances for Supervised Segmentation of White Matter Tractography

abstract: Tractograms are mathematical representations of the main paths of axons within the white matter of the brain, from diffusion MRI data. Such representations are in the form of polylines, called streamlines, and one streamline approximates the common path of tens of thousands of axons. The analysis of tractograms is a task of interest in multiple fields, like neurosurgery and neurology. A basic building block of many pipelines of analysis is the definition of a distance function between streamlines. Multiple distance functions have been proposed in the literature, and different authors use different distances, usually without a specific reason other than invoking the "common practice". To this end, in this work we want to test such common practices, in order to obtain factual reasons for choosing one distance over another. For these reasons, in this work we compare many streamline distance functions available in the literature. We focus on the common task of automatic bundle segmentation and we adopt the recent approach of supervised segmentation from expert-based examples. Using the HCP dataset, we compare several distances obtaining guidelines on the choice of which distance function one should use for supervised bundle segmentation.


title: Weak lensing deflection of three-point correlation functions

abstract: Weak gravitational lensing alters the apparent separations between observed sources, potentially affecting clustering statistics. We derive a general expression for the lensing deflection which is valid for any three-point statistic, and investigate its effect on the three-point clustering correlation function. We find that deflection of the clustering correlation function is greatest at around $z=2$. It is most prominent in regions where the correlation function varies rapidly, in particular at the baryon acoustic oscillation scale where it smooths out the peaks and troughs, reducing the peak-to-trough difference by about 0.1 percent at $z=1$ and around 2.3 percent at $z=10$. The modification due to lensing deflection is typically at the per cent level of the expected errors in a Euclid-like survey and therefore undetectable.


title: Theory and Application on Adaptive-Robust Control of Euler-Lagrange Systems with Linearly Parametrizable Uncertainty Bound

abstract: This work proposes a new adaptive-robust control (ARC) architecture for a class of uncertain Euler-Lagrange (EL) systems where the upper bound of the uncertainty satisfies linear in parameters (LIP) structure. Conventional ARC strategies either require structural knowledge of the system or presume that the overall uncertainties or its time derivative are norm bounded by a constant. Due to unmodelled dynamics and modelling imperfection, true structural knowledge of the system is not always available. Further, for the class of systems under consideration, prior assumption regarding the uncertainties (or its time derivative) being upper bounded by a constant, puts a restriction on states beforehand. Conventional ARC laws invite overestimation-underestimation problem of switching gain. Towards this front, Adaptive Switching-gain based Robust Control (ASRC) is proposed which alleviates the overestimation-underestimation problem of switching gain. Moreover, ASRC avoids any presumption of constant upper bound on the overall uncertainties and can negotiate uncertainties regardless of being linear or nonlinear in parameters. Experimental results of ASRC using a wheeled mobile robot notes improved control performance in comparison to adaptive sliding mode control.


title: Impact of Line-of-Sight and Unequal Spatial Correlation on Uplink MU-MIMO Systems

abstract: Closed-form approximations of the expected per-terminal signal-to-interference-plus-noise-ratio (SINR) and ergodic sum spectral efficiency of a multiuser multiple-input multiple-output system are presented. Our analysis assumes spatially correlated Ricean fading channels with maximum-ratio combining on the uplink. Unlike previous studies, our model accounts for the presence of unequal correlation matrices, unequal Rice factors, as well as unequal link gains to each terminal. The derived approximations lend themselves to useful insights, special cases and demonstrate the aggregate impact of line-of-sight (LoS) and unequal correlation matrices. Numerical results show that while unequal correlation matrices enhance the expected SINR and ergodic sum spectral efficiency, the presence of strong LoS has an opposite effect. Our approximations are general and remain insensitive to changes in the system dimensions, signal-to-noise-ratios, LoS levels and unequal correlation levels.


title: Fast and exact search for the partition with minimal information loss

abstract: In analysis of multi-component complex systems, such as neural systems, identifying groups of units that share similar functionality will aid understanding of the underlying structures of the system. To find such a grouping, it is useful to evaluate to what extent the units of the system are separable. Separability or inseparability can be evaluated by quantifying how much information would be lost if the system were partitioned into subsystems, and the interactions between the subsystems were hypothetically removed. A system of two independent subsystems are completely separable without any loss of information while a system of strongly interacted subsystems cannot be separated without a large loss of information. Among all the possible partitions of a system, the partition that minimizes the loss of information, called the Minimum Information Partition (MIP), can be considered as the optimal partition for characterizing the underlying structures of the system. Although the MIP would reveal novel characteristics of the neural system, an exhaustive search for the MIP is numerically intractable due to the combinatorial explosion of possible partitions. Here, we propose a computationally efficient search to precisely identify the MIP among all possible partitions by exploiting the submodularity of the measure of information loss. Mutual information is one such submodular information loss functions, and is a natural choice for measuring the degree of statistical dependence between paired sets of random variables. By using mutual information as a loss function, we show that the search for MIP can be performed in a practical order of computational time for a reasonably large system. We also demonstrate that MIP search allows for the detection of underlying global structures in a network of nonlinear oscillators.


title: Origin of the abnormal diffusion of transition metal in rutile

abstract: Diffusion of dopants in rutile is the fundamental process that determines the performance of many devices in which rutile is used. The diffusion behavior is known to be highly sample-dependent, but the reasons for this are less well understood. Here, rutile is studied by using first-principles calculations, in order to unravel the microscopic origins of the diverse diffusion behaviors for different doping elements. Anomalous diffusion behavior in the open channel along [001] direction is found: larger atoms include Sc and Zr have lower energy barrier for diffusion via interstitial mechanism, apparently contradicting their known slow diffusion rate. To resolve this, we present an alternate model for the overall diffusion rate of the large-size dopants in rutile, showing that parallel to the [001] channel, it is limited by the formation of the interstitial states, whereas in the direction perpendicular to [001], it proceeds via a kick-out mechanism. By contrast, Co and Ni, prefer to stay in the interstitial site of rutile, and have conventional diffusion with a very small migration barrier in the [001] channel. This leads to highly anisotropic and fast diffusion. The diffusion mechanisms found in the present study can explain the diffusion data measured by experiments, and these findings provide novel understanding for the classic diffusion topic.


title: Summary of the 12th IACHEC Meeting

abstract: We summarize the outcome of the 12th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC), held at the UCLA conference center in Lake Arrowhead (California) in March 2017. 56 scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the X-ray payload inter-calibration, as well as possible ways to improve it. The "Thermal Supernovas Remnant" (SNR) Working Group presented a recently published paper on 1E0102.2-7219 as a calibration standard in the 0.5-1.0 keV band. A new method to measure the high-energy spectrum of the Crab Nebula and pulsar with NuSTAR without using its optics may yield a new absolute flux standard in the 3-7 keV band. A new ACIS contamination model - released with CALDB version 4.7.3 - leads to a significant improvement in modeling the spectral, spatial, and temporal properties of the contaminant. The first calibration results of the scientific payload on board Hitomi confirm the excellent performance of the instruments before the spacecraft operation problems leading to its loss. Finally, the meeting discussed extensively a novel statistic approach to formally identify in which direction the effective areas of different instruments would need to be changed to bring them into concordance. This method could inform future further calibration efforts.


title: Video Salient Object Detection Using Spatiotemporal Deep Features

abstract: This paper presents a method for detecting salient objects in videos where temporal information in addition to spatial information is fully taken into account. Following recent reports on the advantage of deep features over conventional hand-crafted features, we propose a new set of SpatioTemporal Deep (STD) features that utilize local and global contexts over frames. We also propose new SpatioTemporal Conditional Random Field (STCRF) to compute saliency from STD features. STCRF is our extension of CRF to the temporal domain and describes the relationships among neighboring regions both in a frame and over frames. STCRF leads to temporally consistent saliency maps over frames, contributing to the accurate detection of salient objects' boundaries and noise reduction during detection. Our proposed method first segments an input video into multiple scales and then computes a saliency map at each scale level using STD features with STCRF. The final saliency map is computed by fusing saliency maps at different scale levels. Our experiments, using publicly available benchmark datasets, confirm that the proposed method significantly outperforms state-of-the-art methods. We also applied our saliency computation to the video object segmentation task, showing that our method outperforms existing video object segmentation methods.


title: Correlation and Class Based Block Formation for Improved Structured Dictionary Learning

abstract: In recent years, the creation of block-structured dictionary has attracted a lot of interest. Learning such dictionaries involve two step process: block formation and dictionary update. Both these steps are important in producing an effective dictionary. The existing works mostly assume that the block structure is known a priori while learning the dictionary. For finding the unknown block structure given a dictionary commonly sparse agglomerative clustering (SAC) is used. It groups atoms based on their consistency in sparse coding with respect to the unstructured dictionary. This paper explores two innovations towards improving the reconstruction as well as the classification ability achieved with the block-structured dictionary. First, we propose a novel block structuring approach that makes use of the correlation among dictionary atoms. Unlike the SAC approach, which groups diverse atoms, in the proposed approach the blocks are formed by grouping the top most correlated atoms in the dictionary. The proposed block clustering approach is noted to yield significant reductions in redundancy as well as provides a direct control on the block size when compared with the existing SAC-based block structuring. Later, motivated by works using supervised \emph{a priori} known block structure, we also explore the incorporation of class information in the proposed block formation approach to further enhance the classification ability of the block dictionary. For assessment of the reconstruction ability with proposed innovations is done on synthetic data while the classification ability has been evaluated in large variability speaker verification task.


title: Thermodynamically consistent formulation of quasiparticle viscous hydrodynamics

abstract: A novel formulation of second-order relativistic viscous fluid dynamics based on the effective Boltzmann equation for quasi-particles with medium-dependent masses is briefly reviewed.~The evolution equations for the shear and bulk dissipative corrections, and the corresponding transport coefficients, are presented. Resulting approach allows for thermodynamically-consistent incorporation of lattice QCD equation of state in the fluid dynamical framework.


title: Sensitivity analysis of grain surface chemistry to binding energies of ice species

abstract: Advanced telescopes, such as ALMA and JWST, are likely to show that the chemical universe may be even more complex than currently observed, requiring astrochemical modelers to improve their models to account for the impact of new data. However, essential input information for gas-grain models, such as binding energies of molecules to the surface, have been derived experimentally only for a handful of species, leaving hundreds of species with highly uncertain estimates. We present in this paper a systematic study of the effect of uncertainties in the binding energies on an astrochemical two-phase model of a dark molecular cloud, using the rate equations approach. A list of recommended binding energy values based on a literature search of published data is presented. Thousands of simulations of dark cloud models were run, and in each simulation a value for the binding energy of hundreds of species was randomly chosen from a normal distribution. Our results show that the binding energy of H$_{2}$ is critical for the surface chemistry. For high binding energy, H$_{2}$ freezes out on the grain forming an H$_{2}$ ice. This is not physically realistic and we suggest a change in the rate equations. The abundance ranges found are in reasonable agreement with astronomical ice observations. Pearson correlation coefficients revealed that the binding energy of HCO, HNO, CH$_{2}$, and C correlate most strongly with the abundance of dominant ice species. Finally, the formation route of complex organic molecules was found to be sensitive to the branching ratios of H$_{2}$CO hydrogenation.


title: Uniform cohomological expansion of uniformly quasiregular mappings

abstract: Let $f\colon M \to M$ be a uniformly quasiregular self-mapping of a compact, connected, and oriented Riemannian $n$-manifold $M$ without boundary, $n\ge 2$. We show that, for $k \in \{0,\ldots, n\}$, the induced homomorphism $f^* \colon H^k(M;\mathbb{R}) \to H^k(M;\mathbb{R})$, where $H^k(M;\mathbb{R})$ is the $k$:th singular cohomology of $M$, is complex diagonalizable and the eigenvalues of $f^*$ have modulus $(\mathrm{deg}\ f)^{k/n}$. As an application, we obtain a degree restriction for uniformly quasiregular self-mappings of closed manifolds. In the proof of the main theorem, we use a Sobolev--de Rham cohomology based on conformally invariant differential forms and an induced push-forward operator.


title: Reliable uncertainties in indirect measurements

abstract: In this article we present very intuitive, easy to follow, yet mathematically rigorous, approach to the so called data fitting process. Rather than minimizing the distance between measured and simulated data points, we prefer to find such an area in searched parameters' space that generates simulated curve crossing as many acquired experimental points as possible, but at least half of them. Such a task is pretty easy to attack with interval calculations. The problem is, however, that interval calculations operate on guaranteed intervals, that is on pairs of numbers determining minimal and maximal values of measured quantity while in vast majority of cases our measured quantities are expressed rather as a pair of two other numbers: the average value and its standard deviation. Here we propose the combination of interval calculus with basic notions from probability and statistics. This approach makes possible to obtain the results in familiar form as reliable values of searched parameters, their standard deviations, and their correlations as well. There are no assumptions concerning the probability density distributions of experimental values besides the obvious one that their variances are finite. Neither the symmetry of uncertainties of experimental distributions is required (assumed) nor those uncertainties have to be `small.' As a side effect, outliers are quietly and safely ignored, even if numerous.


title: Canonical structure and orthogonality of forces and currents in irreversible Markov chains

abstract: We discuss a canonical structure that provides a unifying description of dynamical large deviations for irreversible finite state Markov chains (continuous time), Onsager theory, and Macroscopic Fluctuation Theory. For Markov chains, this theory involves a non-linear relation between probability currents and their conjugate forces. Within this framework, we show how the forces can be split into two components, which are orthogonal to each other, in a generalised sense. This splitting allows a decomposition of the pathwise rate function into three terms, which have physical interpretations in terms of dissipation and convergence to equilibrium. Similar decompositions hold for rate functions at level 2 and level 2.5. These results clarify how bounds on entropy production and fluctuation theorems emerge from the underlying dynamical rules. We discuss how these results for Markov chains are related to similar structures within Macroscopic Fluctuation Theory, which describes hydrodynamic limits of such microscopic models.c Fluctuation Theory, which describes hydrodynamic limits of such microscopic models.


title: A >99.9%-fidelity quantum-dot spin qubit with coherence limited by charge noise

abstract: Recent advances towards spin-based quantum computation have been primarily fuelled by elaborate isolation from noise sources, such as surrounding nuclear spins and spin-electric susceptibility, to extend spin coherence. In the meanwhile, addressable single-spin and spin-spin manipulations in multiple-qubit systems will necessitate sizable spin-electric coupling. Given background charge fluctuation in nanostructures, however, its compatibility with enhanced coherence should be crucially questioned. Here we realise a single-electron spin qubit with isotopically-enriched phase coherence time (20 microseconds) and fast electrical control speed (up to 30 MHz) mediated by extrinsic spin-electric coupling. Using rapid spin rotations, we reveal that the free-evolution dephasing is caused by charge (instead of conventional magnetic) noise featured by a 1/f spectrum over seven decades of frequency. The qubit nevertheless exhibits superior performance with single-qubit gate fidelities exceeding 99.9% on average. Our work strongly suggests that designing artificial spin-electric coupling with account taken of charge noise is a promising route to large-scale spin-qubit systems having fault-tolerant controllability.


title: A globally convergent filter-trust-region method for large deformation contact problems

abstract: We present a globally convergent method for the solution of frictionless large deformation contact problems for hyperelastic materials. The discretisation uses the mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter-trust-region scheme, and we prove global convergence towards first-order optimal points. The constrained Newton problems are solved robustly and efficiently using a Truncated Non-smooth Newton Multigrid (TNNMG) method with a Monotone Multigrid (MMG) linear correction step. For this we introduce a cheap basis transformation that decouples the contact constraints. Numerical experiments confirm the stability and efficiency of our approach.


title: Observation of 1-D time dependent non-propagating laser plasma structures using Fluid and PIC codes

abstract: The manuscript reports the observation of time dependent localized and non-propagating structures in the coupled laser plasma system through 1-D fluid and PIC simulations. It is reported that such structures form spontaneously as a result of collision amongst certain exact solitonic solutions. They are seen to survive as coherent entities for a long time up to several hundreds of plasma periods. Furthermore, it is shown that such time dependence can also be artificially recreated by significantly disturbing the delicate balance between the radiation and the density fields required for the exact non-propagating solution obtained by Esirkepov et al. [1]. The ensuing time evolution is an interesting interplay between kinetic and field energies of the system. The electrostatic plasma oscillations are coupled with oscillations in the electromagnetic field. The inhomogeneity of the background and the relativistic nature, however, invariably produce large amplitude density perturbations leading to its wave breaking. In the fluid simulations, the signature of wave breaking can be discerned by a drop in the total energy which evidently gets lost to the grid. The PIC simulations are observed to closely follow the fluid simulations till the point of wave breaking. However, the total energy in the case of PIC simulations is seen to remain conserved all throughout the simulations. At the wave breaking the particles are observed to acquire thermal kinetic energy in the case of PIC. Interestingly, even after wave breaking, compact coherent structures with trapped radiation inside high-density peaks, continue to exist both in PIC and fluid simulations. Though the time evolution does not exactly match in the two simulations as it does prior to the process of wave breaking, the time-dependent features exhibited by the remnant structures are characteristically similar.


title: Geometric Embedding of Path and Cycle Graphs in Pseudo-convex Polygons

abstract: Given a graph $ G $ with $ n $ vertices and a set $ S $ of $ n $ points in the plane, a point-set embedding of $ G $ on $ S $ is a planar drawing such that each vertex of $ G $ is mapped to a distinct point of $ S $. A straight-line point-set embedding is a point-set embedding with no edge bends or curves. The point-set embeddability problem is NP-complete, even when $ G $ is $ 2 $-connected and $ 2 $-outerplanar. It has been solved polynomially only for a few classes of planar graphs. Suppose that $ S $ is the set of vertices of a simple polygon. A straight-line polygon embedding of a graph is a straight-line point-set embedding of the graph onto the vertices of the polygon with no crossing between edges of graph and the edges of polygon. In this paper, we present $ O(n) $-time algorithms for polygon embedding of path and cycle graphs in simple convex polygon and same time algorithms for polygon embedding of path and cycle graphs in a large type of simple polygons where $n$ is the number of vertices of the polygon.


title: On dimension and regularity of bundle measures

abstract: In this paper we quantify the notion of antisymmetry of the Fourier transform of certain vector valued measures. The introduced scale is related to the condition appearing in Uchiyama's theorem and is used to give a lower bound for the rectifiable dimension of those measures. Moreover, we obtain an estimate of the lower Hausdorff dimension assuming certain more restrictive version of the $2$-wave cone condition. Results of our considerations can be viewed as an uncertainty-type principle in the following way: it is impossible to simultaneously localize a (bundle) measure and a direction of its Fourier transform on small sets. The investigated class is modeled on the example of gradients of BV functions. The article contains also a theorem concerning regularity: we prove that elements of considered class vanish on 1-purely unrectifiable sets. Our results can be applied to studying the properties of PDE-constrainted measures.


title: A simple approach to distributed observer design for linear systems

abstract: This note investigates the distributed estimation problem for continuous-time linear time-invariant (LTI) systems observed by a network of observers. Each observer in the network has access to only part of the output of the observed system, and communicates with its neighbors according to a given network graph. In this note we recover the known result that if the observed system is observable and the network graph is a strongly connected digraph, then a distributed observer exists. Moreover, the estimation error can be made to converge to zero at any a priori given decay rate. Our approach leads to a relatively straightforward proof of this result, using the mirror of the balanced graph associated with the original network graph. The numerical design of our distributed observer is reduced to solving linear matrix inequalities (LMI's). Each observer in the network has state dimension equal to that of the observed plant.


title: Stacking Characteristics of Close Packed Materials

abstract: It is shown that the enthalpy of any close packed structure for a given element can be characterised as a linear expansion in a set of continuous variables $\alpha_n$ which describe the stacking configuration. This enables us to represent the infinite, discrete set of stacking sequences within a finite, continuous space of the expansion parameters $H_n$. These $H_n$ determine the stable structure and vary continuously in the thermodynamic space of pressure, temperature or composition. The continuity of both spaces means that only transformations between stable structures adjacent in the $H_n$ space are possible, giving the model predictive and well as descriptive ability. We calculate the $H_n$ using density functional theory and interatomic potentials for a range of materials. Some striking results are found: e.g. the Lennard-Jones potential model has 11 possible stable structures and over 50 phase transitions as a function of cutoff range. The very different phase diagrams of Sc, Tl, Y and the lanthanides are understood within a single theory. We find that the widely-reported 9R-fcc transition is not allowed in equilibrium thermodynamics, and in cases where it has been reported in experiments (Li, Na), we show that DFT theory is also unable to predict it.


title: A Linear-time Algorithm for Orthogonal Watchman Route Problem with Minimum Bends

abstract: Given an orthogonal polygon $ P $ with $ n $ vertices, the goal of the watchman route problem is finding a path $ S $ of the minimum length in $ P $ such that every point of the polygon $ P $ is visible from at least one of the point of $ S $. In the other words, in the watchman route problem we must compute a shortest watchman route inside a simple polygon of $ n $ vertices such that all the points interior to the polygon and on its boundary are visible to at least one point on the route. If route and polygon be orthogonal, it is called orthogonal watchman route problem. One of the targets of this problem is finding the orthogonal path with the minimum number of bends as possible. We present a linear-time algorithm for the orthogonal watchman route problem, in which the given polygon is monotone. Our algorithm can be used also for the problem on simple orthogonal polygons $ P $ for which the dual graph induced by the vertical decomposition of $ P $ is a path, which is called path polygon.


title: How Amdahl's low restricts supercomputer applications and building ever bigger supercomputers

abstract: This paper reinterprets Amdahl's law in terms of execution time and applies this simple model to supercomputing. The systematic discussion results in practical formulas enabling to calculate expected running time using large number of processors from experimental runs using low number of processors, delivers a quantitative measure of computational efficiency of supercomputing applications. Through separating non-parallelizable contribution to fractions according to their origin, Amdahl's law enables to derive a timeline for supercomputers (quite similar to Moore's law) and describes why Amdahl's law limits the size of supercomputers. The paper validates that Amdahl's 50-years old model (with slight extension) correctly describes the performance limitations of the present supercomputers. Using some simple and reasonable assumptions, the absolute performance bound of supercomputers is concluded, furthermore that serious enhancements are still necessary to achieve the exaFLOPS dream value.


title: Detection of thermal bridges from thermographic images for the analysis of buildings energy performance

abstract: In this paper, we develop a procedure for the detection of the contours of thermal bridges from thermographic images, in order to study the energetic performance of buildings. Two main steps of the above method are: the enhancement of the thermographic images by an optimized version of the mathematical algorithm for digital image processing based on the theory of sampling Kantorovich operators, and the application of a suitable thresholding based on the analysis of the histogram of the enhanced thermographic images. Finally, an accuracy improvement of the parameter that defines the thermal bridge is obtained.


title: Massively Multilingual Neural Grapheme-to-Phoneme Conversion

abstract: Grapheme-to-phoneme conversion (g2p) is necessary for text-to-speech and automatic speech recognition systems. Most g2p systems are monolingual: they require language-specific data or handcrafting of rules. Such systems are difficult to extend to low resource languages, for which data and handcrafted rules are not available. As an alternative, we present a neural sequence-to-sequence approach to g2p which is trained on spelling--pronunciation pairs in hundreds of languages. The system shares a single encoder and decoder across all languages, allowing it to utilize the intrinsic similarities between different writing systems. We show an 11% improvement in phoneme error rate over an approach based on adapting high-resource monolingual g2p models to low-resource languages. Our model is also much more compact relative to previous approaches.


title: Brain Responses During Robot-Error Observation

abstract: Brain-controlled robots are a promising new type of assistive device for severely impaired persons. Little is however known about how to optimize the interaction of humans and brain-controlled robots. Information about the human's perceived correctness of robot performance might provide a useful teaching signal for adaptive control algorithms and thus help enhancing robot control. Here, we studied whether watching robots perform erroneous vs. correct action elicits differential brain responses that can be decoded from single trials of electroencephalographic (EEG) recordings, and whether brain activity during human-robot interaction is modulated by the robot's visual similarity to a human. To address these topics, we designed two experiments. In experiment I, participants watched a robot arm pour liquid into a cup. The robot performed the action either erroneously or correctly, i.e. it either spilled some liquid or not. In experiment II, participants observed two different types of robots, humanoid and non-humanoid, grabbing a ball. The robots either managed to grab the ball or not. We recorded high-resolution EEG during the observation tasks in both experiments to train a Filter Bank Common Spatial Pattern (FBCSP) pipeline on the multivariate EEG signal and decode for the correctness of the observed action, and for the type of the observed robot. Our findings show that it was possible to decode both correctness and robot type for the majority of participants significantly, although often just slightly, above chance level. Our findings suggest that non-invasive recordings of brain responses elicited when observing robots indeed contain decodable information about the correctness of the robot's action and the type of observed robot.


title: SNR Estimation in Linear Systems with Gaussian Matrices

abstract: This paper proposes a highly accurate algorithm to estimate the signal-to-noise ratio (SNR) for a linear system from a single realization of the received signal. We assume that the linear system has a Gaussian matrix with one sided left correlation. The unknown entries of the signal and the noise are assumed to be independent and identically distributed with zero mean and can be drawn from any distribution. We use the ridge regression function of this linear model in company with tools and techniques adapted from random matrix theory to achieve, in closed form, accurate estimation of the SNR without prior statistical knowledge on the signal or the noise. Simulation results are provided, and show that the proposed method is very accurate.


title: Single-photon-emitting optical centers in diamond fabricated upon Sn implantation

abstract: The fabrication of luminescent defects in single-crystal diamond upon Sn implantation and annealing is reported. The relevant spectral features of the optical centers (emission peaks at 593.5 nm, 620.3 nm, 630.7 nm and 646.7 nm) are attributed to Sn-related defects through the correlation of their photoluminescence (PL) intensity with the implantation fluence. Single Sn-related defects were identified and characterized through the acquisition of their second-order auto-correlation emission functions, by means of Hanbury-Brown-Twiss interferometry. The investigation of their single-photon emission regime as a function of excitation laser power revealed that Sn-related defects are based on three-level systems with a 6 ns radiative decay lifetime. In a fraction of the studied centers, the observation of a blinking PL emission is indicative of the existence of a dark state. Furthermore, absorption dependence from the polarization of the excitation radiation with about 45 percent contrast was measured. This work shed light on the existence of a new optical center associated with a group-IV impurity in diamond, with similar photo-physical properties to the already well-known Si-V and Ge-V emitters, thus providing results of interest from both the fundamental and applicative points of view.


title: Generation of soliton bubbles in a sine-Gordon system with localised inhomogeneities

abstract: Nonlinear wave propagation plays a crucial role in the functioning of many physical and biophysical systems. In the propagation regime, disturbances due to the presence of local external perturbations, such as localised defects or boundary interphase walls have gained great attention. In this article, the complex phenomena that occur when sine-Gordon line solitons collide with localised inhomogeneities are investigated. By a one-dimensional theory, it is shown that internal modes of two-dimensional sine-Gordon solitons can be activated depending on the topological properties of the inhomogeneities. Shape mode instabilities cause the formation of bubble-like and drop-like structures for both stationary and travelling line solitons. It is shown that such structures are formed and stabilised by arrays of localised inhomogeneities distributed in space. Implications of the observed phenomena in physical and biological systems are discussed.


title: Multisymplectic geometry and covariant formalism for mechanical systems with a Lie group as configuration space: application to the Reissner beam

abstract: Many physically important mechanical systems may be described with a Lie group $G$ as configuration space. According to the well-known Noether's theorem, underlying symmetries of the Lie group may be used to considerably reduce the complexity of the problems. However, these reduction techniques, used without care for general problems (waves, field theory), may lead to uncomfortable infinite dimensional spaces. As an alternative, the \emph{covariant} formulation allows to consider a finite dimensional configuration space by increasing the number of independent variables. But the geometric elements needed for reduction, adapted to the specificity of covariant problems which admit Lie groups as configuration space, are difficult to apprehend in the literature (some are even missing to our knowledge). To fill this gap, this article reconsiders the historical geometric construction made by E. Cartan in this particular "covariant Lie group" context. Thus, and it is the main interest of this work, the Poincar\'e-Cartan and multi-symplectic forms are obtained for a principal $G$ bundle. It allows to formulate the \emph{Euler-Poincar\'e equations} of motion and leads to a Noether's current form defined in the dual Lie algebra.


title: Theory of Out-of-Equilibrium Ultrafast Relaxation Dynamics in Metals

abstract: Ultrafast laser excitation of a metal causes correlated, highly nonequilibrium dynamics of electronic and ionic degrees of freedom, which are however only poorly captured by the widely-used two-temperature model. Here we develop an out-of-equilibrium theory that captures the full dynamic evolution of the electronic and phononic populations and provides a microscopic description of the transfer of energy delivered optically into electrons to the lattice. All essential nonequilibrium energy processes, such as electron-phonon and phonon-phonon interactions are taken into account. Moreover, as all required quantities are obtained from first-principles calculations, the model gives an exact description of the relaxation dynamics without the need for fitted parameters. We apply the model to FePt and show that the detailed relaxation is out-of-equilibrium for picoseconds.


title: Multi-modal Factorized Bilinear Pooling with Co-Attention Learning for Visual Question Answering

abstract: Visual question answering (VQA) is challenging because it requires a simultaneous understanding of both the visual content of images and the textual content of questions. The approaches used to represent the images and questions in a fine-grained manner and questions and to fuse these multi-modal features play key roles in performance. Bilinear pooling based models have been shown to outperform traditional linear models for VQA, but their high-dimensional representations and high computational complexity may seriously limit their applicability in practice. For multi-modal feature fusion, here we develop a Multi-modal Factorized Bilinear (MFB) pooling approach to efficiently and effectively combine multi-modal features, which results in superior performance for VQA compared with other bilinear pooling approaches. For fine-grained image and question representation, we develop a co-attention mechanism using an end-to-end deep network architecture to jointly learn both the image and question attentions. Combining the proposed MFB approach with co-attention learning in a new network architecture provides a unified model for VQA. Our experimental results demonstrate that the single MFB with co-attention model achieves new state-of-the-art performance on the real-world VQA dataset. Code available at https://github.com/yuzcccc/mfb.


title: Discrete Time-Crystalline Order in Cavity and Circuit QED Systems

abstract: Discrete time crystals are a recently proposed and experimentally observed out-of-equilibrium dynamical phase of Floquet systems, where the stroboscopic evolution of a local observable repeats itself at an integer multiple of the driving period. We address this issue in a driven-dissipative setup, focusing on the modulated open Dicke model, which can be implemented by cavity or circuit QED systems. In the thermodynamic limit, we employ semiclassical approaches and find rich dynamical phases on top of the discrete time-crystalline order. In a deep quantum regime with few qubits, we find clear signatures of a transient discrete time-crystalline behavior, which is absent in the isolated counterpart. We establish a phenomenology of dissipative discrete time crystals by generalizing the Landau theory of phase transitions to Floquet open systems.

categories: cond-mat.stat-mech cond-mat.quant-gas physics.atom-ph quant-ph


title: Predicate Pairing for Program Verification

abstract: It is well-known that the verification of partial correctness properties of imperative programs can be reduced to the satisfiability problem for constrained Horn clauses (CHCs). However, state-of-the-art solvers for CHCs (CHC solvers) based on predicate abstraction are sometimes unable to verify satisfiability because they look for models that are definable in a given class A of constraints, called A-definable models. We introduce a transformation technique, called Predicate Pairing (PP), which is able, in many interesting cases, to transform a set of clauses into an equisatisfiable set whose satisfiability can be proved by finding an A-definable model, and hence can be effectively verified by CHC solvers. We prove that, under very general conditions on A, the unfold/fold transformation rules preserve the existence of an A-definable model, i.e., if the original clauses have an A-definable model, then the transformed clauses have an A-definable model. The converse does not hold in general, and we provide suitable conditions under which the transformed clauses have an A-definable model iff the original ones have an A-definable model. Then, we present the PP strategy which guides the application of the transformation rules with the objective of deriving a set of clauses whose satisfiability can be proved by looking for A-definable models. PP introduces a new predicate defined by the conjunction of two predicates together with some constraints. We show through some examples that an A-definable model may exist for the new predicate even if it does not exist for its defining atomic conjuncts. We also present some case studies showing that PP plays a crucial role in the verification of relational properties of programs (e.g., program equivalence and non-interference). Finally, we perform an experimental evaluation to assess the effectiveness of PP in increasing the power of CHC solving.


title: Interaction effects in assembly of magnetic nanoparticles

abstract: A specific absorption rate of a dilute assembly of various random clusters of iron oxide nanoparticles in alternating magnetic field has been calculated using Landau- Lifshitz stochastic equation. This approach simultaneously takes into account both the presence of thermal fluctuations of the nanoparticle magnetic moments, and magneto-dipole interaction between the nanoparticles of the clusters. It is shown that for usual 3D clusters the intensity of magneto- dipole interaction is determined mainly by the cluster packing density eta = Np*V/Vcl, where Np is the average number of the particles in the cluster, V is the nanoparticle volume, and Vcl is the cluster volume. The area of the low frequency hysteresis loop and the assembly specific absorption rate have been found to be considerably reduced when the packing density of the clusters increases in the range of 0.005 < eta < 0.4. The dependence of the specific absorption rate on the mean nanoparticle diameter is retained with increase of eta, but becomes less pronounced. For fractal clusters of nanoparticles, which arise in biological media, in addition to considerable reduction of the absorption rate, the absorption maximum is shifted to smaller particle diameters. It is found also that the specific absorption rate of fractal clusters increases appreciably with increase of the thickness of nonmagnetic shells at the nanoparticle surfaces.


title: The Coburn-Simonenko theorem for Toeplitz operators acting between Hardy type subspaces of different Banach function spaces

abstract: Let $\Gamma$ be a rectifiable Jordan curve, let $X$ and $Y$ be two reflexive Banach function spaces over $\Gamma$ such that the Cauchy singular integral operator $S$ is bounded on each of them, and let $M(X,Y)$ denote the space of pointwise multipliers from $X$ to $Y$. Consider the Riesz projection $P=(I+S)/2$, the corresponding Hardy type subspaces $PX$ and $PY$, and the Toeplitz operator $T(a):PX\to PY$ defined by $T(a)f=P(af)$ for a symbol $a\in M(X,Y)$. We show that if $X\hookrightarrow Y$ and $a\in M(X,Y)\setminus\{0\}$, then $T(a)\in\mathcal{L}(PX,PY)$ has a trivial kernel in $PX$ or a dense image in $PY$. In particular, if $1<q\le p<\infty$, $1/r=1/q-1/p$, and $a\in L^{r}\equiv M(L^p,L^q)$ is a nonzero function, then the Toeplitz operator $T(a)$, acting from the Hardy space $H^p$ to the Hardy space $H^q$, has a trivial kernel in $H^p$ or a dense image in $H^q$.


title: LIKWID Monitoring Stack: A flexible framework enabling job specific performance monitoring for the masses

abstract: System monitoring is an established tool to measure the utilization and health of HPC systems. Usually system monitoring infrastructures make no connection to job information and do not utilize hardware performance monitoring (HPM) data. To increase the efficient use of HPC systems automatic and continuous performance monitoring of jobs is an essential component. It can help to identify pathological cases, provides instant performance feedback to the users, offers initial data to judge on the optimization potential of applications and helps to build a statistical foundation about application specific system usage. The LIKWID monitoring stack is a modular framework build on top of the LIKWID tools library. It aims on enabling job specific performance monitoring using HPM data, system metrics and application-level data for small to medium sized commodity clusters. Moreover, it is designed to integrate in existing monitoring infrastructures to speed up the change from pure system monitoring to job-aware monitoring.


title: Diffusion, Influence and Best-Response Dynamics in Networks: An Action Model Approach

abstract: Threshold models and their dynamics may be used to model the spread of `behaviors' in social networks. Regarding such from a modal logical perspective, it is shown how standard update mechanisms may be emulated using action models -- graphs encoding agents' decision rules. A small class of action models capturing the possible sets of decision rules suitable for threshold models is identified, and shown to include models characterizing best-response dynamics of both coordination and anti-coordination games played on graphs.


title: Dilation-commuting operators on power-weighted Orlicz classes

abstract: Let $\Phi_1$ and $\Phi_2$ be nondecreasing functions from $\mathbb{R_+}=(0,\infty)$ onto itself. For $i=1,2$ and $\gamma \in \mathbb{R}$, define the Orlicz class $L_{\Phi_{i}}(\mathbb{R_+})$ to be the set of Lebesgue-measurable functions $f$ on $\mathbb{R_+}$ such that \begin{equation*} \int_{\mathbb{R_+}} \Phi_{i} \left( k|(Tf)(t)| \right) t^{\gamma}dt < \infty \end{equation*} for some $k>0$. Our goal in this paper is to find conditions on $\Phi_1$, $\Phi_2$, $\gamma$ and an operator $T$ so that the assertions \begin{equation} T : L_{\Phi_2,t^{\gamma}}(\mathbb{R_+}) \rightarrow L_{\Phi_1,t^{\gamma}}(\mathbb{R_+}), \tag{I} \end{equation} and \begin{equation}\label{modularA} \int_{\mathbb{R_+}} \Phi_1 \left( |(Tf)(t)| \right)t^{\gamma}dt \leq K \int_{\mathbb{R_+}} \Phi_2 \left( K|f(s)| \right)s^{\gamma}ds, \tag{M} \end{equation} in which $K>0$ is independent of $f$, say, simple on $\mathbb{R_+}$, are equivalent and to then find necessary and sufficient conditions in order that (\ref{modularA}) holds.


title: Convergence analysis of domain decomposition based time integrators for degenerate parabolic equations

abstract: Domain decomposition based time integrators allow the usage of parallel and distributed hardware, making them well-suited for the temporal discretization of parabolic systems, in general, and degenerate parabolic problems, in particular. The latter is due to the degenerate equations' finite speed of propagation. In this study, a rigours convergence analysis is given for such integrators without assuming any restrictive regularity on the solutions or the domains. The analysis is conducted by first deriving a new variational framework for the domain decomposition, which is applicable to the two standard degenerate examples. That is, the $p$-Laplace and the porous medium type vector fields. Secondly, the decomposed vector fields are restricted to the underlying pivot space and the time integration of the parabolic problem can then be interpreted as an operators splitting applied to a dissipative evolution equation. The convergence results then follow by employing elements of the approximation theory for nonlinear semigroups.


title: Valley polarization generated in 3-dimensional group-IV monochalcogenids

abstract: Valleytronics is one of the breaking-through to the technology of electronics, which provides a new degree of freedom to manipulate the properties of electrons. Combining DFT calculations, optical absorption analysis and the linear polarization-resolved transmission measurement together, we report that three pairs of valleys, which feature opposite optical absorption, existing in the 3-dimensional (3D) group-IV monochalcogenids. By applying the linearly-polarized light, valley polarization is successfully generated for the first time in a 3D system, which opens a new direction for the exploration of the valley materials and provides a good platform for the photodetector and valleytronic devices. Valley modulation versus the in-plane strain in GeSe is also studied, suggesting an effective way to get the optimized valleytronic properties.


title: Multivariate Design of Experiments for Engineering Dimensional Analysis

abstract: We consider the design of dimensional analysis experiments when there is more than a single response. We first give a brief overview of dimensional analysis experiments and the dimensional analysis (DA) procedure. The validity of the DA method for univariate responses was established by the Buckingham $\Pi$-Theorem in the early 20th century. We extend the theorem to the multivariate case, develop basic criteria for multivariate design of DA and give guidelines for design construction. Finally, we illustrate the construction of designs for DA experiments for an example involving the design of a heat exchanger.


title: Short-range test of the universality of gravitational constant $G$ at the millimeter scale using a digital image sensor

abstract: The composition dependence of gravitational constant $G$ is measured at the millimeter scale to test the weak equivalence principle, which may be violated at short range through new Yukawa interactions such as the dilaton exchange force. A torsion balance on a turning table with two identical tungsten targets surrounded by two different attractor materials (copper and aluminum) is used to measure gravitational torque by means of digital measurements of a position sensor. Values of the ratios $\tilde{G}_{Al-W}/\tilde{G}_{Cu-W} -1$ and $\tilde{G}_{Cu-W}/G_{N} -1$ were $(0.9 \pm 1.1_{\mathrm{sta}} \pm 4.8_{\mathrm{sys}}) \times 10^{-2}$ and $ (0.2 \pm 0.9_{\mathrm{sta}} \pm 2.1_{\mathrm{sys}}) \times 10^{-2}$ , respectively; these were obtained at a center to center separation of 1.7 cm and surface to surface separation of 4.5 mm between target and attractor, which is consistent with the universality of $G$. A weak equivalence principle (WEP) violation parameter of $\eta_{Al-Cu}(r\sim 1\: \mathrm{cm})=(0.9 \pm 1.1_{\mathrm{sta}} \pm 4.9_{\mathrm{sys}}) \times 10^{-2} $ at the shortest range of around 1 cm was also obtained.


title: Magnetized, Relativistic Jets

abstract: Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydrodynamics. As a first step, we have studied the structure of steady-state configurations of jets by using numerical simulations. We have used a helical field configuration and have changed different relevant parameters that control the way in which the energy flux is distributed in jets (namely, the proportion of the energy flux carried by internal, kinetic or magnetic energy). Our results show significant differences among the different kinds of jets. Finally, we also report on results based on synthetic maps of our simulated jets.


title: Probing sub-GeV mass SIMP dark matter with a low-threshold surface experiment

abstract: Using data from the $\nu$-cleus detector, based on the surface of the Earth, we place constraints on dark matter in the form of Strongly Interacting Massive Particles (SIMPs) which interact with nucleons via nuclear-scale cross sections. For large SIMP-nucleon cross sections the sensitivity of traditional direct dark matter searches using underground experiments is limited by the energy loss experienced by SIMPs, due to scattering with the rock overburden and experimental shielding on their way to the detector apparatus. Hence a surface-based experiment is ideal for a SIMP search, despite the much larger background, resulting from the lack of shielding. We show using data from a recent surface run of a low-threshold cryogenic detector that values of the SIMP-nucleon cross section up to approximately $10^{-27}$ cm$^2$ can be excluded for SIMPs with masses above 100 MeV.


title: Graviton propagator, renormalization scale and black-hole like states

abstract: We study the analytic structure of the resummed graviton propagator, inspired by the possible existence of black hole precursors in its spectrum. We find an infinite number of poles with positive mass, but both positive and negative effective width, and studied their asymptotic behaviour in the infinite sheet Riemann surface. We find that the stability of these precursors depend crucially on the value of the normalisation point scale.


title: A small-scale dynamo in feedback-dominated galaxies - III. Cosmological simulations

abstract: Magnetic fields are widely observed in the Universe in virtually all astrophysical objects, from individual stars to entire galaxies, even in the intergalactic medium, but their specific generation has long been debated. Due to the development of more realistic models of galaxy formation, viable scenarios are emerging to explain cosmic magnetism, thanks to both deeper observations and more efficient and accurate computer simulations. We present here a new cosmological high-resolution zoom-in magnetohydrodynamic (MHD) simulation, using the adaptive mesh refinement (AMR) technique, of a dwarf galaxy with an initially weak and uniform magnetic seed field that is amplified by a small-scale dynamo driven by supernova-induced turbulence. As first structures form from the gravitational collapse of small density fluctuations, the frozen-in magnetic field separates from the cosmic expansion and grows through compression. In a second step, star formation sets in and establishes a strong galactic fountain, self-regulated by supernova explosions. Inside the galaxy, the interstellar medium becomes highly turbulent, dominated by strong supersonic shocks, as demonstrated by the spectral analysis of the gas kinetic energy. In this turbulent environment, the magnetic field is quickly amplified via a small-scale dynamo process and is finally carried out into the circumgalactic medium by a galactic wind. This realistic cosmological simulation explains how initially weak magnetic seed fields can be amplified quickly in early, feedback-dominated galaxies, and predicts, as a consequence of the small scale dynamo process, that high-redshift magnetic fields are likely to be dominated by their small scale components.


title: Origin of superconductivity and latent charge density wave in NbS$_2$

abstract: We elucidate the origin of the phonon-mediated superconductivity in 2$H$-NbS$_2$ using the ab initio anisotropic Migdal-Eliashberg theory including Coulomb interactions. We demonstrate that superconductivity is associated with Fermi surface hot spots exhibiting an unusually strong electron-phonon interaction. The electron-lattice coupling is dominated by low-energy anharmonic phonons, which place the system on the verge of a charge density wave instability. We also provide definitive evidence for two-gap superconductivity in 2$H$-NbS$_2$, and show that the low- and high-energy peaks observed in tunneling spectra correspond to the $\Gamma$- and $K$-centered Fermi surface pockets, respectively. The present findings call for further efforts to determine whether our proposed mechanism underpins superconductivity in the whole family of metallic transition metal dichalcogenides.


title: Focusing RKKY interaction by graphene P-N junction

abstract: The carrier-mediated RKKY interaction between local spins plays an important role for the application of magnetically doped graphene in spintronics and quantum computation. Previous studies largely concentrate on the influence of electronic states of uniform systems on the RKKY interaction. Here we reveal a very different way to manipulate the RKKY interaction by showing that the anomalous focusing - a well-known electron optics phenomenon in graphene P-N junctions - can be utilized to refocus the massless Dirac electrons emanating from one local spin to the other local spin. This gives rise to rich spatial interference patterns and symmetry-protected non-oscillatory RKKY interaction with a strongly enhanced magnitude. It may provide a new way to engineer the long-range spin-spin interaction in graphene.


title: Spectral backtests of forecast distributions with application to risk management

abstract: We study a class of backtests for forecast distributions in which the test statistic depends on a spectral transformation that weights exceedance events by a function of the modeled probability level. The weighting scheme is specified by a kernel measure which makes explicit the user's priorities for model performance. The class of spectral backtests includes tests of unconditional coverage and tests of conditional coverage. We show how the class embeds a wide variety of backtests in the existing literature, and further propose novel variants which are easily implemented, well-sized and have good power. In an empirical application, we backtest forecast distributions for the overnight P&L of ten bank trading portfolios. For some portfolios, test results depend materially on the choice of kernel.


title: Quantitative stability of the free boundary in the obstacle problem

abstract: We prove some detailed quantitative stability results for the contact set and the solution of the classical obstacle problem in $\mathbb{R}^n$ ($n \ge 2$) under perturbations of the obstacle function, which is also equivalent to studying the variation of the equilibrium measure in classical potential theory under a perturbation of the external field. To do so, working in the setting of the whole space, we examine the evolution of the free boundary $\Gamma^t$ corresponding to the boundary of the contact set for a family of obstacle functions $h^t$. Assuming that $h=h^t (x) = h(t,x)$ is $C^{k+1,\alpha}$ in $[-1,1]\times \mathbb{R}^n$ and that the initial free boundary $\Gamma^0$ is regular, we prove that $\Gamma^t$ is twice differentiable in $t$ in a small neighborhood of $t=0$. Moreover, we show that the "normal velocity" and the "normal acceleration" of $\Gamma^t$ are respectively $C^{k-1,\alpha}$ and $C^{k-2,\alpha}$ scalar fields on $\Gamma^t$. This is accomplished by deriving equations for these velocity and acceleration and studying the regularity of their solutions via single and double layers estimates from potential theory.


title: Improving the Simulation of Quark and Gluon Jets with Herwig 7

abstract: The properties of quark and gluon jets, and the differences between them, are increasingly important at the LHC. However, Monte Carlo event generators are normally tuned to data from $e^+e^-$ collisions which are primarily sensitive to quark-initiated jets. In order to improve the description of gluon jets we make improvements to the perturbative and the non-perturbative modelling of gluon jets and include data with gluon-initiated jets in the tuning for the first time. The resultant tunes significantly improve the description of gluon jets and are now the default in Herwig 7.1.


title: Minimap2: pairwise alignment for nucleotide sequences

abstract: Motivation: Recent advances in sequencing technologies promise ultra-long reads of $\sim$100 kilo bases (kb) in average, full-length mRNA or cDNA reads in high throughput and genomic contigs over 100 mega bases (Mb) in length. Existing alignment programs are unable or inefficient to process such data at scale, which presses for the development of new alignment algorithms. Results: Minimap2 is a general-purpose alignment program to map DNA or long mRNA sequences against a large reference database. It works with accurate short reads of $\ge$100bp in length, $\ge$1kb genomic reads at error rate $\sim$15%, full-length noisy Direct RNA or cDNA reads, and assembly contigs or closely related full chromosomes of hundreds of megabases in length. Minimap2 does split-read alignment, employs concave gap cost for long insertions and deletions (INDELs) and introduces new heuristics to reduce spurious alignments. It is 3-4 times faster than mainstream short-read mappers at comparable accuracy and $\ge$30 times faster at higher accuracy for both genomic and mRNA reads, surpassing most aligners specialized in one type of alignment. Availability and implementation: https://github.com/lh3/minimap2 Contact: hengli@broadinstitute.org


title: Simulations of submonolayer Xe on Pt$(111)$: the case for a chaotic low temperature phase

abstract: Molecular Dynamics simulations are reported for the structural and thermodynamic properties of submonolayer xenon adsorbed on the $(111)$ surface of platinum for temperatures up to the (apparently incipient) triple point and beyond. While the motion of the atoms in the surface plane is treated with a standard two-dimensional molecular dynamics simulation, the model takes into consideration the thermal excitation of quantum states associated with surface-normal dynamics in an attempt to describe the apparent smoothing of the corrugation with increasing temperature. We examine the importance of this thermal smoothing to the relative stability of several observed and proposed low-temperature structures. Structure factor calculations are compared to experimental results in an attempt to determine the low temperature structure of this system. These calculations provide strong evidence that, at very low temperatures, the domain wall structure of a xenon monolayer adsorbed on a Pt$(111)$ substrate possesses a chaotic-like nature, exhibiting long-lived meta-stable states with pinned domain walls, these walls having narrow widths and irregular shapes. This result is contrary to the standard wisdom regarding this system, namely that the very low temperature phase of this system is a striped incommensurate phase. We present the case for further experimental investigation of this and similar systems as possible examples of chaotic low temperature phases in two dimensions.


title: Hierarchical Metric Learning for Optical Remote Sensing Scene Categorization

abstract: We address the problem of scene classification from optical remote sensing (RS) images based on the paradigm of hierarchical metric learning. Ideally, supervised metric learning strategies learn a projection from a set of training data points so as to minimize intra-class variance while maximizing inter-class separability to the class label space. However, standard metric learning techniques do not incorporate the class interaction information in learning the transformation matrix, which is often considered to be a bottleneck while dealing with fine-grained visual categories. As a remedy, we propose to organize the classes in a hierarchical fashion by exploring their visual similarities and subsequently learn separate distance metric transformations for the classes present at the non-leaf nodes of the tree. We employ an iterative max-margin clustering strategy to obtain the hierarchical organization of the classes. Experiment results obtained on the large-scale NWPU-RESISC45 and the popular UC-Merced datasets demonstrate the efficacy of the proposed hierarchical metric learning based RS scene recognition strategy in comparison to the standard approaches.


title: Verification of the Quantum Nonequilibrium Work Relation in the Presence of Decoherence

abstract: Although nonequilibrium work and fluctuation relations have been studied in detail within classical statistical physics, extending these results to open quantum systems has proven to be conceptually difficult. For systems that undergo decoherence but not dissipation, we argue that it is natural to define quantum work exactly as for isolated quantum systems, using the two-point measurement protocol. Complementing previous theoretical analysis using quantum channels, we show that the nonequilibrium work relation remains valid in this situation, and we test this assertion experimentally using a system engineered from an optically trapped ion. Our experimental results reveal the work relation's validity over a variety of driving speeds, decoherence rates, and effective temperatures and represent the first confirmation of the work relation for non-unitary dynamics.


title: Force and Mass Dynamics in Non-Newtonian Suspensions

abstract: Above a certain solid fraction, dense granular suspensions in water exhibit non-Newtonian behavior, including impact-activated solidification. Although it has been suggested that solidification depends on boundary interactions, quantitative experiments on the boundary forces have not been reported. Using high-speed video, tracer particles, and photoelastic boundaries, we determine the impactor kinematics and the magnitude and timings of impactor-driven events in the body and at the boundaries of cornstarch suspensions. We observe mass shocks in the suspension during impact. The shockfront dynamics are strongly correlated to those of the intruder. However, the total momentum associated with this shock never approaches the initial impactor momentum. We also observe a faster second front, associated with the propagation of pressure to the boundaries of the suspension. The two fronts depend differently on the initial impactor speed, $v_0$, and the suspension packing fraction. The speed of the pressure wave is at least an order of magnitude smaller than (linear) ultrasound speeds obtained for much higher frequencies, pointing to complex amplitude and frequency response of cornstarch suspensions to compressive strains.


title: Identification of spin wave resonances and crystal field transitions in simple chromites RCrO3 (R=Pr, Sm, Er) at ultralow temperatures in the THz spectral region

abstract: We report on THz absorption spectroscopy combined with high magnetic fields of polycrystalline RCrO3 (R= Pr, Sm, Er) aiming understanding spin wave resonances at their low temperature magnetic phases. Our measurements show that the temperature, and the implicit anisotropies at which the Cr3+ spin reorientation at TSR takes place, are determinant on the ferromagnetic-like (FM) and the antiferromagnetic-like (AFM) spin modes being optically active. It is found that they are dependent on Rare Earth 4f moment and ion size. We also studied temperature and field dependence of crystal field levels in the same spectroscopic region. Pr3+ non-Kramers emerges at 100 K and Zeeman splits. An observed absence of spin wave resonances in PrCrO3 is attributed to Pr3+ remaining paramagnetic. In SmCrO3 near cancelation of the spin and orbital moments is proposed as the possible reason for not detecting Sm3+ ground state transitions. Here, the FM and AFM resonant modes harden when the temperature decreases and split linearly on applied fields at 5 K and below. In ErCrO3 the Er3+ Kramers doublet becomes active at about the TSR onset. Each line further experiences Zeeman splitting under magnetic fields while an spin reversal induced by a 2.5 T field, back to the Gamma4 (Fz) from the Gamma1 phase at 2 K, produces a secondary splitting. The 5 K AFM and FM excitations in ErCrO3 have a concerted frequency-intensity temperature dependence and a shoulder pointing to the Er3+ smaller ion size also disrupting the two magnetic sublattice approximation . Both resonances reduce to one when the temperature is lowered to 2 K in the Gamma1 representation.


title: Diffractive quarkonium photoproduction in $pp$ and $pA$ collisions at the LHC: Predictions of the Resolved Pomeron model for the Run 2 energies

abstract: The inclusive diffractive quarkonium photoproduction in $pp$ and $pA$ collisions is investigated considering the Resolved Pomeron Model to describe the diffractive interaction. We estimate the rapidity and transverse momentum distributions for the $J/\Psi$, $\Psi(2S)$ and $\Upsilon$ photoproduction in hadronic collisions at the LHC and present our estimate for the total cross sections at the Run 2 energies. A comparison with the predictions associated to the exclusive production also is presented. Our results indicate that the inclusive diffractive production is a factor $\gtrsim 10$ smaller than the exclusive one in the kinematical range probed by the LHC.


title: Converting of algebraic Diophantine equations to a diagonal form with the help of an integer non-orthogonal transformation, maintaining the asymptotic behavior of the number of its integer solutions

abstract: The author showed that any homogeneous algebraic Diophantine equation of the second order can be converted to a diagonal form using an integer non-orthogonal transformation maintaining asymptotic behavior of the number of its integer solutions. In this paper, we consider the transformation to the diagonal form of a wider class of algebraic second-order Diophantine equations, and also we consider the conditions for converting higher order algebraic Diophantine equations to this form. The author found an asymptotic estimate for the number of integer solutions of the diagonal Thue equation of odd degree with an amount of variables greater than two, and also he got and asymptotic estimates of the number of integer solutions of other types of diagonal algebraic Diophantine equations.


title: A Multi-band Study of the remarkable Jet in Quasar 4C+19.44

abstract: We present arc-second-resolution data in the radio, IR, optical and X-ray for 4C+19.44 (=PKS 1354+195), the longest and straightest quasar jet with deep X-ray observations. We report results from radio images with half to one arc-second angular resolution at three frequencies, plus HST and Spitzer data. The Chandra data allow us to measure the X-ray spectral index in 10 distinct regions along the 18 arcsec jet and compare with the radio index. The radio and X-ray spectral indices of the jet regions are consistent with a value of $\alpha =0.80$ throughout the jet, to within 2 sigma uncertainties. The X-ray jet structure to the south extends beyond the prominent radio jet and connects to the southern radio lobe, and there is extended X-ray emission in the direction of the unseen counter jet and coincident with the northern radio lobe. This jet is remarkable since its straight appearance over a large distance allows the geometry factors to be taken as fixed along the jet. Using the model of inverse Compton scattering of the cosmic microwave background (iC/CMB) by relativistic electrons, we find that the magnetic field strengths and Doppler factors are relatively constant along the jet. If instead the X-rays are synchrotron emission, they must arise from a population of electrons distinct from the particles producing the radio synchrotron spectrum.


title: Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

abstract: Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 \eV. Other less probable conformations having different solvation pattern are also investigated.


title: Predicting potential treatments for complex diseases based on miRNA and tissue specificity

abstract: Drug repositioning, that is finding new uses for existing drugs to treat more patients. Cumulative studies demonstrate that the mature miRNAs as well as their precursors can be targeted by small molecular drugs. At the same time, human diseases result from the disordered interplay of tissue- and cell lineage-specific processes. However, few computational researches predict drug-disease potential relationships based on miRNA data and tissue specificity. Therefore, based on miRNA data and the tissue specificity of diseases, we propose a new method named as miTS to predict the potential treatments for diseases. Firstly, based on miRNAs data, target genes and information of FDA approved drugs, we evaluate the relationships between miRNAs and drugs in the tissue-specific PPI network. Then, we construct a tripartite network: drug-miRNA-disease Finally, we obtain the potential drug-disease associations based on the tripartite network. In this paper, we take breast cancer as case study and focus on the top-30 predicted drugs. 25 of them (83.3%) are found having known connections with breast cancer in CTD benchmark and the other 5 drugs are potential drugs for breast cancer. We further evaluate the 5 newly predicted drugs from clinical records, literature mining, KEGG pathways enrichment analysis and overlapping genes between enriched pathways. For each of the 5 new drugs, strongly supported evidences can be found in three or more aspects. In particular, Regorafenib has 15 overlapping KEGG pathways with breast cancer and their p-values are all very small. In addition, whether in the literature curation or clinical validation, Regorafenib has a strong correlation with breast cancer. All the facts show that Regorafenib is likely to be a truly effective drug, worthy of our further study. It further follows that our method miTS is effective and practical for predicting new drug indications.


title: Maximum genus of the Jenga like configurations

abstract: We treat the boundary of the union of blocks in the Jenga game as a surface with a polyhedral structure and consider its genus. We generalize the game and determine the maximum genus of the generalized game.


title: Comparison of performance of van der Waals-corrected exchange-correlation functionals for interlayer interaction in graphene and hexagonal boron nitride

abstract: Exchange-correlation functionals with corrections for van der Waals interactions (PBE-D2, PBE-D3, PBE-D3(BJ), PBE-TS, optPBE-vdW and vdW-DF2) are tested for graphene and hexagonal boron nitride, both in the form of bulk and bilayer. The characteristics of the potential energy surface, such as the barrier to relative sliding of the layers and magnitude of corrugation, and physically measurable properties associated with relative in-plane and out-of-plane motion of the layers including the shear modulus and modulus for axial compression, shear mode frequency and frequency of out-of-plane vibrations are considered. The PBE-D3(BJ) functional gives the best results for the stackings of hexagonal boron nitride and graphite that are known to be ground-state from the experimental studies. However, it fails to describe the order of metastable states of boron nitride in energy. The PBE-D3 and vdW-DF2 functionals, which reproduce this order correctly, are identified as the optimal choice for general studies. The vdW-DF2 functional is preferred for evaluation of the modulus for axial compression and frequency of out-of-plane vibrations, while the PBE-D3 functional is somewhat more accurate in calculations of the shear modulus and shear mode frequency. The best description of the latter properties, however, is achieved also using the vdW-DF2 functional combined with consideration of the experimental interlayer distance. In the specific case of graphene, the PBE-D2 functional works very well and can be further improved by adjustment of the parameters.


title: Lasso Guarantees for $ \beta $-Mixing Heavy Tailed Time Series

abstract: Many theoretical results for the lasso require the samples to be iid. Recent work has provided guarantees for the lasso assuming that the time series is generated by a sparse Vector Auto-Regressive (VAR) model with Gaussian innovations. Proofs of these results rely critically on the fact that the true data generating mechanism (DGM) is a finite-order Gaussian VAR. This assumption is quite brittle: linear transformations, including selecting a subset of variables, can lead to the violation of this assumption. In order to break free from such assumptions, we derive non-asymptotic inequalities for estimation error and prediction error of the lasso estimate of the best linear predictor without assuming any special parametric form of the DGM. Instead, we rely only on (strict) stationarity and geometrically decaying \b{eta}-mixing coefficients to establish error bounds for the lasso for subweibull random vectors. The class of subweibull random variables that we introduce includes subgaussian and subexponential random variables but also includes random variables with tails heavier than an exponential. We also show that, for Gaussian processes, the \b{eta}-mixing condition can be relaxed to summability of the {\alpha}-mixing coefficients. Our work provides an alternative proof of the consistency of the lasso for sparse Gaussian VAR models. But the applicability of our results extends to non-Gaussian and non-linear times series models as the examples we provide demonstrate.


title: Edge stacking dislocations in two-dimensional bilayers with a small lattice mismatch

abstract: Incomplete stacking dislocations are predicted to form at edges of the shorter upper layer in two-dimensional hexagonal bilayers upon stretching the longer bottom layer. A concept of the edge Burgers vector is introduced to describe such dislocations by analogy with the Burgers vector of standard bulk dislocations. Analytical solutions for the structure and energy of edge stacking dislocations in bilayer graphene are obtained depending on the magnitude of elongation and angles between the edge Burgers vector, direction of elongation and edge. The barrier for penetration of stacking dislocations inside the bilayer is estimated. The possibilities to measure the barrier to relative motion of graphene layers and strain of graphene on a substrate by observation of edge stacking dislocations are discussed.


title: Automorphisms and deformations of conformally K\"ahler, Einstein-Maxwell metrics

abstract: We obtain a structure theorem for the group of holomorphic automorphisms of a conformally K\"ahler, Einstein-Maxwell metric, extending the classical results of Matsushima, Licherowicz and Calabi in the K\"ahler-Einstein, cscK, and extremal K\"ahler cases. Combined with previous results of LeBrun, Apostolov-Maschler and Futaki-Ono, this completes the classification of the conformally K\"ahler, Einstein--Maxwell metrics on $\mathbb{{CP}}^1 \times \mathbb{{CP}}^1$. We also use our result in order to introduce a (relative) Mabuchi energy in the more general context of $(K, q, a)$-extremal K\"ahler metrics in a given K\"ahler class, and show that the existence of $(K, q, a)$-extremal K\"ahler metrics is stable under small deformation of the K\"ahler class, the Killing vector field $K$ and the normalization constant $a$.


title: PoMSA: An Efficient and Precise Position-based Multiple Sequence Alignment Technique

abstract: Analyzing the relation between a set of biological sequences can help to identify and understand the evolutionary history of these sequences and the functional relations among them. Multiple Sequence Alignment (MSA) is the main obstacle to proper design and develop homology and evolutionary modeling applications since these kinds of applications require an effective MSA technique with high accuracy. This work proposes a novel Position-based Multiple Sequence Alignment (PoMSA) technique -- which depends on generating a position matrix for a given set of biological sequences. This position matrix can be used to reconstruct the given set of sequences in more aligned format. On the contrary of existing techniques, PoMSA uses position matrix instead of distance matrix to correctly adding gaps in sequences which improve the efficiency of the alignment operation. We have evaluated the proposed technique with different datasets benchmarks such as BAliBASE, OXBench, and SMART. The experiments show that PoMSA technique satisfies higher alignment score compared to existing state-of-art algorithms: Clustal-Omega, MAFTT, and MUSCLE.


title: Exact Solution of the Zakharov-Shabat Scattering Problem for Doubly-Truncated Multi-Soliton Potentials

abstract: Recent studies have revealed that multi-soliton solutions of the nonlinear Schr\"odinger equation, as carriers of information, offer a promising solution to the problem of nonlinear signal distortions in fiber optic channels. In any nonlinear Fourier transform based transmission methodology seeking to modulate the discrete spectrum of the multi-solitons, choice of an appropriate windowing function is an important design issue on account of the unbounded support of such signals. Here, we consider the rectangle function as the windowing function for the multi-solitonic signal and provide the exact solution of the associated Zakharov-Shabat scattering problem for the windowed/doubly-truncated multi-soliton potential. This method further allows us to avoid prohibitive numerical computations normally required in order to accurately quantify the effect of time-domain windowing on the nonlinear Fourier spectrum of the multi-solitonic signals. The method devised in this work also applies to general type of signals and may prove to be a useful tool in the theoretical analysis of such systems.


title: Ball characterizations in spaces of constant curvature

abstract: High proved the following theorem. If the intersections of any two congruent copies of a plane convex body are centrally symmetric, then this body is a circle. In our paper we extend the theorem of High to spherical, Euclidean and hyperbolic spaces, under some regularity assumptions. Suppose that in any of these spaces there is a pair of closed convex sets of class $C^2_+$ with interior points, different from the whole space, and the intersections of any congruent copies of these sets are centrally symmetric (provided they have non-empty interiors). Then our sets are congruent balls. Under the same hypotheses, but if we require only central symmetry of small intersections, then our sets are either congruent balls, or paraballs, or have as connected components of their boundaries congruent hyperspheres (and the converse implication also holds). Under the same hypotheses, if we require central symmetry of all compact intersections, then either our sets are congruent balls or paraballs, or have as connected components of their boundaries congruent hyperspheres, and either $d \ge 3$, or $d=2$ and one of the sets is bounded by one hypercycle, or both sets are congruent parallel domains of straight lines, or there are no more compact intersections than those bounded by two finite hypercycle arcs (and the converse implication also holds). We also prove a dual theorem. If in any of these spaces there is a pair of smooth closed convex sets, such that both of them have supporting spheres at any of their boundary points --- for $S^d$ of radius less than $ \pi /2$ --- and the closed convex hulls of any congruent copies of these sets are centrally symmetric, then our sets are congruent balls.


title: Ghost characters and character varieties of 2-fold branched covers

abstract: It is known that for any knot $K$ every (meridionally) trace-free $SL_2(C)$-representation of the knot group $G(K)$ gives an $SL_2(C)$-representation of the fundamental group $\pi_1(\Sigma_2K)$ of the 2-fold branched covering $\Sigma_2K$ of the 3-sphere branched along $K$. In this paper, we show by using a notion called a ghost character of a knot that for the (4,5)-torus knot $T_{4,5}$ the fundamental group $\pi_1(\Sigma_2T_{4,5})$ has an $SL(2,C)$-representation which cannot be given by any trace-free $SL_2(C)$-representation of $G(T_{4,5})$.


title: On Moment Condition and Center Condition for Abel Equation

abstract: In this paper we consider Abel equation $x' = g(t)x^2+f(t)x^3$, where $f$ and $g$ are analytical functions. We proved that if the equation has a center at $x=0$, then the Moment Conditions, i. e., $m_k=\int_{-1}^1f(t)(G(t))^kdt=0,~~k=0,1,2$, is satisfied where $G(t)=\int_{-1}^tg(s)ds$. Besides, we give partial a positive answer to a conjecture proposed by Y. Lijun and T. Yun in 2001.


title: Spatial Mixing and Non-local Markov chains

abstract: We consider spin systems with nearest-neighbor interactions on an $n$-vertex $d$-dimensional cube of the integer lattice graph $\mathbb{Z}^d$. We study the effects that exponential decay with distance of spin correlations, specifically the strong spatial mixing condition (SSM), has on the rate of convergence to equilibrium distribution of non-local Markov chains. We prove that SSM implies $O(\log n)$ mixing of a block dynamics whose steps can be implemented efficiently. We then develop a methodology, consisting of several new comparison inequalities concerning various block dynamics, that allow us to extend this result to other non-local dynamics. As a first application of our method we prove that, if SSM holds, then the relaxation time (i.e., the inverse spectral gap) of general block dynamics is $O(r)$, where $r$ is the number of blocks. A second application of our technology concerns the Swendsen-Wang dynamics for the ferromagnetic Ising and Potts models. We show that SSM implies an $O(1)$ bound for the relaxation time. As a by-product of this implication we observe that the relaxation time of the Swendsen-Wang dynamics in square boxes of $\mathbb{Z}^2$ is $O(1)$ throughout the subcritical regime of the $q$-state Potts model, for all $q \ge 2$. We also prove that for monotone spin systems SSM implies that the mixing time of systematic scan dynamics is $O(\log n (\log \log n)^2)$. Systematic scan dynamics are widely employed in practice but have proved hard to analyze. Our proofs use a variety of techniques for the analysis of Markov chains including coupling, functional analysis and linear algebra.


title: How anisotropy beats fractality in two-dimensional on-lattice DLA growth

abstract: We study the fractal structure of Diffusion-Limited Aggregation (DLA) clusters on the square lattice by extensive numerical simulations (with clusters having up to $10^8$ particles). We observe that DLA clusters undergo strongly anisotropic growth, with the maximal growth rate along the axes. The naive scaling limit of a DLA cluster by its diameter is thus deterministic and one-dimensional. At the same time, on all scales from the particle size to the size of the entire cluster it has non-trivial box-counting fractal dimension which corresponds to the overall growth rate which, in turn, is smaller than the growth rate along the axes. This suggests that the fractal nature of the lattice DLA should be understood in terms of fluctuations around one-dimensional backbone of the cluster.


title: Cramer's rules for the solution to the two-sided restricted quaternion matrix equation

abstract: Weighted singular value decomposition (WSVD) of a quaternion matrix and with its help determinantal representations of the quaternion weighted Moore-Penrose inverse have been derived recently by the author. In this paper, using these determinantal representations, explicit determinantal representation formulas for the solution of the restricted quaternion matrix equations, ${\bf A}{\bf X}{\bf B}={\bf D}$, and consequently, ${\bf A}{\bf X}={\bf D}$ and ${\bf X}{\bf B}={\bf D}$ are obtained within the framework of the theory of column-row determinants. We consider all possible cases depending on weighted matrices.


title: Lepton number violation in $B_s$ meson decays induced by an on-shell Majorana neutrino

abstract: Lepton-number violation can be induced by the exchange of an on-shell Majorana neutrino $N$ in semileptonic $|\Delta L|=2$ decays of the $B_s$ meson, $B_s^0 \to P^-\pi^-\mu^+ \mu^+$ with $P=K, D_s$. We investigate the production of such a heavy sterile neutrino through these four-body $\mu^+\mu^+$ channels and explore the sensitivity that can be reached at the LHCb and CMS experiments. For heavy neutrino lifetimes of $\tau_N$ = [1, 100, 1000] ps and integrated luminosities collected of 10 and 50 fb${}^{-1}$ at the LHCb and 30, 300, and 3000 fb${}^{-1}$ at the CMS, we find a significant sensitivity on branching fractions of the orders ${\rm BR}(B_s^0 \to K^- \pi^-\mu^+ \mu^+) \lesssim \mathcal{O}(10^{-9} - 10^{-8})$ and ${\rm BR}(B_s^0 \to D_s^- \pi^-\mu^+ \mu^+) \lesssim \mathcal{O}(10^{-8} - 10^{-7})$. In the kinematically allowed mass ranges of $m_N \in [0.25,4.77]$ GeV and $m_N \in [0.25,3.29]$ GeV, respectively, we exclude regions on the parameter space ($m_N,|V_{\mu N}|^{2}$) associated with the heavy neutrino, which could slightly improve the limits from $B^- \to \pi^+\mu^-\mu^-$ (LHCb).


title: Time quasi-periodic gravity water waves in finite depth

abstract: We prove the existence and the linear stability of Cantor families of small amplitude time quasi-periodic standing water wave solutions - namely periodic and even in the space variable x - of a bi-dimensional ocean with finite depth under the action of pure gravity. Such a result holds for all the values of the depth parameter in a Borel set of asymptotically full measure. This is a small divisor problem. The main difficulties are the quasi-linear nature of the gravity water waves equations and the fact that the linear frequencies grow just in a sublinear way at infinity. We overcome these problems by first reducing the linearized operators obtained at each approximate quasi-periodic solution along the Nash-Moser iteration to constant coefficients up to smoothing operators, using pseudo-differential changes of variables that are quasi-periodic in time. Then we apply a KAM reducibility scheme which requires very weak Melnikov non-resonance conditions (losing derivatives both in time and space), which we are able to verify for most values of the depth parameter using degenerate KAM theory arguments.


title: Elliptic diffeomorphisms of symplectic 4-manifolds

abstract: We show that symplectically embedded $(-1)$-tori give rise to certain elements in the symplectic mapping class group of $4$-manifolds. An example is given where such elements are proved to be of infinite order.


title: A Latent Variable Model for Two-Dimensional Canonical Correlation Analysis and its Variational Inference

abstract: Describing the dimension reduction (DR) techniques by means of probabilistic models has recently been given special attention. Probabilistic models, in addition to a better interpretability of the DR methods, provide a framework for further extensions of such algorithms. One of the new approaches to the probabilistic DR methods is to preserving the internal structure of data. It is meant that it is not necessary that the data first be converted from the matrix or tensor format to the vector format in the process of dimensionality reduction. In this paper, a latent variable model for matrix-variate data for canonical correlation analysis (CCA) is proposed. Since in general there is not any analytical maximum likelihood solution for this model, we present two approaches for learning the parameters. The proposed methods are evaluated using the synthetic data in terms of convergence and quality of mappings. Also, real data set is employed for assessing the proposed methods with several probabilistic and none-probabilistic CCA based approaches. The results confirm the superiority of the proposed methods with respect to the competing algorithms. Moreover, this model can be considered as a framework for further extensions.