title: Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices

abstract: We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The technique combines time-stepping in canonical coordinates with quadrature in non-canonical coordinates and is applicable in systems where a global transformation to canonical coordinates is known but its inverse is not. A fully implicit class of symplectic Runge-Kutta schemes has recently been introduced and applied to guiding-center motion by [Zhang et al., Phys. Plasmas 21, 32504 (2014)]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance their performance, in particular avoiding evaluation of non-canonical variables at full time steps. For application in toroidal plasma confinement configurations with nested magnetic flux surfaces a global canonicalization of coordinates for the guiding-center Lagrangian by a spatial transform is presented that allows for pre-computation of the required map in a parallel algorithm in the case of time-independent magnetic field geometry. Guiding-center orbits are studied in stationary magnetic equilibrium fields of an axisymmetric tokamak and a realistic three-dimensional stellarator configuration. Superior long-term properties of symplectic methods are demonstrated in comparison to a conventional adaptive Runge-Kutta scheme. Finally statistics of fast fusion alpha particle losses over their slowing-down time are computed in the stellarator field on a representative sample, reaching a speed-up of the symplectic Euler scheme by more than a factor three compared to usual Runge-Kutta schemes while keeping the same statistical accuracy and linear scaling with the number of computing threads.


title: Minimizing Age of Information in Cognitive Radio-based IoT Systems: Underlay or Overlay?

abstract: We consider a cognitive radio-based Internet-of-Things (CR-IoT) network consisting of one primary IoT (PIoT) system and one secondary IoT (SIoT) system. The IoT devices of both the PIoT and the SIoT respectively monitor one physical process and send randomly generated status updates to their associated access points (APs). The timeliness of the status updates is important as the systems are interested in the latest condition (e.g., temperature, speed and position) of the IoT device. In this context, two natural questions arise: (1) How to characterize the timeliness of the status updates in CR-IoT systems? (2) Which scheme, overlay or underlay, is better in terms of the timeliness of the status updates. To answer these two questions, we adopt a new performance metric, named the age of information (AoI). We analyze the average peak AoI of the PIoT and the SIoT for overlay and underlay schemes, respectively. Simple asymptotic expressions of the average peak AoI are also derived when the PIoT operates at high signal-to-noise ratio (SNR). Based on the asymptotic expressions, we characterize a critical generation rate of the PIoT system, which can determine the superiority of overlay and underlay schemes in terms of the average peak AoI of the SIoT. Numerical results validate the theoretical analysis and uncover that the overlay and underlay schemes can outperform each other in terms of the average peak AoI of the SIoT for different system setups.


title: Rodier type theorem for generalized principal series

abstract: Given a regular supercuspidal representation $\rho$ of the Levi subgroup $M$ of a standard parabolic subgroup $P=MN$ in a connected reductive group $G$ defined over a non-archimedean local field $F$, we serve you a Rodier type structure theorem which provides us a geometrical parametrization of the set $JH(Ind^G_P(\rho))$ of Jordan--H{\"o}lder constituents of the Harish-Chandra parabolic induction representation $Ind^G_P(\rho)$, vastly generalizing Rodier structure theorem for $P=B=TU$ Borel subgroup of a connected split reductive group about 40 years ago. Our novel contribution is to overcome the essential difficulty that the relative Weyl group $W_M=N_G(M)/M$ is not a coxeter group in general, as opposed to the well-known fact that the Weyl group $W_T=N_G(T)/T$ is a coxeter group. Indeed, such a beautiful structure theorem also holds for finite central covering groups.


title: Rethinking Uplink Hybrid Processing: When is Pure Analog Processing Suggested?

abstract: In this correspondence, we analytically characterize the benefit of digital processing in uplink massive multiple-input multiple-output (MIMO) with sub-connected hybrid architecture. By assuming that the number of radio frequency (RF) chains is equal to that of users, we characterize achievable rates of both pure analog detection and hybrid detection under the i.i.d. Rayleigh fading channel model. From the derived expressions, we discover that the analog processing can outperform the hybrid processing using the maximal ratio combining (MRC) or zero-forcing (ZF) criterion in cases under some engineering assumptions. Performance comparison of the schemes are presented under tests with various numbers of users and numbers of antennas at the base station.


title: MultiK: A Framework for Orchestrating Multiple Specialized Kernels

abstract: We present, MultiK, a Linux-based framework 1 that reduces the attack surface for operating system kernels by reducing code bloat. MultiK "orchestrates" multiple kernels that are specialized for individual applications in a transparent manner. This framework is flexible to accommodate different kernel code reduction techniques and, most importantly, run the specialized kernels with near-zero additional runtime overheads. MultiK avoids the overheads of virtualization and runs natively on the system. For instance, an Apache instance is shown to run on a kernel that has (a) 93.68% of its code reduced, (b) 19 of 23 known kernel vulnerabilities eliminated and (c) with negligible performance overheads (0.19%). MultiK is a framework that can integrate with existing code reduction and OS security techniques. We demonstrate this by using D-KUT and S-KUT -- two methods to profile and eliminate unwanted kernel code. The whole process is transparent to the user applications because MultiK does not require a recompilation of the application.


title: [Plasma 2020 Decadal] The essential role of multi-point measurements in turbulence investigations: the solar wind beyond single scale and beyond the Taylor Hypothesis

abstract: This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.


title: Description of infinite orbits on multiple projective spaces

abstract: Let $G$ be the general linear group of the degree $n\geq 2$ over the field $\mathbb{K}=\mathbb{R}$ or $\mathbb{C}$. In this article, we give a description of orbit decomposition of the multiple projective space $G^m/P^m$ under the diagonal action of $G$ where $P$ is the maximal parabolic subgroup of $G$ such that $G/P\cong\mathbb{P}^{n-1}\mathbb{K}$. We also construct representatives of orbits. If $m\geq 4$, the number of orbits is infinite, and we give a description of those uncountably many orbits.


title: Critical phenomena and nonlinear dynamics in a spin ensemble strongly coupled to a cavity. I. Semiclassical approach

abstract: We present a theoretical study on the nonlinear dynamics and stationary states of an inhomogeneously broadened spin ensemble coupled to a single-mode cavity driven by an external drive with constant amplitude. Assuming a sizeable number of constituents within the ensemble allows us to use a semiclassical approach and to formally reduce the theoretical description to the Maxwell-Bloch equations for the cavity and spin amplitudes. We explore the critical slowing-down effect, quench dynamics, and asymptotic behavior of the system near a steady-state dissipative phase transition accompanied by a bistability effect. Some of our theoretical findings have recently been successfully verified in a specific experimental realization based on a spin ensemble of negatively charged nitrogen-vacancy centers in diamond strongly coupled to a single-mode microwave cavity (see Science Adv. 3, e1701626 (2017)).

categories: quant-ph cond-mat.mes-hall nlin.CD physics.atom-ph


title: Critical phenomena and nonlinear dynamics in a spin ensemble strongly coupled to a cavity. II. Semiclassical-to-quantum boundary

abstract: We numerically study the dynamics and stationary states of a spin ensemble strongly coupled to a single-mode resonator subjected to loss and external driving. Employing a generalized cumulant expansion approach we analyze finite-size corrections to a semiclassical description of amplitude bistability, which is a paradigm example of a driven-dissipative phase transition. Our theoretical model allows us to include inhomogeneous broadening of the spin ensemble and to capture in which way the quantum corrections approach the semiclassical limit for increasing ensemble size $N$. We set up a criterion for the validity of the Maxwell-Bloch equations and show that close to the critical point of amplitude bistability even very large spin ensembles consisting of up to $10^4$ spins feature significant deviations from the semiclassical theory.


title: Theoretical realization of two-dimensional M3(C6X6)2(M= Co, Cr, Cu, Fe, Mn, Ni, Pd, Rh and X= O, S, Se) metal-organic frameworks

abstract: Most recently, Cu-hexahydroxybenzene MOF was for the time experimentally realized, through a kinetically controlled approach. Cu-HHB belongs to the family of conductive MOFs with a chemical formula of M3(C6X6)2(X=NH, O, S). Motivated by the recent experimental advance in the fabrication of Cu-HHB, we conducted extensive first-principles simulations to explore the thermal stability, mechanical properties and electronic characteristics of M3(C6X6)2(M= Co, Cr, Cu, Fe, Mn, Ni, Pd, Rh and X= O, S, Se) monolayers. First-principles results confirm that all considered 2D porous lattices are thermally stable at high temperatures over 1500 K. It was moreover found that these novel 2D structures can exhibit linear elasticity with considerable tensile strengths, revealing their suitability for practical applications in nanodevices.Depending on the metal and chalcogen atoms in M3(C6X6)2 monolayers, they can yield various electronic and magnetic properties, such as; magnetic semiconducting, perfect half metallic, magnetic and nonmagnetic metallic behaviours. This work highlights the outstanding physics of M3(C6X6)2 2D porous lattices and will hopefully help to expand this conductive MOF family, as promising candidates to design advanced energy storage/conversion, electronics and spintronics systems.


title: Recover and RELAX: Concern-Oriented Software Architecture Recovery for Systems Development and Maintenance

abstract: The stakeholders of a system are legitimately interested in whether and how its architecture reflects their respective concerns at each point of its development and maintenance processes. Having such knowledge available at all times would enable them to continually adjust their systems structure at each juncture and reduce the buildup of technical debt that can be hard to reduce once it has persisted over many iterations. Unfortunately, software systems often lack reliable and current documentation about their architecture. In order to remedy this situation, researchers have conceived a number of architectural recovery methods, some of them concern-oriented. However, the design choices forming the bases of most existing recovery methods make it so none of them have a complete set of desirable qualities for the purpose stated above. Tailoring a recovery to a system is either not possible or only through iterative experiments with numeric parameters. Furthermore, limitations in their scalability make it prohibitive to apply the existing techniques to large systems. Finally, since several current recovery methods employ non-deterministic sampling, their inconsistent results do not lend themselves well to tracking a systems course over several versions, as needed by its stakeholders. RELAX (RELiable Architecture EXtraction), a new concern-based recovery method that uses text classification, addresses these issues efficiently by (1) assembling the overall recovery result from smaller, independent parts, (2) basing it on an algorithm with linear time complexity and (3) being tailorable to the recovery of a single system or a sequence thereof through the selection of meaningfully named, semantic topics. An intuitive, informative architectural visualization rounds out RELAX's contributions. RELAX is illustrated on a number of existing open-source systems and compared to other recovery methods.


title: $S$ =1/2 ferromagnetic Heisenberg chain in a verdazyl-based complex

abstract: We present a model compound for the $S$=1/2 ferromagnetic Heisenberg chain composed of the verdazyl-based complex $[$Zn(hfac)$_2]$$[$4-Cl-$o$-Py-V-(4-F)$_2]$. $Ab$ $initio$ MO calculations indicate a predominant ferromagnetic interaction forming an $S$=1/2 ferromagnetic chain. The magnetic susceptibility and specific heat indicate a phase transition to an AF order owing to the finite interchain couplings. We explain the magnetic susceptibility and magnetization curve above the phase transition temperature based on the $S$=1/2 ferromagnetic Heisenberg chain. The magnetization curve in the ordered phase is described by a conventional AF two-sublattice model. Furthermore, the obtained magnetic specific heat reproduces the almost temperature-independent behavior of the $S$=1/2 ferromagnetic Heisenberg chain. In the low-temperature region, the magnetic specific heat exhibits $\sqrt{T}$ dependence, which is attributed to the energy dispersion in the ferromagnetic chain.


title: Wavelength Dependent Tunneling Delay Time

abstract: A clear consensus on how long it takes a particle to tunnel through a potential barrier has never been so urgently required, since the electron dynamics in strong-field ionization can be resolved on attosecond time-scale in experiment and the exact nature of the tunneling process is the key to trigger subsequent attosecond techniques. Here a general picture of tunneling time is suggested by introducing a quantum travel time, which is defined as the ratio of the travel distance to the expected value of the velocity operator under the barrier. Specially, if applied to rectangular barrier tunneling, it can retrieve the B\"{u}ttiker-Landauer time $\tau_{BL}$ in the case of an opaque barrier, and has a clear physical meaning in the case of a very thin barrier wherein $\tau_{BL}$ can not be well defined. In the case of strong-field tunneling process, with the help of the newly defined time, the tunneling delay time measured by attoclock experiment can be interpreted as a travel time spent by the electron to tunnel from a point under the barrier to the tunnel exit. In addition, a peculiar oscillation structure in the wavelength dependence of tunneling delay time in deep tunneling regime is observed, which is beyond the scope of adiabatic tunneling picture. This oscillation structure can be attributed to the interference between the ground state tunneling channel and the excited states tunneling channels.


title: On-Line Balancing of Random Inputs

abstract: We consider an online vector balancing game where vectors $v_t$, chosen uniformly at random in $\{-1,+1\}^n$, arrive over time and a sign $x_t \in \{-1,+1\}$ must be picked immediately upon the arrival of $v_t$. The goal is to minimize the $L^\infty$ norm of the signed sum $\sum_t x_t v_t$. We give an online strategy for picking the signs $x_t$ that has value $O(n^{1/2})$ with high probability. Up to constants, this is the best possible even when the vectors are given in advance.


title: Magnetically-Sensitive Valley Polarization Reversal and Revival of Defect-Localized Excitons in WSe2-WS2

abstract: Manipulating and reserving the valley pseudospin of excitons is one core aim in the two-dimensional transition metal dichalcogenides (TMDs). However, due to the strong electron-hole exchange and spin-orbit coupling interactions, the exciton recombination lifetime is subject to picosecond timescale intrinsically, and the valley polarization is hardly modulated by a moderate magnetic field. It is fortunate that interlayer and defect-localized excitons promise to overcome these difficulties by suppressing the above interactions. Here we clearly reveal that the valley polarization can be reversed and revived in the defect-localized excitons with a microsecond lifetime in AB-stacked WSe2-WS2 heterobilayer. Specifically, for the interlayer defect-localized exciton, the valley polarization is reversed and can be efficiently enhanced by a weak out-of-plane magnetic field (<0.4 T). In sharp contrast, the valley polarization of the intralayer defect-localized exciton can revive after a fast decay process and follows the direction of the moderate out-of-plane magnetic field (<3 T). We explain the reversed valley polarization with highly magnetic sensitivity by the delocalization of defect-localized holes under a weak magnetic field and the revival of valley polarization by the valley Zeeman effect under a moderate magnetic field. Our results demonstrate that the valley pseudospin of defect-localized excitons can be efficiently modulated by the external magnetic field and enrich both the understanding and the technical approaches on manipulating the valley dynamics in TMDs and their heterostructure.


title: Multi-topological semantics for intuitionistic modal logic

abstract: We present three examples of \textit{multi-topological} semantics for intuitionistic modal logic with one modal operator $\Box$ (which behaves in some sense like necessity). We show that it is possible to treat neighborhood models, introduced earlier, as multi-topological. From the neighborhood point of view, our method is based on differences between properties of minimal and maximal neighborhoods. Also we propose transformation of multi-topological spaces into the neighborhood structures. Although our neighborhoods can be considered as bi-relational models, we believe that reasoning in terms of neighborhoods gives better topological intuitions.


title: Emotion Action Detection and Emotion Inference: the Task and Dataset

abstract: Many Natural Language Processing works on emotion analysis only focus on simple emotion classification without exploring the potentials of putting emotion into "event context", and ignore the analysis of emotion-related events. One main reason is the lack of this kind of corpus. Here we present Cause-Emotion-Action Corpus, which manually annotates not only emotion, but also cause events and action events. We propose two new tasks based on the data-set: emotion causality and emotion inference. The first task is to extract a triple (cause, emotion, action). The second task is to infer the probable emotion. We are currently releasing the data-set with 10,603 samples and 15,892 events, basic statistic analysis and baseline on both emotion causality and emotion inference tasks. Baseline performance demonstrates that there is much room for both tasks to be improved.


title: A Deep Look into Neural Ranking Models for Information Retrieval

abstract: Ranking models lie at the heart of research on information retrieval (IR). During the past decades, different techniques have been proposed for constructing ranking models, from traditional heuristic methods, probabilistic methods, to modern machine learning methods. Recently, with the advance of deep learning technology, we have witnessed a growing body of work in applying shallow or deep neural networks to the ranking problem in IR, referred to as neural ranking models in this paper. The power of neural ranking models lies in the ability to learn from the raw text inputs for the ranking problem to avoid many limitations of hand-crafted features. Neural networks have sufficient capacity to model complicated tasks, which is needed to handle the complexity of relevance estimation in ranking. Since there have been a large variety of neural ranking models proposed, we believe it is the right time to summarize the current status, learn from existing methodologies, and gain some insights for future development. In contrast to existing reviews, in this survey, we will take a deep look into the neural ranking models from different dimensions to analyze their underlying assumptions, major design principles, and learning strategies. We compare these models through benchmark tasks to obtain a comprehensive empirical understanding of the existing techniques. We will also discuss what is missing in the current literature and what are the promising and desired future directions.


title: Modelling of shallow aquifers in interaction with overland water

abstract: In this work, we present a class of new efficient models for water flow in shallow unconfined aquifers, giving an alternative to the classical but less tractable 3D-Richards model. Its derivation is guided by two ambitions: any new model should be low cost in computational time and should still give relevant results at every time scale. We thus keep track of two types of flow occurring in such a context and which are dominant when the ratio thickness over longitudinal length is small: the first one is dominant in a small time scale and is described by a vertical 1D-Richards problem; the second one corresponds to a large time scale, when the evolution of the hydraulic head turns to become independent of the vertical variable. These two types of flow are appropriately modelled by, respectively, a one-dimensional and a two-dimensional system of PDEs boundary value problems. They are coupled along an artificial level below which the Dupuit hypothesis holds true (i.e. the vertical flow is instantaneous) in a way ensuring that the global model is mass conservative. Tuning the artificial level, which even can depend on an unknown of the problem, we browse the new class of models. We prove using asymptotic expansions that the 3D-Richards problem and each model of the class behaves the same at every considered time scale (short, intermediate and large) in thin aquifers. The results are illustrated by numerical simulations, showing especially that the new models results fit well with the ones obtained with the original 3d-Richards problem even in non-thin aquifers.


title: k-Means Clustering of Lines for Big Data

abstract: The input to the $k$-median for lines problem is a set $L$ of $n$ lines in $\mathbb{R}^d$, and the goal is to compute a set of $k$ centers (points) in $\mathbb{R}^d$ that minimizes the sum of squared distances over every line in $L$ and its nearest center. This is a straightforward generalization of the $k$-median problem where the input is a set of $n$ points instead of lines. We suggest the first PTAS that computes a $(1+\epsilon)$-approximation to this problem in time $O(n \log n)$ for any constant approximation error $\epsilon \in (0, 1)$, and constant integers $k, d \geq 1$. This is by proving that there is always a weighted subset (called coreset) of $dk^{O(k)}\log (n)/\epsilon^2$ lines in $L$ that approximates the sum of squared distances from $L$ to any given set of $k$ points. Using traditional merge-and-reduce technique, this coreset implies results for a streaming set (possibly infinite) of lines to $M$ machines in one pass (e.g. cloud) using memory, update time and communication that is near-logarithmic in $n$, as well as deletion of any line but using linear space. These results generalized for other distance functions such as $k$-median (sum of distances) or ignoring farthest $m$ lines from the given centers to handle outliers. Experimental results on 10 machines on Amazon EC2 cloud show that the algorithm performs well in practice. Open source code for all the algorithms and experiments is also provided. This thesis is an extension of the following accepted paper: "$k$-Means Clustering of Lines for Big Data", by Yair Marom & Dan Feldman, Proceedings of NeurIPS 2019 conference, to appear on December 2019.


title: Quantum sensing of curvature

abstract: We address the problem of sensing the curvature of a manifold by performing measurements on a particle constrained to the manifold itself. In particular, we consider situations where the dynamics of the particle is quantum mechanical and the manifold is a surface embedded in the three-dimensional Euclidean space. We exploit ideas and tools from quantum estimation theory to quantify the amount of information encoded into a state of the particle, and to seek for optimal probing schemes, able to actually extract this information. Explicit results are found for a free probing particle and in the presence of a magnetic field. We also address precision achievable by position measurement, and show that it provides a nearly optimal detection scheme, at least to estimate the radius of a sphere or a cylinder.


title: Microdroplet self-propulsion during dropwise condensation on lubricant-infused surfaces

abstract: Water vapor condensation is common in nature and widely used in industrial applications, including water harvesting, power generation, and desalination. As compared to traditional filmwise condensation, dropwise condensation on lubricant-infused surfaces (LIS) can lead to an order-of-magnitude increase in heat transfer rates. Small droplets (with the diameter below 100 $\mu$m) account for nearly 85 percent of the total heat transfer and droplet sweeping plays a crucial role in clearing nucleation sites, allowing for frequent re-nucleation. Here, we focus on the dynamic interplay of microdroplets with the thin lubricant film during water vapor condensation on LIS. Coupling high-speed imaging, optical microscopy, and interferometry, we show that the initially uniform lubricant film re-distributes during condensation. Governed by lubricant height gradients, microdroplets as small as 2 $\mu$m in diameter undergo rigorous and gravity-independent self-propulsion, travelling distances multiples of their diameters at velocities up to 1100 $\mu$m/s. Although macroscopically the movement appears to be random, we show that on a microscopic level capillary attraction due to asymmetrical lubricant menisci causes this gravity-independent droplet motion. Based on a lateral force balance analysis, we quantitatively find that the sliding velocity initially increases during movement, but decreases sharply at shorter inter-droplet spacing. The maximum sliding velocity is inversely proportional to the oil viscosity and is strongly dependent of the droplet size, which is in excellent agreement with the experimental observations. This novel and non-traditional droplet movement is expected to significantly enhance the sweeping efficiency during dropwise condensation, leading to higher nucleation and heat transfer rates.


title: Charges of the gravitational field and (3+1) decomposition of CYK tensors part 2

abstract: The work describes the method of construction of charges (conserved quantities) for the gravity field in the (3 + 1) decomposition. The presented construction uses tensors of the electrical and magnetic parts of the Weyl tensor and conformal Killing vectors. In the case of conformally flat spatial hypersurfaces, we get twenty local charges, which can be expressed in terms of the initial data (three-dimensional metric and extrinsic curvature tensor). The work shows the relationships between charges which are obtained by this method and the usual ADM approach. In traditional ADM approach Killing vectors are used to construct corresponding charges e.g. time translation corresponds to ADM mass, spatial translations give linear momentum and rotations correspond to angular momentum. Gravito-electric and gravito-magnetic charges need conformal Killing vectors e.g. mass corresponds to dilation, linear momentum is related to rotation and angular momentum needs conformal acceleration. The analyzed example of the Schwarzschild--de Sitter spacetime suggests that in some cases the mass calculated by the presented method has better properties than the traditional ADM mass. Next we discuss asymptotic charges which are no longer rigidly conserved but rather approach finite value at spatial infinity.


title: Non-intrusive speech quality assessment using neural networks

abstract: Estimating the perceived quality of an audio signal is critical for many multimedia and audio processing systems. Providers strive to offer optimal and reliable services in order to increase the user quality of experience (QoE). In this work, we present an investigation of the applicability of neural networks for non-intrusive audio quality assessment. We propose three neural network-based approaches for mean opinion score (MOS) estimation. We compare our results to three instrumental measures: the perceptual evaluation of speech quality (PESQ), the ITU-T Recommendation P.563, and the speech-to-reverberation energy ratio. Our evaluation uses a speech dataset contaminated with convolutive and additive noise, labeled using a crowd-based QoE evaluation, evaluated with Pearson correlation with MOS labels, and mean-squared-error of the estimated MOS. Our proposed approaches outperform the aforementioned instrumental measures, with a fully connected deep neural network using Mel-frequency features providing the best correlation (0.87) and the lowest mean squared error (0.15)


title: Classification of dry age-related macular degeneration and diabetic macular edema from optical coherence tomography images using dictionary learning

abstract: Age-related Macular Degeneration (AMD) and Diabetic Macular Edema (DME) are the major causes of vision loss in developed countries. Alteration of retinal layer structure and appearance of exudate are the most significant signs of these diseases. With the aim of automatic classification of DME, AMD and normal subjects from Optical Coherence Tomography (OCT) images, we proposed a classification algorithm. The two important issues intended in this approach are, not utilizing retinal layer segmentation which by itself is a challenging task and attempting to identify diseases in their early stages, where the signs of diseases appear in a small fraction of B-Scans. We used a histogram of oriented gradients (HOG) feature descriptor to well characterize the distribution of local intensity gradients and edge directions. In order to capture the structure of extracted features, we employed different dictionary learning-based classifiers. Our dataset consists of 45 subjects: 15 patients with AMD, 15 patients with DME and 15 normal subjects. The proposed classifier leads to an accuracy of 95.13%, 100.00%, and 100.00% for DME, AMD, and normal OCT images, respectively, only by considering the 4% of all B-Scans of a volume which outperforms the state of the art methods.


title: On the Dynamics of Comets in Extrasolar Planetary Systems

abstract: Since very recently, we acquired knowledge on the existence of comets in extrasolar planetary systems. The formation of comets together with planets around host stars now seems evident. As stars are often born in clusters of interstellar clouds, the interaction between the systems will lead to the exchange of material at the edge of the clouds. Therefore, almost every planetary system should have leftover remnants as a result of planetary formation in form of comets at the edges of those systems. These Oort clouds around stars are often disturbed by different processes (e.g., galactic tides, passing stars, etc.), which consequently scatter bodies from the distant clouds into the system close to the host star. Regarding the Solar System, we observe this outcome in the form of cometary families. This knowledge supports the assumption of the existence of comets around other stars. In the present work, we study the orbital dynamics of hypothetical exocomets, based on detailed computer simulations, in three star-planet systems, which are: HD~10180, 47~UMa, and HD~141399. These systems host one or more Jupiter-like planets, which change the orbits of the incoming comets in characteristic ways.


title: Adaptive image processing: a bilevel structure learning approach for mixed-order total variation regularizers

abstract: A class of mixed-order \emph{PDE}-constraint regularizer for image processing problem is proposed, generalizing the standard first order total variation $(TV)$. A semi-supervised (bilevel) training scheme, which provides a simultaneous optimization with respect to parameters and the new class of regularizers, is studied. Also, A finite approximation method, which used to solve the global optimization solutions of such training scheme, is introduced and analyzed.


title: Semimartingale theory of monotone mean--variance portfolio allocation

abstract: We study dynamic optimal portfolio allocation for monotone mean--variance preferences in a general semimartingale model. Armed with new results in this area we revisit the work of Cui, Li, Wang and Zhu (2012, MAFI) and fully characterize the circumstances under which one can set aside a non-negative cash flow while simultaneously improving the mean--variance efficiency of the left-over wealth. The paper analyzes, for the first time, the monotone hull of the Sharpe ratio and highlights its relevance to the problem at hand.


title: Dzyaloshinskii-Moriya interaction and spin-orbit torque at the Ir/Co interface

abstract: We studied the spin torque efficiency and the Dzyaloshinskii-Moriya interaction (DMI) of heterostructures that contain interface(s) of Ir and Co. The current-induced shifts of the anomalous Hall loops were used to determine the spin torque efficiency and DMI of [Pt/Co/X] multilayers (X=Ir, Cu) as well as Ir/Co and Pt/Ir/Co reference films. We find the effective spin Hall angle and the spin diffusion length of Ir to be ~0.01 and less than ~1 nm, respectively. The short spin diffusion length and the high conductivity make Ir an efficient spin sink layer. Such spin sink layer can be used to control the flow of spin current in heterostructures and to induce sufficient spin-orbit torque on the magnetic layer. The DMI of Ir and Co interface is found to be in the range of ~1.4 to ~2.2 mJ/m2, similar in magnitude to that of the Pt and Co interface. The Ir/Co and Pt/Co interfaces possess the same sign of DMI, resulting in a reduced DMI for the [Pt/Co/Ir] multilayers compared to that of the [Pt/Co/Cu] multilayers. These results show the unique role the Ir layer plays in defining spin-orbit torque and chiral magnetism in thin film heterostructures.


title: An Accurate Bilinear Cavern Model for Compressed Air Energy Storage

abstract: Compressed air energy storage is suitable for large-scale electrical energy storage, which is important for integrating renewable energy sources into electric power systems. A typical compressed air energy storage plant consists of compressors, expanders, caverns, and a motor/generator set. Current cavern models used for compressed air energy storage are either accurate but highly nonlinear or linear but inaccurate. The application of highly nonlinear cavern models in power system optimization problems renders them computationally challenging to solve. In this regard, an accurate bilinear cavern model for compressed air energy storage is proposed in this paper. The charging and discharging processes in a cavern are divided into several real/virtual states. The first law of thermodynamics and ideal gas law are then utilized to derive a cavern model, i.e., a model for the variation of temperature and pressure in these processes. Thereafter, the heat transfer between the air in the cavern and the cavern wall is considered and integrated into the cavern model. By subsequently eliminating several negligible terms, the cavern model reduces to a bilinear model. The accuracy of the bilinear cavern model is verified via comparison with both an accurate nonlinear model and two sets of field-measured data. The bilinear cavern model can be easily linearized and is then suitable for integration into optimization problems considering compressed air energy storage. This is verified via comparatively solving a self-scheduling problem of compressed air energy storage using different cavern models.


title: The isomorphism problem for universal enveloping algebras of four-dimensional solvable Lie algebras

abstract: This paper is a contribution to the isomorphism problem for universal enveloping algebras of finite-dimensional Lie algebras. We focus on solvable Lie algebras of small dimensions over fields of arbitrary characteristic. We prove, over an arbitrary field, that the isomorphism type of a metabelian Lie algebra whose derived subalgebra has codimension one is determined by its universal enveloping algebra. As an application of the results in this paper, we solve the isomorphism problem for solvable Lie algebras of dimension four over fields of characteristic zero and also point out the problems that occur in prime characteristic.


title: A Cross-Season Correspondence Dataset for Robust Semantic Segmentation

abstract: In this paper, we present a method to utilize 2D-2D point matches between images taken during different image conditions to train a convolutional neural network for semantic segmentation. Enforcing label consistency across the matches makes the final segmentation algorithm robust to seasonal changes. We describe how these 2D-2D matches can be generated with little human interaction by geometrically matching points from 3D models built from images. Two cross-season correspondence datasets are created providing 2D-2D matches across seasonal changes as well as from day to night. The datasets are made publicly available to facilitate further research. We show that adding the correspondences as extra supervision during training improves the segmentation performance of the convolutional neural network, making it more robust to seasonal changes and weather conditions.


title: Molecular Polar Belief Propagation Decoder and Successive Cancellation Decoder

abstract: By constructing chemical reaction networks (CRNs), this paper proposes a method of synthesizing polar decoder using belief propagation (BP) algorithm and successive cancellation (SC) algorithm, respectively. Theoretical analysis and simulation results have validated the feasibility of the method. Reactions in the proposed design could be experimentally implemented with DNA strand displacement reactions, making the proposed polar decoders promising for wide application in nanoscale devices.


title: Master equation for operational amplifiers: stability of negative differential converters, crossover frequency and pass-bandwidth

abstract: The time dependent master equation from the seminal article by Ragazzini, Randall and Russell [J. R. Ragazzini, R. H. Randall and F. A. Russell, "Analysis of Problems in Dynamics by Electronic Circuits", Proc. of the I.R.E., Vol. 35, pp. 444--452, (1947)] is recovered as necessary tool for the analysis of contemporary circuits with operational amplifiers. This equation gives the relation between time dependent the output voltage $U_0(t)$ and the difference between the input voltages ($U_{+}(t)$ and $U_{-}(t)$). The crossover frequency $f_0$ is represented the time constant $\tau_{_0}$ of this equation. The work of the master equation is illustrated by two typical examples: a) the stability criterion of the devices with negative impedance converters, which we consider as a new result b) the frequency dependence of the amplifiers with operational amplifiers given in the technical specifications without citations of time dependent equation. A simple circuit for determination of $f_0$ is suggested and the method is illustrated by determination of crossover frequency for the low-noise and high speed ADA4898 operational amplifier. It is concluded that for an exact calculation of the pass bandwidth of amplifiers with active filters the 70 years old master equation is a useful technique implicitly included in the contemporary software. The frequency dependent formulae for the amplification coefficient of inverting and non-inverting amplifiers are given for the case of non-zero conductivity between the inputs of the operational amplifiers.


title: An unstructured finite element model for incompressible two-phase flow based on a monolithic conservative level set method

abstract: We present a robust numerical method for solving incompressible, immiscible two-phase flows. The method extends the monolithic phase conservative level set method with embedded redistancing by Quezada de Luna et al. [38] and a semi-implicit high-order projection scheme for variable-density flows by Guermond and Salgado [17]. The level set method can be initialized conveniently via a simple phase indicator field, which is pre-processed to obtain an approximate signed distance function. To do this, we propose a new PDE-based redistancing method. We also improve the scheme in [38] to provide more accuracy and robustness in full two-phase flow simulations. Specifically, we perform an extra step to ensure convergence to the signed distance level set function and simplify other aspects of the original scheme. Lastly, we introduce consistent artificial viscosity to stabilize the momentum equations in the context of the projection scheme. This stabilization is algebraic, has no tunable parameters and is suitable for unstructured meshes and arbitrary refinement levels. The overall methodology includes few numerical tuning parameters; however, for the wide range of problems that we solve, we identify only one parameter that strongly affects performance of the computational model and provide a value that provides accurate results across all the benchmarks presented. The result is a robust, accurate, and efficient two-phase flow model, which is mass- and volume-conserving on unstructured meshes and has low user input requirements for real applications.


title: Robust Super-Resolution GAN, with Manifold-based and Perception Loss

abstract: Super-resolution using deep neural networks typically relies on highly curated training sets that are often unavailable in clinical deployment scenarios. Using loss functions that assume Gaussian-distributed residuals makes the learning sensitive to corruptions in clinical training sets. We propose novel loss functions that are robust to corruptions in training sets by modeling heavy-tailed non-Gaussian distributions on the residuals. We propose a loss based on an autoencoder-based manifold-distance between the super-resolved and high-resolution images, to reproduce realistic textural content in super-resolved images. We propose to learn to super-resolve images to match human perceptions of structure, luminance, and contrast. Results on a large clinical dataset shows the advantages of each of our contributions, where our framework improves over the state of the art.


title: A new Capacity-Achieving Private Information Retrieval Scheme with (Almost) Optimal File Length for Coded Servers

abstract: In a distributed storage system, private information retrieval (PIR) guarantees that a user retrieves one file from the system without revealing any information about the identity of its interested file to any individual server. In this paper, we investigate $(N,K,M)$ coded sever model of PIR, where each of $M$ files is distributed to the $N$ servers in the form of $(N,K)$ maximum distance separable (MDS) code for some $N>K$ and $M>1$. As a result, we propose a new capacity-achieving $(N,K,M)$ coded linear PIR scheme such that it can be implemented with file length $\frac{K(N-K)}{\gcd(N,K)}$, which is much smaller than the previous best result $K\big(\frac{N}{\gcd(N,K)}\big)^{M-1}$. Notably, among all the capacity-achieving coded linear PIR schemes, we show that the file length is optimal if $M>\big\lfloor \frac{K}{\gcd(N,K)}-\frac{K}{N-K}\big\rfloor+1$, and within a multiplicative gap $\frac{K}{\gcd(N,K)}$ of a lower bound on the minimum file length otherwise.


title: Real time backbone for semantic segmentation

abstract: The rapid development of autonomous driving in recent years presents lots of challenges for scene understanding. As an essential step towards scene understanding, semantic segmentation thus received lots of attention in past few years. Although deep learning based state-of-the-arts have achieved great success in improving the segmentation accuracy, most of them suffer from an inefficiency problem and can hardly applied to practical applications. In this paper, we systematically analyze the computation cost of Convolutional Neural Network(CNN) and found that the inefficiency of CNN is mainly caused by its wide structure rather than the deep structure. In addition, the success of pruning based model compression methods proved that there are many redundant channels in CNN. Thus, we designed a very narrow while deep backbone network to improve the efficiency of semantic segmentation. By casting our network to FCN32 segmentation architecture, the basic structure of most segmentation methods, we achieved 60.6\% mIoU on Cityscape val dataset with 54 frame per seconds(FPS) on $1024\times2048$ inputs, which already outperforms one of the earliest real time deep learning based segmentation methods: ENet.


title: Spatiotemporal Feature Learning for Event-Based Vision

abstract: Unlike conventional frame-based sensors, event-based visual sensors output information through spikes at a high temporal resolution. By only encoding changes in pixel intensity, they showcase a low-power consuming, low-latency approach to visual information sensing. To use this information for higher sensory tasks like object recognition and tracking, an essential simplification step is the extraction and learning of features. An ideal feature descriptor must be robust to changes involving (i) local transformations and (ii) re-appearances of a local event pattern. To that end, we propose a novel spatiotemporal feature representation learning algorithm based on slow feature analysis (SFA). Using SFA, smoothly changing linear projections are learnt which are robust to local visual transformations. In order to determine if the features can learn to be invariant to various visual transformations, feature point tracking tasks are used for evaluation. Extensive experiments across two datasets demonstrate the adaptability of the spatiotemporal feature learner to translation, scaling and rotational transformations of the feature points. More importantly, we find that the obtained feature representations are able to exploit the high temporal resolution of such event-based cameras in generating better feature tracks.


title: Zero Forcing Detection For Short Packet Transmission Under Channel Estimation Errors

abstract: A multiuser multiple-input multiple-output wireless communication system is analytically studied under the short-packet transmission regime. The practical scenario of channel estimation errors is adopted when the signals undergo Rayleigh channel fading conditions. Unlike most previous works, the channel estimation error term is treated as a signal rather than interference or noise, which may further enhance the overall system performance. The spatial multiplexing mode of transmit operation is considered, while the zero-forcing detection is applied at the receiver. New and exact closed-form expressions for the outage probability and total system goodput are derived. Capitalizing on the analytical results, some new engineering insights are also presented; such as the impact of channel estimation errors with respect to the number of antennas, transmit power, and/or coding rate.


title: Pythia: a Framework for the Automated Analysis of Web Hosting Environments

abstract: A common approach when setting up a website is to utilize third party Web hosting and content delivery networks. Without taking this trend into account, any measurement study inspecting the deployment and operation of websites can be heavily skewed. Unfortunately, the research community lacks generalizable tools that can be used to identify how and where a given website is hosted. Instead, a number of ad hoc techniques have emerged, e.g., using Autonomous System databases, domain prefixes for CNAME records. In this work we propose Pythia, a novel lightweight approach for identifying Web content hosted on third-party infrastructures, including both traditional Web hosts and content delivery networks. Our framework identifies the organization to which a given Web page belongs, and it detects which Web servers are self-hosted and which ones leverage third-party services to provide contents. To test our framework we run it on 40,000 URLs and evaluate its accuracy, both by comparing the results with similar services and with a manually validated groundtruth. Our tool achieves an accuracy of 90% and detects that under 11% of popular domains are self-hosted. We publicly release our tool to allow other researchers to reproduce our findings, and to apply it to their own studies.


title: $X(3872)$ as virtual companion pole of the charm-anticharm state $\chi_{c1}(2P)$

abstract: We study the spectral function of the axial-vector charmonium state $\chi_{c1}(2P)$ coupled to $DD^{\ast}$ mesons, by employing a quantum field theoretical approach: a pronounced enhancement close to the $D^{0}% D^{\ast0}$ threshold, to be identified with the $X(3872)$, emerges. In the complex plane we find two poles: a pole for the broad seed state $\chi_{c1}(2P)$, and -- in the easiest scenario -- a virtual pole for the $X(3872)$. Thus, our approach describes both the seed state and the dynamically generated $X(3872)$ simultaneously. In particular, it explains the most prominent, both molecular-like and quarkonium-like, features of the $X(3872)$: its very small width (the decay into $D^{0}D^{\ast0}$ is predicted to be about 0.5 MeV), the enhanced radiative decay into $\psi(2S)\gamma$ w.r.t. $\psi(1S)\gamma$, and the isospin breaking decay into $J/\psi\rho$ (thanks to $DD^{\ast}$ loops mediating this decay channel). At the same time, we aim to determine the pole position and the properties of the charmonium seed state: quite interestingly, even if a pole is always present, it is possible that there is no peak corresponding to this state in the spectral function, thus potentially explaining why the corresponding resonance could not yet be seen in experiments.


title: $B_{s(d)}-\bar{B}_{s(d)}$ Mixing and $B_s\to\mu^+\mu^-$ Decay in the NMSSM with the Flavour Expansion Theorem

abstract: In this paper, motivated by the observation that the Standard Model predictions are now above the experimental data for the mass difference $\Delta M_{s(d)}$, we perform a detailed study of $B_{s(d)}-\bar{B}_{s(d)}$ mixing and $B_s\to\mu^+\mu^-$ decay in the $\mathbb{Z}_3$-invariant NMSSM with non-minimal flavour violation, using the recently developed procedure based on the Flavour Expansion Theorem, with which one can perform a purely algebraic mass-insertion expansion of an amplitude written in the mass eigenstate basis without performing any diagrammatic calculations in the interaction/flavour basis. Specifically, we consider the finite orders of mass insertions for neutralinos but the general orders for squarks and charginos, under two sets of assumptions for the squark flavour structures (\textit{i.e.}, while the flavour-conserving off-diagonal element $\delta_{33}^\text{LR}$ is kept in both of these two sectors, only the flavour-violating off-diagonal elements $\delta_{23}^\text{LL}$ and $\delta_{i3}^\text{RR}$ ($i=1,2$) are kept in the \text{LL} and \text{RR} sectors, respectively). Our analytic results are then expressed directly in terms of the initial Lagrangian parameters in the interaction/flavour basis, making it easy to impose the experimental bounds on them. It is found numerically that the NMSSM effects with the above two assumptions for the squark flavour structures can accommodate the observed deviation for $\Delta M_{s(d)}$, while complying with the experimental constraints from the branching ratios of $B_s\to \mu^+ \mu^-$ and $B\to X_s\gamma$ decays.


title: Active and Passive Portfolio Management with Latent Factors

abstract: We address a portfolio selection problem that combines active (outperformance) and passive (tracking) objectives using techniques from convex analysis. We assume a general semimartingale market model where the assets' growth rate processes are driven by a latent factor. Using techniques from convex analysis we obtain a closed-form solution for the optimal portfolio and provide a theorem establishing its uniqueness. The motivation for incorporating latent factors is to achieve improved growth rate estimation, an otherwise notoriously difficult task. To this end, we focus on a model where growth rates are driven by an unobservable Markov chain. The solution in this case requires a filtering step to obtain posterior probabilities for the state of the Markov chain from asset price information, which are subsequently used to find the optimal allocation. We show the optimal strategy is the posterior average of the optimal strategies the investor would have held in each state assuming the Markov chain remains in that state. Finally, we implement a number of historical backtests to demonstrate the performance of the optimal portfolio.


title: Unimodular multipliers on $\alpha$-modulation spaces: A revisit with new method under weaker conditions

abstract: By a new method derived from Nicola--Primo--Tabacco[24], we study the boundedness on $\alpha$-modulation spaces of unimodular multipliers with symbol $e^{i\mu(\xi)}$. Comparing with the previous results, the boundedness result is established for a larger family of unimodular multipliers under weaker assumptions.


title: Architecture for Integrating Learning Platforms Using Adapter

abstract: The advantage of the electronic and mobile learning platforms is the dissemination of learning contents with ease, but these operate differently to exchange the learning the learning contents from the server to the client. integrating these learning platforms to operate as a single platform and exchange the contents based on learners' request could improve the learning efficiency and reduce the operational cost. this work introduces a Web Services approach based on the client-server model to develop an integrated architecture that join the two learning platforms. in this paper the architecture of the learning platforms is presented and explained. furthermore, an adapter in a form of Web services is developed and as a middleware for the client-server communication.


title: Ordering properties of the smallest order statistic from Weibull G random variables

abstract: In this paper we compare the minimums of two heterogeneous samples each following Weibull-G distribution under three scenarios. In the Fifirst scenario, the units of the samples are assumed to be independently distributed and the comparisons are carried out through vector majorization. The minimums of the samples are compared in the second scenario when the independent units of the samples also experience random shocks. The last scenario describes the comparison when the units have a dependent structure sharing Archimedean copula.


title: On the distinguishable cluster approximation for triple excitations

abstract: The distinguishable cluster approximation applied to coupled cluster doubles equations greatly improves absolute and relative energies. We apply the same approximation to the triples equations and demonstrate that it can also improve results of the coupled cluster with singles, doubles and triples. The resulting method has a nominal computational scaling of ${\cal O}(N^7)$ in the real-space representation, is orbital invariant, size extensive, and exact for three electrons.


title: Orientation and structure of triple step staircase on vicinal Si(1 1 1) surfaces

abstract: The vicinal Si(111) surface, inclined towards the [-1-12] direction, was investigated by scanning tunnelling microscopy and spot profile analysing low energy electron diffraction. It has been established that the surface, consisting of regularly spaced triple steps and (111) terraces with a width equal to that of a single unit cell of the Si(111)-7x7 surface structure, has the (7 7 10) orientation. An atomic model of the triple step is proposed.


title: swCaffe: a Parallel Framework for Accelerating Deep Learning Applications on Sunway TaihuLight

abstract: This paper reports our efforts on swCaffe, a highly efficient parallel framework for accelerating deep neural networks (DNNs) training on Sunway TaihuLight, the current fastest supercomputer in the world that adopts a unique many-core heterogeneous architecture, with 40,960 SW26010 processors connected through a customized communication network. First, we point out some insightful principles to fully exploit the performance of the innovative many-core architecture. Second, we propose a set of optimization strategies for redesigning a variety of neural network layers based on Caffe. Third, we put forward a topology-aware parameter synchronization scheme to scale the synchronous Stochastic Gradient Descent (SGD) method to multiple processors efficiently. We evaluate our framework by training a variety of widely used neural networks with the ImageNet dataset. On a single node, swCaffe can achieve 23\%\~{}119\% overall performance compared with Caffe running on K40m GPU. As compared with the Caffe on CPU, swCaffe runs 3.04\~{}7.84x faster on all the networks. Finally, we present the scalability of swCaffe for the training of ResNet-50 and AlexNet on the scale of 1024 nodes.


title: Thermoelectric conversion at 30K in InAs/InP nanowire quantum dots

abstract: We demonstrate high-temperature thermoelectric conversion in InAs/InP nanowire quantum dots by taking advantage of their strong electronic confinement. The electrical conductance G and the thermopower S are obtained from charge transport measurements and accurately reproduced with a theoretical model accounting for the multi-level structure of the quantum dot. Notably, our analysis does not rely on the estimate of co-tunnelling contributions since electronic thermal transport is dominated by multi-level heat transport. By taking into account two spin-degenerate energy levels we are able to evaluate the electronic thermal conductance K and investigate the evolution of the electronic figure of merit ZT as a function of the quantum dot configuration and demonstrate ZT ~ 35 at 30 K, corresponding to an electronic effciency at maximum power close to the Curzon- Ahlborn limit.


title: EM-Based Smooth Graphon Estimation Using Bayesian and Spline-Based Approaches

abstract: This paper proposes the estimation of a smooth graphon model for network data analysis using principles of the EM algorithm. The approach considers both variability with respect to ordering the nodes of a network and smooth estimation of the graphon by nonparametric regression. To do so, (linear) B-splines are used, which allow for smooth estimation of the graphon, conditional on the node ordering. This provides the M-step. The true ordering of the nodes arising from the graphon model remains unobserved and Bayesian ideas are employed to obtain posterior samples given the network data. This yields the E-step. Combining both steps gives an EM-based approach for smooth graphon estimation. Unlike common other methods, this procedure does not require the restriction of a monotonic marginal function. The proposed graphon estimate allows to explore node-ordering strategies and therefore to compare the common degree-based node ranking with the ordering conditional on the network. Variability and uncertainty are taken into account using MCMC techniques. Examples and simulation studies support the applicability of the approach.


title: Besov-ish spaces through atomic decomposition

abstract: We use the method of atomic decomposition to build new families of function spaces, similar to Besov spaces, in measure spaces with grids, a very mild assumption. Besov spaces with low regularity are considered in measure spaces with good grids, and results on multipliers and left compositions are obtained.


title: Trans-Planckian Effects in Warm Inflation

abstract: We study the effect of a non-trivial vacuum prescription on warm inflation observables, namely the power spectrum of the comoving curvature perturbation. Non-trivial choice of vacuum, can provide information about trans-Planckian physics. Traditionally, the initial condition for inflation is chosen to be the usual Bunch-Davies vacuum. Another more reasonable choice of vacuum is the so-called $\alpha$-vacua. Because the duration of inflation is not infinite, as it is assumed in the Bunch-Davies case, imposing the initial condition at infinite past is not sensible and one must utilize another vacuum prescription. In this paper, working in the slow-roll regime during warm inflation, the initial condition for inflaton fluctuations is imposed at finite past, i.e. the $\alpha$-vacua. We show that this non-trivial vacuum prescription results in oscillatory correction to the comoving curvature power spectrum, which is scale dependent both in amplitude and frequency. Having obtained this scale dependent power spectrum, we consider its late time footprints and compare our results with observational data and other proposed models for the comoving curvature power spectrum.


title: Improving Lemmatization of Non-Standard Languages with Joint Learning

abstract: Lemmatization of standard languages is concerned with (i) abstracting over morphological differences and (ii) resolving token-lemma ambiguities of inflected words in order to map them to a dictionary headword. In the present paper we aim to improve lemmatization performance on a set of non-standard historical languages in which the difficulty is increased by an additional aspect (iii): spelling variation due to lacking orthographic standards. We approach lemmatization as a string-transduction task with an encoder-decoder architecture which we enrich with sentence context information using a hierarchical sentence encoder. We show significant improvements over the state-of-the-art when training the sentence encoder jointly for lemmatization and language modeling. Crucially, our architecture does not require POS or morphological annotations, which are not always available for historical corpora. Additionally, we also test the proposed model on a set of typologically diverse standard languages showing results on par or better than a model without enhanced sentence representations and previous state-of-the-art systems. Finally, to encourage future work on processing of non-standard varieties, we release the dataset of non-standard languages underlying the present study, based on openly accessible sources.


title: Predicted electric field induced surface magnetization of vortex phase of superconductors

abstract: We predict a new effect in condensed matter surface magnetization of the vortex phase of a superconductor induced by electric field. The magnetized superconductor should be one of the plates of a plane capacitor on which a voltage has to be applied. Applying alternating voltage to the capacitor, electrostatic induction leads to oscillations of the magnetic moment which has to be measured by electromotive force inducted in a detector coil. The derived explicit formula for the magnetization contains the effective mass of Cooper pairs and a systematic investigation of the predicted magnetization will lead to a creation of an effective Cooper pair mass spectroscopy. For cleaved superconductors this effective mass is a property of the bulk material.


title: Classic and exotic Besov spaces induced by good grids

abstract: In a previous work we introduced Besov spaces $\mathcal{B}^s_{p,q}$ defined on a measure spaces with a good grid, with $p\in [1,\infty)$, $q\in [1,\infty]$ and $0< s< 1/p$. Here we show that classical Besov spaces on compact homogeneous spaces are examples of such Besov spaces. On the other hand we show that even Besov spaces defined by a good grid made of partitions by intervals may differ from a classical Besov space, giving birth to exotic Besov spaces.


title: Lights Out on graphs

abstract: We model the Lights Out game on general simple graphs in the framework of linear algebra over the field $\mathbb F_2$. Based upon a version of the Fredholm alternative, we introduce a separating invariant of the game, i.e., an initial state can be transformed into a final state if and only if the invariant of both states agrees. We also investigate certain states with particularly interesting properties. Apart from the classical version of the game, we propose several variants, in particular a version with more than only two states (light on, light off), where the analysis resides on systems of linear equations over the ring $\mathbb Z_n$. Although it is easy to find a concrete solution of the Lights Out problem, we show that it is NP-hard to find a minimal solution. We also propose electric circuit diagrams to actually realize the Lights Out game.


title: Transfer operators and atomic decomposition

abstract: We use the method of atomic decomposition and a new family of Banach spaces to study the action of transfer operators associated to piecewise-defined maps. It turns out that these transfer operators are quasi-compact even when the associated potential, the dynamics and the underlying phase space have very low regularity. In particular it is often possible to obtain exponential decay of correlations, the Central Limit Theorem and almost sure invariance principle for fairly general observables, including unbounded ones.


title: Astro2020 Science White Paper: Synoptic Studies of the Sun as a Key to Understanding Stellar Astrospheres

abstract: Ground-based solar observations provide key contextual data (i.e., the 'big picture') to produce a complete description of the only astrosphere we can study in situ: our Sun's heliosphere. The next decade will see the beginning of operations of the Daniel K. Inouye Solar Telescope (DKIST). DKIST will join NASA's Parker Solar Probe and the NASA/ESA Solar Orbital mission, which together will study our Sun's atmosphere with unprecedented detail. This white paper outlines the current paradigm for ground-based solar synoptic observations, and indicates those areas that will benefit from focused attention.


title: Thermodynamic Properties of Graphene in Magnetic Field and Rashba Coupling

abstract: We study the thermodynamic properties of massless Dirac fermions in graphene subjected to a uniform magnetic field $B$ together with Rashba coupling parameter $\lambda_R$. The thermodynamic functions such as the Helmholtz free energy, total energy, entropy and heat capacity are obtained in the high temperature regime using an approach based on the zeta function. These functions will be numerically examined by considering two cases related to $\lambda_R$ smaller or greater than $B$. In particular, we show that the Dulong-Petit law is verified for both cases.

categories: cond-mat.mes-hall math-ph math.MP quant-ph


title: Unsupervised Part-Based Disentangling of Object Shape and Appearance

abstract: Large intra-class variation is the result of changes in multiple object characteristics. Images, however, only show the superposition of different variable factors such as appearance or shape. Therefore, learning to disentangle and represent these different characteristics poses a great challenge, especially in the unsupervised case. Moreover, large object articulation calls for a flexible part-based model. We present an unsupervised approach for disentangling appearance and shape by learning parts consistently over all instances of a category. Our model for learning an object representation is trained by simultaneously exploiting invariance and equivariance constraints between synthetically transformed images. Since no part annotation or prior information on an object class is required, the approach is applicable to arbitrary classes. We evaluate our approach on a wide range of object categories and diverse tasks including pose prediction, disentangled image synthesis, and video-to-video translation. The approach outperforms the state-of-the-art on unsupervised keypoint prediction and compares favorably even against supervised approaches on the task of shape and appearance transfer.


title: An energy-based discontinuous Galerkin method for the wave equation with advection

abstract: An energy-based discontinuous Galerkin method for the advective wave equation is proposed and analyzed. Energy-conserving or energy-dissipating methods follow from simple, mesh-independent choices of the inter-element fluxes, and both subsonic and supersonic advection is allowed. Error estimates in the energy norm are established, and numerical experiments on structured grids display optimal convergence in the $L^2$ norm for upwind fluxes. The method generalizes earlier work on energy-based discontinuous Galerkin methods for second order wave equations which was restricted to energy forms written as a simple sum of kinetic and potential energy.


title: Coding in graphs and linear orderings

abstract: There is a Turing computable embedding $\Phi$ of directed graphs $A$ in undirected graphs. Moreover, there is a fixed tuple of formulas that give a uniform interpretation; i.e., for all directed graphs $A$, these formulas interpret $A$ in $\Phi(G)$. It follows that A is Medvedev reducible to $\Phi(A)$ uniformly; i.e., there is a fixed Turing operator that serves for all $A$. We observe that there is a graph $G$ that is not Medvedev reducible to any linear ordering. Hence, $G$ is not effectively interpreted in any linear ordering. Similarly, there is a graph that is not interpreted in any linear ordering using computable $\Sigma_2$ formulas. Any graph can be interpreted in a linear ordering using computable $\Sigma_3$ formulas. Friedman and Stanley gave a Turing computable embedding L of directed graphs in linear orderings. We show that there is no fixed tuple of $L_{\omega_1,\omega}$ formulas that, for all $G$, interpret the input graph $G$ in the output linear ordering $L(G)$. Harrison-Trainor and Montalb\'an have also shown this, by a quite different proof.


title: Hand range of motion evaluation for Rheumatoid Arthritis patients

abstract: We introduce a framework for dynamic evaluation of the fingers movements: flexion, extension, abduction and adduction. This framework estimates angle measurements from joints computed by a hand pose estimation algorithm using a depth sensor (Realsense SR300). Given depth maps as input, our framework uses Pose-REN, which is a state-of-art hand pose estimation method that estimates 3D hand joint positions using a deep convolutional neural network. The pose estimation algorithm runs in real-time, allowing users to visualise 3D skeleton tracking results at the same time as the depth images are acquired. Once 3D joint poses are obtained, our framework estimates a plane containing the wrist and MCP joints and measures flexion/extension and abduction/aduction angles by applying computational geometry operations with respect to this plane. We analysed flexion and abduction movement patterns using real data, extracting the movement trajectories. Our preliminary results show that this method allows an automatic discrimination of hands with Rheumatoid Arthritis (RA) and healthy patients. The angle between joints can be used as an indicative of current movement capabilities and function. Although the measurements can be noisy and less accurate than those obtained statically through goniometry, the acquisition is much easier, non-invasive and patient-friendly, which shows the potential of our approach. The system can be used with and without orthosis. Our framework allows the acquisition of measurements with minimal intervention and significantly reduces the evaluation time.


title: A Partition-centric Distributed Algorithm for Identifying Euler Circuits in Large Graphs

abstract: Finding the Eulerian circuit in graphs is a classic problem, but inadequately explored for parallel computation. With such cycles finding use in neuroscience and Internet of Things for large graphs, designing a distributed algorithm for finding the Euler circuit is important. Existing parallel algorithms are impractical for commodity clusters and Clouds. We propose a novel partition-centric algorithm to find the Euler circuit, over large graphs partitioned across distributed machines and executed iteratively using a Bulk Synchronous Parallel (BSP) model. The algorithm finds partial paths and cycles within each partition, and refines these into longer paths by recursively merging the partitions. We describe the algorithm, analyze its complexity, validate it on Apache Spark for large graphs, and offer experimental results. We also identify memory bottlenecks in the algorithm and propose an enhanced design to address it.


title: Weak gravitational lensing by two-power-law densities using the Gauss-Bonnet theorem

abstract: We study the weak-field deflection of light by mass distributions described by two-power-law densities $\rho(R)=\rho_0 R^{-\alpha}(R+1)^{\beta-\alpha}$, where $\alpha$ and $\beta$ are non-negative integers. New analytic expressions of deflection angles are obtained via the application of the Gauss-Bonnet theorem to a chosen surface on the optical manifold. Some of the well-known models of this two-power law form are the Navarro-Frenk-White (NFW) model $(\alpha,\beta)=(1,3)$, Hernquist $(1,4)$, Jaffe $(2,4)$, and the singular isothermal sphere $(2,2)$. The calculated deflection angles for Hernquist and NFW agrees with that of Keeton and Bartelmann, respectively. The limiting values of these deflection angles (at zero or infinite impact parameter) are either vanishing or similar to the deflection due to a singular isothermal sphere. We show that these behaviors can be attributed to the topological properties of the optical manifold, thus extending the pioneering insight of Werner and Gibbons to a broader class of mass densities.


title: Correlated states in twisted double bilayer graphene

abstract: Electron-electron interactions play an important role in graphene and related systems and can induce exotic quantum states, especially in a stacked bilayer with a small twist angle. For bilayer graphene where the two layers are twisted by a "magic angle", flat band and strong many-body effects lead to correlated insulating states and superconductivity. In contrast to monolayer graphene, the band structure of untwisted bilayer graphene can be further tuned by a displacement field, providing an extra degree of freedom to control the flat band that should appear when two bilayers are stacked on top of each other. Here, we report the discovery and characterization of such displacement-field tunable electronic phases in twisted double bilayer graphene. We observe insulating states at a half-filled conduction band in an intermediate range of displacement fields. Furthermore, the resistance gap in the correlated insulator increases with respect to the in-plane magnetic fields and we find that the g factor according to spin Zeeman effect is ~2, indicating spin polarization at half filling. These results establish the twisted double bilayer graphene as an easily tunable platform for exploring quantum many-body states.


title: State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics

abstract: Recently proposed non-adiabatic ring polymer molecular dynamics (NRPMD) approach has shown to provide accurate quantum dynamics by incorporating explicit electronic state descriptions and nuclear quantizations. Here, we present a rigorous derivation of the NRPMD Hamiltonian and investigate its performance on simulating excited state non-adiabatic dynamics. Our derivation is based on the Meyer-Miller-Stock-Thoss (MMST) mapping representation for electronic states and the ring-polymer path-integral description for nuclei, resulting in the same Hamiltonian proposed in the original NRPMD approach. In addition, we investigate the accuracy of using NRPMD to simulate photoinduced non-adiabatic dynamics in simple model systems. These model calculations suggest that NRPMD can alleviate the zero-point energy leakage problem that is commonly encountered in the classical Wigner dynamics, and provide accurate excited states non-adiabatic dynamics. This work provides a solid theoretical foundation of the promising NRPMD Hamiltonian and demonstrates the possibility of using state-dependent RPMD approach to accurately simulate electronic non-adiabatic dynamics while explicitly quantize nuclei.


title: Genuine time-bin-encoded quantum key distribution over a turbulent depolarizing free-space channel

abstract: Despite its widespread use in fiber optics, encoding quantum information in photonic time-bin states is usually considered impractical for free-space quantum communication as turbulence-induced spatial distortion impedes the analysis of time-bin states at the receiver. Here, we demonstrate quantum key distribution using time-bin photonic states distorted by turbulence and depolarization during free-space transmission. Utilizing a novel analyzer apparatus, we observe stable quantum bit error ratios of 5.32 %, suitable for generating secure keys, despite significant wavefront distortions and polarization fluctuations across a 1.2 km channel. This shows the viability of time-bin quantum communication over long-distance free-space channels, which will simplify direct fiber/free-space interfaces and enable new approaches for practical free-space quantum communication over multi-mode, turbulent, or depolarizing channels.


title: Homotopy Reconstruction via the Cech Complex and the Vietoris-Rips Complex

abstract: We derive conditions under which the reconstruction of a target space is topologically correct via the \v{C}ech complex or the Vietoris-Rips complex obtained from possibly noisy point cloud data. We provide two novel theoretical results. First, we describe sufficient conditions under which any non-empty intersection of finitely many Euclidean balls intersected with a positive reach set is contractible, so that the Nerve theorem applies for the restricted \v{C}ech complex. Second, we demonstrate the homotopy equivalence of a positive $\mu$-reach set and its offsets. Applying these results to the restricted \v{C}ech complex and using the interleaving relations with the \v{C}ech complex (or the Vietoris-Rips complex), we formulate conditions guaranteeing that the target space is homotopy equivalent to the \v{C}ech complex (or the Vietoris-Rips complex), in terms of the $\mu$-reach. Our results sharpen existing results.


title: Spontaneous photon-pair generation at the nanoscale

abstract: Optical nanoantennas have shown a great capacity for efficient extraction of photons from the near to the far-field, enabling directional emission from nanoscale single-photon sources. However, their potential for the generation and extraction of multi-photon quantum states remains unexplored. Here we demonstrate experimentally the nanoscale generation of two-photon quantum states at telecommunication wavelengths based on spontaneous parametric down-conversion in an optical nanoantenna. The antenna is a crystalline AlGaAs nanocylinder, possessing Mie-type resonances at both the pump and the bi-photon wavelengths and when excited by a pump beam generates photonpairs with a rate of 35 Hz. Normalized to the pump energy stored by the nanoantenna, this rate corresponds to 1.4 GHz/Wm, being one order of magnitude higher than conventional on-chip or bulk photon-pair sources. Our experiments open the way for multiplexing several antennas for coherent generation of multi-photon quantum states with complex spatial-mode entanglement and applications in free-space quantum communications and sensing.


title: The Second AGILE Catalog of Gamma-Ray Sources

abstract: Aims. We present the second AGILE-GRID Catalog (2AGL) of {\gamma}-ray sources in the 100 MeV-10 GeV energy range. Methods. With respect to previous AGILE-GRID catalogs, the current 2AGL Catalog is based on the first 2.3 years of science data from the AGILE mission (the so called 'pointing mode') and incorporates more data and several analysis improvements, including better calibrations at the event reconstruction level, an updated model for the Galactic diffuse gamma-ray emission, a refined procedure for point-like source detection, and the inclusion of a search for extended {\gamma}-ray sources. Results. The 2AGL Catalog includes 175 high-confidence sources (above 4{\sigma} significance) with their location regions and spectral properties, and a variability analysis with 4-day light curves for the most significant ones. Relying on the error region of each source position, including systematic uncertainties, 121 sources are considered as positionally associated with known couterparts at different wavelengths or detected by other {\gamma}-ray instruments. Among the identified or associated sources, 62 are Active Galactic Nuclei (AGNs) of the blazar class. Pulsars represent the largest Galactic source class, with 40 associated pulsars, 7 of them with detected pulsation; 8 Supernova Remnants and 4 high-mass X-ray binaries have also been identified. A substantial number of 2AGL sources are unidentified: for 54 sources no known counterpart is found at different wavelengths. Among these sources, we discuss a sub-class of 29 AGILE-GRID-only {\gamma}-ray sources that are not present in 1FGL, 2FGL or 3FGL catalogs; the remaining sources are unidentified in both 2AGL and 3FGL Catalogs. We also present an extension of the analysis of 2AGL sources detected in the 50-100 MeV energy range.


title: Early-career setback and future career impact

abstract: Setbacks are an integral part of a scientific career, yet little is known about whether an early-career setback may augment or hamper an individual's future career impact. Here we examine junior scientists applying for U.S. National Institutes of Health (NIH) R01 grants. By focusing on grant proposals that fell just below and just above the funding threshold, we compare "near-miss" with "near-win" individuals to examine longer-term career outcomes. Our analyses reveal that an early-career near miss has powerful, opposing effects. On one hand, it significantly increases attrition, with one near miss predicting more than a 10% chance of disappearing permanently from the NIH system. Yet, despite an early setback, individuals with near misses systematically outperformed those with near wins in the longer run, as their publications in the next ten years garnered substantially higher impact. We further find that this performance advantage seems to go beyond a screening mechanism, whereby a more selected fraction of near-miss applicants remained than the near winners, suggesting that early-career setback appears to cause a performance improvement among those who persevere. Overall, the findings are consistent with the concept that "what doesn't kill me makes me stronger." Whereas science is often viewed as a setting where early success begets future success, our findings unveil an intimate yet previously unknown relationship where early-career setback can become a marker for future achievement, which may have broad implications for identifying, training and nurturing junior scientists whose career will have lasting impact.


title: Induced exchange and spin-orbit effects by proximity in graphene on Ni and Co

abstract: The induced-proximity effects of nearly commensurate lattice structure of a graphene layer on Ni(111) and Co(0001) substrates in the AC stacking configuration are addressed through an analytical tight-binding approach within the Slater-Koster method. A minimal Hamiltonian is constructed by considering the hybridizations of the magnetic $3d$-orbitals of Ni(Co) atoms with the $p_z$-orbitals of graphene, in addition to the atomic spin-orbit coupling and the magnetization of the Ni(Co) atoms. A low-energy effective Hamiltonian for graphene/Ni(Co) describing the perturbed $\pi$-bands in the vicinity of the Dirac points is derived which enable us to get further insight on the physical nature of the induced-effective couplings to the graphene layer. It is shown that a magneto-spin-orbit type effect may emerge through two competing mechanisms simultaneously present, namely the proximity induced exchange and Rashba spin-orbit interaction. Such effects results in giant exchange splittings and robust Rashba spin-orbit coupling transferred to the graphene layer in agreement with recent density functional theory calculations and experimental observations. We further analyze the physical conditions for the appearance of intact Dirac cones in the minority spin bands as observed by recent photoemission measurements with spin resolution.


title: Reversible and non-reversible Markov Chain Monte Carlo algorithms for reservoir simulation problems

abstract: We compare numerically the performance of reversible and non-reversible Markov Chain Monte Carlo algorithms for high dimensional oil reservoir problems; because of the nature of the problem at hand, the target measures from which we sample are supported on bounded domains. We compare two strategies to deal with bounded domains, namely reflecting proposals off the boundary and rejecting them when they fall outside of the domain. We observe that for complex high dimensional problems reflection mechanisms outperform rejection approaches and that the advantage of introducing non-reversibility in the Markov Chain employed for sampling is more and more visible as the dimension of the parameter space increases.


title: Entropy modulo a prime

abstract: Building on work of Kontsevich, we introduce a definition of the entropy of a finite probability distribution in which the "probabilities" are integers modulo a prime p. The entropy, too, is an integer mod p. Entropy mod p is shown to be uniquely characterized by a functional equation identical to the one that characterizes ordinary Shannon entropy. We also establish a sense in which certain real entropies have residues mod p, connecting the concepts of entropy over R and over Z/pZ. Finally, entropy mod p is expressed as a polynomial which is shown to satisfy several identities, linking into work of Cathelineau, Elbaz-Vincent and Gangl on polylogarithms.


title: Evolution of the Gas Mass Fraction of Progenitors to Today's Massive Galaxies: ALMA Observations in the CANDELS GOODS-S Field

abstract: We present an ALMA survey of dust continuum emission in a sample of 70 galaxies in the redshift range z=2-5 selected from the CANDELS GOODS-S field. Multi-Epoch Abundance Matching (MEAM) is used to define potential progenitors of a z = 0 galaxy of stellar mass 1.5 10^11 M_sun. Gas masses are derived from the 850um luminosity. Ancillary data from the CANDELS GOODS-S survey are used to derive the gas mass fractions. The results at z<=3 are mostly in accord with expectations: The detection rates are 75% for the z=2 redshift bin, 50% for the z=3 bin and 0% for z>=4. The average gas mass fraction for the detected z=2 galaxies is f_gas = 0.55+/-0.12 and f_gas = 0.62+/-0.15 for the z=3 sample. This agrees with expectations for galaxies on the star-forming main sequence, and shows that gas fractions have decreased at a roughly constant rate from z=3 to z=0. Stacked images of the galaxies not detected with ALMA give upper limits to f_gas of <0.08 and <0.15, for the z=2 and z=3 redshift bins. None of our galaxies in the z=4 and z=5 sample are detected and the upper limit from stacked images, corrected for low metallicity, is f_gas<0.66. We do not think that lower gas-phase metallicities can entirely explain the lower dust luminosities. We briefly consider the possibility of accretion of very low-metallicity gas to explain the absence of detectable dust emission in our galaxies at z>4.


title: Imbalanced multi-label classification using multi-task learning with extractive summarization

abstract: Extractive summarization and imbalanced multi-label classification often require vast amounts of training data to avoid overfitting. In situations where training data is expensive to generate, leveraging information between tasks is an attractive approach to increasing the amount of available information. This paper employs multi-task training of an extractive summarizer and an RNN-based classifier to improve summarization and classification accuracy by 50% and 75%, respectively, relative to RNN baselines. We hypothesize that concatenating sentence encodings based on document and class context increases generalizability for highly variable corpuses.


title: Fast Markov chain Monte Carlo for high dimensional Bayesian regression models with shrinkage priors

abstract: In the past decade, many Bayesian shrinkage models have been developed for linear regression problems where the number of covariates, $p$, is large. Computing the intractable posterior are often done with three-block Gibbs samplers (3BG), based on representing the shrinkage priors as scale mixtures of Normal distributions. An alternative computing tool is a state of the art Hamiltonian Monte Carlo (HMC) method, which can be easily implemented in the Stan software. However, we found both existing methods to be inefficient and often impractical for large $p$ problems. Following the general idea of Rajaratnam et al. (2018), we propose two-block Gibbs samplers (2BG) for three commonly used shrinkage models, namely, the Bayesian group lasso, the Bayesian sparse group lasso and the Bayesian fused lasso models. We demonstrate with simulated and real data examples that the Markov chains underlying 2BG's converge much faster than that of 3BG's, and no worse than that of HMC. At the same time, the computing costs of 2BG's per iteration are as low as that of 3BG's, and can be several orders of magnitude lower than that of HMC. As a result, the newly proposed 2BG is the only practical computing solution to do Bayesian shrinkage analysis for datasets with large $p$. Further, we provide theoretical justifications for the superior performance of 2BG's. We establish geometric ergodicity (GE) of Markov chains associated with the 2BG for each of the three Bayesian shrinkage models. We also prove, for most cases of the Bayesian group lasso and the Bayesian sparse group lasso model, the Markov operators for the 2BG chains are trace-class. Whereas for all cases of all three Bayesian shrinkage models, the Markov operator for the 3BG chains are not even Hilbert-Schmidt.


title: Feather: A Feature Model Transformation Language

abstract: Feature modeling has been a very popular approach for variability management in software product lines. Building a feature model requires substantial domain expertise, however, even experts cannot foresee all future possibilities. Changing requirements can force a feature model to evolve in order to adapt to the new conditions. Feather is a language to describe model transformations that will evolve a feature model. This article presents the structure and foundations of Feather. First, the language elements, which consist of declarations to characterize the model to evolve and commands to manipulate its structure, are introduced. Then, semantics grounding in feature model properties are given for the commands in order to provide precise command definitions. Next, an interpreter that can realize the transformations described by the commands in a Feather script is presented. Finally, effectiveness of the language is discussed using two realistic examples, where one of the examples includes a system from a dynamic environment and the other employs a system that has a large feature model containing 1,227 features.


title: Invariant Manifolds and Rate Constants in Driven Chemical Reactions

abstract: Reaction rates of chemical reactions under nonequilibrium conditions can be determined through the construction of the normally hyperbolic invariant manifold (NHIM) [and moving dividing surface (DS)] associated with the transition state trajectory. Here, we extend our recent methods by constructing points on the NHIM accurately even for multidimensional cases. We also advance the implementation of machine learning approaches to construct smooth versions of the NHIM from a known high-accuracy set of its points. That is, we expand on our earlier use of neural nets, and introduce the use of Gaussian process regression for the determination of the NHIM. Finally, we compare and contrast all of these methods for a challenging two-dimensional model barrier case so as to illustrate their accuracy and general applicability.


title: On the hypersurfaces contained in their Hessian

abstract: This article presents the theory of focal locus applied to the hypersurfaces in the projective space which are (finitely) covered by linear spaces and such that the tangent space is constant along these spaces.


title: Creation of discontinuities in circle maps

abstract: Circle maps frequently arise in mathematical models of physical or biological systems. Motivated by Cherry flows and `threshold' systems such as integrate and fire neuronal models, models of cardiac arrhythmias, and models of sleep/wake regulation, we consider how structural transitions in circle maps occur. In particular we describe how maps evolve near the creation of a discontinuity. We show that the natural way to create discontinuities in the maps associated with both threshold systems and Cherry flows results in a square root singularity in the map as the discontinuity is approached from either one or both sides. We analyse the generic properties of maps with gaps displaying this local behaviour, showing how border collisions and saddle-node bifurcations are interspersed. This highlights how the Arnold tongue picture for tongues bordered by saddle-node bifurcations is amended once gaps are present. For the threshold models we also show that a loss of injectivity naturally results in the creation of multiple gaps giving rise to a novel codimension two bifurcation.


title: Domain adaptation for holistic skin detection

abstract: Human skin detection in images is a widely studied topic of Computer Vision for which it is commonly accepted that analysis of pixel color or local patches may suffice. This is because skin regions appear to be relatively uniform and many argue that there is a small chromatic variation among different samples. However, we found that there are strong biases in the datasets commonly used to train or tune skin detection methods. Furthermore, the lack of contextual information may hinder the performance of local approaches. In this paper we present a comprehensive evaluation of holistic and local Convolutional Neural Network (CNN) approaches on in-domain and cross-domain experiments and compare with state-of-the-art pixel-based approaches. We also propose a combination of inductive transfer learning and unsupervised domain adaptation methods, which are evaluated on different domains under several amounts of labelled data availability. We show a clear superiority of CNN over pixel-based approaches even without labelled training samples on the target domain. Furthermore, we provide experimental support for the counter-intuitive superiority of holistic over local approaches for human skin detection.


title: On the convergence of stochastic MPC to terminal modes of operation

abstract: The stability of stochastic Model Predictive Control (MPC) subject to additive disturbances is often demonstrated in the literature by constructing Lyapunov-like inequalities that guarantee closed-loop performance bounds and boundedness of the state, but convergence to a terminal control law is typically not shown. In this work we use results on general state space Markov chains to find conditions that guarantee convergence of disturbed nonlinear systems to terminal modes of operation, so that they converge in probability to a priori known terminal linear feedback laws and achieve time-average performance equal to that of the terminal control law. We discuss implications for the convergence of control laws in stochastic MPC formulations, in particular we prove convergence for two formulations of stochastic MPC.


title: Sums and Products of Regular Polytopes' Squared Chord Lengths

abstract: Although previous research has found several facts concerning chord lengths of regular polytopes, none of these investigations has considered whether any of these facts define relationships that might generalize to the chord lengths of all regular polytopes. Consequently, this paper explores whether four findings of previous studies-- viz., the four facts relating to the sums and products of squared chord lengths of regular polygons inscribed in unit circles-- can be generalized to all regular ($n$-dimensional) polytopes (inscribed in unit $n$-spheres). We show that (a) one of these four facts actually does generalize to all regular polytopes (of dimension $n \geq 2$), (b) one generalizes to all regular polytopes except most simplices, (c) one generalizes only to the family of crosspolytopes, and (d) one generalizes only to the crosspolytopes and 24-cell. We also discover several corollaries (due to reciprocation) and some theorems specific to the three-dimensional regular polytopes along the way.