title: Tensor calculus in spherical coordinates using Jacobi polynomials. Part-I: Mathematical analysis and derivations

abstract: This paper presents a method for the accurate and efficient computations on scalar, vector and tensor fields in three-dimensional spherical polar coordinates. The methods uses spin-weighted spherical harmonics in the angular directions and rescaled Jacobi polynomials in the radial direction. For the 2-sphere, spin-weighted harmonics allow for automating calculations in a fashion as similar to Fourier series as possible. Derivative operators act as wavenumber multiplication on a set of spectral coefficients. After transforming the angular directions, a set of orthogonal tensor rotations put the radially dependent spectral coefficients into individual spaces each obeying a particular regularity condition at the origin. These regularity spaces have remarkably simple properties under standard vector-calculus operations, such as \textit{grad} and \textit{div}. We use a hierarchy of rescaled Jacobi polynomials for a basis on these regularity spaces. It is possible to select the Jacobi-polynomial parameters such that all relevant operators act in a minimally banded way. Altogether, the geometric structure allows for the accurate and efficient solution of general partial differential equations in the unit ball.

categories: math.NA astro-ph.IM math.CA physics.flu-dyn


title: Solar Kinetic Energy and Cross Helicity Spectra

abstract: We develop a formalism that treats the calculation of solar kinetic energy and cross helicity spectra in an equal manner to that of magnetic energy and helicity spectra. The magnetic helicity spectrum is shown to be equal to the vertical part of the current helicity spectrum divided by the square of the wavenumber. For the cross helicity, we apply the recently developed two-scale approach globally over an entire active region to account for the sign change between the two polarities. Using vector magnetograms and Dopplergrams of NOAA 11158 and 12266, we show that kinetic and magnetic energy spectra have similar slopes at intermediate wavenumbers, where the contribution from the granulation velocity has been removed. At wavenumbers around 0.3 Mm$^{-1}$, the magnetic helicity is found to be close to its maximal value. The cross helicity spectra are found to be within about 10% of the maximum possible value. Using the two-scale method for NOAA 12266, the global cross helicity spectrum is found to be particularly steep, similarly to what has previously been found in theoretical models of spot generation. In the quiet Sun, by comparison, the cross helicity spectrum is found to be small.


title: Classification of auditory stimuli from EEG signals with a regulated recurrent neural network reservoir

abstract: The use of electroencephalogram (EEG) as the main input signal in brain-machine interfaces has been widely proposed due to the non-invasive nature of the EEG. Here we are specifically interested in interfaces that extract information from the auditory system and more specifically in the task of classifying heard speech from EEGs. To do so, we propose to limit the preprocessing of the EEGs and use machine learning approaches to automatically extract their meaningful characteristics. More specifically, we use a regulated recurrent neural network (RNN) reservoir, which has been shown to outperform classic machine learning approaches when applied to several different bio-signals, and we compare it with a deep neural network approach. Moreover, we also investigate the classification performance as a function of the number of EEG electrodes. A set of 8 subjects were presented randomly with 3 different auditory stimuli (English vowels a, i and u). We obtained an excellent classification rate of 83.2% with the RNN when considering all 64 electrodes. A rate of 81.7% was achieved with only 10 electrodes.


title: Adversarial Training of Variational Auto-encoders for High Fidelity Image Generation

abstract: Variational auto-encoders (VAEs) provide an attractive solution to image generation problem. However, they tend to produce blurred and over-smoothed images due to their dependence on pixel-wise reconstruction loss. This paper introduces a new approach to alleviate this problem in the VAE based generative models. Our model simultaneously learns to match the data, reconstruction loss and the latent distributions of real and fake images to improve the quality of generated samples. To compute the loss distributions, we introduce an auto-encoder based discriminator model which allows an adversarial learning procedure. The discriminator in our model also provides perceptual guidance to the VAE by matching the learned similarity metric of the real and fake samples in the latent space. To stabilize the overall training process, our model uses an error feedback approach to maintain the equilibrium between competing networks in the model. Our experiments show that the generated samples from our proposed model exhibit a diverse set of attributes and facial expressions and scale up to high-resolution images very well.


title: Fall to the Centre in Atom Traps and Point-Particle EFT for Absorptive Systems

abstract: Polarizable atoms interacting with a charged wire do so through an inverse-square potential, $V = - g/r^2$. This system is known to realize scale invariance in a nontrivial way and to be subject to ambiguities associated with the choice of boundary condition at the origin, often termed the problem of `fall to the center'. Point-particle effective field theory (PPEFT) provides a systematic framework for determining the boundary condition in terms of the properties of the source residing at the origin. We apply this formalism to the charged-wire/polarizable-atom problem, finding a result that is not a self-adjoint extension because of absorption of atoms by the wire. We explore the RG flow of the complex coupling constant for the dominant low-energy effective interactions, finding flows whose character is qualitatively different when $g$ is above or below a critical value, $g_c$. Unlike the self-adjoint case, (complex) fixed points exist when $g> g_c$, which we show correspond to perfect absorber (or perfect emitter) boundary conditions. We describe experimental consequences for wire-atom interactions and the possibility of observing the anomalous breaking of scale invariance.

categories: hep-ph cond-mat.quant-gas physics.atom-ph quant-ph


title: Simulating dysarthric speech for training data augmentation in clinical speech applications

abstract: Training machine learning algorithms for speech applications requires large, labeled training data sets. This is problematic for clinical applications where obtaining such data is prohibitively expensive because of privacy concerns or lack of access. As a result, clinical speech applications are typically developed using small data sets with only tens of speakers. In this paper, we propose a method for simulating training data for clinical applications by transforming healthy speech to dysarthric speech using adversarial training. We evaluate the efficacy of our approach using both objective and subjective criteria. We present the transformed samples to five experienced speech-language pathologists (SLPs) and ask them to identify the samples as healthy or dysarthric. The results reveal that the SLPs identify the transformed speech as dysarthric 65% of the time. In a pilot classification experiment, we show that by using the simulated speech samples to balance an existing dataset, the classification accuracy improves by about 10% after data augmentation.


title: Accretivity of the general second order linear differential operator

abstract: For the general second order linear differential operator $$\mathcal L_0 = \sum_{j, \, k=1}^n \, a_{jk} \, \partial_j \partial_k + \sum_{j=1}^n \, b_{j} \, \partial_j + c$$ with complex-valued distributional coefficients $a_{jk}$, $b_{j}$, and $c$ in an open set $\Omega \subseteq \mathbf{R}^n$ ($n \ge 1$), we present conditions which ensure that $-\mathcal L_0$ is accretive, i.e., ${\rm Re} \, \langle -\mathcal L_0 \phi, \phi\rangle \ge 0$ for all $\phi \in C^\infty_0(\Omega).$


title: Unusual layered order and charge disproportionation in double perovskite Ca2FeMnO6

abstract: While double perovskites A2BB'O6, if ordered, usually form a rock-salt-type structure with a checkerboard B/B' ordering, it is surprising that Ca2FeMnO6 has alternate FeO2 and MnO2 layers in its perovskite structure. Here we demonstrate, using density functional calculations, that this unusual layered ordering facilitates, and is largely helped by, the Fe3+-Fe5+ charge disproportionation (CD) of the formal Fe4+ ions, which would otherwise be frustrated in the common rock salt structure. To further verify the important role of the CD for stabilization of this layered ordering, we carry out a comparative study for the isostructural Ca2TiMnO6 which has a simple Ti4+-Mn4+ state free of the CD. Our calculations indicate that Ca2TiMnO6 instead prefers the standard rock salt structure to the layered one. Thus our study shows a nontrivial interplay between the CD and the type of ion ordering, and proves that the CD is strongly involved in stabilizing the unusual layered order of Ca2FeMnO6.


title: Scalable Bilinear $\pi$ Learning Using State and Action Features

abstract: Approximate linear programming (ALP) represents one of the major algorithmic families to solve large-scale Markov decision processes (MDP). In this work, we study a primal-dual formulation of the ALP, and develop a scalable, model-free algorithm called bilinear $\pi$ learning for reinforcement learning when a sampling oracle is provided. This algorithm enjoys a number of advantages. First, it adopts (bi)linear models to represent the high-dimensional value function and state-action distributions, using given state and action features. Its run-time complexity depends on the number of features, not the size of the underlying MDPs. Second, it operates in a fully online fashion without having to store any sample, thus having minimal memory footprint. Third, we prove that it is sample-efficient, solving for the optimal policy to high precision with a sample complexity linear in the dimension of the parameter space.


title: A high-order HDG method for the Biot's consolidation model

abstract: We propose a novel high-order HDG method for the Biot's consolidation model in poroelasticity. We present optimal error analysis for both the semi-discrete and full-discrete (combined with temporal backward differentiation formula) schemes. Numerical tests are provided to demonstrate the performance of the method.


title: A short proof of Bresinski's Theorem on Gorenstein semigroup rings generated by 4 elements

abstract: Let $H=\langle n_1, \ldots ,n_4\rangle$ be a numerical semigroup generated by $4$ elements, which is symmetric and let $k[H]$ be the semigroup ring of $H$ over a field $k$. H. Bresinski proved that the defining ideal of $k[H]$ is minimally generated by $3$ or $5$ elements. We give a new short proof of Bresinski's Theorem using the structure theorem of Buchsbaum and Eisenbud on the minimal free resolution of Gorenstein rings of embedding codimension $3$.


title: Rateless Codes for Near-Perfect Load Balancing in Distributed Matrix-Vector Multiplication

abstract: Large-scale machine learning and data mining applications require computer systems to perform massive matrix-vector and matrix-matrix multiplication operations that need to be parallelized across multiple nodes. The presence of straggling nodes -- computing nodes that unpredictably slowdown or fail -- is a major bottleneck in such distributed computations. Ideal load balancing strategies that dynamically allocate more tasks to faster nodes require knowledge or monitoring of node speeds as well as the ability to quickly move data. Recently proposed fixed-rate erasure coding strategies can handle unpredictable node slowdown, but they ignore partial work done by straggling nodes thus resulting in a lot of redundant computation. We propose a \emph{rateless fountain coding} strategy that achieves the best of both worlds -- we prove that its latency is asymptotically equal to ideal load balancing, and it performs asymptotically zero redundant computations. Our idea is to create linear combinations of the $m$ rows of the matrix and assign these encoded rows to different worker nodes. The original matrix-vector product can be decoded as soon as slightly more than $m$ row-vector products are collectively finished by the nodes. We conduct experiments in three computing environments: local parallel computing, Amazon EC2, and Amazon Lambda, which show that rateless coding gives as much as $3\times$ speed-up over uncoded schemes.


title: Sim-to-Real: Learning Agile Locomotion For Quadruped Robots

abstract: Designing agile locomotion for quadruped robots often requires extensive expertise and tedious manual tuning. In this paper, we present a system to automate this process by leveraging deep reinforcement learning techniques. Our system can learn quadruped locomotion from scratch using simple reward signals. In addition, users can provide an open loop reference to guide the learning process when more control over the learned gait is needed. The control policies are learned in a physics simulator and then deployed on real robots. In robotics, policies trained in simulation often do not transfer to the real world. We narrow this reality gap by improving the physics simulator and learning robust policies. We improve the simulation using system identification, developing an accurate actuator model and simulating latency. We learn robust controllers by randomizing the physical environments, adding perturbations and designing a compact observation space. We evaluate our system on two agile locomotion gaits: trotting and galloping. After learning in simulation, a quadruped robot can successfully perform both gaits in the real world.


title: Emergent antiferromagnetism of YTiO3 in YTiO3-CaTiO3 superlattices

abstract: Transport and magnetoresistance measurements are performed on metallic, high-carrier density YTiO3-CaTiO3 superlattices as a probe towards the investigation of an emergent magnetic order of YTiO3. On varying the thickness of YTiO3 while keeping the CaTiO3 layer thickness constant in the superlattices, a low-temperature upturn in sheet-resistance, a non-Fermi liquid-like charge transport and positive magnetoresistance are observed. Analyses of the origin of such effects suggest that a unique antiferromagnetic order is realized in the ultra-thin, epitaxially strained YTiO3 layers, which corroborates well with some recent theoretical predictions in this regard.


title: Deep Learning Coordinated Beamforming for Highly-Mobile Millimeter Wave Systems

abstract: Supporting high mobility in millimeter wave (mmWave) systems enables a wide range of important applications such as vehicular communications and wireless virtual/augmented reality. Realizing this in practice, though, requires overcoming several challenges. First, the use of narrow beams and the sensitivity of mmWave signals to blockage greatly impact the coverage and reliability of highly-mobile links. Second, highly-mobile users in dense mmWave deployments need to frequently hand-off between base stations (BSs), which is associated with critical control and latency overhead. Further, identifying the optimal beamforming vectors in large antenna array mmWave systems requires considerable training overhead, which significantly affects the efficiency of these mobile systems. In this paper, a novel integrated machine learning and coordinated beamforming solution is developed to overcome these challenges and enable highly-mobile mmWave applications. In the proposed solution, a number of distributed yet coordinating BSs simultaneously serve a mobile user. This user ideally needs to transmit only one uplink training pilot sequence that will be jointly received at the coordinating BSs using omni or quasi-omni beam patterns. These received signals draw a defining signature not only for the user location, but also for its interaction with the surrounding environment. The developed solution then leverages a deep learning model that learns how to use these signatures to predict the beamforming vectors at the BSs. This renders a comprehensive solution that supports highly-mobile mmWave applications with reliable coverage, low latency, and negligible training overhead. Simulation results show that the proposed deep-learning coordinated beamforming strategy approaches the achievable rate of the genie-aided solution that knows the optimal beamforming vectors with no training overhead.


title: Communication, Computing and Caching for Mobile VR Delivery: Modeling and Trade-off

abstract: Mobile virtual reality (VR) delivery is gaining increasing attention from both industry and academia due to its ability to provide an immersive experience. However, achieving mobile VR delivery requires ultra-high transmission rate, deemed as a first killer application for 5G wireless networks. In this paper, in order to alleviate the traffic burden over wireless networks, we develop an implementation framework for mobile VR delivery by utilizing caching and computing capabilities of mobile VR device. We then jointly optimize the caching and computation offloading policy for minimizing the required average transmission rate under the latency and local average energy consumption constraints. In a symmetric scenario, we obtain the optimal joint policy and the closed-form expression of the minimum average transmission rate. Accordingly, we analyze the tradeoff among communication, computing and caching, and then reveal analytically the fact that the communication overhead can be traded by the computing and caching capabilities of mobile VR device, and also what conditions must be met for it to happen. Finally, we discuss the optimization problem in a heterogeneous scenario, and propose an efficient suboptimal algorithm with low computation complexity, which is shown to achieve good performance in the numerical results.


title: There are no Cubic Graphs on 26 Vertices with Crossing Number 10 or 11

abstract: We show that no cubic graphs of order 26 have crossing number larger than 9, which proves a conjecture of Ed Pegg Jr and Geoffrey Exoo that the smallest cubic graphs with crossing number 11 have 28 vertices. This result is achieved by first eliminating all girth 3 graphs from consideration, and then using the recently developed QuickCross heuristic to find good embeddings of each remaining graph. In the cases where the embedding found has 10 or more crossings, the heuristic is re-run with a different settings of parameters until an embedding with fewer than 10 crossings is found. We provide a minimal example of a cubic graph on 28 vertices with crossing number 10, and also exhibit for the first time a cubic graph on 30 vertices with crossing number 12, which we conjecture is minimal.


title: Latent Fingerprint Recognition: Role of Texture Template

abstract: We propose a texture template approach, consisting of a set of virtual minutiae, to improve the overall latent fingerprint recognition accuracy. To compensate for the lack of sufficient number of minutiae in poor quality latent prints, we generate a set of virtual minutiae. However, due to a large number of these regularly placed virtual minutiae, texture based template matching has a large computational requirement compared to matching true minutiae templates. To improve both the accuracy and efficiency of the texture template matching, we investigate: i) both original and enhanced fingerprint patches for training convolutional neural networks (ConvNets) to improve the distinctiveness of descriptors associated with each virtual minutiae, ii) smaller patches around virtual minutiae and a fast ConvNet architecture to speed up descriptor extraction, iii) reduce the descriptor length, iv) a modified hierarchical graph matching strategy to improve the matching speed, and v) extraction of multiple texture templates to boost the performance. Experiments on NIST SD27 latent database show that the above strategies can improve the matching speed from 11 ms (24 threads) per comparison (between a latent and a reference print) to only 7.7 ms (single thread) per comparison while improving the rank-1 accuracy by 8.9% against 10K gallery.


title: Conditional quantum nonlocality in dimeric and trimeric arrays of organic molecules

abstract: Arrays of covalently bound organic molecules possess potential for light-harvesting and energy transfer applications due to the strong coherent dipole-dipole coupling between the transition dipole moments of the molecules involved. Here, we show that such molecular systems, based on perylene-molecules, can be considered as arrays of qubits that are amenable for laser-driven quantum coherent control. The perylene monomers exhibit dephasing times longer than four orders of magnitude a typical gating time, thus allowing for the execution of a large number of gate operations on the sub-picosecond timescale. Specifically, we demonstrate quantum logic gates and entanglement in bipartite (dimer) and tripartite (trimer) systems of perylene-based arrays. In dimers, naturally entangled states with a tailored degree of entanglement can be produced. The nonlocality of the molecular trimer entanglement is demonstrated by testing Mermin's (Bell-like) inequality violation.


title: Scattering of a dark-bright soliton by an impurity

abstract: We study the dynamics of a dark-bright soliton interacting with a fixed impurity using a mean-field approach. The system is described by a vector nonlinear Schrodinger equation (NLSE) appropriate to multicomponent Bose-Einstein condensates. We use the variational approximation, based on hyperbolic functions, where we have the center of mass of the two components to describe the propagation of the dark and bright components independently. Therefore, it allows the dark-bright soliton to oscillate. The fixed local impurity is modeled by a delta function. Also, we use perturbation methods to derive the equations of motion for the center of mass of the two components. The interaction of the dark-bright soliton with a delta function potential excites different modes in the system. The analytical model capture two of these modes: the relative oscillation between the two components and the oscillation in the widths. The numerical simulations show additional internal modes play an important role in the interaction problem. The excitation of internal modes corresponds to inelastic scattering. In addition, we calculate the maximum velocity for a dark-bright soliton and find it is limited to a value below the sound speed, depending on the relative number of atoms present in the bright soliton component and excavated by the dark soliton component, respectively. Above a critical value of the maximum velocity, the two components are no longer described by one center of mass variable and develop internal oscillations, eventually breaking apart when pushed to higher velocities. This effect limits the incident kinetic energy in scattering studies and presents a smoking gun experimental signal.


title: The Circumstellar Disk and Asymmetric outflow of the EX Lup Outburst System

abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations at 0.3 arcsec-resolution of EX Lup, the prototype of the EXor class of outbursting pre-main sequence stars. The circumstellar disk of EX Lup is resolved for the first time in 1.3mm continuum emission and in the $J$=2--1 spectral line of three isotopologues of CO. At the spatial resolution and sensitivity achieved, the compact dust continuum disk shows no indications of clumps, fragments, or asymmetries above 5-sigma level. Radiative transfer modeling constrains the characteristic radius of the dust disk to 23 au and a total dust mass of 1.0$\times$10$^{-4}$ M$_\odot$ (33 M_earth), similar to other EXor sources. The $^{13}$CO and C$^{18}$O line emission trace the disk rotation and are used to constrain the disk geometry, kinematics, and a total gas disk mass of 5.1$\times$10$^{-4}$ M$_\odot$. The $^{12}$CO emission extends out to a radius of 200 au and is asymmetric, with one side deviating from Keplerian rotation. We detect blue-shifted, $^{12}$CO arc-like emission located 0.8 arcsec to the north-west, and spatially disconnected from the disk emission. We interpret this extended structure as the brightened walls of a cavity excavated by an outflow, which are more commonly seen in FUor sources. Such outflows have also been seen in the borderline FU/EXor object V1647 Ori, but not towards EXor objects. Our detection provides evidence that the outflow phenomenon persists into the EXor phase, suggesting that FUor and EXor objects are a continuous population in which outflow activity declines with age, with transitional objects such as EX Lup and V1647 Ori.


title: On Some Applications of Sakai's Geometric Theory of Discrete Painlev\'e Equations

abstract: Although the theory of discrete Painlev\'e (dP) equations is rather young, more and more examples of such equations appear in interesting and important applications. Thus, it is essential to be able to recognize these equations, to be able to identify their type, and to see where they belong in the classification scheme. The definite classification scheme for dP equations was proposed by H. Sakai, who used geometric ideas to identify 22 different classes of these equations. However, in a major contrast with the theory of ordinary differential Painlev\'e equations, there are infinitely many non-equivalent discrete equations in each class. Thus, there is no general form for a dP equation in each class, although some nice canonical examples in each equation class are known. The main objective of this paper is to illustrate that, in addition to providing the classification scheme, the geometric ideas of Sakai give us a powerful tool to study dP equations. We consider a very complicated example of a dP equation that describes a simple Schlesinger transformation of a Fuchsian system and we show how this equation can be identified with a much simpler canonical example of the dP equation of the same type and moreover, we give an explicit change of coordinates transforming one equation into the other. Among our main tools are the birational representation of the affine Weyl symmetry group of the equation and the period map. Even though we focus on a concrete example, the techniques that we use are general and can be easily adapted to other examples.


title: Effective field theory of the Higgs mode in a two-dimensional dilute Bose gas

abstract: We investigate the spectral function of the Higgs mode in a two dimensional Bose gas, by using the effective field theory in the zero temperature limit. Our approach explains the experimental feature that the peak of the spectral function is a soft continuum rather than a sharp peak, broadened and vanishing in the superfluid phase, which cannot be explained in terms of the $O(2)$ model. We also find that the scalar susceptibility is the same as the longitudinal susceptibility.


title: Stacked U-Nets: A No-Frills Approach to Natural Image Segmentation

abstract: Many imaging tasks require global information about all pixels in an image. Conventional bottom-up classification networks globalize information by decreasing resolution; features are pooled and downsampled into a single output. But for semantic segmentation and object detection tasks, a network must provide higher-resolution pixel-level outputs. To globalize information while preserving resolution, many researchers propose the inclusion of sophisticated auxiliary blocks, but these come at the cost of a considerable increase in network size and computational cost. This paper proposes stacked u-nets (SUNets), which iteratively combine features from different resolution scales while maintaining resolution. SUNets leverage the information globalization power of u-nets in a deeper network architectures that is capable of handling the complexity of natural images. SUNets perform extremely well on semantic segmentation tasks using a small number of parameters.


title: Semiconductor quantum well irradiated by a two-mode electromagnetic field as a terahertz emitter

abstract: We study theoretically the nonlinear optical properties of a semiconductor quantum well (QW) irradiated by a two-mode electromagnetic wave consisting of a strong resonant dressing field and a weak off-resonant driving field. In the considered strongly coupled electron-field system, the dressing field opens dynamic Stark gaps in the electron energy spectrum of the QW, whereas the driving field induces electron oscillations in the QW plane. Since the gapped electron spectrum restricts the amplitude of the oscillations, the emission of a frequency comb from the QW appears. Therefore, the doubly-driven QW operates as a nonlinear optical element which can be used, particularly, for optically controlled generation of terahertz radiation.


title: A note on the centers of a closed chain of circles

abstract: In this note we prove that the centers of a closed chain of circles for which every two consecutive members meet in the points of two given circles form a tangent polygon of a conic.


title: Nonlinear System Identification using Neural Networks and Trajectory-Based Optimization

abstract: In this paper, we study the identification of two challenging benchmark problems using neural networks. Two different global optimization approaches are used to train a recurrent neural network to identify two challenging nonlinear models, the cascaded tanks, and the Bouc-Wen system. The first approach, quotient gradient system (QGS), uses the trajectories of the nonlinear dynamical system to find the local minima of the optimization problem. The second approach, dynamical trajectory-based methodology, uses two different nonlinear dynamical systems to find the connected components of the feasible region and then searches the regions for local minima of the optimization problem. Simulation results show that both approaches effectively identify the model of the cascade tanks and the Bouc-Wen model.


title: Observation of exciton redshift-blueshift crossover in monolayer WS2

abstract: We report a rare atom-like interaction between excitons in monolayer WS2, measured using ultrafast absorption spectroscopy. At increasing excitation density, the exciton resonance energy exhibits a pronounced redshift followed by an anomalous blueshift. Using both material-realistic computation and phenomenological modeling, we attribute this observation to plasma effects and an attraction-repulsion crossover of the exciton-exciton interaction that mimics the Lennard-Jones potential between atoms. Our experiment demonstrates a strong analogy between excitons and atoms with respect to inter-particle interaction, which holds promise to pursue the predicted liquid and crystalline phases of excitons in two-dimensional materials.


title: An introduction to the two-dimensional Schr\"odinger equation with nonlinear point interactions

abstract: We present an introduction to the nonlinear Schr\"odinger equation (NLSE) with concentrated nonlinearities in $\mathbb{R}^2$. Precisely, taking a cue from the linear problem, we sketch the main challenges and the typical difficulties that arise in the two-dimensional case, and mention some recent results obtained by the authors on local and global well-posedness.


title: Measures of Noncompactness in $\left(\bar{N}_{\Delta^{-}}^{q}\right)$ Summable Difference Sequence Spaces

abstract: In the given paper we first introduce $\bar{N}_{\Delta^{-}}^{q}$ summable difference sequence spaces and prove some properties of these spaces. We then obtain the necessary and sufficient conditions for infinite matrices $A$ to map these sequence spaces on the spaces $c_0, c$ and $\ell_\infty$, the Hausdorff measure of noncompactness is then used to obtain the necessary and sufficient conditions for the compactness of the linear operators defined on these spaces.


title: Weak$^*$-sequential properties of Johnson-Lindenstrauss spaces

abstract: A Banach space $X$ is said to have Efremov's property ($\mathcal{E}$) if every element of the weak$^*$-closure of a convex bounded set $C \subseteq X^*$ is the weak$^*$-limit of a sequence in $C$. By assuming the Continuum Hypothesis, we prove that there exist maximal almost disjoint families of infinite subsets of $\mathbb{N}$ for which the corresponding Johnson-Lindenstrauss spaces enjoy (resp. fail) property ($\mathcal{E}$). This is related to a gap in [A. Plichko, Three sequential properties of dual Banach spaces in the weak$^*$ topology, Topology Appl. 190 (2015), 93--98] and allows to answer (consistently) questions of Plichko and Yost.


title: Applying the weighted horizontal magnetic gradient method to a simulated flaring Active Region

abstract: Here, we test the weighted horizontal magnetic gradient ($WG_M$) as a flare precursor, introduced by Korsos et. al. (2015, ApJ, 802, L21), by applying it to a magneto-hydrodynamic (MHD) simulation of solar-like flares Chatterjee et. al. (2016, Physical Review Letters, 116, 10, 101101). The pre-flare evolution of the $WG_M$ and the behavior of the distance parameter between the area-weighted barycenters of opposite polarity sunspots at various heights is investigated in the simulated $\delta$-type sunspot. Four flares emanated from this sunspot. We found the optimum heights above the photosphere where the flare precursors of the $WG_M$ method are identifiable prior to each flare. These optimum heights agree reasonably well with the heights of the occurrence of flares identified from the analysis of their thermal and Ohmic heating signatures in the simulation. We also estimated the expected time of the flare onsets from the duration of the approaching-receding motion of the barycenters of opposite polarities before each single flare. The estimated onset time and the actual time of occurrence of each flare are in good agreement at the corresponding optimum heights. This numerical experiment further supports the use of flare precursors based on the $WG_M$ method.


title: Recurrence Relations of the Multi-Indexed Orthogonal Polynomials V : Racah and $q$-Racah types

abstract: In previous papers, we discussed the recurrence relations of the multi-indexed orthogonal polynomials of the Laguerre, Jacobi, Wilson and Askey-Wilson types. In this paper we explore those of the Racah and $q$-Racah types. For the $M$-indexed ($q$-)Racah polynomials, we derive $3+2M$ term recurrence relations with variable dependent coefficients and $1+2L$ term ($L\geq M+1$) recurrence relations with constant coefficients. Based on the latter, the generalized closure relations and the creation and annihilation operators of the quantum mechanical systems described by the multi-indexed ($q$-)Racah polynomials are obtained. In appendix we present a proof and some data of the recurrence relations with constant coefficients for the multi-indexed Wilson and Askey-Wilson polynomials.


title: Subgame Perfect Equilibria of Sequential Matching Games

abstract: We study a decentralized matching market in which firms sequentially make offers to potential workers. For each offer, the worker can choose "accept" or "reject," but the decision is irrevocable. The acceptance of an offer guarantees her job at the firm, but it may also eliminate chances of better offers from other firms in the future. We formulate this market as a perfect-information extensive-form game played by the workers. Each instance of this game has a unique subgame perfect equilibrium (SPE), which does not necessarily lead to a stable matching and has some perplexing properties. We show a dichotomy result that characterizes the complexity of computing the SPE. The computation is tractable if each firm makes offers to at most two workers or each worker receives offers from at most two firms. In contrast, it is PSPACE-hard even if both firms and workers are related to at most three offers. We also study engineering aspects of this matching market. It is shown that, for any preference profile, we can design an offering schedule of firms so that the worker-optimal stable matching is realized in the SPE.


title: Attacks and Defenses in Mobile IP: Modeling with Stochastic Game Petri Net

abstract: The urging need for seamless connectivity in mobile environment has contributed to the rapid expansion of Mobile IP. Mobile IP uses wireless transmission medium, thereby making it subject to many security threats during various phases of route optimization. Modeling Mobile IP attacks reasonably and efficiently is the basis for defending against those attacks, which requires quantitative analysis and modeling approaches for expressing threat propagation in Mobile IP. In this Paper, we present four well-known Mobile IP attacks, such as Denial-of-Service (DoS) attack, bombing attack, redirection attack and replay attack and model them with Stochastic Game Petri Net (SGPN). Furthermore, we propose mixed strategy based defense strategies for the aforementioned attacks and model them with SGPN. Finally, we calculate the Nash Equilibrium of the attacker-defender game and thereby obtain the steady state probability of the vulnerable attack states. We show that, under the optimal strategy, an IDS needs to remain active 72.4%, 70%, 68.4% and 66.6% of the time to restrict the attacker's success rate to 8.5%, 6.4%, 7.2% and 8.3% respectively for the aforementioned attacks, thus performing better than the state-of-the-art approach.


title: History-Preserving Bisimulations on Reversible Calculus of Communicating Systems

abstract: History-and hereditary history-preserving bisimulation (HPB and HHPB) are equivalences relations for denotational models of concurrency. Finding their counterpart in process algebras is an open problem, with some partial successes: there exists in calculus of communicating systems (CCS) an equivalence based on causal trees that corresponds to HPB. In Reversible CSS (RCCS), there is a bisimulation that corresponds to HHPB, but it considers only processes without auto-concurrency. We propose equivalences on CCS with auto-concurrency that correspond to HPB and HHPB, and their so-called "weak" variants. The equivalences exploit not only reversibility but also the memory mechanism of RCCS.


title: Principle of Minimum Energy in Magnetic Reconnection in a Self-Organized Critical Model for Solar Flares

abstract: Solar flares are an abrupt release of magnetic energy in the Sun's atmosphere due to reconnection of the coronal magnetic field. This occurs in response to turbulent flows at the photosphere which twist the coronal field. Similar to earthquakes, solar flares represent the behavior of a complex system, and expectedly their energy distribution follows a power law. We present a statistical model based on the principle of minimum energy in a coronal loop undergoing magnetic reconnection, which is described as an avalanche process. We show the distribution of peaks for the flaring events in this self-organized critical system is scale-free. The obtained power law index 1.84 $ \pm $ 0.02 for the peaks is in good agreement with satellite observations of soft X-ray flares. The principle of minimum energy can be applied for general avalanche models to describe many other phenomena.


title: Multi frequency matching for voltage waveform tailoring

abstract: Customized voltage waveforms composed of a number of frequencies and used as the excitation of radio-frequency plasmas can control various plasma parameters such as energy distribution functions, homogeneity of the ionflux or ionization dynamics. So far this technology, while being extensively studied in academia, has yet to be established in applications. One reason for this is the lack of a suitable multi-frequency matching network that allows for maximum power absorption for each excitation frequency that is generated and transmitted via a single broadband amplifier. In this work, a method is introduced for designing such a network based on network theory and synthesis. Using this method, a circuit simulation is established that connects an exemplary matching network to an equivalent circuit plasma model of a capacitive radio-frequency discharge. It is found that for a range of gas pressures and number of excitation frequencies the matching conditions can be satisfied, which proves the functionality and feasibility of the proposed concept. Based on the proposed multi-frequency impedance matching, tailored voltage waveforms can be used at an industrial level.


title: Measuring Hopf links and Hopf invariants in a quenched topological Raman lattice

abstract: In a recent experimental work [Z. Wu et al., Science 354, 83 (2016)], the PKU-USTC group realized a two-dimensional two-band quantum anomalous Hall model on a square Raman lattice. By quenching the atom-laser detuning of such a Raman lattice, the time-dependent Bloch vectors for each quasimomentum points define a Hopf mapping from quasimomentum-time space $(k_x,k_y,t)\in T^3$ to the Bloch sphere $S^2$. The Hopf links between the preimages of any two Bloch vectors on $S^2$ can be measured experimentally through Bloch state tomography using spin-resolved time-of-flight measurement together with suitable radio-frequency manipulations. The dynamical Hopf invariants, which are quasimomentum-time integrations of Chern-Simons densities, can also be extracted experimentally through certain quench processes. As the Hopf invariant equals the Chern number of the post-quench Hamiltonian, the measurements of Hopf links and Hopf invariants provide an alternative way of understanding the topological phase diagram of the equilibrium system.


title: The edge spectrum of $K_4^-$-saturated graphs

abstract: Given graphs $G$ and $H$, $G$ is $H$-saturated if $G$ does not contain a copy of $H$ but the addition of any edge $e\notin E(G)$ creates at least one copy of $H$ within $G$. The edge spectrum of $H$ is the set of all possible sizes of an $H$-saturated graph on $n$ vertices. Let $K_4^-$ be a graph obtained from $K_4$ by deleting an edge. In this note, we show that (a) if $G$ is a $K_4^-$-saturated graph with $|V(G)|=n$ and $|E(G)|>\lfloor \frac{n-1}{2} \rfloor \lceil \frac{n-1}{2} \rceil +2$, then $G$ must be a bipartite graph; (b) there exists a $K_4^-$-saturated non-bipartite graph on $n\ge 10$ vertices with size being in the interval $\left[3n-11, \left\lfloor \frac{n-1}{2} \right\rfloor \left\lceil \frac{n-1}{2} \right\rceil +2\right]$. Together with a result of Fuller and Gould in [{\it On ($\hbox{K}_t-e$)-Saturated Graphs. Graphs Combin., 2018}], we determine the edge spectrum of $K_4^-$ completely, and a conjecture proposed by Fuller and Gould in the same paper also has been resolved.


title: The Category TOF

abstract: We provide a complete set of identities for the symmetric monoidal category, TOF, generated by the Toffoli gate and computational ancillary bits. We do so by demonstrating that the functor which evaluates circuits on total points, is an equivalence into the full subcategory of sets and partial isomorphisms with objects finite powers of the two element set. The structure of the proof builds -- and follows the proof of Cockett et al.-- which provided a full set of identities for the cnot gate with computational ancillary bits. Thus, first it is shown that TOF is a discrete inverse category in which all of the identities for the cnot gate hold; and then a normal form for the restriction idempotents is constructed which corresponds precisely to subobjects of the total points of TOF. This is then used to show that TOF is equivalent to FPinj2, the full subcategory of sets and partial isomorphisms in which objects have cardinality a power of 2.


title: An Element Sensitive Saliency Model with Position Prior Learning for Web Pages

abstract: Understanding human visual attention is important for multimedia applications. Many studies have attempted to learn from eye-tracking data and build computational saliency prediction models. However, limited efforts have been devoted to saliency prediction for Web pages, which are characterized by more diverse content elements and spatial layouts. In this paper, we propose a novel end-to-end deep generative saliency model for Web pages. To capture position biases introduced by page layouts, a Position Prior Learning sub-network is proposed, which models position biases as multivariate Gaussian distribution using variational auto-encoder. To model different elements of a Web page, a Multi Discriminative Region Detection (MDRD) branch and a Text Region Detection(TRD) branch are introduced, which target to extract discriminative localizations and "prominent" text regions likely to correspond to human attention, respectively. We validate the proposed model with FiWI, a public Web-page dataset, and shows that the proposed model outperforms the state-of-art models for Web-page saliency prediction.


title: Thermodynamic uncertainty relation of interacting oscillators in synchrony

abstract: The thermodynamic uncertainty relation sets the minimal bound of the cost-precision trade-off relation for dissipative processes. Examining the dynamics of an internally coupled system that is driven by a constant thermodynamic force, we however find that the trade-off relation of a sub-system is not constrained by the minimal bound of conventional uncertainty relation. We made our point explicit by using an exactly solvable model of interacting oscillators. As the number (N) of interacting oscillators increases, the uncertainty bound of individual oscillators is reduced to 2k_B T/N upon full synchronization under strong coupling. The cost-precision trade-off for the sub-system is particularly relevant for sub-cellular processes where collective dynamics emerges from multiple energy-expending components interacting with each other.


title: SaC2Vec: Information Network Representation with Structure and Content

abstract: Network representation learning (also known as information network embedding) has been the central piece of research in social and information network analysis for the last couple of years. An information network can be viewed as a linked structure of a set of entities. A set of linked web pages and documents, a set of users in a social network are common examples of information network. Network embedding learns low dimensional representations of the nodes, which can further be used for downstream network mining applications such as community detection or node clustering. Information network representation techniques traditionally use only the link structure of the network. But in real world networks, nodes come with additional content such as textual descriptions or associated images. This content is semantically correlated with the network structure and hence using the content along with the topological structure of the network can facilitate the overall network representation. In this paper, we propose Sac2Vec, a network representation technique that exploits both the structure and content. We convert the network into a multi-layered graph and use random walk and language modeling technique to generate the embedding of the nodes. Our approach is simple and computationally fast, yet able to use the content as a complement to structure and vice-versa. We also generalize the approach for networks having multiple types of content in each node. Experimental evaluations on four real world publicly available datasets show the merit of our approach compared to state-of-the-art algorithms in the domain.


title: Bayesian error regions in quantum estimation I: analytical reasonings

abstract: Results concerning the construction of quantum Bayesian error regions as a means to certify the quality of parameter point estimators have been reported in recent years. This task remains numerically formidable in practice for large dimensions and so far, no analytical expressions of the region size and credibility (probability of any given true parameter residing in the region) are known, which form the two principal region properties to be reported alongside a point estimator obtained from collected data. We first establish analytical formulas for the size and credibility that are valid for a uniform prior distribution over parameters, sufficiently large data samples and general constrained convex parameter-estimation settings. These formulas provide a means to an efficient asymptotic error certification for parameters of arbitrary dimensions. Next, we demonstrate the accuracies of these analytical formulas as compared to numerically computed region quantities with simulated examples in qubit and qutrit quantum-state tomography where computations of the latter are feasible.


title: Bayesian error regions in quantum estimation II: region accuracy and adaptive methods

abstract: Bayesian error analysis paves the way to the construction of credible and plausible error regions for a point estimator obtained from a given dataset. We introduce the concept of region accuracy for error regions (a generalization of the point-estimator mean squared-error) to quantify the average statistical accuracy of all region points with respect to the unknown true parameter. We show that the increase in region accuracy is closely related to the Bayesian-region dual operations in [1]. Next with only the given dataset as viable evidence, we establish various adaptive methods to maximize the region accuracy relative to the true parameter subject to the type of reported Bayesian region for a given point estimator. We highlight the performance of these adaptive methods by comparing them with nonadaptive procedures in three quantum-parameter estimation examples. The results of and mechanisms behind the adaptive schemes can be understood as the region analog of adaptive approaches to achieving the quantum Cramer--Rao bound for point estimators.


title: Online Convolutional Sparse Coding with Sample-Dependent Dictionary

abstract: Convolutional sparse coding (CSC) has been popularly used for the learning of shift-invariant dictionaries in image and signal processing. However, existing methods have limited scalability. In this paper, instead of convolving with a dictionary shared by all samples, we propose the use of a sample-dependent dictionary in which filters are obtained as linear combinations of a small set of base filters learned from the data. This added flexibility allows a large number of sample-dependent patterns to be captured, while the resultant model can still be efficiently learned by online learning. Extensive experimental results show that the proposed method outperforms existing CSC algorithms with significantly reduced time and space requirements.


title: Autonomous Vehicles for Smart and Sustainable Cities: An In-Depth Exploration of Privacy and Cybersecurity Implications

abstract: Amidst rapid urban development, sustainable transportation solutions are required to meet the increasing demands for mobility whilst mitigating the potentially negative social, economic, and environmental impacts. This study analyses autonomous vehicles (AVs) as a potential transportation solution for smart and sustainable development. We identified privacy and cybersecurity risks of AVs as crucial to the development of smart and sustainable cities and examined the steps taken by governments around the world to address these risks. We highlight the literature that supports why AVs are essential for smart and sustainable development. We then identify the aspects of privacy and cybersecurity in AVs that are important for smart and sustainable development. Lastly, we review the efforts taken by federal governments in the US, the UK, China, Australia, Japan, Singapore, South Korea, Germany, France, and the EU, and by US state governments to address AV-related privacy and cybersecurity risks in-depth. Overall, the actions taken by governments to address privacy risks are mainly in the form of regulations or voluntary guidelines. To address cybersecurity risks, governments have mostly resorted to regulations that are not specific to AVs and are conducting research and fostering research collaborations with the private sector.


title: Explicit lower bounds on strong quantum simulation

abstract: We consider the problem of strong (amplitude-wise) simulation of $n$-qubit quantum circuits, and identify a subclass of simulators we call monotone. This subclass encompasses almost all prominent simulation techniques. We prove an unconditional (i.e. without relying on any complexity theoretic assumptions) and explicit $(n-2)(2^{n-3}-1)$ lower bound on the running time of simulators within this subclass. Assuming the Strong Exponential Time Hypothesis (SETH), we further remark that a universal simulator computing any amplitude to precision $2^{-n}/2$ must take at least $2^{n - o(n)}$ time. Finally, we compare strong simulators to existing SAT solvers, and identify the time-complexity below which a strong simulator would improve on state-of-the-art SAT solving.


title: Three-Parameter Solutions of the PV Schlesinger-Type Equation near the Point at Infinity and the Monodromy Data

abstract: For the Schlesinger-type equation related to the fifth Painlev\'e equation (V) via isomonodromy deformation, we present a three-parameter family of matrix solutions along the imaginary axis near the point at infinity, and also the corresponding monodromy data. Two-parameter solutions of (V) with their monodromy data immediately follow from our results. Under certain conditions, these solutions of (V) admit sequences of zeros and of poles along the imaginary axis. The monodromy data are obtained by matching techniques for a perturbed linear system.


title: Refined enumerations of alternating sign triangles

abstract: This article introduces and investigates a refinement of alternating sign trapezoids by means of Catalan objects and Motzkin paths. Alternating sign trapezoids are a generalisation of alternating sign triangles, which were recently introduced by Ayyer, Behrend and Fischer. We show that the number of alternating sign trapezoids associated with a Catalan object (resp. Motzkin path) is a polynomial function in the length of the shorter base of the trapezoid. We also study the rational roots of the polynomials and formulate several conjectures and derive some partial results. Lastly, we deduce a constant term identity for the refined counting of alternating sign trapezoids.


title: dhSegment: A generic deep-learning approach for document segmentation

abstract: In recent years there have been multiple successful attempts tackling document processing problems separately by designing task specific hand-tuned strategies. We argue that the diversity of historical document processing tasks prohibits to solve them one at a time and shows a need for designing generic approaches in order to handle the variability of historical series. In this paper, we address multiple tasks simultaneously such as page extraction, baseline extraction, layout analysis or multiple typologies of illustrations and photograph extraction. We propose an open-source implementation of a CNN-based pixel-wise predictor coupled with task dependent post-processing blocks. We show that a single CNN-architecture can be used across tasks with competitive results. Moreover most of the task-specific post-precessing steps can be decomposed in a small number of simple and standard reusable operations, adding to the flexibility of our approach.


title: Splitting gradient algorithms for solving monotone equilibrium problems

abstract: It is well known that the projection method is not convergent for monotone equilibrium problems. Recently Sosa \textit{et al.} in \cite{SS2011} proposed a projection algorithm ensuring convergence for paramonotone equilibrium problems. In this paper we modify this algorithm to obtain a splitting convergent one for the case when the bifunction is the sum of the two ones. At each iteration, two strongly convex subprograms are required to solve separately, one for each component bifunction. We show that the algorithm is convergent for paramonotone bifunction without any Lipschitz type condition as well as H\"older continuity of the involved bifunctions. Furthermore, we show that the ergodic sequence defined by the algorithm's iterates converges to a solution without paramonotonicity property. We use the proposed algorithm to solve a jointly constrained Cournot-Nash model. The computational results show that this algorithm is efficient for the model with a restart strategy.


title: Models of Total, Elastic and Diffractive Cross Sections

abstract: The LHC has brought much new information on total, elastic and diffractive cross sections, which is not always in agreement with extrapolations from lower energies. The default framework in the Pythia event generator is one case in point. In this article we study and implement two recent models, as more realistic alternatives. Both describe total and elastic cross sections, whereas one also includes single diffraction. Noting some issues at high energies, a variant of the latter is proposed, and extended also to double and central diffraction. Further, the experimental definition of diffraction is based on the presence of rapidity gaps, which however also could be caused by colour reconnection in nondiffractive events, a phenomenon that is studied in the context of a specific model. Throughout comparisons with LHC and other data are presented.


title: Constraints on tidal charge of the supermassive black hole at the Galactic Center with trajectories of bright stars

abstract: As it was pointed out recently in Hees et al. (2017), observations of stars near the Galactic Center with current and future facilities provide an unique tool to test general relativity (GR) and alternative theories of gravity in a strong gravitational field regime. In particular, the authors showed that the Yukawa gravity could be constrained with Keck and TMT observations. Some time ago, Dadhich et al. (2001) showed that the Reissner -- Nordstr\"om metric with a tidal charge is naturally appeared in the framework of Randall -- Sundrum model with an extra dimension ($Q^2$ is called tidal charge and it could be negative in such an approach). Astrophysical consequences of of presence of black holes with a tidal charge are considerered, in particular, geodesics and shadows in Kerr -- Newman braneworld metric are analyzed in (Schee and Stuchlik, 2009a), while profiles of emission lines generated by rings orbiting braneworld Kerr black hole are considered in (Schee and Stuchlik, 2009b). Possible observational signatures of gravitational lensing in a presence of the Reissner -- Nordstr\"om black hole with a tidal charge at the Galactic Center are discussed in papers by Bin-Nun (2010a, 2010b, 2011). Here we are following such an approach and we obtain analytical expressions for orbital precession for Reissner -- Nordstrom -- de-Sitter solution in post-Newtonian approximation and discuss opportunities to constrain parameters of the metric from observations of bright stars with current and future astrometric observational facilities such as VLT, Keck, GRAVITY, E-ELT and TMT.


title: On interrelations between divergence-free and Hamiltonian dynamics

abstract: A mathematically correct description is presented on the interrelations between the dynamics of divergence free vector fields on an oriented 3-dimensional manifold $M$ and the dynamics of Hamiltonian systems. It is shown that for a given divergence free vector field $X$ with a global cross-section there exist some 4-dimensional symplectic manifold $\tilde{M}\supset M$ and a smooth Hamilton function $H: \tilde{M}\to \mathbb R$ such that for some $c\in \mathbb R$ one gets $M = \{H=c\}$ and the Hamiltonian vector field $X_H$ restricted on this level coincides with $X$. For divergence free vector fields with singular points such the extension is impossible but the existence of local cross-section allows one to reduce the dynamics to the study of symplectic diffeomorphisms in some sub-domains of $M$. We also consider the case of a divergence free vector field $X$ with a smooth integral having only finite number of critical levels. It is shown that such a noncritical level is always a 2-torus and restriction of $X$ on it possesses a smooth invariant 2-form. The linearization of the flow on such a torus (i.e. the reduction to the constant vector field) is not always possible in contrast to the case of an integrable Hamiltonian system but in the analytic case ($M$ and $X$ are real analytic), due to the Kolmogorov theorem, such the linearization is possible for tori with Diophantine rotation numbers.


title: Compact gravimeter with an ensemble of ultracold atoms in spin-dependent optical lattices

abstract: Atomic interferometry in optical lattices is a new trend of developing practical quantum gravimeter. Here, we propose a compact and portable gravimetry scheme with an ensemble of ultracold atoms in gravitationally tilted spin-dependent optical lattices. The fast, coherent separation and recombination of atoms can be realized via polarization-synthesized optical lattices. The input atomic wavepacket is coherently split into two parts by a spin-dependent shift and a subsequent $\frac{\pi}{2}$ pulse. Then the two parts are held for accumulating a relative phase related to the gravity. Lastly the two parts are recombined for interference by a $\frac{\pi}{2}$ pulse and a subsequent spin-dependent shift. The $\frac{\pi}{2}$ pulses not only preclude the spin-dependent energies in the accumulated phase, but also avoid the error sources such as dislocation of optical lattices in the holding process. In addition, we develop an analytical method for the sensitivity in multi-path interferometry.


title: Differentiability Properties of a Parametric Consumer Problem

abstract: We study the budget map and the indirect utility function of a parametric consumer problem in a Banach space setting by some advanced tools from set-valued and variational analysis. The Lipschitz-likeness and differentiability properties of the budget map, as well as formulas for finding subdifferentials of the infimal nuisance function, which is obtained from the indirect utility function by changing its sign, are established. Our investigation is mainly based on the paper by Mordukhovich [J. Global Optim. 28 (2004), 347--362] on coderivative analysis of variational systems and the paper of Mordukhovich, Nam, and Yen [Math. Program. 116 (2009), 369--396] on subgradients of marginal functions. Economic meanings of the obtained subdifferential estimates are explained in details.


title: Axion-Photon Conversion and Effects on 21cm Observation

abstract: Recently the EDGES experiment reported an enhanced 21cm absorption signal in the radio wave observation, which may be interpreted as either anomalous cooling of baryons or heating of cosmic microwave background photons. In this paper, we pursue the latter possibility. We point out that dark radiation consisting of axion-like particles can resonantly convert into photons under the intergalactic magnetic field, which can effectively heat up the radiation in the frequency range relevant for the EDGES experiment. This may explain the EDGES anomaly.


title: Riemannian optimal identification method for linear systems with symmetric positive-definite matrix

abstract: This study develops identification methods for linear continuous-time symmetric systems, such as electrical network systems, multi-agent network systems, and temperature dynamics in buildings. To this end, we formulate three system identification problems for the corresponding discrete-time systems. The first is a least-squares problem in which we wish to minimize the sum of squared errors between the true and model outputs on the product manifold of the manifold of symmetric positive-definite matrices and two Euclidean spaces. In the second problem, to reduce the search dimensions, the product manifold is replaced with the quotient set under a specified group action by the orthogonal group. In the third problem, the manifold of symmetric positive-definite matrices in the first problem is replaced by the manifold of matrices with only positive diagonal elements. In particular, we examine the quotient geometry in the second problem. We propose Riemannian conjugate gradient methods for the three problems, and select initial points using a popular subspace method. The effectiveness of our proposed methods is demonstrated through numerical simulations and comparisons with the Gauss--Newton method, which is one of the most popular approach for solving least-squares problems.


title: Testing the Kerr nature of the supermassive black hole in Ark 564

abstract: Einstein's theory of general relativity has been extensively tested in weak gravitational fields, mainly with experiments in the Solar System and observations of radio pulsars, and current data agree well with the theoretical predictions. Nevertheless, there are a number of scenarios beyond Einstein's gravity that have the same predictions for weak fields and present deviations only when gravity becomes strong. Here we try to test general relativity in the strong field regime. We fit the X-ray spectrum of the supermassive black hole in Ark 564 with a disk reflection model beyond Einstein's gravity, and we are able to constrain the black hole spin $a_*$ and the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$ separately. For $\alpha_{22} = 0$, we find $a_* > 0.96$ and $-1.0 < \alpha_{13} < 0.2$ with a 99% confidence level. For $\alpha_{13} = 0$, we get $a_* > 0.96$ and $-0.1 < \alpha_{22} < 0.9$ with a 99% confidence level. Our measurements are thus consistent with the hypothesis that the supermassive compact object in Ark 564 can be described by the Kerr metric.


title: Development of Rehabilitation System (ReHabgame) through Monte-Carlo Tree Search Algorithm

abstract: Computational Intelligence (CI) in computer games plays an important role that could simulate various aspects of real-life problems. CI in real-time decision-making games can provide a platform for the examination of tree search algorithms. In this paper, we present a rehabilitation serious game (ReHabgame) in which the Monte-Carlo Tree Search (MCTS) algorithm is utilized. The game is designed to combat the physical impairment of post-stroke/brain injury casualties in order to improve upper limb movement. Through the process of ReHabgame the player chooses paths via upper limb according to his/her movement ability to reach virtual goal objects. The system adjusts the difficulty level of the game based on the player's quality of activity through MCTS. It learns from the movements made by a player and generates further subsequent objects for collection. The system collects orientation, muscle and joint activity data and utilizes them to make decisions. Players data are collected through Kinect Xbox One and Myo Armband. The results show the effectiveness of the MCTS in the ReHabgame that progresses from highly achievable paths to the less achievable ones, thus configuring and personalizing the rehabilitation process.


title: E and B families of the Stokes parameters in the polarized synchrotron and thermal dust foregrounds

abstract: Better understanding of Galactic foregrounds is one of the main obstacles to detection of primordial gravitational waves through measurement of the B mode in the polarized microwave sky. We generalize the method proposed in \cite{2011PhRvD..83h3003B} and decompose the polarization signals into the E and B families directly in the domain of the Stokes $Q$, $U$ parameters as $(Q,U)\equiv(Q_E, U_E)+(Q_B,U_B)$. This also enables an investigation of the morphology and the frequency dependence of these two families, which has been done in the WMAP K, Ka (tracing synchrotron emission) and Planck 2015 HFI maps (tracing thermal dust). The results reveal significant differences in spectra between the E and B families. The spectral index of the E family fluctuates less across the sky than that of the B family, and the same tendency occurs for the polarization angles of the dust and synchrotron channels. The new insight from WMAP and Planck data on the North Polar Spur and BICEP2 zones through our method clearly indicates that these zones are characterized by very low polarization intensity of the B family compared to the E family. We have detected global structure of the B family polarization angles at high Galactic latitudes which cannot be attributed to the cosmic microwave background or instrumental noise. However, we cannot exclude instrumental systematics as a partial contributor to these anomalies.


title: Games and Ramsey-like cardinals

abstract: We generalise the $\alpha$-Ramsey cardinals introduced in Holy and Schlicht (2018) for cardinals $\alpha$ to arbitrary ordinals $\alpha$, and answer several questions posed in that paper. In particular, we show that $\alpha$-Ramseys are downwards absolute to the core model $K$ for all $\alpha$ of uncountable cofinality, that strategic $\omega$-Ramsey cardinals are equiconsistent with remarkable cardinals and that strategic $\alpha$-Ramsey cardinals are equiconsistent with measurable cardinals for all $\alpha>\omega$. We also show that the $n$-Ramseys satisfy indescribability properties and use them to provide a game-theoretic characterisation of completely ineffable cardinals, as well as establishing further connections between the $\alpha$-Ramsey cardinals and the Ramsey-like cardinals introduced in Gitman (2011), Feng (1990) and Sharpe and Welch (2011).


title: Unique excitonic effects in tungsten disulphide monolayers on two-layer graphene

abstract: Light emission in atomically thin heterostructures is known to depend on the type of materials, number and stacking sequence of the constituent layers. Here we show that the thickness of a two-dimensional substrate can be crucial in modulating the light emission. We study the layer-dependent charge transfer in vertical heterostructures built from monolayer tungsten disulphide (WS2) on one- and two-layer epitaxial graphene, unravelling the effect that the interlayer electronic coupling has on the excitonic properties of such heterostructures. We bring evidence that the excitonic properties of WS2 can be effectively tuned by the number of supporting graphene layers. Integrating WS2 monolayers with two-layer graphene leads to a significant enhancement of the photoluminescence response, up to one order of magnitude higher compared to WS2 supported on one-layer graphene. Our findings highlight the importance of substrate engineering when constructing atomically thin layered heterostructures.


title: Self-interacting dark matter with a stable vector mediator

abstract: Light vector mediators can naturally induce velocity-dependent dark matter self-interactions while at the same time allowing for the correct dark matter relic abundance via thermal freeze-out. If these mediators subsequently decay into Standard Model states such as electrons or photons however, this is robustly excluded by constraints from the Cosmic Microwave Background. We study to what extent this conclusion can be circumvented if the vector mediator is stable and hence contributes to the dark matter density while annihilating into lighter degrees of freedom. We find viable parts of parameter space which lead to the desired self-interaction cross section of dark matter to address the small-scale problems of the collisionless cold dark matter paradigm while being compatible with bounds from the Cosmic Microwave Background and Big Bang Nucleosynthesis observations.


title: Trudinger-Moser inequalities on a closed Riemannian surface with the action of a finite isometric group

abstract: Let $(\Sigma,g)$ be a closed Riemannian surface, $W^{1,2}(\Sigma,g)$ be the usual Sobolev space, $\textbf{G}$ be a finite isometric group acting on $(\Sigma,g)$, and $\mathscr{H}_\textbf{G}$ be a function space including all functions $u\in W^{1,2}(\Sigma,g)$ with $\int_\Sigma udv_g=0$ and $u(\sigma(x))=u(x)$ for all $\sigma\in \textbf{G}$ and all $x\in\Sigma$. Denote the number of distinct points of the set $\{\sigma(x): \sigma\in \textbf{G}\}$ by $I(x)$ and $\ell=\inf_{x\in \Sigma}I(x)$. Let $\lambda_1^\textbf{G}$ be the first eigenvalue of the Laplace-Beltrami operator on the space $\mathscr{H}_\textbf{G}$. Using blow-up analysis, we prove that if $\alpha<\lambda_1^\textbf{G}$ and $\beta\leq 4\pi\ell$, then there holds $$\sup_{u\in\mathscr{H}_\textbf{G},\,\int_\Sigma|\nabla_gu|^2dv_g-\alpha \int_\Sigma u^2dv_g\leq 1}\int_\Sigma e^{\beta u^2}dv_g<\infty;$$ if $\alpha<\lambda_1^\textbf{G}$ and $\beta>4\pi\ell$, or $\alpha\geq \lambda_1^\textbf{G}$ and $\beta>0$, then the above supremum is infinity; if $\alpha<\lambda_1^\textbf{G}$ and $\beta\leq 4\pi\ell$, then the above supremum can be attained. Moreover, similar inequalities involving higher order eigenvalues are obtained. Our results partially improve original inequalities of J. Moser \cite{Moser}, L. Fontana \cite{Fontana} and W. Chen \cite{Chen-90}.


title: Cohomology and Deformations of $n$-Lie algebra morphisms

abstract: The study of $n$-Lie algebras which are natural generalization of Lie algebras is motivated by Nambu Mechanics and recent developments in String Theory and M-branes. The purpose of this paper is to define cohomology complexes and study deformation theory of $n$-Lie algebra morphisms. We discuss infinitesimal deformations, equivalent deformations and obstructions. Moreover, we provide various examples.


title: Event Forecasting with Pattern Markov Chains

abstract: We present a system for online probabilistic event forecasting. We assume that a user is interested in detecting and forecasting event patterns, given in the form of regular expressions. Our system can consume streams of events and forecast when the pattern is expected to be fully matched. As more events are consumed, the system revises its forecasts to reflect possible changes in the state of the pattern. The framework of Pattern Markov Chains is used in order to learn a probabilistic model for the pattern, with which forecasts with guaranteed precision may be produced, in the form of intervals within which a full match is expected. Experimental results from real-world datasets are shown and the quality of the produced forecasts is explored, using both precision scores and two other metrics: spread, which refers to the "focusing resolution" of a forecast (interval length), and distance, which captures how early a forecast is reported.


title: Local module identification in dynamic networks: do more inputs guarantee smaller variance?

abstract: Recent developments in science and engineering have motivated control systems to be considered as interconnected and networked systems. From a system identification point of view, modelling of a local module in such a structured system is a relevant and interesting problem. This work focuses on the quality, in terms of variance, of an estimate of a local module. We analyse which predictor input signals are relevant and contribute to variance reduction, while still guaranteeing the consistency of the estimate. For a targeted local module, a comparison of its estimate variance is made between a full-MISO approach and an immersed network setting, where a reduced number of inputs is used, while still guaranteeing consistency. A case study of a four-node network is considered and it is shown that a smaller set of predictor inputs can, under some conditions, result in a smaller variance compared to the full-MISO approach.


title: Automatic classification of trees using a UAV onboard camera and deep learning

abstract: Automatic classification of trees using remotely sensed data has been a dream of many scientists and land use managers. Recently, Unmanned aerial vehicles (UAV) has been expected to be an easy-to-use, cost-effective tool for remote sensing of forests, and deep learning has attracted attention for its ability concerning machine vision. In this study, using a commercially available UAV and a publicly available package for deep learning, we constructed a machine vision system for the automatic classification of trees. In our method, we segmented a UAV photography image of forest into individual tree crowns and carried out object-based deep learning. As a result, the system was able to classify 7 tree types at 89.0% accuracy. This performance is notable because we only used basic RGB images from a standard UAV. In contrast, most of previous studies used expensive hardware such as multispectral imagers to improve the performance. This result means that our method has the potential to classify individual trees in a cost-effective manner. This can be a usable tool for many forest researchers and managements.


title: Kernels of block Hankel operators and independency of vector-valued functions modulo Nevanlinna class

abstract: For a matrix-valued function $\Phi\in L^2_{M_{n\times m}}$, it is well-known that the kernel of a block Hankel operator $H_\Phi$ is an invariant subspace for the shift operator. Thus, if the kernel is nontrivial, then $\ker H_\Phi= \Theta H^2_{\mathbb C^r}$ for a natural number $r$ and an $m\times r$ matrix inner function $\Theta$ by Beurling-Lax-Halmos Theorem. It will be shown that the size of the matrix inner function $\Theta$ associated with the kernel of a block Hankel operator $H_\Phi$ is closely related with a certain independency of the columns of $\Phi$, which is defined in this paper. As an important application of this result, the shape of shift invariant, or, backward shift invariant subspaces of $H^2_{\mathbb C^n}$ generated by finite elements will be studied.


title: An adaptive self-organizing fuzzy logic controller in a serious game for motor impairment rehabilitation

abstract: Rehabilitation robotics combined with video game technology provides a means of assisting in the rehabilitation of patients with neuromuscular disorders by performing various facilitation movements. The current work presents ReHabGame, a serious game using a fusion of implemented technologies that can be easily used by patients and therapists to assess and enhance sensorimotor performance and also increase the activities in the daily lives of patients. The game allows a player to control avatar movements through a Kinect Xbox, Myo armband and rudder foot pedal, and involves a series of reach-grasp-collect tasks whose difficulty levels are learnt by a fuzzy interface. The orientation, angular velocity, head and spine tilts and other data generated by the player are monitored and saved, whilst the task completion is calculated by solving an inverse kinematics algorithm which orientates the upper limb joints of the avatar. The different values in upper body quantities of movement provide fuzzy input from which crisp output is determined and used to generate an appropriate subsequent rehabilitation game level. The system can thus provide personalised, autonomously-learnt rehabilitation programmes for patients with neuromuscular disorders with superior predictions to guide the development of improved clinical protocols compared to traditional theraputic activities.


title: Quantum interference among heavy NMSSM Higgs bosons

abstract: In the Next-to-Minimal Supersymmetric Standard Model (NMSSM), it is possible to have strong mass degeneracies between the new singlet-like scalar and the heavy doublet-like scalar, as well as between the singlet-like and doublet-like pseudoscalar Higgs states. When the difference in the masses of such states is comparable with the sum of their widths, the quantum mechanical interference between their propagators can become significant. We study these effects by taking into account the full Higgs boson propagator matrix in the calculation of the production process of $\tau^+\tau^-$ pairs in gluon fusion at the Large Hadron Collider (LHC). We find that, while these interference effects are sizeable, they are not resolvable in terms of the distributions of differential cross sections, owing to the poor detector resolution of the $\tau^+\tau^-$ invariant mass. They are, however, identifiable via the inclusive cross sections, which are subject to significant variations with respect to the standard approaches, wherein the propagating Higgs bosons are treated independently from one another. We quantify these effects for several representative benchmark points, extracted from a large set of points, obtained by numerical scanning of the NMSSM parameter space, that satisfy the most important experimental constraints currently available.


title: Why are classical bulges more common in S0 galaxies than in spiral galaxies?

abstract: In this paper, we try to understand why the classical bulge fraction observed in S0 galaxies is significantly higher than that in spiral galaxies. We carry out a comparative study of the bulge and global properties of a sample of spiral and S0 galaxies in a fixed environment. Our sample is flux limited and contains 262 spiral and 155 S0 galaxies drawn from the Sloan Digital Sky Survey. We have classified bulges into classical and pseudobulge categories based on their position on the Kormendy diagram. Dividing our sample into bins of galaxy stellar mass, we find that the fraction of S0 galaxies hosting a classical bulge is significantly higher than the classical bulge fraction seen in spirals even at fixed stellar mass. We have compared the bulge and the global properties of spirals and S0 galaxies in our sample and find indications that spiral galaxies which host a classical bulge, preferentially get converted into S0 population as compared to pseudobulge hosting spirals. By studying the star formation properties of our galaxies in the NUV-r color-mass diagram, we find that the pseudobulge hosting spirals are mostly star forming while the majority of classical bulge host spirals are in the green valley or in the passive sequence. We suggest that some internal process, such as AGN feedback or morphological quenching due to the massive bulge, quenches these classical bulge hosting spirals and transforms them into S0 galaxies, thus resulting in the observed predominance of the classical bulge in S0 galaxies.


title: Lethal radiation from nearby supernovae helps to explain the small cosmological constant

abstract: The observed value $\Lambda_{\rm obs}$ of the cosmological constant $\Lambda$ is extremely smaller than theoretical expectations, and the anthropic argument has been proposed as a solution to this problem because galaxies do not form when $\Lambda \gg \Lambda_{\rm obs}$. However, the contemporary galaxy formation theory predicts that stars form even with a high value of $\Lambda / \Lambda_{\rm obs} \sim$ 50, which makes the anthropic argument less persuasive. Here we calculate the probability distribution of $\Lambda$ using a model of cosmological galaxy formation, considering extinction of observers caused by radiation from nearby supernovae. The life survival probability decreases in a large $\Lambda$ universe because of higher stellar density. Using a reasonable rate of lethal supernovae, we find that the mean expectation value of $\Lambda$ can be close to $\Lambda_{\rm obs}$, and hence this effect may be essential to understand the small but nonzero value of $\Lambda$. It is predicted that we are located on the edge of habitable regions about stellar density in the Galaxy, which may be tested by future exoplanet studies.

categories: astro-ph.CO astro-ph.EP astro-ph.GA astro-ph.HE


title: Average Case Analysis of Leaf-Centric Binary Tree Sources

abstract: We study the average number of distinct fringe subtrees in random trees generated by leaf-centric binary tree sources as introduced by Zhang, Yang and Kieffer. A leaf-centric binary tree source induces for every $n \geq 2$ a probability distribution on the set of binary trees with $n$ leaves. We generalize a result by Flajolet, Gourdon, Martinez and Devroye, according to which the average number of distinct fringe subtrees in a random binary search tree of size $n$ is in $\Theta(n/\log n)$, as well as a result by Flajolet, Sipala and Steayert, according to which the number of distinct fringe subtrees in a uniformly random binary tree of size $n$ is in $\Theta(n/\sqrt{\log n})$.


title: Implicit Copulas from Bayesian Regularized Regression Smoothers

abstract: We show how to extract the implicit copula of a response vector from a Bayesian regularized regression smoother with Gaussian disturbances. The copula can be used to compare smoothers that employ different shrinkage priors and function bases. We illustrate with three popular choices of shrinkage priors --- a pairwise prior, the horseshoe prior and a g prior augmented with a point mass as employed for Bayesian variable selection --- and both univariate and multivariate function bases. The implicit copulas are high-dimensional, have flexible dependence structures that are far from that of a Gaussian copula, and are unavailable in closed form. However, we show how they can be evaluated by first constructing a Gaussian copula conditional on the regularization parameters, and then integrating over these. Combined with non-parametric margins the regularized smoothers can be used to model the distribution of non-Gaussian univariate responses conditional on the covariates. Efficient Markov chain Monte Carlo schemes for evaluating the copula are given for this case. Using both simulated and real data, we show how such copula smoothing models can improve the quality of resulting function estimates and predictive distributions.


title: Re-examining the quadratic approximation in theory of a weakly interacting Bose gas with condensate: the role of nonlocal interaction potentials

abstract: We derive and analyze the coupled equations of quadratic approximation of the Bogoliubov model for a weakly interacting Bose gas. The first equation determines the condensate density as a variational parameter and ensures the minimum of the grand thermodynamic potential. The second one provides a relation between the total number of particles and chemical potential. Their consistent theoretical analysis is performed for a number of model interaction potentials including contact (local) and nonlocal interactions, where the latter provide nontrivial dependencies in momentum space. We demonstrate that the derived equations have no solutions for the local potential, although they formally reproduce the well-known results of the Bogoliubov approach. At the same time, it is shown that these equations have the solutions for physically relevant nonlocal potentials. We show that in the regimes close to experimental realizations with ultracold atoms, the contribution of the terms originating from the quadratic part of the truncated Hamiltonian to the chemical potential can be of the same order of magnitude as from its $c$-number part. Due to this fact, in particular, the spectrum of single-particle excitations in the quadratic approximation acquires a gap. The issue of the gap is also discussed.


title: Tunable Electronic Structure and Magnetic Coupling in Strained Two-Dimensional Semiconductor MnPSe3

abstract: The electronic structures and magnetic properties of strained monolayer MnPSe3 are investigated systematically by first-principles calculations. It is found that the magnetic ground state (GS) of monolayer MnPSe3 can be significantly affected by biaxial strain engineering, while the semiconducting characteristics are well preserved. Owing to the sensitivity of the magnetic coupling towards the structural deformation, a biaxial tensile strain about 13% can lead to an antiferromagnetic-ferromagnetic (AFM-FM) transition. The underlying physical mechanism of strain-dependent magnetic stability is mainly attributed to the competition effect of direct AFM interaction and indirect FM superexchange interaction between the nearest-neighbor (NN) two Mn atoms. In addition, we find that FM MnPSe3 is an intrinsic half semiconductor with a large spin exchange splitting in conduction bands, which is crucial for the spin-polarized carrier injection and detection. The sensitive interdependence among external stimuli, electronic structure and magnetic coupling suggests that monolayer MnPSe3 can be a promising candidate in spintronics.


title: Phase-Matched Second Harmonic Generation in a flux grown KTiOPO4 crystal ridge optical waveguide

abstract: Type II second-harmonic generation was performed in a 15.8-mm-long KTiOPO4 micrometric ridge waveguide with an average transversal section of 38 micrometer^2. Theoretical predictions are compared with experiments. Strong agreements are obtained for both phase-matching wavelengths and second-harmonic intensity. This work opens wide perspectives for integrated parametric optics.


title: A Technique for Estimating the Absolute Gain of a Photomultiplier Tube

abstract: Detection of low-intensity light relies on the conversion of photons to photoelectrons, which are then multiplied and detected as an electrical signal. To measure the actual intensity of the light, one must know the factor by which the photoelectrons have been multiplied. To obtain this amplification factor, we have developed a procedure for estimating precisely the signal caused by a single photoelectron. The method utilizes the fact that the photoelectrons conform to a Poisson distribution. The average signal produced by a single photoelectron can then be estimated from the number of noise events, without requiring analysis of the distribution of the signal produced by a single photoelectron. The signal produced by one or more photoelectrons can be estimated experimentally without any assumptions. This technique, and an example of the analysis of a signal from a photomultiplier tube, are described in this study.


title: Information and Statistical Efficiency When Quantizing Noisy DC Values

abstract: This paper considers estimation of a quantized constant in noise when using uniform and nonuniform quantizers. Estimators based on simple arithmetic averages, on sample statistical moments and on the maximum-likelihood procedure are considered. It provides expressions for the statistical efficiency of the arithmetic mean by comparing its variance to the proper Cram\'er-Rao lower bound. It is conjectured that the arithmetic mean is optimal among all estimators with an exactly known bias. Conditions under which its statistical performance are improved by the other estimation procedures when the exact bias is not known are found and analyzed. Using simulations and analysis of experimental data, it is shown that both moment-based and maximum-likelihood-based estimators provide better results, when the noise standard deviation is comparable with the quantization step and the noise model of quantization can not be applied.


title: Well-posedness and stability results for some periodic Muskat problems

abstract: We study the two-dimensional Muskat problem in a horizontally periodic setting and for fluids with arbitrary densities and viscosities. We show that in the presence of surface tension effects the Muskat problem is a quasilinear parabolic problem which is well-posed in the Sobolev space $H^r(\mathbb{S})$ for each $r\in(2,3)$. When neglecting surface tension effects, the Muskat problem is a fully nonlinear evolution equation and of parabolic type in the regime where the Rayleigh-Taylor condition is satisfied. We then establish the well-posedness of the Muskat problem in the open subset of $H^2(\mathbb{S})$ defined by the Rayleigh-Taylor condition. Besides, we identify all equilibrium solutions and study the stability properties of trivial and of small finger-shaped equilibria. Also other qualitative properties of solutions such as parabolic smoothing, blow-up behavior, and criteria for global existence are outlined.


title: Noncommutative fluid and Growing Modes of Inhomogeneity in (Newtonian) Cosmology

abstract: Time evolution of modes of density contrast, in particular the growing modes, dictate the structure formation in Universe. In this paper we explicitly show how (spatial) Non-Commutativity (NC) can affect the behavior of the modes, that is we compute NC corrected power law profiles of the density contrast modes. We develop a generalized fluid model that lives in NC space. The dynamical equations of fluid, namely the continuity and Euler equations receive NC contributions. When mapped to comoving coordinates these generate the NC extended versions of continuity and Friedmann equations for cosmology. Introducing cosmological perturbations finally yield the NC corrected evolution of density contrast modes. The construction of the NC fluid model from first principles and development of the formal aspects of its Hamiltonian formulation have been presented in the first part.


title: Hausdorff dimension of limsup sets of rectangles in the Heisenberg group

abstract: The almost sure value of the Hausdorff dimension of limsup sets generated by randomly distributed rectangles in the Heisenberg group is computed in terms of directed singular value functions.


title: A rational parametrization of B\'ezier like curves

abstract: In this paper, we construct a family of Bernstein functions using a class of rational parametrization. The new family of rational Bernstein basis on an index $\alpha \in {\left(-\infty \, , \, 0 \right)}\cup {\left(1 \, , \, +\infty\right)}$, and for a given degree $k\in \mathbb{N}^*$, these basis functions are rational with a numerator and a denominator are polynomials of degree k. All of the classical properties as positivity, partition of unity are hold for these rational Bernstein basis and they constitute approximation basis functions for continuous functions spaces. The B\'ezier curves obtained verify the classical properties and we have the classical computational algorithms like the deCasteljau Algorithm and the algorithm of subdivision with the similar accuracy. Given a degree k and a control polygon points all of these algorithms converge to the same B\'ezier curve as the classical case. That means the B\'ezier curve is independent of the index $\alpha$. The classical polynomial Bernstein basis seems a asymptotic case of our new class of rational Bernstein basis.


title: Characterization of half-radial matrices

abstract: Numerical radius $r(A)$ is the radius of the smallest ball with the center at zero containing the field of values of a given square matrix $A$. It is well known that $r(A)\leq \|A\| \leq 2r(A)$, where $\| \cdot \|$ is the matrix 2-norm. Matrices attaining the lower bound are called radial, and have been analyzed thoroughly. This is not the case for matrices attaining the upper bound where only partial results are available. In this paper we consider matrices satisfying $r(A)=\|A\|/2$, and call them half-radial. We summarize the existing results and formulate new ones. In particular, we investigate their singular value decomposition and algebraic structure, and provide other necessary and sufficient conditions for a matrix to be half-radial. Based on that, we study the extreme case of the attainable constant~$2$ in Crouzeix's conjecture. The presented results support the conjecture of Greenbaum and Overton, that the Crabb-Choi-Crouzeix matrix always plays an important role in this extreme case.


title: Random constructions for translates of non-negative functions

abstract: Suppose $\Lambda$ is a discrete infinite set of nonnegative real numbers. We say that $ {\Lambda}$ is type $2$ if the series $s(x)=\sum_{\lambda\in\Lambda}f(x+\lambda)$ does not satisfy a zero-one law. This means that we can find a non-negative measurable "witness function" $f: {\mathbb R}\to [0,+ {\infty})$ such that both the convergence set $C(f, {\Lambda})=\{x: s(x)<+ {\infty} \}$ and its complement the divergence set $D(f, {\Lambda})=\{x: s(x)=+ {\infty} \}$ are of positive Lebesgue measure. If $ {\Lambda}$ is not type $2$ we say that $ {\Lambda}$ is type $1$. The main result of our paper answers a question raised by Z. Buczolich, J-P. Kahane, and D. Mauldin. By a random construction we show that one can always choose a witness function which is the characteristic function of a measurable set. We also consider the effect on the type of a set $ {\Lambda}$ if we randomly delete its elements. Motivated by results concerning weighted sums $\sum c_n f(nx)$ and the Khinchin conjecture, we also discuss some results about weighted sums $\sum_{n=1}^{\infty}c_n f(x+\lambda_n)$.


title: A construction of free digroup

abstract: We give a construction of a free digroup $F(X)$ on a set $X$ and formulate the halo and the group parts of $F(X)$. We prove that $F(X)$ is isomorphic to $F(Y)$ if and only if $card(X)=card(Y)$.


title: Two-Dimensional Computation of Pulsed Magnetic Field Diffusion Dynamics in Gold Cone with Consideration of Inductive Heating and Temperature Dependence of Electrical Conductivity

abstract: Application of an external kilo-tesla-level magnetic field, which can be generated using high-intensity laser, to a target is a promising scheme to reduce spray angle of a laser-driven relativistic electron beam (REB) for enhancing the isochoric heating of a dense plasma with the laser-driven REB. Here we have developed a two-dimensional electro-magnetic dynamics (2D-EMD) simulation code to solve Maxwell equations with considerations of the inductive heating and temperature-dependence of electrical conductivity of a material for calculating temporally and spatially resolved two-dimensional profile of the externally applied magnetic field in a gold-cone-attached target.


title: Automation Of Road Intersections Using Consensus-based Auction Algorithms

abstract: This paper investigates a consensus-based auction algorithm in the context of decentralized traffic control. In particular, we study the automation of a road intersection, where a set of vehicles is required to cross without collisions. The crossing order will be negotiated in a decentralized fashion. An on-board model predictive controller (MPC) will compute an optimal trajectory which avoids collisions with higher priority vehicles, thus retaining convex safety constraints. Simulations are then performed in a time-variant traffic environment.